Distribution of T1, Q, Pegasus and mariner transposable elements on the polytene chromosomes of PEST, a standard strain of Anopheles gambiae

Gypsy, RTE and Mariner transposable elements populate Eyprepocnemis plorans genome

We analyze here the presence and abundance of three types of transposable elements (TEs), i.e. Gypsy, RTE and Mariner, in the genome of the grasshopper Eyprepocnemis plorans. PCR experiments allowed amplification, cloning and sequencing of these elements (EploGypI, EploRTE5, EploMar20) from the E. plorans genome. Fluorescent in situ hybridization (FISH) showed that all three elements are restricted to euchromatic regions, thus being absent from the pericentromeric region of all A chromosomes, which contain a satellite DNA (satDNA) and ribosomal DNA (rDNA), and being very scarce in B chromosomes mostly made up of these two types of repetitive DNA. FISH suggested that EploGypI is the most abundant and EploMar20 is the least abundant, with EploRTE5 showing intermediate abundance. An estimation of copy number, by means of quantitative PCR, showed that EploGypI is, by far, the most abundant element, followed by EploRTE5 and EploMar20, in consistency with FISH results. RNA isolation and PCR experiments on complementary DNA (cDNA) showed the presence of transcripts for the three TE elements. The implications of the preferential location of these TE elements into euchromatin, the significance of TE abundance in the giant genome of this species, and a possible relationship between TEs and B chromosome mutability, are discussed.

Breakpoint structure of the Anopheles gambiae 2Rb chromosomal inversion

Background

Alternative arrangements of chromosome 2 inversions in Anopheles gambiae are important sources of population structure, and are associated with adaptation to environmental heterogeneity. The forces responsible for their origin and maintenance are incompletely understood. Molecular characterization of inversion breakpoints provides insight into how they arose, and provides the basis for development of molecular karyotyping methods useful in future studies.

Methods

Sequence comparison of regions near the cytological breakpoints of 2Rb allowed the molecular delineation of breakpoint boundaries. Comparisons were made between the standard 2R+^b arrangement in the An. gambiae PEST reference genome and the inverted 2Rb arrangements in the An. gambiae M and S genome assemblies. Sequence differences between alternative 2Rb arrangements were exploited in the design of a PCR diagnostic assay, which was evaluated against the known chromosomal banding pattern of laboratory colonies and field-collected samples from Mali and Cameroon.

Results

The breakpoints of the 7.55 Mb 2Rb inversion are flanked by extensive runs of the same short (72 bp) tandemly organized sequence, which was likely responsible for chromosomal breakage and rearrangement. Application of the molecular diagnostic assay suggested that 2Rb has a single common origin in An. gambiae and its sibling species, Anopheles arabiensis, and also that the standard arrangement (2R+^b) may have arisen twice through breakpoint reuse. The molecular diagnostic was reliable when applied to laboratory colonies, but its accuracy was lower in natural populations.

Conclusions

The complex repetitive sequence flanking the 2Rb breakpoint region may be prone to structural and sequence-level instability. The 2Rb molecular diagnostic has immediate application in studies based on laboratory colonies, but its usefulness in natural populations awaits development of complementary molecular tools.

Chromosomal plasticity and evolutionary potential in the malaria vector Anopheles gambiae sensu stricto: insights from three decades of rare paracentric inversions

Background

In the Anopheles gambiae complex, paracentric chromosomal inversions are non-randomly distributed along the complement: 18/31 (58%) of common polymorphic inversions are on chromosome arm 2R, which represents only ~30% of the complement. Moreover, in An. gambiae sensu stricto, 6/7 common polymorphic inversions occur on 2R. Most of these inversions are considered markers of ecological adaptation that increase the fitness of the carriers of alternative karyotypes in contrasting habitats. However, little is known about the evolutionary forces responsible for their origin and subsequent establishment in field populations.

Results

