Genome sequence of Anopheles sinensis provides insight into genetics basis of mosquito competence for malaria parasites

Uncovering nutrients and energy related gene functions of black soldier fly Hermetia illucens strain KUP

The black soldier fly (Hermetia illucens) has emerged as a significant insect species in the decomposition of organic waste for sustainable agricultural practices. Due to its remarkable characteristics and performance, H. illucens is increasingly utilised for insect farming, particularly for industrial-scale rearing throughout the world. In this study, we employed whole-genome sequencing to annotate the gene and protein functions of H. illucens and to explore the functional genomics related to nutrients and energy. As a result, a genome size of H. illucens strain KUP 1.68 Gb with a GC content of 42.13 % was achieved. Of the 14,036 coding sequences, we determined the function of 12,046 protein-coding genes. Based on metabolic functional assignment, we classified 4,218 protein-coding genes; the main category was metabolism (32.86 %). Comparative genomic analysis across the other H. illucens strain and insect species revealed that the major metabolic gene functions and pathways related to nutrient and energy sources of H. illucens KUP are involved in key amino acid metabolism (e.g., cysteine and methionine) as well as fatty acid biosynthesis and glycerolipid metabolism. These findings underscore the metabolic capability and versatility of H. illucens, which is regarded as a potential source of proteins and lipids. Our study contributes to the knowledge regarding the feed utilisation of H. illucens and offers insights into transforming waste into valuable products. H. illucens has the potential to create globally sustainable nutrients and environmentally friendly solutions, aligning with the goal of responsible resource utilisation.

An insight into the female and male Sabethes cyaneus mosquito salivary glands transcriptome

Mosquitoes are responsible for the death and debilitation of millions of people every year due to the pathogens they can transmit while blood feeding. While a handful of mosquitoes, namely those in the Aedes, Anopheles, and Culex genus, are the dominant vectors, many other species belonging to different genus are also involved in various pathogen cycles. Sabethes cyaneus is one of the many poorly understood mosquito species involved in the sylvatic cycle of Yellow Fever Virus. Here, we report the expression profile differences between male and female of Sa.cyaneus salivary glands (SGs). We find that female Sa.cyaneus SGs have 165 up-regulated and 18 down-regulated genes compared to male SGs. Most of the up-regulated genes have unknown functions, however, odorant binding proteins, such as those in the D7 protein family, and mucins were among the top 30 genes. We also performed various in vitro activity assays of female SGs. In the activity analysis we found that female SG extracts inhibit coagulation by blocking factor Xa and has endonuclease activity. Knowledge about mosquitoes and their physiology are important for understanding how different species differ in their ability to feed on and transmits pathogens to humans. These results provide us with an insight into the Sabethes SG activity and gene expression that expands our understanding of mosquito salivary glands.

Dipteran endoparasitoid Exorista bombycis utilizes antihemocyte components against host defense of silkworm Bombyx mori

Dipteran endoparasitoids avoid host immune response; however, antidefense components from the Dipterans are unknown. Infestation of commercial silkworm Bombyx mori Linnaeus (Lepidoptera: Bombycidae) by endoparasitoid Exorista bombycis Louis (Diptera: Tachinidae) induced immune reactions, cytotoxicity, granulation, degranulation, and augmented release of cytotoxic marker enzyme lactate dehydrogenase (LDH), and degranulation-mediator enzyme β-hexosaminidase in hemocytes. In this study, by reverse phase high-performance liquid chromatography, fractions of E. bombycis larval tissue protein with antihemocytic activity are separated. From the fraction, peptides of hemocyte aggregation inhibitor protein (HAIP) and pyridoxamine phosphate oxidase (PNPO) are identified by mass spectrometry. Interacting partners of HAIP and PNPO are retrieved that further enhance the virulence of the parasitoid. PNPO and HAIP genes showed a four- to seven fold increase in expression in the integument of the parasitoid larva. Together, the dipteran endoparasitoid E. bombycis exploit antihemocyte activity to inhibit host defense reactions in addition to defense evasion contemplated.

CRISPR/Cas9 mediates efficient site-specific mutagenesis of the odorant receptor co-receptor (Orco) in the malaria vector Anopheles sinensis

BACKGROUND

Anopheles sinensis is the most widely distributed mosquito species and is the main transmitter of Plasmodium vivax malaria in China. Most previous research has focused on the mechanistic understanding of biological processes in An. sinensis and novel ways of interrupting malaria transmission. However, the development of functional genomics and genetics-based vector control strategies against An. sinensis remain limited because of insufficient site-specific genome editing tools.

RESULTS

We report the first successful application of the CRISPR/Cas9 mediated knock-in for highly efficient, site-specific mutagenesis in An. sinensis. The EGFP marker gene driven by the 3 × P3 promoter was precisely integrated into the odorant receptor co-receptor (Orco) by direct injections of Cas9 protein, double-stranded DNA donor, and Orco-gRNA. We achieved a mutation rate of 3.77%, similar to rates in other mosquito species. Precise knock-in at the intended locus was confirmed by polymerase chain reaction (PCR) amplification and sequencing. The Orco mutation severely impaired mosquito sensitivity to some odors and their ability to locate and discriminate a human host.

CONCLUSION

Orco was confirmed as a key mediator of multiple olfactory-driven behaviors in the An. sinensis life cycle, highlighting the importance of Orco as a key molecular target for malaria control. The results also demonstrated that CRISPR/Cas9 was a simple and highly efficient genome editing technique for An. sinensis and could be used to develop genetic control tools for this vector. © 2022 Society of Chemical Industry.

Genome-wide identification and expression profiling of odorant receptor genes in the malaria vector Anopheles sinensis

Background

The olfactory system plays a crucial role in regulating insect behaviors. The detection of odorants is mainly mediated by various odorant receptors (ORs) that are expressed in the dendrites of olfactory neurons of chemosensilla. Anophelessinensis is a major malaria vector in Eastern Asia and its genome has recently been successfully sequenced and annotated. In this study, we present genome-wide identification and expression profiling of OR genes in different chemosensory tissues of An.sinensis.

Methods

The OR genes were identified using the available genome sequences of An.sinensis. A series of bioinformatics analyses were conducted to investigate the structure, genome distribution, selective pressure and phylogenetic relationships of OR genes, the conserved domains and specific functional sites in the OR amino acid sequences. The expression levels of OR genes were analyzed from transcriptomic data from An.sinensis antennae, proboscis and maxillary palps of both sexes.

Results

