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

Identification and functional validation of P450 genes associated with pyrethroid resistance in the malaria vector Anopheles sinensis (Diptera Culicidae)

Cytochrome P450 monooxygenases (P450s), a multifunctional protein superfamily, are one of three major classes of detoxification enzymes. However, the diversity and functions of P450 genes from pyrethroid-resistant populations of Anopheles sinensis have not been fully explored. In this study, P450 genes associated with pyrethroid resistance were systematically screened using RNA-seq in three field pyrethroid-resistant populations (AH-FR, CQ-FR, YN-FR) and one laboratory resistant strain (WX-LR) at developmental stages, tissues, and post blood-meal in comparison to the laboratory susceptible strain (WX-LS) in An. sinensis. Importantly, the expression of significantly upregulated P450s was verified using RT-qPCR, and the function of selected P450s in pyrethroid detoxification was determined with RNA interference using four laboratory pyrethroid-resistant strains (WX-LR, AH-LR, CQ-LR, YN-LR). Sixteen P450 genes were significantly upregulated in at least one field-resistant population, and 44 were significantly upregulated in different developmental stages, tissues or post blood-meal. A total of 19 P450s were selected to verify their association with pyrethroid resistance, and four of them (AsCYP6P3v1, AsCYP6P3v2, AsCYP9J10, and AsCYP9K1) demonstrated significant upregulation in laboratory pyrethroid-resistant strains using RT-qPCR. Knockdown of these four genes all significantly reduced pyrethroid resistance and increased the mortality by 57.19% (AsCYP6P3v1 and AsCYP6P3v2 knockdown group), 38.39% (AsCYP9K1 knockdown group) and 48.87% (AsCYP9J10 knockdown group) in An. sinensis by RNAi, which determined the pyrethroid detoxification function of these four genes. This study revealed the diversity of P450 genes and provided functional evidence for four P450s in pyrethroid detoxification in An. sinensis for the first time, which increases our understanding of the pyrethroid resistance mechanism, and is of potential value for pyrethroid resistance detection and surveillance.

Characteristics of Trypsin genes and their roles in insecticide resistance based on omics and functional analyses in the malaria vector Anopheles sinensis

Trypsin is one of the most diverse and widely studied protease hydrolases. However, the diversity and characteristics of the Trypsin superfamily of genes have not been well understood, and their role in insecticide resistance is yet to be investigated. In this study, a total of 342 Trypsin genes were identified and classified into seven families based on homology, characteristic domains and phylogenetics in Anopheles sinensis, and the LY-Domain and CLECT-Domain families are specific to the species. Four Trypsin genes, (Astry2b, Astry43a, Astry90, Astry113c) were identified to be associated with pyrethroid resistance based on transcriptome analyses of three field resistant populations and qRT-PCR validation, and the knock-down of these genes significantly decrease the pyrethroid resistance of Anopheles sinensis based on RNAi. The activity of Astry43a can be reduced by five selected insecticides (indoxacarb, DDT, temephos, imidacloprid and deltamethrin); and however, the Astry43a could not directly metabolize these five insecticides, like the trypsin NYD-Tr did in earlier reports. This study provides the overall information frame of Trypsin genes, and proposes the role of Trypsin genes to insecticide resistance. Further researches are necessary to investigate the metabolism function of these trypsins to insecticides.

Phylogenetic Analysis of Some Species of the Anopheles hyrcanus Group (Diptera: Culicidae) in China Based on Complete Mitochondrial Genomes

