Codon usage patterns in chromosomal and retrotransposon genes of the mosquito Anopheles gambiae

On the Base Composition of Transposable Elements

Transposable elements exhibit a base composition that is often different from the genomic average and from hosts’ genes. The most common compositional bias is towards Adenosine and Thymine, although this bias is not universal, and elements with drastically different base composition can coexist within the same genome. The AT-richness of transposable elements is apparently maladaptive because it results in poor transcription and sub-optimal translation of proteins encoded by the elements. The cause(s) of this unusual base composition remain unclear and have yet to be investigated. Here, I review what is known about the nucleotide content of transposable elements and how this content can affect the genome of their host as well as their own replication. The compositional bias of transposable elements could result from several non-exclusive processes including horizontal transfer, mutational bias, and selection. It appears that mutation alone cannot explain the high AT-content of transposons and that selection plays a major role in the evolution of the compositional bias. The reason why selection would favor a maladaptive nucleotide content remains however unexplained and is an area of investigation that clearly deserves attention.

Next-generation site-directed transgenesis in the malaria vector mosquito Anopheles gambiae: self-docking strains expressing germline-specific phiC31 integrase

Diseases transmitted by mosquitoes have a devastating impact on global health and the situation is complicated due to difficulties with both existing control measures and the impact of climate change. Genetically modified mosquitoes that are refractory to disease transmission are seen as having great potential in the delivery of novel control strategies. The Streptomyces phage phiC31 integrase system has been successfully adapted for site-directed transgene integration in a range of insects, thus overcoming many limitations due to size constraints and random integration associated with transposon-mediated transformation. Using this technology, we previously published the first site-directed transformation of Anopheles gambiae, the principal vector of human malaria. Mosquitoes were initially engineered to incorporate the phiC31 docking site at a defined genomic location. A second phase of genetic modification then achieved site-directed integration of an anti-malarial effector gene. In the current publication we report improved efficiency and utility of the phiC31 integrase system following the generation of Anopheles gambiae self-docking strains. Four independent strains, with docking sites at known locations on three different chromosome arms, were engineered to express integrase under control of the regulatory regions of the nanos gene from Anopheles gambiae. The resulting protein accumulates in the posterior oocyte to provide integrase activity at the site of germline development. Two self-docking strains, exhibiting significantly different levels of integrase expression, were assessed for site-directed transgene integration and found to demonstrate greatly improved survival and efficiency of transformation. In the fight against malaria, it is imperative to establish a broad repertoire of both anti-malarial effector genes and tissue-specific promoters to regulate their expression, enabling those offering maximum effect with minimum fitness cost to be identified. The improved technology we describe here will facilitate comparative studies of effector transgenes, allowing informed choices to be made that potentially lead to transmission blockade.

Mosquito genomics: progress and challenges

The whole-genome sequencing of mosquitoes has facilitated our understanding of fundamental biological processes at their basic molecular levels and holds potential for application to mosquito control and prevention of mosquito-borne disease transmission. Draft genome sequences are available for Anopheles gambiae, Aedes aegypti, and Culex quinquefasciatus. Collectively, these represent the major vectors of African malaria, dengue fever and yellow fever viruses, and lymphatic filariasis, respectively. Rapid advances in genome technologies have revealed detailed information on genome architecture as well as phenotype-specific transcriptomics and proteomics. These resources allow for detailed comparative analyses within and across populations as well as species. Next-generation sequencing technologies will likely promote a proliferation of genome sequences for additional mosquito species as well as for individual insects. Here we review the current status of genome research in mosquitoes and identify potential areas for further investigations.

Coadaptation of isoacceptor tRNA genes and codon usage bias for translation efficiency in Aedes aegypti and Anopheles gambiae

The transfer RNAs (tRNAs) are essential components of translational machinery. We determined that tRNA isoacceptors (tRNAs with different anticodons but incorporating the same amino acid in protein synthesis) show differential copy number abundance, genomic distribution patterns and sequence evolution between Aedes aegypti and Anopheles gambiae mosquitoes. The tRNA-Ala genes are present in unusually high copy number in the Ae. aegypti genome but not in An. gambiae. Many of the tRNA-Ala genes of Ae. aegypti are flanked by a highly conserved sequence that is not observed in An. gambiae. The relative abundance of tRNA isoacceptor genes is correlated with preferred (or optimal) and nonpreferred (or rare) codons for ∼2-4% of the predicted protein coding genes in both species. The majority (∼74-85%) of these genes are related to pathways involved with translation, energy metabolism and carbohydrate metabolism. Our results suggest that these genes and the related pathways may be under translational selection in these mosquitoes.