A total of 59 putative OR genes have been identified and characterized in An.sinensis. This number is significantly less than that in An.gambiae. Whether this difference is caused by the contraction or expansion of OR genes after divergence of the two species remains unknown. The RNA-seq analysis showed that AsORs have obvious tissue- and sex-specific expression patterns. Most AsORs are highly expressed in the antennae and the expression pattern and number of AsORs expressed in antennae are similar in males and females. However, the relative levels of AsOR transcripts are much higher in female antennae than in male antennae, which indicates that the odor sensitivity is likely to be increased in female mosquitoes. Based on the expression patterns and previous studies, we have speculated on the functions of some OR genes but this needs to be validated by further behavioral, molecular and electrophysiological studies. Further studies are necessary to compare the olfactory-driven behaviors and identify receptors that respond strongly to components of human odors that may act in the process of human recognition.

Conclusions

This is the first genome-wide analysis of the entire repertoire of OR genes in An.sinensis. Characterized features and profiled expression patterns of ORs suggest their involvement in the odorous reception of this species. Our findings provide a basis for further research on the functions of OR genes and additional genetic and behavioral targets for more sustainable management of An.sinensis in the future.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-022-05259-x.

Whole Genome Sequencing Contributions and Challenges in Disease Reduction Focused on Malaria

Malaria elimination remains an important goal that requires the adoption of sophisticated science and management strategies in the era of the COVID-19 pandemic. The advent of next generation sequencing (NGS) is making whole genome sequencing (WGS) a standard today in the field of life sciences, as PCR genotyping and targeted sequencing provide insufficient information compared to the whole genome. Thus, adapting WGS approaches to malaria parasites is pertinent to studying the epidemiology of the disease, as different regions are at different phases in their malaria elimination agenda. Therefore, this review highlights the applications of WGS in disease management, challenges of WGS in controlling malaria parasites, and in furtherance, provides the roles of WGS in pursuit of malaria reduction and elimination. WGS has invaluable impacts in malaria research and has helped countries to reach elimination phase rapidly by providing required information needed to thwart transmission, pathology, and drug resistance. However, to eliminate malaria in sub-Saharan Africa (SSA), with high malaria transmission, we recommend that WGS machines should be readily available and affordable in the region.

Parallel Evolution of C-Type Lectin Domain Gene Family Sizes in Insect-Vectored Nematodes

The dispersal stage of pathogens is crucial for the successful spread and infection of their hosts. Some plant-parasitic nematodes (PPNs) have evolved specialized dispersal stages to reach healthy hosts by being carried out by insect vectors. Because gene gain and loss is a major factor contributing to the evolution of novel characteristics, it is essential to clarify the gene family characteristics among nematodes with different dispersal modes to disentangle the evolution of insect-mediated dispersal. Here, the size of the C-type lectin (CTL) family genes of insect-vectored nematodes was found to be drastically reduced compared with those of self-dispersing nematodes, whereas the diversity of their functional domains was significantly higher. The gene family sizes of vector-dispersed nematodes were only a twentieth of the size of that of a self-dispersing (i.e., without a biotic vector) nematode model Caenorhabditis elegans, and these genes were inactive during the dispersal stage. Phylogenetic analysis showed that some CTL genes of vector-borne PPNs shared higher homology to the animal parasitic nematodes compared with other PPNs. Moreover, homology modeling predicted that the CTLs of insect-vectored nematodes bear remarkable structural similarity to the lectin genes of their vector's immune system. Because CTL genes are important sugar-binding proteins for the innate immune response of C. elegans, the loss of some CTL genes of vector-transmitted PPNs might be responsible for their parallel adaptations to a mutualistic relationship with their vector. These results expand our understanding of the evolutionary benefits of vector-mediated transmission for the nematode and vector-nematode co-evolution.

Improving mosquito control strategies with population genomics

Mosquito control strategies increasingly apply knowledge from population genomics research. This review highlights recent applications to three research domains: mosquito invasions, insecticide resistance evolution, and rear and release programs. Current research trends follow developments in reference assemblies, either as improvements to existing assemblies (particularly Aedes) or assemblies for new taxa (particularly Anopheles). With improved assemblies, studies of invasive and rear and release target populations are better able to incorporate adaptive as well as demographic hypotheses. New reference assemblies are aiding comparisons of insecticide resistance across sister taxa while helping resolve taxon boundaries amidst frequent introgression. Anopheles gene drive deployments and improved Aedes genome assemblies should lead to a convergence in research aims for Anopheles and Aedes in the coming years.

Electroporation-mediated nucleic acid delivery during non-embryonic stages for gene-function analysis in Anopheles sinensis

The delivery of exogenous nucleic acids to eggs or non-embryonic individuals by microinjection is a vital reverse genetics technique used to determine gene function in mosquitoes. However, DNA delivery to eggs is complex and time-consuming, and conventional, non-embryonic-injection techniques may result in unobvious phenotypes caused by insufficient absorption of nucleic acid fragments by cells at target body parts or tissues. In this study, we developed a set of electroporation-mediated non-embryonic microinjections for the delivery of exogenous nucleic acids in Anopheles sinensis. Gene silencing using this method led to down-regulation of target gene expression (AsCPR128) by 77% in targeted body parts, compared with only 10% in non-targeted body parts, thus increasing the defect-phenotype rate in the target area by 5.3-fold, compared with non-shock injected controls. Electroporation-mediated somatic transgenesis resulted in stable phenotypic characteristics of the reporter gene at the shocked body parts during the pupal-adult stages in about 69% of individuals. Furthermore, injecting plasmid DNA near the ovaries of female mosquitoes after a blood meal followed by electric shock produced three germline G₁ transgenic lines, with a transformation rate of about 11.1% (calculated from ovulatory G₀ females). Among the positive G₁ lines, 42%, 40%, and 31% of individuals emitted red fluorescence in the larval stage. When the red fluorescent larvae developed into adults, green fluorescence was emitted from the ovaries of the females upon feeding. These results suggest that electroporation-mediated non-embryonic microinjection can be an efficient, rapid, and simple technique for analyzing gene function in non-model mosquitoes or other small insects.

The algal selenoproteomes

Background

Selenium is an essential trace element, and selenocysteine (Sec, U) is its predominant form in vivo. Proteins that contain Sec are selenoproteins, whose special structural features include not only the TGA codon encoding Sec but also the SECIS element in mRNA and the conservation of the Sec-flanking region. These unique features have led to the development of a series of bioinformatics methods to predict and research selenoprotein genes. There have been some studies and reports on the evolution and distribution of selenoprotein genes in prokaryotes and multicellular eukaryotes, but the systematic analysis of single-cell eukaryotes, especially algae, has been very limited.

Results

In this study, we predicted selenoprotein genes in 137 species of algae by using a program we previously developed. More than 1000 selenoprotein genes were obtained. A database website was built to record these algae selenoprotein genes (www.selenoprotein.com). These genes belong to 42 selenoprotein families, including three novel selenoprotein gene families.

Conclusions