Some species of the Hyrcanus group are vectors of malaria in China. However, the member species are difficult to identify accurately by morphology. The development of sequencing technologies offers the possibility of further studies based on the complete mitochondrial genome. In this study, samples of mosquitoes of the Hyrcanus group were collected in China between 1997 and 2015. The mitochondrial genomes of ten species of the Hyrcanus group were analyzed, including the structure and base composition, codon usage, secondary structure of tRNA, and base difference sites in protein coding regions. Phylogenetic analyses using maximum-likelihood and Bayesian inference were performed based on mitochondrial genes and complete mitochondrial genomes The mitochondrial genome of 10 Hyrcanus group members ranged from 15,403 bp to 15,475 bp, with an average 78.23% (A + T) content, comprising of 13 PCGs (protein coding genes), 22 tRNAs, and 2 rRNAs. Site differences between some closely related species in the PCGs were small. There were only 36 variable sites between Anopheles sinensis and Anopheles belenrae for a variation ratio of 0.32% in all PCGs. The pairwise interspecies distance based on 13 PCGs was low, with an average of 0.04. A phylogenetic tree constructed with the 13 PCGs was consistent with the known evolutionary relationships. Some phylogenetic trees constructed by single coding regions (such as COI or ND4) or combined coding regions (COI + ND2 + ND4 + ND5 or ND2 + ND4) were consistent with the phylogenetic tree constructed using the 13 PCGs. The phylogenetic trees constructed using some coding genes (COII, ND5, tRNAs, 12S rRNA, and 16S rRNA) differed from the phylogenetic tree constructed using PCGs. The difference in mitochondrial genome sequences between An. sinensis and An. belenrae was very small, corresponding to intraspecies difference, suggesting that the species was in the process of differentiation. The combination of all 13 PCG sequences was demonstrated to be optimal for phylogenetic analysis in closely related species.

Glutathione S-transferase (GST) genes and their function associated with pyrethroid resistance in the malaria vector Anopheles sinensis

BACKGROUND

Glutathione S-transferases (GSTs), a multifunctional protein family, are involved in insecticide resistance. However, a systematic analysis of GSTs in Anopheles sinensis, an important vector for malaria transmission, is lacking. In this study, we investigated the diversity and characteristics of GST genes, and analyzed their expression patterns and functions associated with insecticide resistance in this species.

RESULTS

We identified 32 putative cytosolic and three putative microsomal GST genes in the An. sinensis genome. Transcriptome analysis showed that GSTs were highly expressed in larvae, and mainly expressed in the antennae, midgut and Malpighian tubules of adults. In addition, we found that GSTd2 and GSTe2 were significantly upregulated in four An. sinensis pyrethroid-resistant field populations. Furthermore, silencing of GSTd2 and GSTe2 significantly increased the susceptibility of An. sinensis to deltamethrin, and recombinant GSTd2 and GSTe2 exhibited high enzymatic activity in the metabolism of 1-chloro-2,4-dinitrobenzene and dichlorodiphenyltrichloroethane (DDT).

CONCLUSION

These results showed that GSTs are involved in the development of insecticide resistance in An. sinensis through transcriptional overexpression and enzymatic metabolization, facilitating our understanding of insecticide resistance in insects. © 2022 Society of Chemical Industry.

Proteomic analysis suggests that monoterpenes in lemongrass disrupt Ca2+ homeostasis in Haemaphysalis longicornis leading to mitochondrial depolarization and cytotoxicity

Complex mixtures of bioactive ingredients in plant essential oils present complex chemistries which involve different modes of action. An increasing body of scientific reports has recently focused on the acaricidal activities of plant essential oils attributed to their monoterpene components, but information about their underlying molecular mechanism of action is scarce. Here, after the chemical analysis of lemongrass oil, a proteomic analysis of the ovary, salivary gland, and midgut of Haemaphysalis longicornis exposed to Cymbopogon citratus (lemongrass) essential oil was performed via data-independent acquisition mass spectrometry (DIA-MS) technology to further elucidate the molecular mechanisms involved. Pathway analysis reveals the activation of metabolic pathways mediated by oxidoreductases and transferases. Furthermore, the upregulation of various calcium-associated proteins and the upregulation of cytochrome c1, cytochrome c oxidase polypeptide IV, and programmed cell death protein 6-like isoform X1 suggest a cytotoxic mode of action via the formation of reactive oxygen species (ROS), mitochondrial Ca²⁺ overload, mitochondrial uncoupling, and depolarization, and ATP depletion leading to either apoptotic or necrotic death. Morphological alterations observed after the RNAi of a major detoxification enzyme (glutathione S-transferase) merit further investigation. Hence, the cytotoxic mode of action exhibited by C. citratus oil could be vital for the development of eco-friendly acaricide.