Here, we present data on 82 previously undescribed rare chromosomal inversions (RCIs) recorded during extensive field sampling in 16 African countries over a 30 year period, which may shed light on the dynamics of chromosomal plasticity in An. gambiae. We analyzed breakpoint distribution, length, and geographic distribution of RCIs, and compared these measures to those of the common inversions. We found that RCIs, like common inversions, are disproportionately clustered on 2R, which may indicate that this arm is especially prone to breakages. However, contrasting patterns were observed between the geographic distribution of common inversions and RCIs. RCIs were equally frequent across biomes and on both sides of the Great Rift Valley (GRV), whereas common inversions predominated in arid ecological settings and west of the GRV. Moreover, the distribution of RCI lengths followed a random pattern while common inversions were significantly less frequent at shorter lengths.

Conclusion

Because 17/82 (21%) RCIs were found repeatedly at very low frequencies – at the same sampling location in different years and/or in different sampling locations – we suggest that RCIs are subject mainly to drift under unperturbed ecological conditions. Nevertheless, RCIs may represent an important reservoir of genetic variation for An. gambiae in response to environmental changes, further testifying to the considerable evolutionary potential hidden within this pan-African malaria vector.

SNP discovery and molecular evolution in Anopheles gambiae, with special emphasis on innate immune system

Background

Anopheles innate immunity affects Plasmodium development and is a potential target of innovative malaria control strategies. The extent and distribution of nucleotide diversity in immunity genes might provide insights into the evolutionary forces that condition pathogen-vector interactions. The discovery of polymorphisms is an essential step towards association studies of susceptibility to infection.

Results

We sequenced coding fragments of 72 immune related genes in natural populations of Anopheles gambiae and of 37 randomly chosen genes to provide a background measure of genetic diversity across the genome. Mean nucleotide diversity (π) was 0.0092 in the A. gambiae S form, 0.0076 in the M form and 0.0064 in A. arabiensis. Within each species, no statistically significant differences in mean nucleotide diversity were detected between immune related and non immune related genes. Strong purifying selection was detected in genes of both categories, presumably reflecting strong functional constraints.

Conclusion

Our results suggest similar patterns and rates of molecular evolution in immune and non-immune genes in A. gambiae. The 3,214 Single Nucleotide Polymorphisms (SNPs) that we identified are the first large set of Anopheles SNPs from fresh, field-collected material and are relevant markers for future phenotype-association studies.

Segmental duplication implicated in the genesis of inversion 2Rj of Anopheles gambiae

The malaria vector Anopheles gambiae maintains high levels of inversion polymorphism that facilitate its exploitation of diverse ecological settings across tropical Africa. Molecular characterization of inversion breakpoints is a first step toward understanding the processes that generate and maintain inversions. Here we focused on inversion 2Rj because of its association with the assortatively mating Bamako chromosomal form of An. gambiae, whose distinctive breeding sites are rock pools beside the Niger River in Mali and Guinea. Sequence and computational analysis of 2Rj revealed the same 14.6 kb insertion between both breakpoints, which occurred near but not within predicted genes. Each insertion consists of 5.3 kb terminal inverted repeat arms separated by a 4 kb spacer. The insertions lack coding capacity, and are comprised of degraded remnants of repetitive sequences including class I and II transposable elements. Because of their large size and patchwork composition, and as no other instances of these insertions were identified in the An. gambiae genome, they do not appear to be transposable elements. The 14.6 kb modules inserted at both 2Rj breakpoint junctions represent low copy repeats (LCRs, also called segmental duplications) that are strongly implicated in the recent (∼0.4N_e generations) origin of 2Rj. The LCRs contribute to further genome instability, as demonstrated by an imprecise excision event at the proximal breakpoint of 2Rj in field isolates.

Comparative analysis of BAC and whole genome shotgun sequences from an Anopheles gambiae region related to Plasmodium encapsulation

