Gambiae Genome Research Articles

In-depth transcriptomic analysis of Anopheles gambiae hemocytes uncovers novel genes and the oenocytoid developmental lineage

BackgroundHemocytes are immune cells that patrol the mosquito hemocoel and mediate critical cellular defense responses against pathogens. However, despite their importance, a comprehensive transcriptome of these cells was lacking because they constitute a very small fraction of the total cells in the insect, limiting the study of hemocyte differentiation and immune function.ResultsIn this study, an in-depth hemocyte transcriptome was built by extensive bulk RNA sequencing and assembly of hemocyte RNAs from adult A. gambiae female mosquitoes, based on approximately 2.4 billion short Illumina and about 9.4 million long PacBio high-quality reads that mapped to the A. gambiae PEST genome (P4.14 version). A total of 34,939 transcripts were annotated including 4,020 transcripts from novel genes and 20,008 novel isoforms that result from extensive differential splicing of transcripts from previously annotated genes. Most hemocyte transcripts identified (89.8%) are protein-coding while 10.2% are non-coding RNAs. The number of transcripts identified in the novel hemocyte transcriptome is twice the number in the current annotation of the A. gambiae genome (P4.14 version). Furthermore, we were able to refine the analysis of a previously published single-cell transcriptome (scRNAseq) data set by using the novel hemocyte transcriptome as a reference to re-define the hemocyte clusters and determine the path of hemocyte differentiation. Unsupervised pseudo-temporal ordering using the Tools for Single Cell Analysis software uncovered a novel putative prohemocyte precursor cell type that gives rise to prohemocytes. Pseudo-temporal ordering with the Monocle 3 software, which analyses changes in gene expression during dynamic biological processes, determined that oenocytoids derive from prohemocytes, a cell population that also gives rise to the granulocyte lineage.ConclusionA high number of mRNA splice variants are expressed in hemocytes, and they may account for the plasticity required to mount efficient responses to many different pathogens. This study highlights the importance of a comprehensive set of reference transcripts to perform robust single-cell transcriptomic data analysis of cells present in low abundance. The detailed annotation of the hemocyte transcriptome will uncover new facets of hemocyte development and function in adult dipterans and is a valuable community resource for future studies on mosquito cellular immunity.

Gene drive designs for efficient and localisable population suppression using Y-linked editors.

The sterile insect technique (SIT) has been successful in controlling some pest species but is not practicable for many others due to the large number of individuals that need to be reared and released. Previous computer modelling has demonstrated that the release of males carrying a Y-linked editor that kills or sterilises female descendants could be orders of magnitude more efficient than SIT while still remaining spatially restricted, particularly if combined with an autosomal sex distorter. In principle, further gains in efficiency could be achieved by using a self-propagating double drive design, in which each of the two components (the Y-linked editor and the sex ratio distorter) boosted the transmission of the other. To better understand the expected dynamics and impact of releasing constructs of this new design we have analysed a deterministic population genetic and population dynamic model. Our modelling demonstrates that this design can suppress a population from very low release rates, with no invasion threshold. Importantly, the design can work even if homing rates are low and sex chromosomes are silenced at meiosis, potentially expanding the range of species amenable to such control. Moreover, the predicted dynamics and impacts can be exquisitely sensitive to relatively small (e.g., 25%) changes in allele frequencies in the target population, which could be exploited for sequence-based population targeting. Analysis of published Anopheles gambiae genome sequences indicates that even for weakly differentiated populations with an FST of 0.02 there may be thousands of suitably differentiated genomic sites that could be used to restrict the spread and impact of a release. Our proposed design, which extends an already promising development pathway based on Y-linked editors, is therefore a potentially useful addition to the menu of options for genetic biocontrol.

Annotations of novel antennae-expressed genes in male Glossina morsitans morsitans tsetse flies.