Mitochondrial genes associated with pyrethroid resistance revealed by mitochondrial genome and transcriptome analyses in the malaria vector Anopheles sinensis (Diptera: Culicidae)

BACKGROUND

Insecticides are still the main method of mosquito control, but mosquito resistance presents a large obstacle. The function of mitochondrial genes in the evolution of insecticide resistance is still poorly understood. Pyrethroid is the most commonly used insecticide, and Anopheles sinensis is an important malaria vector in China and Southeast Asia. In this study, we investigated the mitochondrial genes associated with pyrethroid resistance through their genetic and expression variation based on analyses of transcriptomes and 36 individuals with resequencing in three geographical populations in China.

RESULTS

The nucleotide diversity (Pi) in 18 resistant individuals was much lower than that in 18 susceptible individuals, which suggests that some sites experienced purifying selection subject to pyrethroid stress. Ka/Ks and amino acid analyses showed that ND4 experienced positive selection and had 23 amino acid mutations due to pyrethroid stress. These mutations might change the ND4 structure and function and thus alter the efficiency of the respiratory chain. ND5 was significantly upregulated, and ATP8 was significantly downregulated in these three pyrethroid resistant populations, which suggests that these two genes function in the production and maintenance of pyrethroid resistance. There are differences in mitochondrial genes involved in pyrethroid resistance among these three populations.

CONCLUSION

This is the first study to reveal the association of mitochondrial genes in the evolution of insecticide resistance through amino acid mutation and expression patterns and can help us further understand insecticide resistance mechanisms. © 2019 Society of Chemical Industry.

Comparative analyses of simple sequence repeats (SSRs) in 23 mosquito species genomes: Identification, characterization and distribution (Diptera: Culicidae)

Simple sequence repeats (SSRs) exist in both eukaryotic and prokaryotic genomes and are the most popular genetic markers, but the SSRs of mosquito genomes are still not well understood. In this study, we identified and analyzed the SSRs in 23 mosquito species using Drosophila melanogaster as reference at the whole-genome level. The results show that SSR numbers (33 076-560 175/genome) and genome sizes (574.57-1342.21 Mb) are significantly positively correlated (R² = 0.8992, P < 0.01), but the correlation in individual species varies in these mosquito species. In six types of SSR, mono- to trinucleotide SSRs are dominant with cumulative percentages of 95.14%-99.00% and densities of 195.65/Mb-787.51/Mb, whereas tetra- to hexanucleotide SSRs are rare with 1.12%-4.22% and 3.76/Mb-40.23/Mb. The (A/T)n, (AC/GT)n and (AGC/GCT)n are the most frequent motifs in mononucleotide, dinucleotide and trinucleotide SSRs, respectively, and the motif frequencies of tetra- to hexanucleotide SSRs appear to be species-specific. The 10-20 bp length of SSRs are dominant with the number of 110 561 ± 93 482 and the frequency of 87.25% ± 5.73% on average, and the number and frequency decline with the increase of length. Most SSRs (83.34% ± 7.72%) are located in intergenic regions, followed by intron regions (11.59% ± 5.59%), exon regions (3.74% ± 1.95%), and untranslated regions (1.32% ± 1.39%). The mono-, di- and trinucleotide SSRs are the main SSRs in both gene regions (98.55% ± 0.85%) and exon regions (99.27% ± 0.52%). An average of 42.52% of total genes contains SSRs, and the preference for SSR occurrence in different gene subcategories are species-specific. The study provides useful insights into the SSR diversity, characteristics and distribution in 23 mosquito species of genomes.

HSP superfamily of genes in the malaria vector Anopheles sinensis: diversity, phylogenetics and association with pyrethroid resistance

BACKGROUND