The only natural mechanism of malaria transmission in sub-Saharan Africa is the mosquito, generally Anopheles gambiae. Blocking malaria parasite transmission by stopping the development of Plasmodium in the insect vector would provide a useful alternative to the current methods of malaria control. Toward this end, it is important to understand the molecular basis of the malaria parasite refractory phenotype in An. gambiae mosquito strains. We have selected and sequenced six bacterial artificial chromosome (BAC) clones from the Pen-1 region that is the major quantitative trait locus involved in Plasmodium encapsulation. The sequence and the annotation of five overlapping BAC clones plus one adjacent, but not contiguous clone, totaling 585kb of genomic sequence from the centromeric end of the Pen-1 region of the PEST strain were compared to that of the genome sequence of the same strain produced by the whole genome shotgun technique. This project identified 23 putative mosquito genes plus putative copies of the retrotransposable elements BEL12 and TRANSIBN1_AG in the six BAC clones. Nineteen of the predicted genes are most similar to their Drosophila melanogaster homologs while one is more closely related to vertebrate genes. Comparison of these new BAC sequences plus previously published BAC sequences to the cognate region of the assembled genome sequence identified three retrotransposons present in one sequence version but not the other. One of these elements, Indy, has not been previously described. These observations provide evidence for the recent active transposition of these elements and demonstrate the plasticity of the Anopheles genome. The BAC sequences strongly support the public whole genome shotgun assembly and automatic annotation while also demonstrating the benefit of complementary genome sequences and of human curation. Importantly, the data demonstrate the differences in the genome sequence of an individual mosquito compared to that of a hypothetical, average genome sequence generated by whole genome shotgun assembly.

Transposable elements in mosquitoes

We describe the current state of knowledge about transposable elements (TEs) in different mosquito species. DNA-based elements (class II elements), non-LTR retrotransposons (class I elements), and MITEs (Miniature Inverted Repeat Transposable Elements) are found in the three genera, Anopheles, Aedes and Culex, whereas LTR retrotransposons (class I elements) are found only in Anopheles and Aedes. Mosquitoes were the first insects in which MITEs were reported; they have several LTR retrotransposons belonging to the Pao family, which is distinct from the Gypsy-Ty3 and Copia-Ty1 families. The number of TE copies shows huge variations between classes of TEs within a given species (from 1 to 1000), in sharp contrast to Drosophila, which shows only relatively minor differences in copy number between elements (from 1 to 100). The genomes of these insects therefore display major differences in the amount of TEs and therefore in their structure and global composition. We emphasize the need for more population genetic data about the activity of TEs, their distribution over chromosomes and their frequencies in natural populations of mosquitoes, to further the current attempts to develop a transgenic mosquito unable to transmit malaria that is intended to replace the natural populations.

Evolution of the mdg1 lineage of the Ty3/gypsy group of LTR retrotransposons in Anopheles gambiae

So far, only a few retrovirus-like transposable elements (TEs) have been reported in Anopheles mosquitoes, although a large fraction of their genomes is made up of these middle repetitive sequences. By screening the A. gambiae genome databases, we have found 10 element families belonging to the mdg1 lineage of the Ty3/gypsy group of long terminal repeat (LTR) retrotransposons. These Anopheles families constitute a sister clade of the Drosophila representatives of this same lineage. According to the phylogenetic reconstruction of their open reading frame (ORF)2 enzymatic domains, the analysis of patterns of nucleotide substitution therein, and the estimation of the age of particular insertions, all these elements must have been active until quite recently, and some of them must be very young. On the other hand, the fact that all these element families are primarily composed of fragmentary copies (mostly solos) or full-length copies with inactivating mutations indicates that their turnover rate has been probably very low. Finally, incongruent phylogenies obtained from different regions of the elements strongly suggest that recombination has played a significant role in their evolutionary history.

Heterologous expression of four glutathione transferase genes genetically linked to a major insecticide-resistance locus from the malaria vector Anopheles gambiae

A cluster of eight genes encoding glutathione transferases (GSTs) are located on division 33B of polytene chromosome arm 3R of the African malaria mosquito, Anopheles gambiae. This region of the genome contains a major 1,1,1-trichloro-2,2-bis-( p -chlorophenyl)ethane (DDT)-resistance locus, rtd1. These GSTs belong to the insect-specific Epsilon class and share between 22.6 and 65.2% identity at the amino acid level. Two distinct allelic variants of the Epsilon GST, GSTe1, differing at 12 out of 224 amino acid residues, are present in laboratory and field populations of A. gambiae. To investigate the possible role of these GSTs in conferring resistance to the insecticide DDT, both GSTe1 alleles, plus three additional members of this gene cluster, were expressed in Escherichia coli and the recombinant proteins biochemically characterized. The five putative glutathione transferases encoded catalytically active subunits with variable biochemical properties. For example, the two allelic variants of GSTE1-1 encoded proteins with over 100-fold variation in peroxidase activity, while the three remaining GSTs had no detectable peroxidase activity. Only GSTE2-2 was able to metabolize DDT. Western blots using antibodies raised against these GSTs indicated that the expression of GSTE2-2 is elevated in a DDT-resistant strain of A. gambiae.