Tsetse flies use antennal expressed genes to navigate their environment. While most canonical genes associated with chemoreception are annotated, potential gaps with important antennal genes are uncharacterized in Glossina morsitans morsitans. We generated antennae-specific transcriptomes from adult male G. m. morsitans flies fed/unfed on bloodmeal and/or exposed to an attractant (ε-nonalactone), a repellant (δ-nonalactone) or paraffin diluent. Using bioinformatics approach, we mapped raw reads onto G. m. morsitans gene-set from VectorBase and collected un-mapped reads (constituting the gaps in annotation). We de novo assembled these reads (un-mapped) into transcript and identified corresponding genes of the transcripts in G. m. morsitans gene-set and protein homologs in UniProt protein database to further annotate the gaps. We predicted potential protein-coding gene regions associated with these transcripts in G. m. morsitans genome, annotated/curated these genes and identified their putative annotated orthologs/homologs in Drosophila melanogaster, Musca domestica or Anopheles gambiae genomes. We finally evaluated differential expression of the novel genes in relation to odor exposures relative to no-odor control (unfed flies). About 45.21% of the sequenced reads had no corresponding transcripts within G. m. morsitans gene-set, corresponding to the gap in existing annotation of the tsetse fly genome. The total reads assembled into 72,428 unique transcripts, most (74.43%) of which had no corresponding genes in the UniProt database. We annotated/curated 592 genes from these transcripts, among which 202 were novel while 390 were improvements of existing genes in the G. m. morsitans genome. Among the novel genes, 94 had orthologs in D. melanogaster, M. domestica or An. gambiae while 88 had homologs in UniProt. These orthologs were putatively associated with oxidative regulation, protein synthesis, transcriptional and/or translational regulation, detoxification and metal ion binding, thus providing insight into their specific roles in antennal physiological processes in male G. m. morsitans. A novel gene (GMOY014237.R1396) was differentially expressed in response to the attractant. We thus established significant gaps in G. m. morsitans genome annotation and identified novel male antennae-expressed genes in the genome, among which > 53% (108) are potentially G. m. morsitans specific.

Establishment and a comparative transcriptomic analysis of a male-specific cell line from the African malaria mosquito Anopheles gambiae

Cell lines allow studying various biological processes that may not be easily tractable in whole organisms. Here, we have established the first male-specific cell line from the African malaria mosquito, Anopheles gambiae. The cells, named AgMM and derived from the sex-sorted neonate larvae, were able to undergo spontaneous contractions for a number of passages following establishment, indicating their myoblast origin. Comparison of their transcriptome to the transcriptome of an A. gambiae-derived Sua5.1 hemocyte cells revealed distinguishing molecular signatures of each cell line, including numerous muscle-related genes that were highly and uniquely expressed in the AgMM cells. Moreover, the AgMM cells express the primary sex determiner gene Yob and support male sex determination and dosage compensation pathways. Therefore, the AgMM cell line represents a valuable tool for molecular and biochemical in vitro studies of these male-specific processes. In a broader context, a rich transcriptomic data set generated in this study contributes to a better understanding of transcribed regions of the A. gambiae genome and sheds light on the biology of both cell types, facilitating their anticipated use for various cell-based assays.

Identification of new Anopheles gambiae transcriptional enhancers using a cross-species prediction approach.

The success of transgenic mosquito vector control approaches relies on well-targeted gene expression, requiring the identification and characterization of a diverse set of mosquito promoters and transcriptional enhancers. However, few enhancers have been characterized in Anopheles gambiae to date. Here, we employ the SCRMshaw method we previously developed to predict enhancers in the A. gambiae genome, preferentially targeting vector-relevant tissues such as the salivary glands, midgut and nervous system. We demonstrate a high overall success rate, with at least 8 of 11 (73%) tested sequences validating as enhancers in an in vivo xenotransgenic assay. Four tested sequences drive expression in either the salivary gland or the midgut, making them directly useful for probing the biology of these infection-relevant tissues. The success of our study suggests that computational enhancer prediction should serve as an effective means for identifying A. gambiae enhancers with activity in tissues involved in malaria propagation and transmission.

Physical Mapping of the Anopheles (Nyssorhynchus) darlingi Genomic Scaffolds.

