The impact of agrochemical pollutant mixtures on the selection of insecticide resistance in the malaria vector Anopheles gambiae: insights from experimental evolution and transcriptomics.

BACKGROUND: There are several indications that pesticides used in agriculture contribute to the emergence and spread of resistance of mosquitoes to vector control insecticides. However, the impact of such an indirect selection pressure has rarely been quantified and the molecular mechanisms involved are still poorly characterized. In this context, experimental selection with different agrochemical mixtures was conducted in Anopheles gambiae. The multi-generational impact of agrochemicals on insecticide resistance was evaluated by phenotypic and molecular approaches. METHODS: Mosquito larvae were selected for 30 generations with three different agrochemical mixtures containing (i) insecticides, (ii) non-insecticides compounds, and (iii) both insecticide and non-insecticide compounds. Every five generations, the resistance of adults to deltamethrin and bendiocarb was monitored using bioassays. The frequencies of the kdr (L995F) and ace1 (G119S) target-site mutations were monitored every 10 generations. RNAseq was performed on all lines at generation 30 in order to identify gene transcription level variations and polymorphisms associated with each selection regime. RESULTS: Larval selection with agrochemical mixtures did not affect bendiocarb resistance and did not select for ace1 mutation. Contrastingly, an increased deltamethrin resistance was observed in the three selected lines. Such increased resistance was not majorly associated with the presence of kdr L995F mutation in selected lines. RNA-seq identified 63 candidate resistance genes over-transcribed in at least one selected line. These include genes coding for detoxification enzymes or cuticular proteins previously associated with insecticide resistance, and other genes potentially associated with chemical stress response. Combining an allele frequency filtering with a Bayesian FST-based genome scan allowed to identify genes under selection across multiple genomic loci, supporting a multigenic adaptive response to agrochemical mixtures. CONCLUSION: This study supports the role of agrochemical contaminants as a significant larval selection pressure favouring insecticide resistance in malaria vectors. Such selection pressures likely impact kdr mutations and detoxification enzymes, but also more generalist mechanisms such as cuticle resistance, which could potentially lead to cross-tolerance to unrelated insecticide compounds. Such indirect effect of global landscape pollution on mosquito resistance to public health insecticides deserves further attention since it can affect the nature and dynamics of resistance alleles circulating in malaria vectors and impact the efficacy of control vector strategies.

Deltamethrin and transfluthrin select for distinct transcriptomic responses in the malaria vector Anopheles gambiae.

BACKGROUND: The widespread use of pyrethroid insecticides in Africa has led to the development of strong resistance in Anopheles mosquitoes. Introducing new active ingredients can contribute to overcome this phenomenon and ensure the effectiveness of vector control strategies. Transfluthrin is a polyfluorinated pyrethroid whose structural conformation was thought to prevent its metabolism by cytochrome P450 monooxygenases in malaria vectors, thus representing a potential alternative for managing P450-mediated resistance occurring in the field. In this study, a controlled selection was used to compare the dynamics of resistance between transfluthrin and the widely used pyrethroid deltamethrin in the mosquito Anopheles gambiae. Then, the associated molecular mechanisms were investigated using target-site mutation genotyping and RNA-seq. METHODS: A field-derived line of An. gambiae carrying resistance alleles at low frequencies was used as starting material for a controlled selection experiment. Adult females were selected across 33 generations with deltamethrin or transfluthrin, resulting in three distinct lines: the Delta-R line (selected with deltamethrin), the Transflu-R line (selected with transfluthrin) and the Tiassale-S line (maintained without selection). Deltamethrin and transfluthrin resistance levels were monitored in each selected line throughout the selection process, as well as the frequency of the L1014F kdr mutation. At generation 17, cross-resistance to other public health insecticides was investigated and transcriptomes were sequenced to compare gene transcription variations and polymorphisms associated with adaptation to each insecticide. RESULTS: A rapid increase in resistance to deltamethrin and transfluthrin was observed throughout the selection process in each selected line in association with an increased frequency of the L1014F kdr mutation. Transcriptomic data support a broader response to transfluthrin selection as compared to deltamethrin selection. For instance, multiple detoxification enzymes and cuticle proteins were specifically over-transcribed in the Transflu-R line including the known pyrethroid metabolizers CYP6M2, CYP9K1 and CYP6AA1 together with other genes previously associated with resistance in An. gambiae. CONCLUSION: This study confirms that recurrent exposure of adult mosquitoes to pyrethroids in a public health context can rapidly select for various resistance mechanisms. In particular, it indicates that in addition to target site mutations, the polyfluorinated pyrethroid transfluthrin can select for a broad metabolic response, which includes some P450s previously associated to resistance to classical pyrethroids. This unexpected finding highlights the need for an in-depth study on the adaptive response of mosquitoes to newly introduced active ingredients in order to effectively guide and support decision-making programmes in malaria control.

