PCR cloning of a histone H1 gene from Anopheles stephensi mosquito cells: comparison of the protein sequence with histone H1-like, C-terminal extensions on mosquito ribosomal protein S6
© Zhai and Fallon; licensee BioMed Central Ltd. 2005
Received: 29 July 2004
Accepted: 24 January 2005
Published: 24 January 2005
In Aedes and Anopheles mosquitoes, ribosomal protein RPS6 has an unusual C-terminal extension that resembles histone H1 proteins. To explore homology between a mosquito H1 histone and the RPS6 tail, we took advantage of the Anopheles gambiae genome database to clone a histone H1 gene from an Anopheles stephensi mosquito cell line.
We designed specific primers based on RPS6 and histone H1 alignments to recover an Anopheles stephensi histone H1 corresponding to a conceptual An. gambiae protein, with 92% identity. Southern blots suggested that Anopheles stephensi histone H1 gene has multiple variants, as is also the case for histone H1 proteins in Chironomid flies.
Histone H1 proteins from Anopheles stephensi and Anopheles gambiae mosquitoes share 92% identity to each other, but only 50% identity to a Drosophila homolog. In a phylogenetic analysis, Anopheles, Chironomus and Drosophila histone H1 proteins cluster separately from the histone H1-like, C-terminal tails on RPS6 in Aedes and Anopheles mosquitoes. These observations suggest that the resemblance between histone H1 and the C-terminal extensions on mosquito RPS6 has been maintained by convergent evolution.
Ribosomal protein (RP) S6 is a phosphorylated protein that resides on the small subunit of eukaryotic ribosomes. Phosphorylation occurs on a cluster of five serine residues near the C-terminal end of the protein. Although details remain unclear, the phosphorylation state of RPS6 is believed to influence translational efficiency of some mRNAs , possibly mediated by direct contact between RPS6 and the 28S rRNA in the large subunit. RPS6 has also been implicated in ribosome biogenesis, and is thought to play a conserved role in the initiation of protein synthesis .
In Aedes aegypti and Aedes albopictus mosquitoes, the RPS6 protein is ~17 kDa larger than its Drosophila homolog, and on polyacrylamide gels, it migrates as the largest protein from the small ribosomal subunit. Ae. aegypti and Ae. albopictus RPS6 cDNAs encode an approximately 100 amino acid extension at the C-terminal end of the protein. The extension is particularly rich in lysine, alanine and glutamic acid, and most closely resembles the sequence of histone H1 proteins from diverse sources .
Because RPS6 is thought to have regulatory function(s) in a variety of cell signaling pathways , we were surprised to uncover this difference between mosquito and Drosophila RPS6 proteins. We have recently shown that RPS6 protein isolated from ribosomal subunits retains its histone H1-like tail . Thus, unlike the case with the ubiquitinated ribosomal protein S27a in the rat , the histone tail is not removed from the mosquito ribosomal protein prior to ribosome assembly.
RpS6 cDNA from an Anopheles stephensi cell line encodes an approximately 170 amino acid histone H1-like C-terminal extension, and in silico analysis reveals a similar modification encoded by the rpS6 gene in Anopheles gambiae. In both Aedes and Anopheles mosquitoes, the C-terminal extension was completely encoded in Exon 3, directly contiguous with upstream open reading frame encoding the series of serines that may be phosphorylated . Anopheline mosquitoes are believed to be ancestral to the Culicidae, which includes the genera Aedes and Culex . Thus, to a first approximation, we infer that the longer tail in Anopheles mosquitoes represents the ancestral state, and that the RPS6 tail has been lost in the higher Diptera, which include D. melanogaster.
Although mosquito RPS6 tails in general resemble histone H1 proteins, their divergence between Aedes and Anopheles mosquitoes was high, relative to the conventional portion of the RPS6 coding sequence. Because histone H1 is the most variable of the histone proteins, and functions as a linker, rather than as a component of the histone octamer, we set out to clone a cDNA encoding a bona fide histone H1 protein from an An. stephensi cell line. In a phylogenetic comparison, the An. stephensi histone H1 protein clusters with homologs from Drosophila and Chironomus, rather than with RPS6 histone H1-like tails from mosquitoes. These results indicate that the histone H1-like tails on mosquito RPS6 proteins are evolving independently of conspecific histone H1 proteins.
Design of PCR primers
The gene encoding Drosophila melanogaster histone H1 spans 1204 nucleotides, and encodes a 256 amino acid protein in a single exon . There is a single recorded His1 allele in Drosophila , while multiple histone H1 variants have been described in Chironomid flies [9–11]. When the deduced sequence of the Drosophila histone H1 protein (Accession NM_058232) was compared to the Anopheles gambiae genome using the program BLAST  on the NCBI website (National Center for Biotechnology Information; http://www.ncbi.nlm.nih.gov/), we obtained 5 accessions with E values ranging from 3e-35 to 8e-43, distributed on mosquito chromosomes 2 and 3. Upon further examination, we noted that XP_314184 and XP_314186 (chromosome 2) corresponded to the same protein. Two additional histone H1 candidates (XP_311486 and XP_309451) were encoded on chromosome 3. These three conceptual Anopheles proteins shared 70–80% identity to one another, and about 50% identity to the Drosophila H1 protein sequence. In the EST-other database, we found a single uninformative match to an unidentified An. gambiae entry (dbEST id = 11236311), with the relatively modest E value of 0.055. Histone H1 sequences from Aedes mosquitoes are not yet in existing databases.
The 50% identity between Drosophila and Anopheles histone H1 proteins was relatively low, compared to approximately 80% amino acid identity between Drosophila and Anopheles RPS6, exclusive of the histone-H1-like tail in the mosquito protein. The Drosophila H1 histone was also ~50% identical to that from Chironomus thummi, a fly closely related to mosquitoes in the infraorder/superfamily Culicomorpha .
Recovery of An. stephensi histone H1 gene
We used F1 and R1 primers with Hin dIII-digested genomic DNA from An. stephensi cells to obtain an approximately 450 bp PCR product, which was sequenced and verified to encode a histone H1 protein. The 5-end of the gene, which extended 81 nucleotides upstream of the ATG start codon, was obtained using primer R1 with the GeneRacer kit (Invitrogen, Carlsbad, CA), with total RNA as the template. The absence of a poly(A) tail on histone mRNAs required an unconventional strategy to obtain the 3'-end of the coding sequence. First, we used Hin dIII-digested genomic DNA template, with a primer based entirely on the 3'-UTR of An. gambiae XP_314184, without success. When we designed a second primer (R2, in Fig. 2) extending from the 3'-UTR through the TAA stop codon and into the coding region, we obtained the 3'-end of the coding sequence. Finally, primers F2 and R2 (Fig. 2) were used to verify the complete nucleotide sequence.
Southern blots with An. stephensi genomic DNA
Comparisons of histone H1 proteins with mosquito RPS6 C-terminal extensions
The identity between Drosophila and Anopheles (or Drosophila and Chironomus) histone H1 proteins was only 50%. This divergence undoubtedly reflects the ~250 million years  separating Nematoceran from Cyclorrhaphan diptera. In this study, we were interested in comparing mosquito histone H1 proteins to the histone H1-like tails of mosquito RPS6. Fig. 4B shows a neighbor-joining analysis in which we compared protein sequences from Aedes and Anopheles RPS6 histone H1-like tails, exclusive of the conventional RPS6 protein sequence, with histone H1 proteins from the nematode Caenorhabditis elegans (AAM44399), the closely-related flies Chironomus thummi (Q07134) and Chironomus tentans (AAB62239), Drosophila, and the Anopheles gambiae and Anopheles stephensi homologs (Fig. 4A). With the C. elegans sequence designated as the outgroup, the phylogram shows that the RPS6 tails cluster into a distinct group relative to the Dipteran histone H1 proteins. Circled values indicate bootstrap values based on 1000 replicates. When the analysis was repeated with the optimality criterion set to parsimony, we obtained a tree with the same topology, with the 77% value shown in Fig. 4B reduced to 59%, and the 97% value reduced to 94%. The 100% values remained unchanged.
An important rationale for cloning an An. stephensi histone H1 was to compare its sequence to the histone H1-like tails on mosquito RPS6 ribosomal proteins. Our choice of an Anopheles histone H1 was based on the existing database for An. gambiae, the observation that the tail in Anopheles RPS6 is nearly twice as long as that in Aedes RPS6 proteins , and evidence that the genus Anopheles is ancestral to Aedes . Because putative homologies to Drosophila histone H1 protein could be recovered as conceptual translation products from the An. gambiae database, we used these sequences to design primers that would discriminate between an An. stephensi histone H1 gene, and the histone H1-like extension in An. stephensi RPS6. Because the Drosophila gene was encoded in a single exon, and the histone message was unlikely to be polyadenylated , we used genomic DNA from An. stephensi as a template for our PCR reaction.
The gene we recovered had more than 90% identity to XP_314184 in An. gambiae. The proteins differed in length by a single amino acid residue, and showed 92 % identity. When we analyzed RPS6 tails and histone H1 genes, we found that the Dipteran histone H1 proteins and the RPS6 tails each fell into distinct groups, suggesting that in present-day mosquitoes, these proteins are evolving independently. Although these data are consistent with the possibility that present-day histone H1 proteins and the histone H1-like tails on mosquito RPS6 protein share a common ancestral gene, the histone tails seem to be evolving independently in the two mosquito genera, and have changed more rapidly than the conventional portion of mosquito RPS6 proteins.
Because RPS6 is considered an important functional component of the ribosome, it seems surprising that a histone H1-like tail occurs at the C-terminal end of this particular protein. However, histone H1-like tails have been reported at the N-terminus of Drosophila melanogaster ribosomal proteins L22 and L23a . The An. gambiae homolog of D. melanogaster L23a also contains an N-terminal histone-like extension. The N-terminal tails of Drosophila L22 and L23a were found in an effort to identify proteins that interact with poly (ADP-ribose) polymerase (PARP). In future studies, we plan to explore whether the histone H1-like tail undergoes posttranslational modification, and whether it plays a functional role in ribosome biogenesis, perhaps through the activity of PARP.
Mosquito cells and culture conditions
We used the ASE-IV Anopheles stephensi mosquito cell line , which was adapted to Eagle's minimal medium, supplemented with non-essential amino acids, glutamine and 5% heat-inactivated fetal bovine serum . This formulation is called E-5 medium.
Genomic DNA preparation
Cells grown as suspended vesicles for 4 to 5 days in twenty 60 mm plates were collected by centrifugation, and the cell pellet was washed twice with phosphate-buffered saline (PBS; ). The cell pellet was resuspended in 20 ml lysis buffer (10 mM Tris-HCl, pH 7.5, 10 mM EDTA, 200 μg/ml proteinase K), and SDS was added to a final concentration of 0.5%. The lysate was incubated at 37°C overnight. NaCl was added to a final concentration of 0.4 M, and the DNA was extracted once with 20 ml phenol, twice with an equal volume of phenol:chloroform (1:1), and twice with an equal volume of chloroform. Two volumes of ethanol were added, and DNA was spooled onto a clean glass rod. The DNA was dried, and dissolved in 10 ml of TE (10 mM Tris-HCl, pH 8.0, containing 1 mM EDTA) at 37°C. RNase A was added to a final concentration of 200 μg/ml and incubated at 37°C for 4 hours. DNA was phenol extracted, ethanol precipitated and dissolved in TE as described above.
DNA amplification by PCR
Genomic DNA (0.4 mg) was digested with Hin dIII (Promega) at 37°C overnight. Enzyme was removed by phenol:choloroform extraction, and the DNA was recovered by precipitation with ethanol and dissolved in TE. Digested DNA (100 ng) was used as template for the PCR reaction, which contained 1X PCR buffer, 1.5 mM MgCl2, 0.2 mM of each of the four dNTPs, 0.4 μM of primer F1 (5'CCG AAG AAG CCG AAG AAG CCC) and R1 (5'TGC TTT CGG CTT CTT GGC AGC) and 2.5 units of Taq DNA polymerase (Promega, Madison, WI). PCR was performed with an initial denaturation at 94°C for 2 minutes. The next 35 cycles included 94°C denaturation for 45 sec, 55°C annealing for 1 minute, and 72°C extension for 1 minute. The reaction was terminated by a final elongation cycle at 72°C for 2 minutes. The PCR product was recovered from a 0.9% agarose gel, purified using Ultra-Clean 15 (MO Bio Laboratories Inc., Solana Beach, CA) and cloned into PGEM T-Easy vector (Promega). The 3'-end of the gene was obtained in a similar manner, using primers R2 (Fig. 2) and F1.
Amplifying the 5'-end of the cDNA
Total RNA was recovered from ASE-IV cells by guanidine isothiocyanate extraction and cesium chloride centrifugation as described by Davis et al. . The final RNA pellet was dissolved in DEPC-treated water and stored at -70°C. RNA (1 μg) was used with the GeneRacer kit (Invitrogen) to obtain the 5' end of the mRNA, using primer R1 as the reverse primer.
Programs and accession numbers
The analysis in Fig. 4A was produced using the Genetics Computer Group (GCG; Madison, WI) program "gap". The tree in Fig. 4B and the alignment in Fig. 5 were produced by an alignment of amino acid residues using default parameters of Clustal X (version 1.83) . The tree was created in PAUP* , with the C. elegans H1 protein designated as an outgroup. The An. stephensi histone H1 sequence has GenBank accession # AY672907.
We acknowledge support from the National Institutes of Health (AI 20385) and from the University of Minnesota Agricultural Experiment Station, St. Paul, MN.
- Sturgill TW, Wu J: Recent progress in characterization of protein kinase cascades for phosphorylation of ribosomal protein S6. Biochim Biophys Acta. 1991, 1092: 350-357. 10.1016/S0167-4889(97)90012-4.PubMedView ArticleGoogle Scholar
- Fumagalli S, Thomas G: S6 phosphorylation and signal transduction. Translational control of gene expression. Edited by: Sonenberg N, Hershey JWB, Mathews MB. 2000, Cold Spring Harbor, NY: Cold Spring Harbor Press, 695-717.Google Scholar
- Hernandez VP, Fallon AM: Ribosomal protein S6 cDNA from two Aedes mosquitoes encodes a carboxyl-terminal extension that resembles histone H1 proteins. Genetica. 1999, 106: 263-267. 10.1023/A:1003914513726.PubMedView ArticleGoogle Scholar
- Hernandez VP, Higgins LA, Schwientek MS, Fallon AM: The histone-like C-terminal extension in ribosomal protein S6 in Aedes and Anopheles mosquitoes is encoded within the distal portion of exon 3. Insect Biochem Mol Biol. 2003, 33: 901-910. 10.1016/S0965-1748(03)00095-X.PubMedView ArticleGoogle Scholar
- Redman KL, Rechsteiner M: Identification of the long ubiquitin extension as ribosomal protein S27a. Nature. 1989, 338: 438-440. 10.1038/338438a0.PubMedView ArticleGoogle Scholar
- Yeates DK, Weigmann BM: Congruence and controversy: Toward a higher level phylogeny of Diptera. Annu Rev Entomol. 1999, 44: 397-428. 10.1146/annurev.ento.44.1.397.PubMedView ArticleGoogle Scholar
- Murphy TJ, Blumenfeld M: Nucleotide sequence of a Drosophila melanogaster H1 histone gene. Nucleic Acids Res. 1986, 14: 5563-5563.PubMedPubMed CentralView ArticleGoogle Scholar
- The FlyBase Consortium: The FlyBase database of the Drosophila genome projects and community literature. Nucleic Acids Res. 2003, 31: 172-175. 10.1093/nar/gkg094. [http://flybase.org/]View ArticleGoogle Scholar
- Trieschmann L, Schulze E, Schulze B, Grossbach U: The histone H1 genes of the dipteran insect, Chironomus thummi, fall under two divergent classes and encode proteins with distinct intranuclear distribution and potentially different functions. Eur J Biochem. 1997, 250: 184-196. 10.1111/j.1432-1033.1997.00184.x.PubMedView ArticleGoogle Scholar
- Schulz E, Trieschmann L, Schulze B, Schmidt ER, Pitzel S, Zechel K, Grossbach U: Structural and functional differences between histone H1 sequence variants with differential intranuclear distribution. Proc Natl Acad Sci USA. 1993, 90: 2481-2485.View ArticleGoogle Scholar
- Mohr E, Trieschmann L, Grossbach U: Histone H1 in two subspecies of Chironomus thummi with different genome sizes: Homologous chromosome sites differ largely in their content of a specific H1 variant. Proc Natl Acad Sci USA. 1989, 86: 9308-9312.PubMedPubMed CentralView ArticleGoogle Scholar
- Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.PubMedPubMed CentralView ArticleGoogle Scholar
- Oosterbroek P, Courtney G: Phylogeny of the nematocerous families of Diptera (Insecta). Zool J Linnean Soc. 1995, 115: 267-311. 10.1006/zjls.1995.0080.View ArticleGoogle Scholar
- Wang Z-F, Sirotkin AM, Buchold GM, Skoultchi AI, Marzluff WF: The mouse histone H1 genes: Gene organization and differential regulation. J Mol Biol. 1997, 271: 124-138. 10.1006/jmbi.1997.1166.PubMedView ArticleGoogle Scholar
- Cirillo L, Zaret K: A linker histone restricts muscle development. Science. 2004, 304: 1607-1609. 10.1126/science.1099577.PubMedView ArticleGoogle Scholar
- Ashburner M: Drosophila: A Laboratory Handbook. 1989, Cold Spring Harbor: Cold Spring Harbor PressGoogle Scholar
- Koyama Y, Katagiri S, Hanai S, Uchida K, Miwa M: Poly(ADP-ribose) polymerase interacts with novel Drosophila ribosomal proteins, L22 and L23a, with unique histone-like amino-terminal extensions. Gene. 1999, 226: 339-345. 10.1016/S0378-1119(98)00529-0.PubMedView ArticleGoogle Scholar
- Kurtti TJ, Munderloh UG: Advances in the definition of culture media for mosquito cells. Invertebrate Cell Systems Applications. Edited by: Mitsuhashi J. 1989, Boca Raton: CRC Press, 1: 21-29.Google Scholar
- Shih KM, Gerenday A, Fallon AM: Culture of mosquito cells in Eagle's medium. In Vitro Cell Develop Biol – Animal. 1998, 34: 629-630.View ArticleGoogle Scholar
- Dulbecco R, Vogt M: Plaque formation and isolation of pure lines with poliomyelitis virus. J Exp Med. 1954, 99: 167-182. 10.1084/jem.99.2.167.PubMedPubMed CentralView ArticleGoogle Scholar
- Davis LG, Dibner MD, Battey JF: Basic Methods in Molecular Biology. 1986, New York: Elsevier, 130-135.View ArticleGoogle Scholar
- Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG: The ClustalX-Windows interface: Flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucl Acids Res. 1997, 25: 4876-4882. 10.1093/nar/25.24.4876.PubMedPubMed CentralView ArticleGoogle Scholar
- Swofford DL: PAUP*: Phylogenetic analysis using parsimony and other methods (software). 2000, Sunderland, MA: Sinauer AssociatesGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.