Simple SummaryAnopheles darlingi mosquitoes are the main vectors of malaria in the Brazilian Amazon. To assign genomic DNA sequences to chromosomes of this species, we performed fluorescence in situ hybridization of DNA probes with salivary glands polytene chromosomes. We compared the physical locations of the An. darlingi probes with homologous sequences in other Anopheles species, namely Anopheles albimanus and Anopheles gambiae. The results demonstrated that substantial genome rearrangements occurred throughout the evolutionary history of these mosquitoes. The physical mapping data can be useful for improving the structural accuracy of the An. darlingi genome assembly and for understanding the chromosomal evolution of these mosquitoes.The genome assembly of Anopheles darlingi consists of 2221 scaffolds (N50 = 115,072 bp) and has a size spanning 136.94 Mbp. This assembly represents one of the smallest genomes among Anopheles species. Anopheles darlingi genomic DNA fragments of ~37 Kb were cloned, end-sequenced, and used as probes for fluorescence in situ hybridization (FISH) with salivary gland polytene chromosomes. In total, we mapped nine DNA probes to scaffolds and autosomal arms. Comparative analysis of the An. darlingi scaffolds with homologous sequences of the Anopheles albimanus and Anopheles gambiae genomes identified chromosomal rearrangements among these species. Our results confirmed that physical mapping is a useful tool for anchoring genome assemblies to mosquito chromosomes.

Anopheles gambiae Genome Conservation as a Resource for Rational Gene Drive Target Site Selection.

Simple SummaryMalaria is a huge public health burden that affects predominantly sub-Saharan Africa and is transmitted by Anopheles mosquitoes. As a measure for population control, a method called gene drive has been recently developed, which relies on genetic engineering to introduce specific genetic traits into mosquito populations. Gene drives are designed to insert at specific target sites in the mosquito genome. The efficacy of gene drives greatly depends on the selection of appropriate target sites that are functionally or structurally constrained and less likely to tolerate mutations that can hinder the spread of the desired trait in the population. The aim of this study was to perform a genome-wide analysis of highly conserved genomic regions in Anopheles gambiae and introduce a measure of conservation that could indicate sites of functional or structural constraint. The results of this analysis are gathered in a publicly available dataset that can support gene drive target selection and can offer further insights in the nature of conserved genomic regions.The increase in molecular tools for the genetic engineering of insect pests and disease vectors, such as Anopheles mosquitoes that transmit malaria, has led to an unprecedented investigation of the genomic landscape of these organisms. The understanding of genome variability in wild mosquito populations is of primary importance for vector control strategies. This is particularly the case for gene drive systems, which look to introduce genetic traits into a population by targeting specific genomic regions. Gene drive targets with functional or structural constraints are highly desirable as they are less likely to tolerate mutations that prevent targeting by the gene drive and consequent failure of the technology. In this study we describe a bioinformatic pipeline that allows the analysis of whole genome data for the identification of highly conserved regions that can point at potential functional or structural constraints. The analysis was conducted across the genomes of 22 insect species separated by more than hundred million years of evolution and includes the observed genomic variation within field caught samples of Anopheles gambiae and Anopheles coluzzii, the two most dominant malaria vectors. This study offers insight into the level of conservation at a genome-wide scale as well as at per base-pair resolution. The results of this analysis are gathered in a data storage system that allows for flexible extraction and bioinformatic manipulation. Furthermore, it represents a valuable resource that could provide insight into population structure and dynamics of the species in the complex and benefit the development and implementation of genetic strategies to tackle malaria.

In Silico Karyotyping of Chromosomally Polymorphic Malaria Mosquitoes in the Anopheles gambiae Complex.

Chromosomal inversion polymorphisms play an important role in adaptation to environmental heterogeneities. For mosquito species in the Anopheles gambiae complex that are significant vectors of human malaria, paracentric inversion polymorphisms are abundant and are associated with ecologically and epidemiologically important phenotypes. Improved understanding of these traits relies on determining mosquito karyotype, which currently depends upon laborious cytogenetic methods whose application is limited both by the requirement for specialized expertise and for properly preserved adult females at specific gonotrophic stages. To overcome this limitation, we developed sets of tag single nucleotide polymorphisms (SNPs) inside inversions whose biallelic genotype is strongly correlated with inversion genotype. We leveraged 1,347 fully sequenced An. gambiae and Anopheles coluzzii genomes in the Ag1000G database of natural variation. Beginning with principal components analysis (PCA) of population samples, applied to windows of the genome containing individual chromosomal rearrangements, we classified samples into three inversion genotypes, distinguishing homozygous inverted and homozygous uninverted groups by inclusion of the small subset of specimens in Ag1000G that are associated with cytogenetic metadata. We then assessed the correlation between candidate tag SNP genotypes and PCA-based inversion genotypes in our training sets, selecting those candidates with >80% agreement. Our initial tests both in held-back validation samples from Ag1000G and in data independent of Ag1000G suggest that when used for in silico inversion genotyping of sequenced mosquitoes, these tags perform better than traditional cytogenetics, even for specimens where only a small subset of the tag SNPs can be successfully ascertained.