Strong heterogeneity in nucleotidic composition and codon bias in the pea aphid (Acyrthosiphon pisum) shown by EST-based coding genome reconstruction

The aim of this study was to analyze patterns of nucleotidic composition and codon usage in the pea aphid genome (Acyrthosiphon pisum). A collection of 60,000 expressed sequence tags (ESTs) in the pea aphid has been used to automatically reconstruct 5809 coding sequences (CDSs), based on similarity with known proteins and on coding style recognition. Reconstructions were manually checked for ribosomal proteins, leading to tentatively reconstruct the nea-complete set of this category. Pea aphid coding sequences showed a shift toward AT (especially at the third codon position) compared to drosophila homologues. Genes with a putative high level of expression (ribosomal and other genes with high EST support) remained more GC3-rich and had a distinct codon usage from bulk sequences: they exhibited a preference for C-ending codons and CGT (for arginine), which thus appeared optimal for translation. However, the discrimination was not as strong as in drosophila, suggesting a reduced degree of translational selection. The space of variation in codon usage for A. pisum appeared to be larger than in drosophila, with a substantial fraction of genes that remained GC3-rich. Some of those (in particular some structural proteins) also showed high levels of codon bias and a very strong preference for C-ending codons, which could be explained either by strong translational selection or by other mechanisms. Finally, genomic traces were analyzed to build 206 fragments containing a full CDS, which allowed studying the correlations between GC contents of coding and those of noncoding (flanking and introns) sequences.

Large-scale gene discovery in the pea aphid Acyrthosiphon pisum (Hemiptera)

Aphids are the leading pests in agricultural crops. A large-scale sequencing of 40,904 ESTs from the pea aphid Acyrthosiphon pisum was carried out to define a catalog of 12,082 unique transcripts. A strong AT bias was found, indicating a compositional shift between Drosophila melanogaster and A. pisum. An in silico profiling analysis characterized 135 transcripts specific to pea-aphid tissues (relating to bacteriocytes and parthenogenetic embryos). This project is the first to address the genetics of the Hemiptera and of a hemimetabolous insect.

Identification and characterization of gp65, a salivary-gland-specific molecule expressed in the malaria vector Anopheles albimanus

A group of salivary-gland-specific proteins, designated gp65, were identified in the mosquito Anopheles albimanus. Two-dimensional gel electrophoresis resolved this group into at least four molecules with pI 6.4-6.5. The N-terminal amino acid sequence was determined for the major species, gp65-1, and degenerate oligonucleotide primers were used to amplify a specific probe for library screening. A 1312 bp cDNA clone encoding a predicted translation product of 386 amino acids was recovered. gp65-1 is expressed abundantly in the medial and distal-lateral lobes of the adult female glands, and is secreted in the saliva. The amino acid sequence has potential sites for N-glycosylation, phosphorylation and myristylation, and is similar to a number of proteins of unknown function from other mosquito species.

Intraspecific DNA variation in nuclear genes of the mosquito Aedes aegypti

Single nucleotide polymorphisms (SNPs) are an abundant source of genetic variation among individual organisms. To assess the usefulness of SNPs for genome analysis in the yellow fever mosquito, Aedes aegypti, we sequenced 25 nuclear genes in each of three strains and analysed nucleotide diversity. The average frequency of nucleotide variation was 12 SNPs per kilobase, indicating that nucleotide variation in Ae. aegypti is similar to that in other organisms, including Drosophila and the malaria vector Anopheles gambiae. Transition polymorphisms outnumbered transversion polymorphisms, at a ratio of about 2:1. We examined codon usage and confirmed that mutational bias favours G and C ending codons. Codon bias was most pronounced in highly expressed genes. Nucleotide diversity estimates indicated that substitution rates are positively correlated in coding and non-coding regions. Nucleotide diversity varied from one gene to another. The unequal distribution of SNPs among Ae. aegypti nuclear genes suggests that single base variations are non-neutral and are subject to selective constraints. Our analysis showed that ubiquitously expressed genes have lower polymorphism rates and are likely under strong purifying selection, whereas tissue specific genes and genes with a putative role in parasite defence exhibit higher levels of polymorphism that may be associated with diversifying selection.

The evolution of HOM-C homeoboxes in the Dipteran family Drosophilidae