This study reveals the primordial state of the eukaryotic selenoproteome. It is an important clue to explore the significance of selenium for primordial eukaryotes and to determine the complete evolutionary spectrum of selenoproteins in all life forms.

Evolutionary superscaffolding and chromosome anchoring to improve Anopheles genome assemblies

Background

New sequencing technologies have lowered financial barriers to whole genome sequencing, but resulting assemblies are often fragmented and far from ‘finished’. Updating multi-scaffold drafts to chromosome-level status can be achieved through experimental mapping or re-sequencing efforts. Avoiding the costs associated with such approaches, comparative genomic analysis of gene order conservation (synteny) to predict scaffold neighbours (adjacencies) offers a potentially useful complementary method for improving draft assemblies.

Results

We evaluated and employed 3 gene synteny-based methods applied to 21 Anopheles mosquito assemblies to produce consensus sets of scaffold adjacencies. For subsets of the assemblies, we integrated these with additional supporting data to confirm and complement the synteny-based adjacencies: 6 with physical mapping data that anchor scaffolds to chromosome locations, 13 with paired-end RNA sequencing (RNAseq) data, and 3 with new assemblies based on re-scaffolding or long-read data. Our combined analyses produced 20 new superscaffolded assemblies with improved contiguities: 7 for which assignments of non-anchored scaffolds to chromosome arms span more than 75% of the assemblies, and a further 7 with chromosome anchoring including an 88% anchored Anopheles arabiensis assembly and, respectively, 73% and 84% anchored assemblies with comprehensively updated cytogenetic photomaps for Anopheles funestus and Anopheles stephensi.

Conclusions

Experimental data from probe mapping, RNAseq, or long-read technologies, where available, all contribute to successful upgrading of draft assemblies. Our evaluations show that gene synteny-based computational methods represent a valuable alternative or complementary approach. Our improved Anopheles reference assemblies highlight the utility of applying comparative genomics approaches to improve community genomic resources.

Variant Ionotropic Receptors are Expressed in the Antennae of Anopheles sinensis (Diptera: Culicidae)

Mosquitoes transmit many harmful diseases that seriously threaten public health. The mosquito's olfactory system is of great significance for host selection. Inotropic receptors (IRs) and olfactory receptors (ORs) have been demonstrated to be capable of odorant molecular recognition. Analyzing the molecular principles of mosquito olfaction facilitates the development of prevention and therapy techniques. Advances in the understanding of IRs have been seriously inadequate compared to those of ORs. Here, we provide evidence that 35 Anopheles sinensis IR (AsIR) genes are expressed, 7 of which are in the antennae and 2 have expression levels that are upregulated with a blood meal. A homologous analysis of the sequences showed that AsIRs are a subfamily of ionotropic glutamate receptors (iGLURs). This is the first that time IRs have been identified in Anopheles sinensis in vitro. The ultrastructure of the antennae supports the theory that diverse sensilla are distributed in the antennae. The results here may facilitate the revelation of the regulation mechanism in AsIRs, which could mitigate the transmission of diseases by mosquitoes.

An odorant receptor from Anopheles sinensis in China is sensitive to oviposition attractants

Background

Anopheles sinensis is an important vector for the spread of malaria in China. Olfactory-related behaviours, particularly oviposition site seeking, offer opportunities for disrupting the disease-transmission process.

Results

This is the first report of the identification and characterization of AsinOrco and AsinOR10 in An. sinensis. AsinOrco and AsinOR10 share 97.49% and 90.37% amino acid sequence identity, respectively, with related sequences in Anopheles gambiae. A functional analysis demonstrated that AsinOrco- and AsinOR10-coexpressing HEK293 cells were highly sensitive to 3-methylindole, but showed no significant differences in response to other test odorants when compared to DMSO.

Conclusions

AsinOrco was characterized as a new member of the Orco ortholog subfamily. AsinOR10, which appears to be a member of the OR2-10 subfamily, is directly involved in identification of oviposition sites. This finding will help to elucidate the molecular mechanisms underlying olfactory signaling in An. sinensis and provide many more molecular targets for eco-friendly pest control.

Electronic supplementary material

The online version of this article (10.1186/s12936-018-2501-4) contains supplementary material, which is available to authorized users.

Genome-wide identification, characterization and evolution of cuticular protein genes in the malaria vector Anopheles sinensis (Diptera: Culicidae)

Thirteen cuticular protein (CP) families have been recognized in arthropods. In this study, 250 Anopheles sinensis CP genes were identified and named based on genome and transcriptome sequences. They were classified into 10 families based on motifs and phylogenetic analyses. In 11 other insect species, nine had CP numbers > 150 while Apis mellifera and Tribolium castaneum had CP numbers less than 52. The CPs of eight species occupied > 1.4% of the total genomic gene number, whereas in three species the CPs occupied < 1%. The phylogenies for each CP family in An. sinensis were constructed and discussed. The 250 CPs each had 1-8 exons with 144 CPs (57.6%) having two exons. The intron length ranged from 66-3888 bp with 174 introns (54.0%) being 66-100 bp long. Except for two CPs on two contigs, 248 CPs were mapped onto 28 scaffolds with 136 genes (54.4%) restricted to five scaffolds. A total of 107 CPs were clustered and located at 27 loci. The CPR family had the conserved motif GSYSLVEPDGTVRTV. The RR-1 subfamily had an additional 21 amino acid (aa) motifs with the YVADENGF sequence that is common in insects. The RR-2 subfamily had an additional 50 aa motifs with two additional regions RDGDVVKG and G-x(3)-VV. A comparison with 115 orthologous counterparts of An. gambiae CPs suggested purifying selection for all of these genes. This study provides basic information useful for further studies on biological functions of An. sinensis CPs as well as for comparative genomics of insect CPs.

Sialotranscriptome sequencing and analysis of Anopheles sinensis and comparison with Psorophora albipes sialotranscriptome (Diptera: Culicidae)

Most of adult female mosquitoes secrete saliva to facilitate blood sucking, digestion and nutrition, and mosquito-borne disease prevention. The knowledge of classification and characteristics of sialotranscriptome genes are still quite limited. Anopheles sinensis is a major malaria vector in China and southeast Asian countries. In this study, the An. sinensis sialotranscriptome was sequenced using Illumina sequencing technique with a total of 10 907 unigenes to be obtained and annotated in biological functions and pathways, and 10 470 unigenes were mapped to An. sinensis reference genome with 70.46% of genes having 90%-100% genome mapping through bioinformatics analysis. These mapped genes were classified into four categories: housekeeping (6632 genes), secreted (1177), protein-coding genes with function-unknown (2646) and transposable element (15). The housekeeping genes were divided into 27 classes, and the secreted genes were divided into 11 classes and 96 families. The classification, characteristics and evolution of these classes/families of secreted genes are further described and discussed. The comparison of the 1177 secreted genes in An. sinensis in the Anophelinae subfamily with 811 in Psorophora albipes in the Culicinae subfamily show that six classes/subclasses have the gene number more than twice and two classes (uniquely found in anophelines, and Orphan proteins of unique standing) are unique in the former compared with the latter, whereas four classes/subclasses are much expanded and uniquely found in the Aedes class and is unique in the later. The An. sinensis sialotranscriptome sequence data is the most complete in mosquitoes to date, and the analyses provide a comprehensive information frame for further research of mosquito sialotranscriptome.