Introgression of a synthetic sex ratio distortion system from Anopheles gambiae into Anopheles arabiensis

I-PpoI is a homing endonuclease that has a high cleavage activity and specificity for a conserved sequence within the ribosomal rDNA repeats, located in a single cluster on the Anopheles gambiae X chromosome. This property has been exploited to develop a synthetic sex ratio distortion system in this mosquito species. When I-PpoI is expressed from a transgene during spermatogenesis in mosquitoes, the paternal X chromosome is shredded and only Y chromosome-bearing sperm are viable, resulting in a male-biased sex ratio of >95% in the progeny. These distorter male mosquitoes can efficiently suppress caged wild-type populations, providing a powerful tool for vector control strategies. Given that malaria mosquito vectors belong to a species complex comprising at least two major vectors, we investigated whether the sex distorter I-PpoI, originally integrated in the A. gambiae genome, could be transferred via introgression to the sibling vector species Anopheles arabiensis. In compliance with Haldane’s rule, F1 hybrid male sterility is known to occur in all intercrosses among members of the Anopheles gambiae complex. A scheme based on genetic crosses and transgene selection was used to bypass F1 hybrid male sterility and introgress the sex distorter I-PpoI into the A. arabiensis genetic background. Our data suggest that this sex distortion technique can be successfully applied to target A. arabiensis mosquitoes.

UDP-glycosyltransferase genes and their association and mutations associated with pyrethroid resistance in Anopheles sinensis (Diptera: Culicidae)

BackgroundUDP-glycosyltransferase (UGT) is an important biotransformation superfamily of enzymes. They catalyze the transfer of glycosyl residues from activated nucleotide sugars to acceptor hydrophobic molecules, and function in several physiological processes, including detoxification, olfaction, cuticle formation, pigmentation. The diversity, classification, scaffold location, characteristics, phylogenetics, and evolution of the superfamily of genes at whole genome level, and their association and mutations associated with pyrethroid resistance are still little known.MethodsThe present study identified UGT genes in Anopheles sinensis genome, classified UGT genes in An. sinensis, Anopheles gambiae, Aedes aegypti and Drosophila melanogaster genomes, and analysed the scaffold location, characteristics, phylogenetics, and evolution of An. sinensis UGT genes using bioinformatics methods. The present study also identified the UGTs associated with pyrethroid resistance using three field pyrethroid-resistant populations with RNA-seq and RT-qPCR, and the mutations associated with pyrethroid resistance with genome re-sequencing in An. sinensis.ResultsThere are 30 putative UGTs in An. sinensis genome, which are classified into 12 families (UGT301, UGT302, UGT306, UGT308, UGT309, UGT310, UGT313, UGT314, UGT315, UGT36, UGT49, UGT50) and further into 23 sub-families. The UGT308 is significantly expanded in gene number compared with other families. A total of 119 UGTs from An. sinensis, An. gambiae, Aedes aegypti and Drosophila melanogaster genomes are classified into 19 families, of which seven are specific for three mosquito species and seven are specific for Drosophila melanogaster. The UGT308 and UGT302 are proposed to main families involved in pyrethroid resistance. The AsUGT308D3 is proposed to be the essential UGT gene for the participation in biotransformation in pyrethroid detoxification process, which is possibly regulated by eight SNPs in its 3′ flanking region. The UGT302A3 is also associated with pyrethroid resistance, and four amino acid mutations in its coding sequences might enhance its catalytic activity and further result in higher insecticide resistance.ConclusionsThis study provides the diversity, phylogenetics and evolution of UGT genes, and potential UGT members and mutations involved in pyrethroid resistance in An. sinensis, and lays an important basis for the better understanding and further research on UGT function in defense against insecticide stress.

Characterisation of Anopheles gambiae heme oxygenase and metalloporphyrin feeding suggests a potential role in reproduction