Heat shock proteins (HSPs) are molecular chaperones that are involved in many normal cellular processes and various kinds of environmental stress. There is still no report regarding the diversity and phylogenetics research of HSP superfamily of genes at whole genome level in insects, and the HSP gene association with pyrethroid resistance is also not well known. The present study investigated the diversity, classification, scaffold location, characteristics, and phylogenetics of the superfamily of genes in Anopheles sinensis genome, and the HSP genes associated with pyrethroid resistance.

METHODS

The present study identified the HSP genes in the An. sinensis genome, analysed their characteristics, and deduced phylogenetic relationships of all HSPs in An. sinensis, Anopheles gambiae, Culex quinquefasciatus and Aedes aegypti by bioinformatic methods. Importantly, the present study screened the HSPs associated with pyrethroid resistance using three field pyrethroid-resistant populations with RNA-seq and RT-qPCR, and looked over the HSP gene expression pattern for the first time in An. sinensis on the time-scale post insecticide treatment with RT-qPCR.

RESULTS

There are 72 HSP genes in An. sinensis genome, and they are classified into five families and 11 subfamilies based on their molecular weight, homology and phylogenetics. Both RNA-seq and qPCR analysis revealed that the expression of AsHSP90AB, AsHSP70-2 and AsHSP21.7 are significantly upregulated in at least one field pyrethroid-resistant population. Eleven genes are significantly upregulated in different period after pyrethroid exposure. The HSP90, sHSP and HSP70 families are proposed to be involved in pyrethroid stress response based in expression analyses of three field pyrethroid-resistant populations, and expression pattern on the time scale post insecticide treatment. The AsHSP90AB gene is proposed to be the essential HSP gene for pyrethroid stress response in An. sinensis.

CONCLUSIONS

This study provides the information frame for HSP superfamily of genes, and lays an important basis for the better understanding and further research of HSP function in insect adaptability to diverse environments.

ATP-Binding Cassette (ABC) Transporter Genes Involved in Pyrethroid Resistance in the Malaria Vector Anopheles sinensis: Genome-Wide Identification, Characteristics, Phylogenetics, and Expression Profile

background: The ATP-binding cassette (ABC) transporters family is one of the largest families of membrane proteins existing in all living organisms. Pyrethroid resistance has become the largest unique obstacle for mosquito control worldwide. ABC transporters are thought to be associated with pyrethroid resistance in some agricultural pests, but little information is known for mosquitoes. Herein, we investigated the diversity, location, characteristics, phylogenetics, and evolution of ABC transporter family of genes in the Anopheles sinensis genome, and identified the ABC transporter genes associated with pyrethroid resistance through expression profiles using RNA-seq and qPCR. Results: 61 ABC transporter genes are identified and divided into eight subfamilies (ABCA-H), located on 22 different scaffolds. Phylogenetic and evolution analyses with ABC transporters of A. gambiae, Drosophila melanogaster, and Homo sapiens suggest that the ABCD, ABCG, and ABCH subfamilies are monophyly, and that the ABCC and ABCG subfamilies have experienced a gene duplication event. Both RNA-seq and qPCR analyses show that the AsABCG28 gene is uniquely significantly upregulated gene in all three field pyrethroid-resistant populations (Anhui, Chongqing, and Yunnan provinces) in comparison with a laboratory-susceptible strain from Jiangsu province. The AsABCG28 is significantly upregulated at 12-h and 24-h after deltamethrin exposure in three-day-old female adults. Conclusion: This study provides the information frame for ABC transporter subfamily of genes, and lays an important basis for the better understanding and further research of ABC transporter function in insecticide toxification. The AsABCG28 gene is associated with pyrethroid detoxification, and it functions at later period in the detoxification process for xenobiotics transportation.

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

BACKGROUND

UDP-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.

METHODS

The 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.

RESULTS

There 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.

CONCLUSIONS

This 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.

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.

Complete mitogenome of Anopheles sinensis and mitochondrial insertion segments in the nuclear genomes of 19 mosquito species