Forty-five new Homeotic Complex (HOM-C) homeobox sequences from six species of Drosophilidae (Drosophila heteroneura, D. adiastola, Zaprionus vittiger, Chymomyza amoena, Scaptodrosophila pattersoni and Hirtodrosophila pictiventris) were obtained using a PCR-cloning method. These new homeoboxes are from the labial, proboscipedia, Deformed, Sex combs reduced, fushi tarazu, Antennapedia, Ultrabithorax, abdominal-A and Abdominal-B genes. Phylogenetic signal in the homeobox sequences was assessed and several aspects of sequence evolution were examined. In particular, codon bias was examined and found to exist between the drosophilid species examined here and Anopheles gambiae outgroup sequences. In addition, different patterns of codon bias were detected in homeoboxes interrupted with introns when compared to homeoboxes that are uninterrupted.

A family of pupal-specific cuticular protein genes in the mosquito Anopheles gambiae

We have cloned and sequenced members of a cuticular protein multi-gene family from the mosquito Anopheles gambiae. Three genes (agcp2a-c), each approximately 1 kb in length, were found in a 17.4 kb genomic phage clone. Analysis of ten cDNAs revealed that at least four related genes are present. The open reading frame of the genes and cDNAs showed 95% sequence identity. Divergence was observed in the sequence of the 3' ends and the number of copies of two repeated coding sequences. In situ hybridizations with a probe prepared from one of these circular protein genes physically mapped to two loci, 26B on chromosome 2L and 37A on 3R. Transcription of these An. gambiae cuticular protein genes appears to be limited to pharate pupae and the expressed protein(s) is found in early pupae. The deduced amino acid sequence of these proteins contains a hydrophilic region with significant similarity to other cuticular proteins including the pupal-specific cuticular protein, EDG84, of Drosophila melanogaster (Apple and Fristrom).

Molecular cloning and expression of a hexamerin cDNA from the malaria mosquito, Anopheles gambiae

During the last larval instar, dipteran insects synthesize two hexamerins rich in aromatic residues, typified by the larval serum proteins 1 and 2 (LSP-1 and LSP-2) of Drosophila melanogaster. We report here the characterization of a complete cDNA sequence encoding a LSP-1-like protein from a lower dipteran insect, the malaria mosquito Anopheles gambiae. The cDNA encodes the subunit of a homohexamer, A. gambiae hexamerin-1.1 (AgHex-1.1), which is a major pupal protein but only a minor constituent of late larval hemolymph. AgHex-1.1 is moderately rich in methionine (3.9%) and particularly rich in aromatic residues (21% Phe+Tyr). Cytogenetic analysis reveals AgHex-1.1 to be encoded by a single-copy gene localized to division 22F within the proximal 2La inversion breakpoint of chromosome 2 of A. gambiae. The AgHex-1.1 transcript is first detected in fourth-instar larvae (L4) and disappears abruptly in early pupae. In situ hybridization shows accumulation of the transcript uniquely in the larval fat body. AgHex-1.1 mRNA is re-expressed in male and female adults at about 10% of the L4 level, with no effect of bloodfeeding in females. The potential roles of AgHex-1.1 in Anopheles development and reproductive maturation are discussed.

Analysis of a lysozyme gene from the malaria vector mosquito, Anopheles gambiae

A genomic DNA sequence encoding a basic lysozyme was isolated from the malaria vector mosquito Anopheles gambiae by screening a library with a probe prepared by PCR of reverse transcribed adult RNA. The sequence consists of an upstream region of about 2 kb, a coding region containing three exons and two introns, and a short 3' untranslated region. The coding region indicates that this mosquito lysozyme consists of a signal peptide of 20 residues followed by an 120 aa mature protein which is very similar to other basic lysozymes. The two small introns, 67 and 76 bp, are located at evolutionarily conserved sites. RT-PCR indicated that this gene is expressed abundantly in sugar-fed adults, and at considerably lower levels when females have fed on blood. Although it remains to be seen whether this gene is induced by bacterial infection, the surrounding sequence contains six sequence motifs very similar to the consensus binding sites for a transcription factor similar to NF-kappa B that are found associated with most insect immune response genes. This lysozyme gene maps to division 27 on the left arm of polytene chromosome 2L. An ORF unrelated to any animal protein in current data bases was found at the 5' end of the clone.

Pegasus, a small terminal inverted repeat transposable element found in the white gene of Anopheles gambiae