Genome-wide identification, characterization and classification of ionotropic glutamate receptor genes (iGluRs) in the malaria vector Anopheles sinensis (Diptera: Culicidae)

BACKGROUND

Ionotropic glutamate receptors (iGluRs) are conserved ligand-gated ion channel receptors, and ionotropic receptors (IRs) were revealed as a new family of iGluRs. Their subdivision was unsettled, and their characteristics are little known. Anopheles sinensis is a major malaria vector in eastern Asia, and its genome was recently well sequenced and annotated.

METHODS

We identified iGluR genes in the An. sinensis genome, analyzed their characteristics including gene structure, genome distribution, domains and specific sites by bioinformatic methods, and deduced phylogenetic relationships of all iGluRs in An. sinensis, Anopheles gambiae and Drosophila melanogaster. Based on the characteristics and phylogenetics, we generated the classification of iGluRs, and comparatively analyzed the intron number and selective pressure of three iGluRs subdivisions, iGluR group, Antenna IR and Divergent IR subfamily.

RESULTS

A total of 56 iGluR genes were identified and named in the whole-genome of An. sinensis. These genes were located on 18 scaffolds, and 31 of them (29 being IRs) are distributed into 10 clusters that are suggested to form mainly from recent gene duplication. These iGluRs can be divided into four groups: NMDA, non-NMDA, Antenna IR and Divergent IR based on feature comparison and phylogenetic analysis. IR8a and IR25a were suggested to be monophyletic, named as Putative in the study, and moved from the Antenna subfamily in the IR family to the non-NMDA group as a sister of traditional non-NMDA. The generated iGluRs of genes (including NMDA and regenerated non-NMDA) are relatively conserved, and have a more complicated gene structure, smaller ω values and some specific functional sites. The iGluR genes in An. sinensis, An. gambiae and D. melanogaster have amino-terminal domain (ATD), ligand binding domain (LBD) and Lig_Chan domains, except for IR8a that only has the LBD and Lig_Chan domains. However, the new concept IR family of genes (including regenerated Antenna IR, and Divergent IR), especially for Divergent IR are more variable, have a simpler gene structure (intron loss phenomenon) and larger ω values, and lack specific functional sites. These IR genes have no other domains except for Antenna IRs that only have the Lig_Chan domain.

CONCLUSIONS

This study provides a comprehensive information framework for iGluR genes in An. sinensis, and generated the classification of iGluRs by feature and bioinformatics analyses. The work lays the foundation for further functional study of these genes.

ASGDB: a specialised genomic resource for interpreting Anopheles sinensis insecticide resistance

BACKGROUND

Anopheles sinensis is an important malaria vector in Southeast Asia. The widespread emergence of insecticide resistance in this mosquito species poses a serious threat to the efficacy of malaria control measures, particularly in China. Recently, the whole-genome sequencing and de novo assembly of An. sinensis (China strain) has been finished. A series of insecticide-resistant studies in An. sinensis have also been reported. There is a growing need to integrate these valuable data to provide a comprehensive database for further studies on insecticide-resistant management of An. sinensis.

RESULTS

A bioinformatics database named An. sinensis genome database (ASGDB) was built. In addition to being a searchable database of published An. sinensis genome sequences and annotation, ASGDB provides in-depth analytical platforms for further understanding of the genomic and genetic data, including visualization of genomic data, orthologous relationship analysis, GO analysis, pathway analysis, expression analysis and resistance-related gene analysis. Moreover, ASGDB provides a panoramic view of insecticide resistance studies in An. sinensis in China. In total, 551 insecticide-resistant phenotypic and genotypic reports on An. sinensis distributed in Chinese malaria-endemic areas since the mid-1980s have been collected, manually edited in the same format and integrated into OpenLayers map-based interface, which allows the international community to assess and exploit the high volume of scattered data much easier. The database has been given the URL: http://www.asgdb.org /.

CONCLUSIONS

ASGDB was built to help users mine data from the genome sequence of An. sinensis easily and effectively, especially with its advantages in insecticide resistance surveillance and control.

Sequencing and analysis of the complete mitochondrial genome in Anopheles sinensis (Diptera: Culicidae)

Background

Anopheles sinensis (Diptera: Culicidae) is a primary vector of Plasmodium vivax and Brugia malayi in most regions of China. In addition, its phylogenetic relationship with the cryptic species of the Hyrcanus Group is complex and remains unresolved. Mitochondrial genome sequences are widely used as molecular markers for phylogenetic studies of mosquito species complexes, of which mitochondrial genome data of An. sinensis is not available.

Methods

An. sinensis samples was collected from Shandong, China, and identified by molecular marker. Genomic DNA was extracted, followed by the Illumina sequencing. Two complete mitochondrial genomes were assembled and annotated using the mitochondrial genome of An. gambiae as reference. The mitochondrial genomes sequences of the 28 known Anopheles species were aligned and reconstructed phylogenetic tree by Maximum Likelihood (ML) method.

Findings

The length of complete mitochondrial genomes of An. sinensis was 15,076 bp and 15,138 bp, consisting of 13 protein-coding genes, 22 transfer RNA (tRNA) genes, 2 ribosomal RNA (rRNA) genes, and an AT-rich control region. As in other insects, most mitochondrial genes are encoded on the J strand, except for ND5, ND4, ND4L, ND1, two rRNA and eight tRNA genes, which are encoded on the N strand. The bootstrap value was set as 1000 in ML analyses. The topologies restored phylogenetic affinity within subfamily Anophelinae. The ML tree showed four major clades, corresponding to the subgenera Cellia, Anopheles, Nyssorhynchus and Kerteszia of the genus Anopheles.

Conclusions

The complete mitochondrial genomes of An. sinensis were obtained. The number, order and transcription direction of An. sinensis mitochondrial genes were the same as in other species of family Culicidae.

Electronic supplementary material

The online version of this article (10.1186/s40249-017-0362-7) contains supplementary material, which is available to authorized users.

Biology, Bionomics and Molecular Biology of Anopheles sinensis Wiedemann 1828 (Diptera: Culicidae), Main Malaria Vector in China