The mosquito Anopheles gambiae is the principal vector for malaria in sub-Saharan Africa. The ability of A. gambiae to transmit malaria is strictly related to blood feeding and digestion, which releases nutrients for oogenesis, as well as substantial amounts of highly toxic free heme. Heme degradation by heme oxygenase (HO) is a common protective mechanism, and a gene for HO exists in the An. gambiae genome HO (AgHO), although it has yet to be functionally examined. Here, we have cloned and expressed An. gambiae HO (AgHO) in E. coli. Purified recombinant AgHO bound hemin stoichiometrically to form a hemin–enzyme complex similar to other HOs, with a KD of 3.9 ± 0.6 μM; comparable to mammalian and bacterial HOs, but 7-fold lower than that of Drosophila melanogaster HO. AgHO also degraded hemin to biliverdin and released CO and iron in the presence of NADPH cytochrome P450 oxidoreductase (CPR). Optimal AgHO activity was observed at 27.5 °C and pH 7.5. To investigate effects of AgHO inhibition, adult female A. gambiae were fed heme analogues Sn- and Zn-protoporphyrins (SnPP and ZnPP), known to inhibit HO. These led to a dose dependent decrease in oviposition. Cu-protoporphyrin (CuPP), which does not inhibit HO had no effect. These results demonstrate that AgHO is a catalytically active HO and that it may play a key role in egg production in mosquitoes. It also presents a potential target for the development of compounds aimed at sterilising mosquitoes for vector control.

Partial-arm translocations in evolution of malaria mosquitoes revealed by high-coverage physical mapping of the Anopheles atroparvus genome

BackgroundMalaria mosquitoes have had a remarkable stability in the number of chromosomes in their karyotype (2n = 6) during 100 million years of evolution. Moreover, autosomal arms were assumed to maintain their integrity even if their associations with each other changed via whole-arm translocations. Here we use high-coverage comparative physical genome mapping of three Anopheles species to test the extent of evolutionary conservation of chromosomal arms in malaria mosquitoes.ResultsIn this study, we developed a physical genome map for Anopheles atroparvus, one of the dominant malaria vectors in Europe. Using fluorescence in situ hybridization (FISH) of DNA probes with the ovarian nurse cell polytene chromosomes and synteny comparison, we anchored 56 genomic scaffolds to the An. atroparvus chromosomes. The obtained physical map represents 89.6% of the An. atroparvus genome. This genome has the second highest mapping coverage among Anophelinae assemblies after An. albimanus, which has 98.2% of the genome assigned to its chromosomes. A comparison of the An. atroparvus, An. albimanus, and An. gambiae genomes identified partial-arm translocations between the autosomal arms that break down the integrity of chromosome elements in evolution affecting the structure of the genetic material in the pericentromeric regions. Unlike An. atroparvus and An. albimanus, all chromosome elements of An. gambiae are fully syntenic with chromosome elements of the putative ancestral Anopheles karyotype. We also detected nonrandom distribution of large conserved synteny blocks and confirmed a higher rate of inversion fixation in the X chromosome compared with autosomes.ConclusionsOur study demonstrates the power of physical mapping for understanding the genome evolution in malaria mosquitoes. The results indicate that syntenic relationships among chromosome elements of Anopheles species have not been fully preserved because of multiple partial-arm translocations.

Multi-tissue GAL4-mediated gene expression in all Anopheles gambiae life stages using an endogenous polyubiquitin promoter

The ability to manipulate the Anopheles gambiae genome and alter gene expression effectively and reproducibly is a prerequisite for functional genetic analysis and for the development of novel control strategies in this important disease vector. However, in vivo transgenic analysis in mosquitoes is limited by the lack of promoters active ubiquitously. To address this, we used the GAL4/UAS system to investigate the promoter of the An. gambiae Polyubiquitin-c (PUBc) gene and demonstrated its ability to drive expression in mosquito cell culture before incorporation into An. gambiae transgenic driver lines. To generate such lines, piggyBac-mediated insertion was used to identify genomic regions able to sustain widespread expression and to create φC31 docking lines at these permissive sites. Patterns of expression induced by PUBc-GAL4 drivers carrying single intergenic insertions were assessed by crossing with a novel responder UAS-mCD8:mCherry line that was created by φC31-mediated integration. Amongst the drivers created at single, unique chromosomal integration loci, two were isolated that induced differential expression levels in a similar multiple-tissue spatial pattern throughout the mosquito life cycle. This work expands the tools available for An. gambiae functional analysis by providing a novel promoter for investigating phenotypes resulting from widespread multi-tissue expression, as well as identifying and tagging genomic sites that sustain broad transcriptional activity.