Construction of a BAC library and generation of BAC end sequence-tagged connectors for genome sequencing of the African malaria mosquito Anopheles gambiae

A Bacterial Artificial Chromosome (BAC) genomic DNA library of Anopheles gambiae, the major human malaria vector in sub-Saharan Africa, was constructed and characterized. This library (ND-TAM) is composed of 30,720 BAC clones in eighty 384-well plates. The estimated average insert size of the library is 133 kb, with an overall genome coverage of approximately 14-fold. The ends of approximately two-thirds of the clones in the library were sequenced, yielding 32,340 pair-mate ends. A statistical analysis (G-test) of the results of PCR screening of the library indicated a random distribution of BACs in the genome, although one gap encompassing the white locus on the X-chromosome was identified. Furthermore, combined with another previously constructed BAC library (ND-1), ~2,000 BACs have been physically mapped by polytene chromosomal in situ hybridization. These BAC end pair mates and physically mapped BACs have been useful for both the assembly of a fully sequenced A. gambiae genome and for linking the assembled sequence to the three polytene chromosomes. This ND-TAM library is now publicly available at both http://www.malaria.mr4.org/mr4pages/index.html/ and http://hbz.tamu.edu/, providing a valuable resource to the mosquito research community.

The esterase and PHD domains in CR1-like non-LTR retrotransposons

Most active non-LTR (long terminal repeat) retrotransposons carry two open reading frames (ORFs) encoding ORF1p and ORF2p proteins. The ORF2p proteins are relatively well studied and are known to contain endonuclease/reverse transcriptase domains. At the same time, the biological function of ORF1p proteins remains poorly understood, except in that they nonspecifically bind single-stranded mRNA/DNA molecules. CR1-like elements form the most widely distributed clade/superfamily of non-LTR retrotransposons. We found that ORF1p proteins encoded by diverse CR1-like elements contain conserved esterase domain (ES) or plant homeodomain (PHD). This indicates that CR1-like ORF1p proteins are either lipolytic enzymes or are involved in protein-protein interactions related to chromatin remodeling. Sequence conservation of ES suggests that interaction with cellular membranes is an important phase in life circles of CR1-like elements. Presumably such interaction helps in penetrating host cells. As a consequence, the presence of multiple young CR1 families characterized by approximately 10% intrafamily and 40% interfamily identities may be explained by a relatively frequent horizontal transfer of these CR1-like elements. Unexpectedly, ES links together non-LTR retrotransposons and single-stranded RNA viruses like influenza C and coronaviruses, which are known to depend on their own ES.

The genome sequence of the malaria mosquito Anopheles gambiae

Anopheles gambiae is the principal vector of malaria, a disease that afflicts more than 500 million people and causes more than 1 million deaths each year. Tenfold shotgun sequence coverage was obtained from the PEST strain of A. gambiae and assembled into scaffolds that span 278 million base pairs. A total of 91% of the genome was organized in 303 scaffolds; the largest scaffold was 23.1 million base pairs. There was substantial genetic variation within this strain, and the apparent existence of two haplotypes of approximately equal frequency ("dual haplotypes") in a substantial fraction of the genome likely reflects the outbred nature of the PEST strain. The sequence produced a conservative inference of more than 400,000 single-nucleotide polymorphisms that showed a markedly bimodal density distribution. Analysis of the genome sequence revealed strong evidence for about 14,000 protein-encoding transcripts. Prominent expansions in specific families of proteins likely involved in cell adhesion and immunity were noted. An expressed sequence tag analysis of genes regulated by blood feeding provided insights into the physiological adaptations of a hematophagous insect.