Anopheles sinensis is a major malarial vector in China and Southeast Asia. The mitochondria is involved in many important biological functions. Nuclear mitochondrial DNA segments (NUMTs) are common in eukaryotic organisms, but their characteristics are poorly understood. We sequenced and analyzed the complete mitochondrial (mt) genome of An. sinensis. The mt genome is 15,418 bp long and contains 13 protein-coding genes (PCGs), two rRNAs, 22 tRNAs and a large non-coding region. Its gene arrangement is similar to previously published mosquito mt genomes. We identified and analyzed the NUMTs of 19 mosquito species with both nuclear genomes and mt genome sequences. The number, total length and density of NUMTs are significantly correlated with genome size. About half of NUMTs are short (< 200 bp), but larger genomes can house longer NUMTs. NUMTs may help explain genome size expansion in mosquitoes. The expansion due to mitochondrial insertion segments is variable in different insect groups. PCGs are transferred to nuclear genomes at a higher frequency in mosquitoes, but NUMT origination is more different than in mammals. Larger-sized nuclear genomes have longer mt genome sequences in both mosquitoes and mammals. The study provides a foundation for the functional research of mitochondrial genes in An. sinensis and helps us understand the characteristics and origin of NUMTs and the potential contribution to genome expansion.

Genome-wide and expression-profiling analyses suggest the main cytochrome P450 genes related to pyrethroid resistance in the malaria vector, Anopheles sinensis (Diptera Culicidae)

BACKGROUND

Anopheles sinensis is one of the major malaria vectors. However, pyrethroid resistance in An. sinensis is threatening malaria control. Cytochrome P450-mediated detoxification is an important pyrethroid resistance mechanism that has been unexplored in An. sinensis. In this study, we performed a comprehensive analysis of the An. sinensis P450 gene superfamily with special attention to their role in pyrethroid resistance using bioinformatics and molecular approaches.

RESULTS

Our data revealed the presence of 112 individual P450 genes in An. sinensis, which were classified into four major clans (mitochondrial, CYP2, CYP3 and CYP4), 18 families and 50 subfamilies. Sixty-seven genes formed nine gene clusters, and genes within the same cluster and the same gene family had a similar gene structure. Phylogenetic analysis showed that most of An. sinensis P450s (82/112) had very close 1: 1 orthology with Anopheles gambiae P450s. Five genes (AsCYP6Z2, AsCYP6P3v1, AsCYP6P3v2, AsCYP9J5 and AsCYP306A1) were significantly upregulated in three pyrethroid-resistant populations in both RNA-seq and RT-qPCR analyses, suggesting that they could be the most important P450 genes involved in pyrethroid resistance in An. sinensis.

CONCLUSION

Our study provides insight on the diversity of An. sinensis P450 superfamily and basis for further elucidating pyrethroid resistance mechanism in this mosquito species. © 2018 Society of Chemical Industry.

Subcellular stoichiogenomics reveal cell evolution and electrostatic interaction mechanisms in cytoskeleton

Background

Eukaryotic cells contain a huge variety of internally specialized subcellular compartments. Stoichiogenomics aims to reveal patterns of elements usage in biological macromolecules. However, the stoichiogenomic characteristics and how they adapt to various subcellular microenvironments are still unknown.

Results

Here we first updated the definition of stoichiogenomics. Then we applied it to subcellular research, and detected distinctive nitrogen content of nuclear and hydrogen, sulfur content of extracellular proteomes. Specially, we found that acidic amino acids (AAs) content of cytoskeletal proteins is the highest. The increased charged AAs are mainly caused by the eukaryotic originated cytoskeletal proteins. Functional subdivision of the cytoskeleton showed that activation, binding/association, and complexes are the three largest functional categories. Electrostatic interaction analysis showed an increased electrostatic interaction between both primary sequences and PPI interfaces of 3D structures, in the cytoskeleton.

Conclusions

This study creates a blueprint of subcellular stoichiogenomic characteristics, and explains that charged AAs of the cytoskeleton increased greatly in evolution, which offer material basis for the eukaryotic cytoskeletal proteins to act in two ways of electrostatic interactions, and further perform their activation, binding/association and complex formation.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4845-0) contains supplementary material, which is available to authorized users.

De novo RNA-seq and functional annotation of Ornithonyssus bacoti