China has set a goal to eliminate all malaria in the country by 2020, but it is unclear if current understanding of malaria vectors and transmission is sufficient to achieve this objective. Anopheles sinensis is the most widespread malaria vector specie in China, which is also responsible for vivax malaria outbreak in central China. We reviewed literature from 1954 to 2016 on An. sinensis with emphasis on biology, bionomics, and molecular biology. A total of 538 references were relevant and included. An. sienesis occurs in 29 Chinese provinces. Temperature can affect most life-history parameters. Most An. sinensis are zoophilic, but sometimes they are facultatively anthropophilic. Sporozoite analysis demonstrated An. sinensis efficacy on Plasmodium vivax transmission. An. sinensis was not stringently refractory to P. falciparum under experimental conditions, however, sporozoite was not found in salivary glands of field collected An. sinensis. The literature on An. sienesis biology and bionomics was abundant, but molecular studies, such as gene functions and mechanisms, were limited. Only 12 molecules (genes, proteins or enzymes) have been studied. In addition, there were considerable untapped omics resources for potential vector control tools. Existing information on An. sienesis could serve as a baseline for advanced research on biology, bionomics and genetics relevant to vector control strategies.

Complete mitochondrial genomes of Anopheles stephensi and An. dirus and comparative evolutionary mitochondriomics of 50 mosquitoes

To better understand the phylogeny and evolution of mosquitoes, the complete mitochondrial genome (mitogenome) of Anopheles stephensi and An. dirus were sequenced and annotated, and a total of 50 mosquito mitogenomes were comparatively analyzed. The complete mitogenome of An. stephensi and An. dirus is 1,5371 bp and 1,5406 bp long, respectively. The main features of the 50 mosquito mitogenomes are conservative: 13 protein-coding genes (PCGs), two ribosomal RNA genes, 22 transfer RNA genes, positive AT-skew and negative GC-skew. The gene order trnA-trnR in ancestral insects is rearranged. All tRNA genes have the typical clover leaf secondary structure but tRNA Ser . The control regions are highly variable in size. PCGs show signals of purifying selection, but evidence for positive selection in ND2, ND4 and ND6 is found. Bayesian and Maximum Likelihood phylogenetic analyses based on all PCG nucleotides produce an identical tree topology and strongly support the monophyly of subgenera Cellia, Anopheles, Keterszia and Nyssorhynchus, the sister relationship of the subgenera Nyssorhynchus and Keterszia, and Cellia and Anopheles. The most recent ancestor of the genus Anopheles and Culicini + Aedini exited ~145 Mya ago. This is the first comprehensive study of mosquito mitogenomes, which are effective for mosquito phylogeny at various taxonomic levels.

The antenna transcriptome changes in mosquito Anopheles sinensis, pre- and post- blood meal

Antenna is the main chemosensory organ in mosquitoes. Characterization of the transcriptional changes after blood meal, especially those related to chemoreception, may help to explain mosquito blood sucking behavior and to identify novel targets for mosquito control. Anopheles sinensis is an Asiatic mosquito species which transmits malaria and lymphatic filariasis. However, studies on chemosensory biology in female An. sinensis are quite lacking. Here we report a transcriptome analysis of An. sinensis female antennae pre- and post- blood meal. We created six An. sinensis antenna RNA-seq libraries, three from females without blood meal and three from females five hours after a blood meal. Illumina sequencing was conducted to analyze the transcriptome differences between the two groups. In total, the sequenced fragments created 21,643 genes, 1,828 of them were novel. 12,861 of these genes were considered to be expressed (FPKM >1.0) in at least one of the two groups, with 12,159 genes expressed in both groups. 548 genes were differentially expressed in the blood-fed group, with 331 genes up-regulated and 217 genes down-regulated. GO enrichment analysis of the differentially expressed genes suggested that there were no statistically over represented GO terms among down-regulated genes in blood-fed mosquitoes, while the enriched GO terms of the up-regulated genes occurred mainly in metabolic process. For the chemosensory gene families, a subtle distinction in the expression levels can be observed according to our statistical analysis. However, the firstly comprehensive identification of these chemosensory gene families in An. sinensis antennae will help to characterize the precise function of these proteins in odor recognition in mosquitoes. This study provides a first global view in the changes of transcript accumulation elicited by blood meal in An. sinensis female antennae.

Comparative physical genome mapping of malaria vectors Anopheles sinensis and Anopheles gambiae

Background

Anopheles sinensis is a dominant natural vector of Plasmodium vivax in China, Taiwan, Japan, and Korea. Recent genome sequencing of An. sinensis provides important insights into the genomic basis of vectorial capacity. However, the lack of a physical genome map with chromosome assignment and orientation of sequencing scaffolds hinders comparative analyses with other genomes to infer evolutionary changes relevant to the vector capacity.

Results

Here, a physical genome map for An. sinensis was constructed by assigning 52 scaffolds onto the chromosomes using fluorescence in situ hybridization (FISH). This chromosome-based genome assembly composes approximately 36% of the total An. sinensis genome. Comparisons of 3955 orthologous genes between An. sinensis and Anopheles gambiae identified 361 conserved synteny blocks and 267 inversions fixed between these two lineages. The rate of gene order reshuffling on the X chromosome is approximately 3.2 times higher than that on the autosomes.

Conclusions

The physical map will facilitate detailed genomic analysis of An. sinensis and contribute to understanding of the patterns and mechanisms of large-scale genome rearrangements in anopheline mosquitoes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12936-017-1888-7) contains supplementary material, which is available to authorized users.

The B-cell translocation gene 1 (CgBTG1) identified in oyster Crassostrea gigas exhibit multiple functions in immune response

B-cell translocation gene 1 (BTG1) is a member of the anti-proliferative gene family, which plays important roles in regulation of cell cycle. In the present study, a B-cell translocation gene 1 molecule homologue (designed CgBTG1) are identified and characterized in oyster Crassostrea gigas. CgBTG1 contains a conserved BTG domain with Box A and Box B motifs, and it shares high similarities with both BTG1 and BTG2 proteins in vertebrates. CgBTG1 mRNA is predominantly expressed in hemocytes, and its expression level in hemocytes is significantly up-regulated at 6 h (5.40-fold, p < 0.01) post Vibrio splendidus stimulation. The apoptosis rate of oyster hemocytes is significantly decreased (p < 0.05) after CgBTG1 interfered by dsRNA (dsCgBTG1). This is indicated that CgBTG1 participated in the regulation of oyster hemocytes apoptosis. Furthermore, CgBTG1 could also induce the apoptosis of cancer cells (HeLa, A549 and BEL7402) in vitro. Compared with normal oysters, both vessel-like structures and muscle fibers in CgBTG1 interfered oysters are severely damaged after V. splendidus challenge in paraffin section, considering that CgBTG1 possessed an analogous feature of angiogenesis for maintenance of vessel-like structures in adductor muscle of oyster. The results suggests that CgBTG1 is a multi-functional molecule involved in the immune response of C. gigas against pathogen infection, which provides more clues for intensive studies of BTG family proteins in invertebrates.