Serine protease-related proteins in the malaria mosquito, Anopheles gambiae

Insect serine proteases (SPs) and serine protease homologs (SPHs) participate in digestion, defense, development, and other physiological processes. In mosquitoes, some clip-domain SPs and SPHs (i.e. CLIPs) have been investigated for possible roles in antiparasitic responses. In a recent test aimed at improving quality of gene models in the Anopheles gambiae genome using RNA-seq data, we observed various discrepancies between gene models in AgamP4.5 and corresponding sequences selected from those modeled by Cufflinks, Trinity and Bridger. Here we report a comparative analysis of the 337 SP-related proteins in A. gambiae by examining their domain structures, sequence diversity, chromosomal locations, and expression patterns. One hundred and ten CLIPs contain 1 to 5 clip domains in addition to their protease domains (PDs) or non-catalytic, protease-like domains (PLDs). They are divided into five subgroups: CLIPAs (22) are clip1−5-PLD; CLIPBs (29), CLIPCs (12) and CLIPDs (14) are mainly clip-PD; most CLIPEs (33) have a domain structure of PD/PLD-PLD-clip-PLD0−1. While expression of the CLIP genes in group-1 is generally low and detected in various tissue- and stage-specific RNA-seq libraries, some putative GPs/GPHs (i.e. single domain gut SPs/SPHs) in group-2 are highly expressed in midgut, whole larva or whole adult libraries. In comparison, 46 SPs, 26 SPHs, and 37 multi-domain SPs/SPHs (i.e. PD/PLD-PLD≥1) in group-3 do not seem to be specifically expressed in digestive tract. There are 16 SPs and 2 SPH containing other types of putative regulatory domains (e.g. LDLa, CUB, Gd). Of the 337 SP and SPH genes, 159 were sorted into 46 groups (2–8 members/group) based on similar phylogenetic tree position, chromosomal location, and expression profile. This information and analysis, including improved gene models and protein sequences, constitute a solid foundation for functional analysis of the SP-related proteins in A. gambiae.

The Physical Genome Mapping of Anopheles albimanus Corrected Scaffold Misassemblies and Identified Interarm Rearrangements in Genus Anopheles.

The genome of the Neotropical malaria vector Anopheles albimanus was sequenced as part of the 16 Anopheles Genomes Project published in 2015. The draft assembly of this species consisted of 204 scaffolds with an N50 scaffold size of 18.1 Mb and a total assembly size of 170.5 Mb. It was among the smallest genomes with the longest scaffolds in the 16 Anopheles species cluster, making An. albimanus the logical choice for anchoring the genome assembly to chromosomes. In this study, we developed a high-resolution cytogenetic photomap with completely straightened polytene chromosomes from the salivary glands of the mosquito larvae. Based on this photomap, we constructed a chromosome-based genome assembly using fluorescent in situ hybridization of PCR-amplified DNA probes. Our physical mapping, assisted by an ortholog-based bioinformatics approach, identified and corrected nine misassemblies in five large genomic scaffolds. Misassemblies mostly occurred in junctions between contigs. Our comparative analysis of scaffolds with the An. gambiae genome detected multiple genetic exchanges between pericentromeric regions of chromosomal arms caused by partial-arm translocations. The final map consists of 40 ordered genomic scaffolds and corrected fragments of misassembled scaffolds. The An. albimanus physical map comprises 98.2% of the total genome assembly and represents the most complete genome map among mosquito species. This study demonstrates that physical mapping is a powerful tool for correcting errors in draft genome assemblies and for creating chromosome-anchored reference genomes.

Increased production of piRNAs from euchromatic clusters and genes in Anopheles gambiae compared with Drosophila melanogaster.