Ornithonyssus bacoti (Hirst) (Acari: Macronyssidae) is a vector and reservoir of pathogens causing serious infectious diseases, such as epidemic hemorrhagic fever, endemic typhus, tularemia, and leptospirosis. Its genome and transcriptome data are lacking in public databases. In this study, total RNA was extracted from live O. bacoti to conduct RNA-seq, functional annotation, coding domain sequence (CDS) prediction and simple sequence repeats (SSRs) detection. The results showed that 65.8 million clean reads were generated and assembled into 72,185 unigenes, of which 49.4% were annotated by seven functional databases. 23,121 unigenes were annotated and assigned to 457 species by non-redundant protein sequence database. The BLAST top-two hit species were Metaseiulus occidentalis and Ixodes scapularis. The procedure detected 12,426 SSRs, of which tri- and di-nucleotides were the most abundant types and the representative motifs were AAT/ATT and AC/GT. 26,936 CDS were predicted with a mean length of 711 bp. 87 unigenes of 30 functional genes, which are usually involved in stress responses, drug resistance, movement, metabolism and allergy, were further identified by bioinformatics methods. The unigenes putatively encoding cytochrome P450 proteins were further analyzed phylogenetically. In conclusion, this study completed the RNA-seq and functional annotation of O. bacoti successfully, which provides reliable molecular data for its future studies of gene function and molecular markers.

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.

Comparative genomics of chemosensory protein genes (CSPs) in twenty-two mosquito species (Diptera: Culicidae): Identification, characterization, and evolution

Chemosensory proteins (CSP) are soluble carrier proteins that may function in odorant reception in insects. CSPs have not been thoroughly studied at whole-genome level, despite the availability of insect genomes. Here, we identified/reidentified 283 CSP genes in the genomes of 22 mosquitoes. All 283 CSP genes possess a highly conserved OS-D domain. We comprehensively analyzed these CSP genes and determined their conserved domains, structure, genomic distribution, phylogeny, and evolutionary patterns. We found an average of seven CSP genes in each of 19 Anopheles genomes, 27 CSP genes in Cx. quinquefasciatus, 43 in Ae. aegypti, and 83 in Ae. albopictus. The Anopheles CSP genes had a simple genomic organization with a relatively consistent gene distribution, while most of the Culicinae CSP genes were distributed in clusters on the scaffolds. Our phylogenetic analysis clustered the CSPs into two major groups: CSP1-8 and CSE1-3. The CSP1-8 groups were all monophyletic with good bootstrap support. The CSE1-3 groups were an expansion of the CSP family of genes specific to the three Culicinae species. The Ka/Ks ratios indicated that the CSP genes had been subject to purifying selection with relatively slow evolution. Our results provide a comprehensive framework for the study of the CSP gene family in these 22 mosquito species, laying a foundation for future work on CSP function in the detection of chemical cues in the surrounding environment.

Identification of carboxylesterase genes associated with pyrethroid resistance in the malaria vector Anopheles sinensis (Diptera: Culicidae)

BACKGROUND

Carboxylesterases (CCEs) are one of three large detoxification enzyme families. Some CCEs are active on synthetic insecticides with ester structures. Anopheles sinensis is an important malaria vector in eastern Asia. This study identified and characterized the CCE genes in the A. sinensis genome and determined CCE genes associated with pyrethroid resistance using RNA sequencing (RNA-seq) and quantitative reverse transcription - polymerase chain reaction (qRT-PCR), in A. sinensis from Anhui, Chongqing, and Yunnan in China.

RESULTS

Fifty-seven putative CCEs were identified and placed into three classes, 12 subfamilies and 14 clades through phylogenetic and homology analyses. Exon sizes ranged from 31 to 4317 bp, with 49 CCEs having two to five exons and eight having six to 11 exons. A total of 183 introns were recognized with sizes ranging from 31 to 4317 bp. The 57 CCEs were located on 14 scaffolds, with 70% located on four scaffolds. The alpha-esterase subfamily was significantly expanded compared with that of Anopheles gambiae. In a pyrethroid-resistant strain, RNA-seq detected five upregulated CCE genes and qRT-PCR detected 12 upregulated CCE genes. The α-esterase 10 (AsAe10) and acetylcholinesterase 1 (AsAce1) genes were the main CCE genes associated with pyrethroid resistance.