Experimental evolution supports the potential of neonicotinoid-pyrethroid combination for managing insecticide resistance in malaria vectors.

The introduction of neonicotinoids for managing insecticide resistance in mosquitoes is of high interest as they interact with a biochemical target not previously used in public health. In this concern, Bayer developed a combination of the neonicotinoid clothianidin and the pyrethroid deltamethrin (brand name Fludora Fusion) as a new vector control tool. Although this combination proved to be efficient against pyrethroid-resistant mosquitoes, its ability to prevent the selection of pyrethroid and neonicotinoid resistance alleles was not investigated. In this context, the objective of this work was to study the dynamics and the molecular mechanisms of resistance of An. gambiae to the separated or combined components of this combination. A field-derived An. gambiae line carrying resistance alleles to multiple insecticides at low frequencies was used as a starting for 33 successive generations of controlled selection. Resistance levels to each insecticide and target site mutation frequencies were monitored throughout the selection process. Cross resistance to other public health insecticides were also investigated. RNA-seq was used to compare gene transcription variations and polymorphisms across all lines. This study confirmed the potential of this insecticide combination to impair the selection of resistance as compared to its two separated components. Deltamethrin selection led to the rapid enrichment of the kdr L1014F target-site mutation. Clothianidin selection led to the over-transcription of multiple cytochrome P450s including some showing high homology with those conferring neonicotinoid resistance in other insects. A strong selection signature associated with clothianidin selection was also observed on a P450 gene cluster previously associated with resistance. Within this cluster, the gene CYP6M1 showed the highest selection signature together with a transcription profile supporting a role in clothianidin resistance. Modelling the impact of point mutations selected by clothianidin on CYP6M1 protein structure showed that selection retained a protein variant with a modified active site potentially enhancing clothianidin metabolism. In the context of the recent deployment of neonicotinoids for mosquito control and their frequent usage in agriculture, the present study highlights the benefit of combining them with other insecticides for preventing the selection of resistance and sustaining vector control activities.

A genomic amplification affecting a carboxylesterase gene cluster confers organophosphate resistance in the mosquito Aedes aegypti: From genomic characterization to high-throughput field detection.

By altering gene expression and creating paralogs, genomic amplifications represent a key component of short-term adaptive processes. In insects, the use of insecticides can select gene amplifications causing an increased expression of detoxification enzymes, supporting the usefulness of these DNA markers for monitoring the dynamics of resistance alleles in the field. In this context, the present study aims to characterize a genomic amplification event associated with resistance to organophosphate insecticides in the mosquito Aedes aegypti and to develop a molecular assay to monitor the associated resistance alleles in the field. An experimental evolution experiment using a composite population from Laos supported the association between the over-transcription of multiple contiguous carboxylesterase genes on chromosome 2 and resistance to multiple organophosphate insecticides. Combining whole genome sequencing and qPCR on specific genes confirmed the presence of a ~100-Kb amplification spanning at least five carboxylesterase genes at this locus with the co-existence of multiple structural duplication haplotypes. Field data confirmed their circulation in South-East Asia and revealed high copy number polymorphism among and within populations suggesting a trade-off between this resistance mechanism and associated fitness costs. A dual-color multiplex TaqMan assay allowing the rapid detection and copy number quantification of this amplification event in Ae. aegypti was developed and validated on field populations. The routine use of this novel assay will improve the tracking of resistance alleles in this major arbovirus vector.