Discovery of flavivirus-derived endogenous viral elements in Anopheles mosquito genomes supports the existence of Anopheles-associated insect-specific flaviviruses

The Flavivirus genus encompasses several arboviruses of public health significance such as dengue, yellow fever, and Zika viruses. It also includes insect-specific flaviviruses (ISFs) that are only capable of infecting insect hosts. The vast majority of mosquito-infecting flaviviruses have been associated with mosquito species of the Aedes and Culex genera in the Culicinae subfamily, which also includes most arbovirus vectors. Mosquitoes of the Anophelinae subfamily are not considered significant arbovirus vectors; however, flaviviruses have occasionally been detected in field-caught Anopheles specimens. Whether such observations reflect occasional spillover or laboratory contamination or whether Anopheles mosquitoes are natural hosts of flaviviruses is unknown. Here, we provide in silico and in vivo evidence of transcriptionally active, flavivirus-derived endogenous viral elements (EVEs) in the genome of Anopheles minimus and Anopheles sinensis. Such non-retroviral endogenization of RNA viruses is consistent with a shared evolutionary history between flaviviruses and Anopheles mosquitoes. Phylogenetic analyses of the two newly described EVEs support the existence of a distinct clade of Anopheles-associated ISFs.

Imbalanced positive selection maintains the functional divergence of duplicated DIHYDROKAEMPFEROL 4-REDUCTASE genes

Gene duplication could be beneficial by functional division but might increase the risk of genetic load. The dynamics of duplicated paralogs number could involve recombination, positive selection, and functional divergence. Duplication of DIHYDROFLAVONOL 4-REDUCTASE (DFR) has been reported in several organisms and may have been retained by escape from adaptive conflict (EAC). In this study, we screened the angiosperm DFR gene focusing on a diversified genus Scutellaria to investigate how these duplicated genes are retained. We deduced that gene duplication involved multiple independent events in angiosperms, but the duplication of DFR was before the divergence of Scutellaria. Asymmetric positive selective pressures resulted in different evolutionary rates between the duplicates. Different numbers of regulatory elements, differential codon usages, radical amino acid changes, and differential gene expressions provide evidences of functional divergence between the two DFR duplicates in Scutellaria, implying adaptive subfunctionalization between duplicates. The discovery of pseudogenes accompanying a reduced replacement rate in one DFR paralogous gene suggested possibly leading to “loss of function” due to dosage imbalance after the transient adaptive subfunctionalization in the early stage of duplication. Notwithstanding, episodic gene duplication and functional divergence may be relevant to the diversification of ecological function of DFR gene in Scutellaria.

Draft Genomes of Anopheles cracens and Anopheles maculatus: Comparison of Simian Malaria and Human Malaria Vectors in Peninsular Malaysia

Anopheles cracens has been incriminated as the vector of human knowlesi malaria in peninsular Malaysia. Besides, it is a good laboratory vector of Plasmodium falciparum and P. vivax. The distribution of An. cracens overlaps with that of An. maculatus, the human malaria vector in peninsular Malaysia that seems to be refractory to P. knowlesi infection in natural settings. Whole genome sequencing was performed on An. cracens and An. maculatus collected here. The draft genome of An. cracens was 395 Mb in size whereas the size of An. maculatus draft genome was 499 Mb. Comparison with the published Malaysian An. maculatus genome suggested the An. maculatus specimen used in this study as a different geographical race. Comparative analyses highlighted the similarities and differences between An. cracens and An. maculatus, providing new insights into their biological behavior and characteristics.

An Inhibitor of the Alternative Pathway of Complement in Saliva of New World Anopheline Mosquitoes

The complement system present in circulating blood is an effective mechanism of host defense, responsible for the killing of pathogens and the production of potent anaphylatoxins. Inhibitors of the complement system have been described in the saliva of hematophagous arthropods that are involved in the protection of digestive tissues against complement system-mediated damage. In this study, we describe albicin, a novel inhibitor of the alternative pathway of complement from the salivary glands of the malaria vector, Anopheles albimanus The inhibitor was purified from salivary gland homogenates by reverse-phase HPLC and identified by mass spectrometry as a small (13.4-kDa) protein related to the gSG7 protein of Anopheles gambiae and Anopheles stephensi Recombinant albicin was produced in Escherichia coli and found to potently inhibit lysis of rabbit erythrocytes in assays of the alternative pathway while having no inhibitory effect on the classical or lectin pathways. Albicin also inhibited the deposition of complement components on agarose-coated plates, although it could not remove previously bound components. Antisera produced against recombinant albicin recognized both the native and recombinant inhibitors and also blocked their activities in in vitro assays. Using surface plasmon resonance and enzymatic assays, we found that albicin binds and stabilizes the C3-convertase complex (C3bBb) formed on a properdin surface and inhibits the convertase activity of a reconstituted C3bBb complex in solution. The data indicate that albicin specifically recognizes the activated form of the complex, allowing more efficient inhibition by an inhibitor whose quantity is limited.

More than one rabbit out of the hat: Radiation, transgenic and symbiont-based approaches for sustainable management of mosquito and tsetse fly populations

Mosquitoes (Diptera: Culicidae) and tsetse flies (Diptera: Glossinidae) are bloodsucking vectors of human and animal pathogens. Mosquito-borne diseases (malaria, filariasis, dengue, zika, and chikungunya) cause severe mortality and morbidity annually, and tsetse fly-borne diseases (African trypanosomes causing sleeping sickness in humans and nagana in livestock) cost Sub-Saharan Africa an estimated US$ 4750 million annually. Current reliance on insecticides for vector control is unsustainable: due to increasing insecticide resistance and growing concerns about health and environmental impacts of chemical control there is a growing need for novel, effective and safe biologically-based methods that are more sustainable. The integration of the sterile insect technique has proven successful to manage crop pests and disease vectors, particularly tsetse flies, and is likely to prove effective against mosquito vectors, particularly once sex-separation methods are improved. Transgenic and symbiont-based approaches are in development, and more advanced in (particularly Aedes) mosquitoes than in tsetse flies; however, issues around stability, sustainability and biosecurity have to be addressed, especially when considering population replacement approaches. Regulatory issues and those relating to intellectual property and economic cost of application must also be overcome. Standardised methods to assess insect quality are required to compare and predict efficacy of the different approaches. Different combinations of these three approaches could be integrated to maximise their benefits, and all have the potential to be used in tsetse and mosquito area-wide integrated pest management programmes.

Genomic Analysis of Detoxification Supergene Families in the Mosquito Anopheles sinensis