CONCLUSION

This information will be useful for further study of the CCE gene family and pyrethroid resistance mechanisms mediated by CCEs. © 2017 Society of Chemical Industry.

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.

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.

Suppression of Laccase 2 severely impairs cuticle tanning and pathogen resistance during the pupal metamorphosis of Anopheles sinensis (Diptera: Culicidae)

Background

Phenol oxidases (POs) catalyze the oxidation of dopa and dopamine to melanin, which is crucial for cuticle formation and innate immune maintenance in insects. Although, Laccase 2, a member of the PO family, has been reported to be a requirement for melanin-mediated cuticle tanning in the development stages of some insects, whether it participates in cuticle construction and other physiological processes during the metamorphosis of mosquito pupae is unclear.

Methods

The association between the phenotype and the expression profile of Anopheles sinensis Laccase 2 (AsLac2) was assessed from pupation to adult eclosion. Individuals showing an expression deficiency of AsLac2 that was produced by RNAi and their phenotypic defects and physiological characterizations were compared in detail with the controls.

Results

During the dominant expression period, knockdown of AsLac2 in pupae caused the cuticle to be unpigmented, and produced thin and very soft cuticles, which further impeded the eclosion rate of adults as well as their fitness. Moreover, melanization immune responses in the pupae were sharply decreased, leading to poor resistance to microorganism infection. Both the high conservation among Laccase 2 homologs and a very similar genomic synteny of the neighborhood in Anopheles genus implies a conservative function in the pupal stage.

Conclusions

To our knowledge, this is the first study to report the serious phenotypic defects in mosquito pupae caused by the dysfunction of Laccase 2. Our findings strongly suggest that Laccase 2 is crucial for Anopheles cuticle construction and melanization immune responses to pathogen infections during pupal metamorphosis. This irreplaceability provides valuable information on the application of Lacccase 2 and/or other key genes in the melanin metabolism pathway for developing mosquito control strategies.

Electronic supplementary material

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

De novo transcriptome analysis of Spodoptera exigua multiple nucleopolyhedrovirus (SeMNPV) genes in latently infected Se301 cells

Cells of the P8-Se301-C1 strain are Spodoptera exigua cell clones that each harbor a partial version of the S. exigua multiple nucleopolyhedrovirus (SeMNPV) genome and which are resistant to homologous SeMNPV infections. The cells produce no viral progeny, suggesting that the infection is a latent-like viral infection. To investigate the SeMNPV genes harbored in the P8-Se301-C1 cells, the de novo transcriptomes of P8-Se301-C1 cells and S. exigua Se301 cells were analyzed and compared. A total of 54,569,296 reads were obtained from the P8-Se301-C1 cells that yielded 112,565 final unigenes with a mean length of 1,093 nt. A total of 56,865,504 reads were obtained from the Se301 cells that yielded 102,996 final unigenes with a mean length of 1,082 nt. Ten SeMNPV gene transcripts (se5, se7, se8, se12, se43, se45, se89, se90, se124, and se126) were detected in the P8-Se301-C1 cells by RNA-Seq but not in the Se301 cells, which was verified by RTPCR. 5'/3' RACE analyses showed that the 3'- or 5'-end sequences of the viral transcripts are aligned to the host gene sequences in P8-Se301-C1 cells, suggesting that the SeMNPV genes may integrate into and be transcribed with the host genes in the P8-Se301-C1 cells. Furthermore, six additional viral gene transcripts, se11, se42, se44, se88, se91, and se127 (incorporated into chimeric fusion transcripts in the P8-Se301-C1 cells), were detected in the RACE analyses. Taken together, sixteen SeMNPV transcripts were identified in the P8-Se301-C1 cell strain. This study provides information to develop the understanding of baculovirus latent infections and superinfection exclusion.