Pegasus, a novel transposable element, was discovered as a length polymorphism in the white gene of Anopheles gambiae. Sequence analysis revealed that this 535 bp element was flanked by 8 bp target site duplications and 8 bp perfect terminal inverted repeats similar to those found in many members of the Tc1 family. Its small size and lack of long open reading frames preclude protein coding capacity. Southern analysis and in situ hybridization to polytene chromosomes demonstrated that Pegasus occurs in approximately 30 copies in the genomes of An. gambiae and its sibling species and is homogenous in structure but polymorphic in chromosomal location. Characterization of five additional elements by sequencing revealed nucleotide identities of 95% to 99%. Of 30 Pegasus-containing phage clones examined by PCR, only one contained an element exceeding 535 bp in length, due to the insertion of another transposable element-like sequence. Thus, the majority, if not all, extant Pegasus elements may be defective copies of a complete element whose contemporary existence in An. gambiae is uncertain. No Pegasus-hybridizing sequences were detected in nine other anophelines and three culicines examined, suggesting a very limited taxonomic distribution.

Molecular characterization of the Anopheles gambiae 2L telomeric region via an integrated transgene

A Drosophila P-element derivative (pUChsneo) integrated into the telomeric region of the left arm of the second chromosome of Anopheles gambiae was used to clone the proximally flanking An. gambiae sequences. Molecular analyses revealed that the pUChsneo construct was partially duplicated and had integrated into a subterminal minisatellite. This satellite has a repeat unit of 820 bp and is located exclusively at the tip of 2L. No sequence similarity to subterminal minisatellites from other dipterans was detected, but some structural features such as tandem subrepeats are shared. The end of the chromosome was mapped with respect to restriction sites in pUChsneo at approximately generation 100 after the integration event. Considering inevitable terminal nucleotide loss due to incomplete DNA replication, we conclude that the chromosome end must have undergone a dramatic elongation process since it was mapped in generation 23.

Cloning and characterization of the white gene from Anopheles gambiae

A 14 kb region of genomic DNA containing the X-linked Anopheles gambiae eye colour gene, white, was cloned and sequenced. Genomic clones containing distinct white+ alleles were polymorphic for the insertion of a small transposable element in intron 3, and differed at 1% of nucleotide positions compared. Sequence was also determined from a rare 2914 bp cDNA. Comparison of cDNA and genomic sequences established an intron-exon structure distinct from Drosophila white. Despite a common trend in Anopheles and Drosophila of weak codon bias given low levels of gene expression, codon usage by Anopheles gambiae white was strongly biased. Overall amino acid identity between the predicted mosquito and fruitfly proteins was 64%, but dropped to 14% at the amino terminus. To correlate phenotypically white-eyed strains of A. gambiae with structural lesions in white, five available strains were analysed by PCR and Southern blotting. Although these strains carried allelic mutations, independently generated by gamma radiation (three strains) or spontaneous events (two strains), no white lesions were detected. Significantly, another non-allelic X-linked mutation, causing an identical white-eyed phenotype, has been correlated with a structural defect in the cloned white gene (Benedict et al., 1995). Taken together, these observations suggest that the white-eyed mutants analysed in the present study carry mutations in a second eye colour gene and are most likely white+.

A cDNA encoding an ADP/ATP carrier from the mosquito Anopheles gambiae

Two cDNAs are described from Anopheles gambiae that correspond to the ADP/ATP carrier or translocase. The clones are identical except for minor differences in the 5' non-coding region and in the lengths of the poly-A tails. They code for mRNAs of 1261 and 1263 bp and contain one open reading frame of 906 bp. A probe made from the 1263 bp cDNA hybridized to bands of approximately 1260 and 1700 bp on developmental Northern blots. The putative 300 amino acid peptide sequence shows from 53.4-78.5% identity to AAC peptide sequences from a range of organisms from Zea mays to human. Both clones mapped to region 26a on the left arm of chromosome 2 in An. gambiae.

The mitochondrial genome of the mosquito Anopheles gambiae: DNA sequence, genome organization, and comparisons with mitochondrial sequences of other insects

The entire 15,363 bp mitochondrial genome was cloned and sequenced from the mosquito Anopheles gambiae. With respect to the protein-coding genes, rRNA genes and the control region, the gene order was identical to that reported for other insects. There were significant differences, however, in the position and orientation of specific tRNA loci. The overall nucleotide composition was heavily biased towards adenine and thymine, which accounted for 77.6% of all nucleotides. Comparisons were made with the mitochondrial genomes of other insects on the basis genome size and organization, DNA and putative amino acid sequence data, nucleotide substitutions, codon usage and bias, and patterns of AT enrichment.