Absence of known knockdown resistance mutations but fixation of CYP337B3 was detected in field populations of Helicoverpa armigera across China.

The cotton bollworm (Helicoverpa armigera) is a worldwide agricultural pest that infests many important crops. Pyrethroids targeting the voltage-gated sodium channel (VGSC) have been long used in the control of the cotton bollworm. Two amino acid substitutions (D1561V and E1565G) in H. armigera VGSC (HaVGSC) and the presence of a chimeric P450 gene (CYP337B3) have been documented to be associated with pyrethroid resistance. To understand the current occurrence of kdr mutations and the CYP337B3 gene in Chinese H. armigera populations, high-throughput amplicon sequencing was adopted to detect potential nucleotide variations in three fragments of the VGSC gene that cover 10 reported knockdown resistance (kdr) sites in insects, and gene-specific PCR was performed to examine the presence of CYP337B3 gene in H. armigera samples collected across China. The nucleotide variation analysis revealed a wealth of nucleotide variations in not only exons but also introns in the VGSC gene in Chinese H. armigera populations. However, neither previously reported kdr-conferring amino acid replacements nor other non-synonymous mutations were observed in a total of 1439 examined individuals. Population genetic analysis suggested that the H. armigera population in Nanchang, Jiangxi Province (JNC) had a moderate genetic differentiation from other populations, while no significant divergence was observed in other populations in northern and northwestern China. The CYP337B3 was present in all the examined individuals, indicating that CYP337B3 is extensively fixed in H. armigera populations across China. These results support that point mutations in VGSC are not a major factor involved in the current pyrethroid resistance in H. armigera. Instead, CYP337B3 plays a prevalent role in the development of resistance to pyrethroids in H. armigera.

Up-regulation of CYP6G4 mediated by a CncC/maf binding-site-containing insertion confers resistance to multiple classes of insecticides in the house fly Musca domestica.

Populations of many insect species have evolved a variety of resistance mechanisms in response to insecticide selection. Current knowledge about mutations responsible for insecticide resistance is largely achieved from studies on target-site resistance, while much less is known about metabolic resistance. Although it is well known that P450 monooxygenases are one of the major players involved in insecticide metabolism and resistance, understanding mutation(s) responsible for CYP-mediated resistance has been a big challenge. In this study, we used the house fly to pursue a better understanding of P450 mediated insecticide resistance at the molecular level. Metabolism studies illustrated that CYP6G4 had a broad-spectrum metabolic activity in metabolizing insecticides. Population genotyping revealed that the CYP6G4v1 allele harboring a DNA insertion (MdIS1) had been selected in many house fly populations on different continents. Dual luciferase reporter assays identified that the MdIS1 contained a CncC/Maf binding site, and electrophoretic mobility shift assay confirmed that transcription factor CncC was involved in the MdIS1-mediated regulation. This study highlights the common involvement of the CncC pathway in adaptive evolution, and provides an interesting case supportive of parallel evolution in P450-mediated insecticide resistance in insects.

Detection of ryanodine receptor G4911E and I4754M mutation sites and analysis of binding modes of diamide insecticides with RyR on Galeruca daurica (Coleoptera: Chrysomelidae).