Genome-wide identification and characterization of odorant-binding protein (OBP) genes in the malaria vector Anopheles sinensis (Diptera: Culicidae)

Anopheles sinensis is a major malaria vector. Insect odorant-binding proteins (OBPs) may function in the reception of odorants in the olfactory system. The classification and characterization of the An. sinensis OBP genes have not been systematically studied. In this study, 64 putative OBP genes were identified at the whole-genome level of An. sinensis based on the comparison between OBP conserved motifs, PBP_GOBP, and phylogenetic analysis with An. gambiae OBPs. The characterization of An. sinensis OBPs, including the motif's conservation, gene structure, genomic organization and classification, were investigated. A new gene, AsOBP73, belonging to the Plus-C subfamily, was identified with the support of transcript and conservative motifs. These An. sinensis OBP genes were classified into three subfamilies with 37, 15 and 12 genes in the subfamily Classic, Atypical and Plus-C, respectively. The genomic organization of An. sinensis OBPs suggests a clustered distribution across nine different scaffolds. Eight genes (OBP23-28, OBP63-64) might originate from a single gene through a series of historic duplication events at least before divergence of Anopheles, Culex and Aedes. The microsynteny analyses indicate a very high synteny between An. sinensis and An. gambiae OBPs. OBP70 and OBP71 earlier classified under Plus-C in An. gambiae are recognized as belonging to the group Obp59a of the Classic subfamily, and OBP69 earlier classified under Plus-C has been moved to the Atypical subfamily in this study. The study established a basic information frame for further study of the OBP genes in insects as well as in An. sinensis.

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.

Advances in genetics and genomics: use and limitations in achieving malaria elimination goals

Success of the global research agenda towards eradication of malaria will depend on the development of new tools, including drugs, vaccines, insecticides and diagnostics. Genetic and genomic information now available for the malaria parasites, their mosquito vectors and human host, can be harnessed to both develop these tools and monitor their effectiveness. Here we review and provide specific examples of current technological advances and how these genetic and genomic tools have increased our knowledge of host, parasite and vector biology in relation to malaria elimination and in turn enhanced the potential to reach that goal. We then discuss limitations of these tools and future prospects for the successful achievement of global malaria elimination goals.

De novo assembly and sex-specific transcriptome profiling in the sand fly Phlebotomus perniciosus (Diptera, Phlebotominae), a major Old World vector of Leishmania infantum

Background

The phlebotomine sand fly Phlebotomus perniciosus (Diptera: Psychodidae, Phlebotominae) is a major Old World vector of the protozoan Leishmania infantum, the etiological agent of visceral and cutaneous leishmaniases in humans and dogs, a worldwide re-emerging diseases of great public health concern, affecting 101 countries. Despite the growing interest in the study of this sand fly species in the last years, the development of genomic resources has been limited so far. To increase the available sequence data for P. perniciosus and to start studying the molecular basis of the sexual differentiation in sand flies, we performed whole transcriptome Illumina RNA sequencing (RNA-seq) of adult males and females and de novo transcriptome assembly.

Results

We assembled 55,393 high quality transcripts, of which 29,292 were unique, starting from adult whole body male and female pools. 11,736 transcripts had at least one functional annotation, including full-length low abundance salivary transcripts, 981 transcripts were classified as putative long non-coding RNAs and 244 transcripts encoded for putative novel proteins specific of the Phlebotominae sub-family. Differential expression analysis identified 8590 transcripts significantly biased between sexes. Among them, some show relaxation of selective constraints when compared to their orthologs of the New World sand fly species Lutzomyia longipalpis.

Conclusions

In this paper, we present a comprehensive transcriptome resource for the sand fly species P. perniciosus built from short-read RNA-seq and we provide insights into sex-specific gene expression at adult stage. Our analysis represents a first step towards the identification of sex-specific genes and pathways and a foundation for forthcoming investigations into this important vector species, including the study of the evolution of sex-biased genes and of the sexual differentiation in phlebotomine sand flies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2088-x) contains supplementary material, which is available to authorized users.