Anopheles sinensis is an important malaria vector in China and other Southeast Asian countries, and the emergence of insecticide resistance in this mosquito poses a serious threat to the efficacy of malaria control programs. The recently published An. sinensis genome and transcriptome provide an opportunity to understand the molecular mechanisms of insecticide resistance. Analysis of the An. sinensis genome revealed 174 detoxification genes, including 93 cytochrome P450s (P450s), 31 glutathione-S-transferases (GSTs), and 50 choline/carboxylesterases (CCEs). The gene number was similar to that in An. gambiae, but represented a decrease of 29% and 42% compared with Aedes aegypti and Culex quinquefasciatus, respectively. The considerable contraction in gene number in Anopheles mosquitoes mainly occurred in two detoxification supergene families, P450s and CCEs. The available An. sinensis transcriptome was also re-analyzed to further identify key resistance-associated detoxification genes. Among 174 detoxification genes, 124 (71%) were detected. Several candidate genes overexpressed in a deltamethrin-resistant strain (DR-strain) were identified as belonging to the CYP4 or CYP6 family of P450s and the Delta GST class. These generated data provide a basis for identifying the resistance-associated genes of An. sinensis at the molecular level.

A call for benchmarking transposable element annotation methods

DNA derived from transposable elements (TEs) constitutes large parts of the genomes of complex eukaryotes, with major impacts not only on genomic research but also on how organisms evolve and function. Although a variety of methods and tools have been developed to detect and annotate TEs, there are as yet no standard benchmarks-that is, no standard way to measure or compare their accuracy. This lack of accuracy assessment calls into question conclusions from a wide range of research that depends explicitly or implicitly on TE annotation. In the absence of standard benchmarks, toolmakers are impeded in improving their tools, annotators cannot properly assess which tools might best suit their needs, and downstream researchers cannot judge how accuracy limitations might impact their studies. We therefore propose that the TE research community create and adopt standard TE annotation benchmarks, and we call for other researchers to join the authors in making this long-overdue effort a success.

Heterochromatin, histone modifications, and nuclear architecture in disease vectors

Interactions between a pathogen and a vector are plastic and dynamic. Such interactions can be more rapidly accommodated by epigenetic changes than by genetic mutations. Gene expression can be affected by the proximity to the heterochromatin, by local histone modifications, and by the three-dimensional position within the nucleus. Recent studies of disease vectors indicate that gene regulation by these factors can be important for susceptibility to pathogens, reproduction, immunity, development, and longevity. Knowledge about heterochromatin, histone modifications, and nuclear architecture will help our understanding of epigenetic mechanisms that control gene function at traits related to vectorial capacity.

Horizontal transmission of an R4 clade non-long terminal repeat retrotransposon between the divergent Aedes and Anopheles mosquito genera

AaegR4_1 and AgamR4_1 are the sole R4 clade non-long terminal repeat (non-LTR) retrotransposons in Aedes aegypti and Anopheles gambiae, two species that diverged approximately 145-200 million years ago. Twelve full-length copies were found in Ae. aegypti and have less than 1% nucleotide (nt) divergence, suggesting recent activity on an evolutionary time scale. Five of these copies have intact open reading frames and the 3.6 kb open reading frame of AaegR4_1.1 has 78% nt identity to AgamR4_1.1. No intact copies were found in An. gambiae. Searches of 25 genomic databases for 22 mosquito species from three genera revealed R4 clade representatives in Aedes and Anopheles genera but not in Culex. Interestingly, these elements are present in all six species of the An. gambiae species complex that were searched but not in 13 other anopheline species. These results combined with divergence vs. age analysis suggest that horizontal transfer is the most likely explanation for the low divergence between R4 clade retrotransposon sequences of the divergent mosquito species from the Aedes and Anopheles genera. This is the first report of the horizontal transfer of an R4 clade non-LTR retrotransposon and the first report of the horizontal transfer of a non-LTR retrotransposon in mosquitoes.

De novo assembly and annotation of the Asian tiger mosquito (Aedes albopictus) repeatome with dnaPipeTE from raw genomic reads and comparative analysis with the yellow fever mosquito (Aedes aegypti)

Repetitive DNA, including transposable elements (TEs), is found throughout eukaryotic genomes. Annotating and assembling the “repeatome” during genome-wide analysis often poses a challenge. To address this problem, we present dnaPipeTE—a new bioinformatics pipeline that uses a sample of raw genomic reads. It produces precise estimates of repeated DNA content and TE consensus sequences, as well as the relative ages of TE families. We shows that dnaPipeTE performs well using very low coverage sequencing in different genomes, losing accuracy only with old TE families. We applied this pipeline to the genome of the Asian tiger mosquito Aedes albopictus, an invasive species of human health interest, for which the genome size is estimated to be over 1 Gbp. Using dnaPipeTE, we showed that this species harbors a large (50% of the genome) and potentially active repeatome with an overall TE class and order composition similar to that of Aedes aegypti, the yellow fever mosquito. However, intraorder dynamics show clear distinctions between the two species, with differences at the TE family level. Our pipeline’s ability to manage the repeatome annotation problem will make it helpful for new or ongoing assembly projects, and our results will benefit future genomic studies of A. albopictus.

Transposable elements as agents of rapid adaptation may explain the genetic paradox of invasive species

Rapid adaptation of invasive species to novel habitats has puzzled evolutionary biologists for decades, especially as this often occurs in the face of limited genetic variability. Although some ecological traits common to invasive species have been identified, little is known about the possible genomic/genetic mechanisms that may underlie their success. A common scenario in many introductions is that small founder population sizes will often lead to reduced genetic diversity, but that invading populations experience large environmental perturbations, such as changes in habitat and environmental stress. Although sudden and intense stress is usually considered in a negative context, these perturbations may actually facilitate rapid adaptation by affecting genome structure, organization and function via interactions with transposable elements (TEs), especially in populations with low genetic diversity. Stress-induced changes in TE activity can alter gene action and can promote structural variation that may facilitate the rapid adaptation observed in new environments. We focus here on the adaptive potential of TEs in relation to invasive species and highlight their role as powerful mutational forces that can rapidly create genetic diversity. We hypothesize that activity of transposable elements can explain rapid adaptation despite low genetic variation (the genetic paradox of invasive species), and provide a framework under which this hypothesis can be tested using recently developed and emerging genomic technologies.

Genome mining offers a new starting point for parasitology research