BackgroundSpecific genomic loci, termed Piwi-interacting RNA (piRNA) clusters, manufacture piRNAs that serve as guides for the inactivation of complementary transposable elements (TEs). The piRNA pathway has been accurately detailed in Drosophila melanogaster, while it remains poorly examined in other insects. This pathway is increasingly recognized as critical for germline development and reproduction. Understanding of the piRNA functions in mosquitoes could offer an opportunity for disease vector control by the reduction of their reproductive potential.ResultsTo analyze the similarities and differences in this pathway between Drosophila and mosquito, we performed an in-depth analysis of the genomic loci producing piRNAs and their targets in the African malaria vector Anopheles gambiae. We identified 187 piRNA clusters in the An. gambiae genome and 155 piRNA clusters in the D. melanogaster genome. We demonstrate that many more piRNA clusters in the mosquito compared with the fruit fly are uni-directionally transcribed and are located outside pericentromeric heterochromatin. About 11 % of the An. gambiae piRNA population map to gene transcripts. This is a noticeable increase compared with the ~6 % of the piRNA population mapped to genes in D. melanogaster. A subset of the piRNA-enriched genes in An. gambiae has functions related to reproduction and development. At least 24 and 65 % of the mapped piRNAs correspond to genomic TE sequences in An. gambiae and D. melanogaster, respectively. DNA transposons and non-LTR retrotransposons are more abundant in An. gambiae, while LTR retrotransposons are more abundant in D. melanogaster. Yet, piRNAs predominantly target LTR retrotransposons in both species, which may point to a distinct feature of these elements compared to the other classes of TEs concerning their silencing by the piRNA pathway.ConclusionsHere, we demonstrate that piRNA-producing loci have more ubiquitous distribution in the An. gambiae genome than in the genome of D. melanogaster. Also, protein-coding genes have an increased role in production of piRNAs in the germline of this mosquito. Genes involved in germline and embryonic development of An. gambiae generate a substantial portion of piRNAs, suggesting a role of the piRNA pathway in the epigenetic regulation of the reproductive processes in the African malaria vector.Electronic supplementary materialThe online version of this article (doi:10.1186/s13072-015-0041-5) contains supplementary material, which is available to authorized users.

The fucomic potential of mosquitoes: Fucosylated N-glycan epitopes and their cognate fucosyltransferases

Fucoconjugates are key mediators of protein-glycan interactions in prokaryotes and eukaryotes. As examples, N-glycans modified with the non-mammalian core α1,3-linked fucose have been detected in various organisms ranging from plants to insects and are immunogenic in mammals. The rabbit polyclonal antibody raised against plant horseradish peroxidase (anti-HRP) is able to recognize the α1,3-linked fucose epitope and is also known to specifically stain neural tissues in the fruit fly Drosophila melanogaster. In this study, we have detected and localized the anti-HRP cross-reactivity in another insect species, the malaria mosquito vector Anopheles gambiae. We were able to identify and structurally elucidate fucosylated N-glycans including core mono- and difucosylated structures (responsible for anti-HRP cross reactivity) as well as a Lewis-type antennal modification on mosquito anionic N-glycans by applying enzymatic and chemical treatments. The three mosquito fucosyltransferase open reading frames (FucT6, FucTA and FucTC) required for the in vivo biosynthesis of the fucosylated N-glycan epitopes were identified in the Anopheles gambiae genome, cloned and recombinantly expressed in Pichia pastoris. Using a robust MALDI-TOF MS approach, we characterised the activity of the three recombinant fucosyltransferases in vitro and demonstrate that they share similar enzymatic properties as compared to their homologues from D. melanogaster and Apis mellifera. Thus, not only do we confirm the neural reactivity of anti-HRP in a mosquito species, but also demonstrate enzymatic activity for all its α1,3- and α1,6-fucosyltransferase homologues, whose specificity matches the results of glycomic analyses.

Constructing a Genome-Wide LD Map of Wild A. gambiae Using Next-Generation Sequencing.