In recent years, the leaf beetle Galeruca daurica has broken out in the northern grasslands of Inner Mongolia, its management still mainly depends on chemical control using traditional insecticides or with novel action. The study was aim to identify mutation locus associated with resistance to diamide insecticides in field population of G. daurica, to provide a reference for rational selection of insecticides and to avoid the rapid resistance development to diamide insecticides. We cloned the full length of the ryanodine receptor gene of G. daurica (GdRyR), constructed 3D model and transmembrane regions by homologous modeling based on deduced amino acid sequence. Two potential mutation loci (Gly4911Glu and Ile4754Met) and allelic mutation frequencies were detected in individuals of G. daurica. In addition, their binding patterns to two diamide insecticides (chlorantraniliprole, cyantraniliprole) were analyzed separately using a molecular docking method. The full-length cDNA sequence of GdRyR (GenBank accession number: OP828593) was obtained by splicing and assembling, which is 15,399 bp in length and encodes 5,133 amino acids. The amino acid similarity of GdRyR with that of other Coleopteran insects were 86.70%-91.33%, which possessed the typical structural characteristics. An individual resistance allelic mutation frequency test on fifty field leaf beetles has identified 12% and 32% heterozygous individuals at two potential mutation loci Gly4911Glu and Ile4754Met, respectively. The affinity of the I4754M mutant model of GdRyR for chlorantraniliprole and cyantraniliprole was not significantly different from that of the wild type, and all had non-covalent interactions such as hydrogen bonding, hydrophobic interactions and π-cation interactions. However, the G4911E mutant model showed reduced affinity and reduced mode of action with two diamide insecticides, thus affecting the binding stability of the ryanodine receptor to the diamide insecticides. In conclusion, the G4911E mutation in GdRyR may be a potential mechanism for the development of resistance to diamide insecticides on G. daurica and should be a key concern for resistance risk assessment and reasonable applications of diamide insecticides for control in future. Moreover, this study could provide a reference for ryanodine receptor structure-based insecticides design.

CYP4G8 is responsible for the synthesis of methyl-branched hydrocarbons in the polyphagous caterpillar of Helicoverpa armigera.

Insect cuticular hydrocarbons (CHCs) have dual functions as physical barrier and chemical signals. The last step of CHC biosynthesis is known to be catalyzed by cytochrome P450 CYP4G in a number of insects. Until recently, studies on CYP4Gs in the context of functional evolution are rare. In this study, we analyzed sequence similarity and temporal-spatial expression patterns of the five CYP4G genes in the cotton bollworm Helicoverpa armigera, an important agricultural pest and also typical representative of lepidopteran insects. Moreover, the CRISPR/Cas9-induced knockout was used to clarify the roles of the five CYP4Gs in CHC biosynthesis. Temporal-spatial expression patterns revealed that CYP4G8 was highly expressed at all developmental stages and in most tissues examined. Larvae with CYP4G8 knocked out could not produce methyl-branched CHCs and failed to pupate, while larvae with the other four CYP4G genes knocked out (4G1-type-KO) showed no significant changes in their CHC profiles, weight gain and survival. Comparative transcriptomics revealed that knocking out CYP4G8 affected the global gene expression in larvae, especially down-regulated the expression of genes in the fatty acid biosynthetic pathway, while no significant change in 4G1-type-KO transcriptome was observed. These findings indicate that the five members of the CYP4G subfamily have undergone functional divergence: CYP4G8 maintains the essential function in CHC biosynthesis, while the function of the other four CYP4G genes remains unclear. Intriguingly, CYP4G8 has evolved to be a P450 enzyme responsible for the synthesis of larval methyl-branched hydrocarbons. The observation that CYP4G8 knockout is lethal strongly suggest that CYP4G8 may serve as a candidate target for the development of insecticidal agents for the control of cotton bollworms.

Identification and phylogenetic analysis of voltage-gated sodium channel haplotypes in the malaria vector Anopheles sinensis using a high-throughput amplicon sequencing approach.

BACKGROUND: Anopheles sinensis is a dominant vector for malaria transmission in Asian countries. Voltage-gated sodium channel (VGSC) mutation-mediated knock-down resistance (kdr) has developed in many A. sinensis populations because of intensive and long-term use of pyrethroids. Our previous study showed that multiple mutations at position 1014 of the VGSC were heterogeneously distributed in A. sinensis populations across Sichuan, China. METHODS: To understand resistance genotypes at the haplotype level and reconstruct the phylogenetic relationship of VGSC haplotypes, a cost-effective next-generation sequencing (NGS)-based amplicon sequencing approach was established to clarify haplotypes containing codon 1014 of the VGSC gene from a total of 446 adults collected in 12 locations of Sichuan, China. RESULTS: Nineteen (19) haplotypes were identified, including 11 wild 1014L, 6 resistance 1014F, and 2 resistance 1014C haplotypes. We found that resistance haplotypes of A. sinensis VGSC were widely distributed at frequencies ranging from 3.67 to 92.61%. The frequencies of the 1014C haplotype in the southeast of Sichuan (Luzhou, Guangan, and Suining) were relatively higher than those in other sampling locations. Phylogenetic analyses support that kdr-type mutation at position 1014 is not singly originated and resistance 1014C haplotypes evolve from TTT-encoding 1014F. CONCLUSIONS: A cost-effective next-generation sequencing (NGS)-based amplicon sequencing approach has been established in this study. The data revealed the patchy distribution of VGSC resistance haplotypes with overall high frequencies in Sichuan, China. Phylogenetic analyses support multiple origins and sequential evolution (1014L → 1014F → 1014C) for kdr-type mutations in A. sinensis.

