Genome-wide identification, characterization and evolution of cuticular protein genes in 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.

Candidatus Liberibacter asiaticus influences the emergence of the Asian citrus psyllid Diaphorina citri by regulating key cuticular proteins

The Asian citrus psyllid, Diaphorina citri, is the primary vector of the HLB pathogen, Candidatus Liberibacter asiaticus (CLas). The acquisition of CLas shortens the developmental period of nymphs, accelerating the emergence into adulthood and thereby facilitating the spread of CLas. Cuticular proteins (CPs) are involved in insect emergence. In this study, we investigated the molecular mechanisms underlying CLas-promoted emergence in D. citri via CP mediation. Here, a total of 159 CP genes were first identified in the D. citri genome. Chromosomal location analysis revealed an uneven distribution of these CP genes across the 13 D. citri chromosomes. Proteomic analysis identified 54 differentially expressed CPs during D. citri emergence, with 14 CPs exhibiting significant differential expression after CLas acquisition. Five key genes, Dc18aa-1, Dc18aa-2, DcCPR-24, DcCPR-38 and DcCPR-58, were screened from the proteome and CLas acquisition. The silencing of these 5 genes through a modified feeding method significantly reduced the emergence rate and caused various abnormal phenotypes, indicating the crucial role that these genes play in D. citri emergence. This study provides a comprehensive overview of the role of CPs in D. citri and reveals that CLas can influence the emergence process of D. citri by regulating the expression of CPs. These key CPs may serve as potential targets for future research on controlling huanglongbing (HLB) transmission.

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.

Expressions and functions of RR-1 cuticular protein genes in the integument of Mythimna separata

As the most outer layer between itself and the environment, integuments are necessary for insects with various important functions. Cuticular proteins (CPs) are the main components in integuments, while the functions of CP genes remain unknown in Mythimna separata (Walker), which is a devastating agricultural pest. In this study, 79 CP genes were identified from the transcriptomes of larval integuments, 57 of which were from the family containing conserved Rebers & Riddiford (R&R) consensus (CPR family). Amongst these CPRs, 44 genes belonged to the subfamily with RR-1 motif (RR-1 genes) and clustered into three clades, with the top 15 most abundant RR-1 genes identified based on fragments per kilobase per million mapped fragments (FPKM) values. RT-qPCR analysis showed that most of RR-1 genes such as MsCPR1-4 were highly expressed at larval stages and in their integuments. The expression levels of RR-1 genes were generally decreased at the beginning but increased at the late stage of molting process. RNAi was applied for six RR-1 genes, and MsCPR1-4 were knocked down significantly. Silence of MsCPR2 resulted in abnormal integument formed after molting, while knockdown of MsCPR3 and MsCPR4 led to failure of molting, respectively. No phenotype was obtained for the RNAi of MsCPR1. Therefore, the expression of RR-1 genes and their functions were analyzed in the development of integuments in M. separata, providing new insights of RR-1 genes and potential targets for the development of growth regulators and new insecticides for M. separata.

Global Gene Expression Analysis Reveals Complex Cuticle Organization of the Tribolium Compound Eye

The red flour beetle Tribolium castaneum is a resource-rich model for genomic and developmental studies. To extend previous studies on Tribolium eye development, we produced transcriptomes for normal-eyed and eye-depleted heads of pupae and adults to identify differentially transcript-enriched (DE) genes in the visual system. Unexpectedly, cuticle-related genes were the largest functional class in the pupal compound eye DE gene population, indicating differential enrichment in three distinct cuticle components: clear lens facet cuticle, highly melanized cuticle of the ocular diaphragm, which surrounds the Tribolium compound eye for internal fortification, and newly identified facet margins of the tanned cuticle, possibly enhancing external fortification. Phylogenetic, linkage, and high-throughput gene knockdown data suggest that most cuticle proteins (CPs) expressed in the Tribolium compound eye stem from the deployment of ancient CP genes. Consistent with this, TcasCPR15, which we identified as the major lens CP gene in Tribolium, is a beetle-specific but pleiotropic paralog of the ancient CPR RR-2 CP gene family. The less abundant yet most likely even more lens-specific TcasCP63 is a member of a sprawling family of noncanonical CP genes, documenting a role of local gene family expansions in the emergence of the Tribolium compound eye CP repertoire. Comparisons with Drosophila and the mosquito Anopheles gambiae reveal a steady turnover of lens-enriched CP genes during insect evolution.