Anopheles gambiae is the major malaria vector in Africa. Examining the molecular basis of A. gambiae traits requires knowledge of both genetic variation and genome-wide linkage disequilibrium (LD) map of wild A. gambiae populations from malaria-endemic areas. We sequenced the genomes of nine wild A. gambiae mosquitoes individually using next-generation sequencing technologies and detected 2,219,815 common single nucleotide polymorphisms (SNPs), 88% of which are novel. SNPs are not evenly distributed across A. gambiae chromosomes. The low SNP-frequency regions overlay heterochromatin and chromosome inversion domains, consistent with the lower recombinant rates at these regions. Nearly one million SNPs that were genotyped correctly in all individual mosquitoes with 99.6% confidence were extracted from these high-throughput sequencing data. Based on these SNP genotypes, we constructed a genome-wide LD map for wild A. gambiae from malaria-endemic areas in Kenya and made it available through a public Website. The average size of LD blocks is less than 40 bp, and several large LD blocks were also discovered clustered around the para gene, which is consistent with the effect of insecticide selective sweeps. The SNPs and the LD map will be valuable resources for scientific communities to dissect the A. gambiae genome.

In depth annotation of the Anopheles gambiae mosquito midgut transcriptome.

BackgroundGenome sequencing of Anopheles gambiae was completed more than ten years ago and has accelerated research on malaria transmission. However, annotation needs to be refined and verified experimentally, as most predicted transcripts have been identified by comparative analysis with genomes from other species. The mosquito midgut—the first organ to interact with Plasmodium parasites—mounts effective antiplasmodial responses that limit parasite survival and disease transmission. High-throughput Illumina sequencing of the midgut transcriptome was used to identify new genes and transcripts, contributing to the refinement of An. gambiae genome annotation.ResultsWe sequenced ~223 million reads from An. gambiae midgut cDNA libraries generated from susceptible (G3) and refractory (L35) mosquito strains. Mosquitoes were infected with either Plasmodium berghei or Plasmodium falciparum, and midguts were collected after the first or second Plasmodium infection. In total, 22,889 unique midgut transcript models were generated from both An. gambiae strain sequences combined, and 76% are potentially novel. Of these novel transcripts, 49.5% aligned with annotated genes and appear to be isoforms or pre-mRNAs of reference transcripts, while 50.5% mapped to regions between annotated genes and represent novel intergenic transcripts (NITs). Predicted models were validated for midgut expression using qRT-PCR and microarray analysis, and novel isoforms were confirmed by sequencing predicted intron-exon boundaries. Coding potential analysis revealed that 43% of total midgut transcripts appear to be long non-coding RNA (lncRNA), and functional annotation of NITs showed that 68% had no homology to current databases from other species. Reads were also analyzed using de novo assembly and predicted transcripts compared with genome mapping-based models. Finally, variant analysis of G3 and L35 midgut transcripts detected 160,742 variants with respect to the An. gambiae PEST genome, and 74% were new variants. Intergenic transcripts had a higher frequency of variation compared with non-intergenic transcripts.ConclusionThis in-depth Illumina sequencing and assembly of the An. gambiae midgut transcriptome doubled the number of known transcripts and tripled the number of variants known in this mosquito species. It also revealed existence of a large number of lncRNA and opens new possibilities for investigating the biological function of many newly discovered transcripts.Electronic supplementary materialThe online version of this article (doi:10.1186/1471-2164-15-636) contains supplementary material, which is available to authorized users.

ACE for all - a molecular perspective.

Angiotensin-I converting enzyme (ACE, EC 3.4.15.1) is a zinc dependent dipeptidyl carboxypeptidase with an essential role in mammalian blood pressure regulation as part of the renin-angiotensin aldosterone system (RAAS). As such, it has long been targeted in the treatment of hypertension through the use of ACE inhibitors. Although ACE has been studied since the 1950s, only recently have the full range of functions of this enzyme begun to truly be appreciated. ACE homologues have been found in a host of other organisms, and are now known to be conserved in insects. Insect ACE homologues typically share over 30 % amino acid sequence identity with human ACE. Given that insects lack a mammalian type circulatory system, they must have crucial roles in other physiological processes. The first ACE crystal structures were reported during the last decade and have enabled these enzymes to be studied from an entirely different perspective. Here we review many of these key developments and the implications that they have had on our understanding of the diverse functions of these enzymes. Specifically, we consider how structural information is being used in the design of a new generation of ACE inhibitors with increased specificity, and how the structures of ACE homologues are related to their functions. The Anopheles gambiae genome is predicted to code for ten ACE homologues, more than any genome studied so far. We have modelled the active sites of some of these as yet uncharacterised enzymes to try and infer more about their potential roles at the molecular level.