Differences in the fungal communities nursed by two genetic groups of the alpine cushion plant, Silene acaulis.

Foundation plants shape the composition of local biotic communities and abiotic environments, but the impact of a plant's intraspecific variations on these processes is poorly understood. We examined these links in the alpine cushion moss campion (Silene acaulis) on two neighboring mountain ranges in the French Alps. Genotyping of cushion plants revealed two genetic clusters matching known subspecies. The exscapa subspecies was found on both limestone and granite, while the longiscapa one was only found on limestone. Even on similar limestone bedrock, cushion soils from the two S. acaulis subspecies deeply differed in their impact on soil abiotic conditions. They further strikingly differed from each other and from the surrounding bare soils in fungal community composition. Plant genotype variations accounted for a large part of the fungal composition variability in cushion soils, even when considering geography or soil chemistry, and particularly for the dominant molecular operational taxonomic units (MOTUs). Both saprophytic and biotrophic fungal taxa were related to the MOTUs recurrently associated with a single plant genetic cluster. Moreover, the putative phytopathogens were abundant, and within the same genus (Cladosporium) or species (Pyrenopeziza brassicae), MOTUs showing specificity for each plant subspecies were found. Our study highlights the combined influences of bedrock and plant genotype on fungal recruitment into cushion soils and suggests the coexistence of two mechanisms, an indirect selection resulting from the colonization of an engineered soil by free-living saprobes and a direct selection resulting from direct plant-fungi interactions.

Contrasting microbial biogeographical patterns between anthropogenic subalpine grasslands and natural alpine grasslands.

The effect of plant species composition on soil microbial communities was studied at the multiregional level. We compared the soil microbial communities of alpine natural grasslands dominated by Carex curvula and anthropogenic subalpine pastures dominated by Nardus stricta. We conducted paired sampling across the Carpathians and the Alps and used Illumina sequencing to reveal the molecular diversity of soil microbes. We found that bacterial and fungal communities exhibited contrasting regional distributions and that the distribution in each grassland is well discriminated. Beta diversity of microbial communities was much higher in C. curvula grasslands due to a marked regional effect. The composition of grassland-type core microbiomes suggest that C. curvula, and N. stricta to a lesser extent, tend to select a cohort of microbes related to antibiosis/exclusion, pathogenesis and endophytism. We discuss these findings in light of the postglacial history of the studied grasslands, the habitat connectivity and the disturbance regimes. Human-induced disturbance in the subalpine belt of European mountains has led to homogeneous soil microbial communities at large biogeographical scales. Our results confirm the overarching role of the dominant grassland plant species in the distribution of microbial communities and highlight the relevance of biogeographical history.

GeBP and GeBP-like proteins are noncanonical leucine-zipper transcription factors that regulate cytokinin response in Arabidopsis.

Understanding the role of transcription factors (TFs) is essential in reconstructing developmental regulatory networks. The plant-specific GeBP TF family of Arabidopsis thaliana (Arabidopsis) comprises 21 members, all of unknown function. A subset of four members, the founding member GeBP and GeBP-like proteins (GPL) 1, 2, and 3, shares a conserved C-terminal domain. Here we report that GeBP/GPL genes represent a newly defined class of leucine-zipper (Leu-zipper) TFs and that they play a redundant role in cytokinin hormone pathway regulation. Specifically, we demonstrate using yeast, in vitro, and split-yellow fluorescent protein in planta assays that GeBP/GPL proteins form homo- and heterodimers through a noncanonical Leu-zipper motif located in the C-terminal domain. A triple loss-of-function mutant of the three most closely related genes gebp gpl1 gpl2 shows a reduced sensitivity to exogenous cytokinins in a subset of cytokinin responses such as senescence and growth, whereas root inhibition is not affected. We find that transcript levels of type-A cytokinin response genes, which are involved in the negative feedback regulation of cytokinin signaling, are higher in the triple mutant. Using a GPL version that acts as a constitutive transcriptional activator, we show that the regulation of Arabidopsis response regulators (ARRs) is mediated by at least one additional, as yet unknown, repressor acting genetically downstream in the GeBP/GPL pathway. Our results indicate that GeBP/GPL genes encode a new class of unconventional Leu-zipper TF proteins and suggest that their role in the cytokinin pathway is to antagonize the negative feedback regulation on ARR genes to trigger the cytokinin response.