RR-1 cuticular protein TcCPR69 is required for growth and metamorphosis in Tribolium castaneum

Cuticle is not only critical for protecting insects from noxious stimuli but is also involved in a variety of metabolic activities. Cuticular proteins (CPs) affect cuticle structure and mechanical properties during insect growth, reproduction, and environmental adaptation. Here, we describe the identification and characterization of a member of the RR-1 subfamily of CPs with an R&R consensus (CPR) in Tribolium castaneum (TcCPR69). Although it was previously reported to be highly expressed in the wings, we found that knocking down TcCPR69 by RNA interference (RNAi) did not cause obvious wing abnormalities but markedly disrupted the growth and metamorphosis of beetles with 100% cumulative mortality; additionally, the chitin content of the pharate adult was decreased and the new abdominal cuticle was significantly thinner before molting. TcCPR69 showed chitin-binding ability and the expression levels of key genes involved in chitin metabolism (trehalase [TcTRE], chitin synthase [TcCHSA and TcCHSB], and chitinase [TcCHT5 and TcCHT10]) were also decreased by TcCPR69 knockdown. TcCPR69 gene expression peaked shortly after molting and was increased 2.61 fold at 12 h after 20-hydroxyecdysone (20E) injection. This was reversed by RNAi of the ecdysone-related genes ecdysone receptor (TcECR) and fushi tarazu transcription factor 1 (TcFTZ-F1). These results indicate that TcCPR69 is positively regulated by 20E signaling to contribute to cuticle formation and maintain chitin accumulation during the growth and metamorphosis of beetles.

Characterization of cuticular proteins in CPR family in the wolf spider, Pardosa pseudoannulata, and the response of one subfamily genes to environmental stresses

Cuticular protein (CP) plays an essential role in the construction and function of exoskeleton in arthropods. CPR family, CP with Rebers and Riddiford (R&R) Consensus, is the largest CP family in insects, but it lacks systematic research in non-insect arthropods. In this study, we explored CPRs in the wolf spider, Pardosa pseudoannulata, a predator to many insect pests. We totally identified 152 CPRs in P. pseudoannulata genome, which were divided into two subgroups based on R&R Consensus sequences, with 12 CPRs in RR-1 and 140 in RR-2. All RR-2 members presented a novel Consensus with 34 amino acids, G-x(8)-G-x(6)-Y-x-A-x(3)-G-x(7)-N-E-x-G, which was a common characteristic for RR-2 CPRs in chelicerates. Transcriptome data was used to document the expression patterns of CPR genes in different tissues and ecdysis processes. The specific expressions were found for part CPR genes, such as five RR-2 genes that were specifically expressed in male genital bulbs and eleven RR-1 genes that were highly expressed in the integument. Due to the limited number and integument-specific expression of RR-1 genes, we further analyzed their responses to different environmental stresses at the transcriptional level. Except for PapsCPR11, ten RR-1 genes responded to at least one environmental stress, among with the expression of PapsCPR12 was significantly changed by three stresses (dryness, low temperature and imidacloprid treatments). Silencing PapsCPR12 increased the tolerance of P. pseudoannulata to imidacloprid. Overall, the results presented novel Consensus characteristics of CPRs in P. pseudoannulata, which was helpful for the identification and evolution analysis of CPRs in non-insect arthropods.

Genomic and Transcriptomic Analysis Reveals Cuticular Protein Genes Responding to Different Insecticides in Fall Armyworm Spodoptera frugiperda

The fall armyworm (FAW), Spodoptera frugiperda, is a serious pest of crucial crops causing great threats to the food security of the world. It has evolved resistance to various insecticides, while the underlying molecular mechanisms remain largely unknown. Cuticular proteins (CPs), as primary components in cuticle, play an important role in insects' protection against environmental stresses. Few of them have been documented as participating in insecticide resistance in several insect species. In order to explore whether CP genes of the FAW exhibit a functional role in responding to insecticides stress, a total of 206 CPs, classified into eight families, were identified from the genome of the FAW through a homology-based approach coupled with manual efforts. The temporal expression profiles of all identified CP genes across developmental stages and their responses to 23 different insecticides were analyzed using the RNA-seq data. Expression profiling indicated that most of the CP genes displayed stage-specific expression patterns. It was found that the expression of 51 CP genes significantly changed after 48 h exposure to 17 different insecticides. The expression of eight CP genes responding to four insecticides were confirmed by RT-PCR analysis. The results showed that their overall expression profiles were consistent with RNA-seq analysis. The findings provide a basis for further functional investigation of CPs implied in insecticide stress in FAW.