Vegetative development and host immune interaction of Ophiocordyceps sinensis within the hemocoel of the ghost moth larva, Thitarodes xiaojinensis.

Ophiocordyceps sinensis is an entomopathogenic fungus that infects ghost moth larva, forming the most valuable and rare traditional Chinese medicine, Chinese cordyceps. Our knowledge of the basic morphology and developmental biology of Chinese cordyceps is limited. In this study, morphological and ultrastructural observations of O. sinensis development in the hemocoel of Thitarodes xiaojinensis were obtained by multiple light and electron microscopy techniques, and the host immune reaction activities were determined. Our results indicated that fungal cells in the host hemocoel underwent morphotype transformations from blastospores to prehyphae to hyphae in sequence. The fusiform yeast-like blastospores were the initial cell type present in the host hemocoel and remained for 5 months or more; the encapsulation reaction and phenoloxidase activity of T. xiaojinensis hemolymph were inhibited during this period. When larvae entered the last instar, the blastospores switched to prehyphae and expanded throughout the host tissues, and then hyphae germinated from the prehyphae and mycelia formed, which finally led to host death. Considering the distinct differences between blastospores and hyphae, we identified prehyphae, which play important roles in fungal expansion, hyphae germination, and fusion formation among filaments. Notably, the elongation of prehyphae was strongly presumed to occur through fission but without separation of the two sister cells, in contrast to blastospore budding. During the morphotype transformation, the amount and composition of lipid droplets changed greatly, suggesting their important roles in these events. Overall, we provide a morphological and ultrastructural characterization of O. sinensis vegetative development within the hemocoel of T. xiaojinensis, identify and name the prehypha fungal cell type in entomopathogenic fungi for the first time, and conclude that O. sinensis infection causes sustained immunosuppression in T. xiaojinensis.

Comparative analysis of C-type lectin domain proteins in the ghost moth, Thitarodes xiaojinensis (Lepidoptera: Hepialidae).

Insects have a large family of C-type lectins involved in cell adhesion, pathogen recognition and activation of immune responses. In this study, 32 transcripts encoding C-type lectin domain proteins (CTLDPs) were identified from the Thitarodes xiaojinensis transcriptome. According to their domain structures, six CTLDPs with one carbohydrate-recognition domain (CRD) were classified into the CTL-S subfamily. The other 23 CTLDPs with two CRDs were grouped into the immulectin (IML) subfamily. The remaining three with extra regulatory domains were sorted into the CTL-X subfamily. Phylogenetic analysis showed that CTL-S and CTL-X members from different insects could form orthologous groups. In contrast, no T. xiaojinensis IML orthologues were found in other insects. Remarkable lineage-specific expansion in this subfamily was observed reflecting that these CTLDPs, as important receptors, have evolved diversified members in response to a variety of microbes. Prediction of binding ligands revealed that T. xiaojinensis, a cold-adapted species, conserved the ability of CRDs to combine with Ca2+ to keep its receptors from freezing. Comparative analysis of induction of CTLDP genes after different immune challenges indicated that IMLs might play critical roles in immune defenses. This study examined T. xiaojinensis CTLDPs and provides a basis for further studies of their characteristics.