A plant porphyria related to defects in plastid import of protochlorophyllide oxidoreductase A.

The plastid envelope of higher plant chloroplasts is a focal point of plant metabolism. It is involved in numerous pathways, including tetrapyrrole biosynthesis and protein translocation. Chloroplasts need to import a large number of proteins from the cytosol because most are encoded in the nucleus. Here we report that a loss-of-function mutation in the outer plastid envelope 16-kDa protein (oep16) gene causes a conditional seedling lethal phenotype related to defects in import and assembly of NADPH:protochlorophyllide (Pchlide) oxidoreductase A. In the isolated knockout mutant of Arabidopsis thaliana, excess Pchlide accumulated in the dark operated as photosensitizer and provoked cell death during greening. Our results highlight the essential role of the substrate-dependent plastid import pathway of precursor Pchlide oxidoreductase A for seedling survival and the avoidance of developmentally programmed porphyria in higher plants.

The Arabidopsis TUBULIN-FOLDING COFACTOR A gene is involved in the control of the alpha/beta-tubulin monomer balance.

The control of the stoichiometric balance of alpha- and beta-tubulin is important during microtubule biogenesis. This process involves several tubulin-folding cofactors (TFCs), of which only TFC A is not essential in mammalian in vitro systems or in vivo in yeast. Here, we show that the TFC A gene is important in vivo in plants. The Arabidopsis gene KIESEL (KIS) shows sequence similarity to the TFC A gene. Expression of the mouse TFC A gene under the control of the 35S promoter rescues the kis mutation, indicating that KIS is the Arabidopsis ortholog of TFC A. kis plants exhibit a range of defects similar to the phenotypes associated with impaired microtubule function: plants are reduced in size and show meiotic defects, cell division is impaired, and trichomes are bulged and less branched. Microtubule density was indistinguishable from that of the wild type, but microtubule organization was affected in trichomes and hypocotyl cells of dark-grown kis plants. The kis phenotype was rescued by overexpression of an alpha-tubulin, indicating that KIS is involved in the control of the correct balance of alpha- and beta-tubulin monomers.

Functional analysis of the tubulin-folding cofactor C in Arabidopsis thaliana.

The biogenesis of microtubules comprises several steps, including the correct folding of alpha- and beta-tubulin and heterodimer formation. In vitro studies and the genetic analysis in yeast revealed that, after translation, alpha- and beta-tubulin are processed by several chaperonins and microtubule-folding cofactors (TFCs) to produce assembly-competent alpha-/beta-tubulin heterodimers. One of the TFCs, TFC-C, does not exist in yeast, and a potential function of TFC-C is thus based only on the biochemical analysis. In this study and in a very recently published study by Steinborn and coworkers, the analysis of the Arabidopsis porcino (por) mutant has shown that TFC-C is important for microtubule function in vivo. The predicted POR protein shares weak amino acid similarity with the human TFC-C (hTFC-C). Our finding that hTFC-C under the control of the ubiquitously expressed 35S promoter can rescue the por mutant phenotype shows that the POR gene encodes the Arabidopsis ortholog of hTFC-C. The analysis of plants carrying a GFP:POR fusion construct showed that POR protein is localized in the cytoplasm and is not associated with microtubules. While, in por mutants, microtubule density was indistinguishable from wild-type, their organization was affected.