Parasites including helminthes, protozoa, and medical arthropod vectors are a major cause of global infectious diseases, affecting one-sixth of the world's population, which are responsible for enormous levels of morbidity and mortality important and remain impediments to economic development especially in tropical countries. Prevalent drug resistance, lack of highly effective and practical vaccines, as well as specific and sensitive diagnostic markers are proving to be challenging problems in parasitic disease control in most parts of the world. The impressive progress recently made in genome-wide analysis of parasites of medical importance, including trematodes of Clonorchis sinensis, Opisthorchis viverrini, Schistosoma haematobium, S. japonicum, and S. mansoni; nematodes of Brugia malayi, Loa loa, Necator americanus, Trichinella spiralis, and Trichuris suis; cestodes of Echinococcus granulosus, E. multilocularis, and Taenia solium; protozoa of Babesia bovis, B. microti, Cryptosporidium hominis, Eimeria falciformis, E. histolytica, Giardia intestinalis, Leishmania braziliensis, L. donovani, L. major, Plasmodium falciparum, P. vivax, Trichomonas vaginalis, Trypanosoma brucei and T. cruzi; and medical arthropod vectors of Aedes aegypti, Anopheles darlingi, A. sinensis, and Culex quinquefasciatus, have been systematically covered in this review for a comprehensive understanding of the genetic information contained in nuclear, mitochondrial, kinetoplast, plastid, or endosymbiotic bacterial genomes of parasites, further valuable insight into parasite-host interactions and development of promising novel drug and vaccine candidates and preferable diagnostic tools, thereby underpinning the prevention and control of parasitic diseases.

Mosquito genomics. Highly evolvable malaria vectors: the genomes of 16 Anopheles mosquitoes

Variation in vectorial capacity for human malaria among Anopheles mosquito species is determined by many factors, including behavior, immunity, and life history. To investigate the genomic basis of vectorial capacity and explore new avenues for vector control, we sequenced the genomes of 16 anopheline mosquito species from diverse locations spanning ~100 million years of evolution. Comparative analyses show faster rates of gene gain and loss, elevated gene shuffling on the X chromosome, and more intron losses, relative to Drosophila. Some determinants of vectorial capacity, such as chemosensory genes, do not show elevated turnover but instead diversify through protein-sequence changes. This dynamism of anopheline genes and genomes may contribute to their flexible capacity to take advantage of new ecological niches, including adapting to humans as primary hosts.

Genome analysis of a major urban malaria vector mosquito, Anopheles stephensi

Background

Anopheles stephensi is the key vector of malaria throughout the Indian subcontinent and Middle East and an emerging model for molecular and genetic studies of mosquito-parasite interactions. The type form of the species is responsible for the majority of urban malaria transmission across its range.

Results

Here, we report the genome sequence and annotation of the Indian strain of the type form of An. stephensi. The 221 Mb genome assembly represents more than 92% of the entire genome and was produced using a combination of 454, Illumina, and PacBio sequencing. Physical mapping assigned 62% of the genome onto chromosomes, enabling chromosome-based analysis. Comparisons between An. stephensi and An. gambiae reveal that the rate of gene order reshuffling on the X chromosome was three times higher than that on the autosomes. An. stephensi has more heterochromatin in pericentric regions but less repetitive DNA in chromosome arms than An. gambiae. We also identify a number of Y-chromosome contigs and BACs. Interspersed repeats constitute 7.1% of the assembled genome while LTR retrotransposons alone comprise more than 49% of the Y contigs. RNA-seq analyses provide new insights into mosquito innate immunity, development, and sexual dimorphism.

Conclusions

The genome analysis described in this manuscript provides a resource and platform for fundamental and translational research into a major urban malaria vector. Chromosome-based investigations provide unique perspectives on Anopheles chromosome evolution. RNA-seq analysis and studies of immunity genes offer new insights into mosquito biology and mosquito-parasite interactions.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0459-2) contains supplementary material, which is available to authorized users.

De novo transcriptome sequencing and sequence analysis of the malaria vector Anopheles sinensis (Diptera: Culicidae)

Background

Anopheles sinensis is the major malaria vector in China and Southeast Asia. Vector control is one of the most effective measures to prevent malaria transmission. However, there is little transcriptome information available for the malaria vector. To better understand the biological basis of malaria transmission and to develop novel and effective means of vector control, there is a need to build a transcriptome dataset for functional genomics analysis by large-scale RNA sequencing (RNA-seq).

Methods

To provide a more comprehensive and complete transcriptome of An. sinensis, eggs, larvae, pupae, male adults and female adults RNA were pooled together for cDNA preparation, sequenced using the Illumina paired-end sequencing technology and assembled into unigenes. These unigenes were then analyzed in their genome mapping, functional annotation, homology, codon usage bias and simple sequence repeats (SSRs).

Results

Approximately 51.6 million clean reads were obtained, trimmed, and assembled into 38,504 unigenes with an average length of 571 bp, an N50 of 711 bp, and an average GC content 51.26%. Among them, 98.4% of unigenes could be mapped onto the reference genome, and 69% of unigenes could be annotated with known biological functions. Homology analysis identified certain numbers of An. sinensis unigenes that showed homology or being putative 1:1 orthologues with genomes of other Dipteran species. Codon usage bias was analyzed and 1,904 SSRs were detected, which will provide effective molecular markers for the population genetics of this species.

Conclusions

Our data and analysis provide the most comprehensive transcriptomic resource and characteristics currently available for An. sinensis, and will facilitate genetic, genomic studies, and further vector control of An. sinensis.

Transcriptome profiling of pyrethroid resistant and susceptible mosquitoes in the malaria vector, Anopheles sinensis

BACKGROUND

Anopheles sinensis is a major malaria vector in China and other Southeast Asian countries, and it is becoming increasingly resistant to the insecticides used for agriculture, net impregnation, and indoor residual spray. Very limited genomic information on this species is available, which has hindered the development of new tools for resistance surveillance and vector control. We used the 454 GS FLX system and generated expressed sequence tag (EST) databases of various life stages of An. sinensis, and we determined the transcriptional differences between deltamethrin resistant and susceptible mosquitoes.

RESULTS

The 454 GS FLX transcriptome sequencing yielded a total of 624,559 reads (average length of 290 bp) with the pooled An. sinensis mosquitoes across various development stages. The de novo assembly generated 33,411 contigs with average length of 493 bp. A total of 8,057 ESTs were generated with Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) annotation. A total of 2,131 ESTs were differentially expressed between deltamethrin resistant and susceptible mosquitoes collected from the same field site in Jiangsu, China. Among these differentially expressed ESTs, a total of 294 pathways were mapped to the KEGG database, with the predominant ESTs belonging to metabolic pathways. Furthermore, a total of 2,408 microsatellites and 15,496 single nucleotide polymorphisms (SNPs) were identified.

CONCLUSIONS

The annotated EST and transcriptome databases provide a valuable genomic resource for further genetic studies of this important malaria vector species. The differentially expressed ESTs associated with insecticide resistance identified in this study lay an important foundation for further functional analysis. The identified microsatellite and SNP markers will provide useful tools for future population genetic and comparative genomic analyses of malaria vectors.