Comparative genomic analysis in the region of a major Plasmodium-refractoriness locus of Anopheles gambiae

We have sequenced six overlapping clones from a library of bacterial artificial chromosome (BAC) clones derived from a laboratory strain of the mosquito, Anopheles gambiae, the major vector of human malaria in Africa. The resulting uninterrupted 528-kb sequence is from the 8C region of the mosquito 2R chromosome, at or very near the major refractoriness locus associated with melanotic encapsulation of parasites. This sequence represents the first extensive view of the mosquito genome structure encompassing 48 genes. Genomic comparison reveals that the majority of the orthologues are found in six microsyntenic clusters in Drosophila melanogaster. A BAC clone that is wholly contained within this region demonstrates the existence of a remarkable degree of local polymorphism in this species, which may prove important for its population structure and vectorial capacity.

Structure and evolution of mtanga, a retrotransposon actively expressed on the Y chromosome of the African malaria vector Anopheles gambiae

Here we report the discovery of a novel family of long terminal repeat (LTR)-retrotransposons designated mtanga-Y, specific to the Y chromosome of the African malaria vector, Anopheles gambiae. mtanga-Y elements represent the first Y-linked sequences and the first members of the Ty1-copia superfamily of retrotransposons described from this mosquito. Analysis of a full-length 4,284-bp element revealed the presence of two intact overlapping open reading frames bounded by LTRs of 119 bp. Evidence suggests that the elements are capable of retrotransposition, as transcripts and potential replication intermediates (one-LTR circles) were detected. However, the approximately 12 copies of mtanga-Y appear to be clustered rather than dispersed on the Y chromosome. Absent from the Y chromosome of four sibling species (A. arabiensis, A. quadriannulatus, A. melas, and A. merus), similar, but often defective, mtanga elements are present elsewhere in these genomes, as well as in A. gambiae. These data are consistent with a relatively recent invasion of the A. gambiae Y chromosome by an intact element. The presence of functional mtanga-Y elements suggests that the Y chromosome may be a source, not just a sink, for retrotransposons.

Plasmodium, human and Anopheles genomics and malaria

The Plasmodium spp. parasites that cause malaria are transmitted to humans by Anopheles spp. mosquitoes. Scientists have now amassed a great body of knowledge about the parasite, its mosquito vector and human host. Yet this year there will be 300-500 million new malaria infections and 1-3 million deaths caused by the disease. We believe that integrated analyses of genome sequence, DNA polymorphisms, and messenger RNA and protein expression profiles will lead to greater understanding of the molecular basis of vector-human and host-parasite interactions and provide strategies to build upon these insights to develop interventions to mitigate human morbidity and mortality from malaria.

Identification of a novel class of insect glutathione S-transferases involved in resistance to DDT in the malaria vector Anopheles gambiae

The sequence and cytological location of five Anopheles gambiae glutathione S-transferase (GST) genes are described. Three of these genes, aggst1-8, aggst1-9 and aggst1-10, belong to the insect class I family and are located on chromosome 2R, in close proximity to previously described members of this gene family. The remaining two genes, aggst3-1 and aggst3-2, have a low sequence similarity to either of the two previously recognized classes of insect GSTs and this prompted a re-evaluation of the classification of insect GST enzymes. We provide evidence for seven possible classes of insect protein with GST-like subunits. Four of these contain sequences with significant similarities to mammalian GSTs. The largest novel insect GST class, class III, contains functional GST enzymes including two of the A. gambiae GSTs described in this report and GSTs from Drosophila melanogaster, Musca domestica, Manduca sexta and Plutella xylostella. The genes encoding the class III GST of A. gambiae map to a region of the genome on chromosome 3R that contains a major DDT [1,1,1-trichloro-2,2-bis-(p-chlorophenyl)ethane] resistance gene, suggesting that this gene family is involved in GST-based resistance in this important malaria vector. In further support of their role in resistance, we show that the mRNA levels of aggst3-2 are approx. 5-fold higher in a DDT resistant strain than in the susceptible strain and demonstrate that recombinant AgGST3-2 has very high DDT dehydrochlorinase activity.

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.