Identification and expression analysis of cuticular protein genes in the diamondback moth, Plutella xylostella (Lepidoptera: Plutellidae)

Structural cuticular proteins (CPs) are major components of the insect cuticle, and they play critical roles in insect development and insecticide resistance. Here, a total of 196 CP genes were successfully annotated in the Plutella xylostella genome. On the basis of motif analysis, these CPs were classified into 10 different families, including 122 CPR, 12 CPAP1, 8 CPAP3, 9 CPLCP, 2 Tweedle, 1 CPF, 1 CPFL, 1 CPCFC, 17 CPG and 2 18 aa proteins, and the remaining 21 unclassified CPs were classed as cuticular proteins hypothetical (CPH). A phylogenetic analysis of CPs from different insects revealed species-specific clades of RR-1 and RR-2 genes, suggesting that CP gene duplication might occur independently among insect taxa, while we also found that some other CPs (such as CPAP1 and CPAP3) had a closer relationship based on their conserved domain architecture. Using available RNAseq libraries, the expression profiles of the CPs were analyzed over the four developmental stages of the insect (i.e., egg, larva, pupa, and adult), revealing stage-specific expression patterns for the CPs. In a chlorpyrifos resistant strain, 18 CP genes were found to be more than two-fold upregulated compared to the susceptible control strain, and qRT-PCR analysis showed that these CP genes were overexpressed after exposure to chlorpyrifos, suggesting a potential role in the molecular mechanism of insecticide resistance in P. xylostella. This study provides the tools and molecular basis to study the role of CPs in the post-embryonal development and the mechanisms of insecticide resistance of P. xylostella.

Transcriptome profiling of maternal stress-induced wing dimorphism in pea aphids

Wing dimorphism, that is, wingless and winged forms, can be induced by maternal stress signals and is an adaptive response of aphids to environmental changes. Here, we investigated the ecological and molecular effects of three kinds of stress, namely crowding, predation, and aphid alarm pheromone, on wing dimorphism. These three stressors induced high proportion of up to 60% of winged morphs in offspring. Transcriptome analysis of stress-treated female aphids revealed different changes in maternal gene expression induced by the three stressors. Crowding elicited widespread changes in the expression of genes involved in nutrient accumulation and energy mobilization. Distinct from crowding, predation caused dramatic expression changes in cuticle protein (CP) genes. Twenty-three CP genes that belong to CP RR2 subfamily and are highly expressed in legs and embryos were greatly repressed by the presence of ladybird. By contrast, application of alarm pheromone, E-β-farnesene, caused slight changes in gene expression. The three factors shared a responsive gene, cuticle protein 43. This study reveals the adaptive response of aphids to environmental stresses and provides a rich resource on genome-wide expression genes for exploring molecular mechanisms of ecological adaptation in aphids.

OPEN RESEARCH BADGES

This article has earned an Open Data Badge for making publicly available the digitally-shareable data necessary to reproduce the reported results. The data is available at https://doi.org/10.5061/dryad.55b2b15.

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.

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.

3,4-Dihydroxyphenylacetaldehyde synthase and cuticle formation in insects

Cuticle is the most important structure that protects mosquitoes and other insect species from adverse environmental conditions and infections of microorganism. The physiology and biochemistry of insect cuticle formation have been studied for many years and our understanding of cuticle formation and hardening has increased considerably. This is especially true for flexible cuticle. The recent discovery of a novel enzyme that catalyzes the production of 3,4-dihydroxyphenylacetaldehyde (DOPAL) in insects provides intriguing insights concerning the flexible cuticle formation in insects. For convenience, the enzyme that catalyzes the production DOPAL from l-dopa is named DOPAL synthase. In this mini-review, we summarize the biochemical pathways of cuticle formation and hardening in general and discuss DOPAL synthase-mediated protein crosslinking in insect flexible cuticle in particular.