The molecular evolution of cytochrome P450 genes within and between drosophila species

From sub-Saharan Africa to China: Evolutionary history and adaptation of Drosophila melanogaster revealed by population genomics

Drosophila melanogaster is a widely used model organism for studying environmental adaptation. However, the genetic diversity of populations in Asia is poorly understood, leaving a notable gap in our knowledge of the global evolution and adaptation of this species. We sequenced genomes of 292 D. melanogaster strains from various ecological settings in China and analyzed them along with previously published genome sequences. We have identified six global genetic ancestry groups, despite the presence of widespread genetic admixture. The strains from China represent a unique ancestry group, although detectable differentiation exists among populations within China. We deciphered the global migration and demography of D. melanogaster, and identified widespread signals of adaptation, including genetic changes in response to insecticides. We validated the effects of insecticide resistance variants using population cage trials and deep sequencing. This work highlights the importance of population genomics in understanding the genetic underpinnings of adaptation, an effort that is particularly relevant given the deterioration of ecosystems.

Genomes from historical Drosophila melanogaster specimens illuminate adaptive and demographic changes across more than 200 years of evolution

The ability to perform genomic sequencing on long-dead organisms is opening new frontiers in evolutionary research. These opportunities are especially notable in the case of museum collections, from which countless documented specimens may now be suitable for genomic analysis-if data of sufficient quality can be obtained. Here, we report 25 newly sequenced genomes from museum specimens of the model organism Drosophila melanogaster, including the oldest extant specimens of this species. By comparing historical samples ranging from the early 1800s to 1933 against modern-day genomes, we document evolution across thousands of generations, including time periods that encompass the species' initial occupation of northern Europe and an era of rapidly increasing human activity. We also find that the Lund, Sweden population underwent local genetic differentiation during the early 1800s to 1933 interval (potentially due to drift in a small population) but then became more similar to other European populations thereafter (potentially due to increased migration). Within each century-scale time period, our temporal sampling allows us to document compelling candidates for recent natural selection. In some cases, we gain insights regarding previously implicated selection candidates, such as ChKov1, for which our inferred timing of selection favors the hypothesis of antiviral resistance over insecticide resistance. Other candidates are novel, such as the circadian-related gene Ahcy, which yields a selection signal that rivals that of the DDT resistance gene Cyp6g1. These insights deepen our understanding of recent evolution in a model system, and highlight the potential of future museomic studies.

Evolution of five environmentally responsive gene families in a pine-feeding sawfly, Neodiprion lecontei (Hymenoptera: Diprionidae)

A central goal in evolutionary biology is to determine the predictability of adaptive genetic changes. Despite many documented cases of convergent evolution at individual loci, little is known about the repeatability of gene family expansions and contractions. To address this void, we examined gene family evolution in the redheaded pine sawfly Neodiprion lecontei, a noneusocial hymenopteran and exemplar of a pine-specialized lineage evolved from angiosperm-feeding ancestors. After assembling and annotating a draft genome, we manually annotated multiple gene families with chemosensory, detoxification, or immunity functions before characterizing their genomic distributions and molecular evolution. We find evidence of recent expansions of bitter gustatory receptor, clan 3 cytochrome P450, olfactory receptor, and antimicrobial peptide subfamilies, with strong evidence of positive selection among paralogs in a clade of gustatory receptors possibly involved in the detection of bitter compounds. In contrast, these gene families had little evidence of recent contraction via pseudogenization. Overall, our results are consistent with the hypothesis that in response to novel selection pressures, gene families that mediate ecological interactions may expand and contract predictably. Testing this hypothesis will require the comparative analysis of high-quality annotation data from phylogenetically and ecologically diverse insect species and functionally diverse gene families. To this end, increasing sampling in under-sampled hymenopteran lineages and environmentally responsive gene families and standardizing manual annotation methods should be prioritized.

Divergent amplifications of CYP9A cytochrome P450 genes provide two noctuid pests with differential protection against xenobiotics

Here, we provide mechanistic support for the involvement of the CYP9A subfamily of cytochrome P450 monooxygenases in the detoxification of host plant defense compounds and chemical insecticides in Spodoptera exigua and Spodoptera frugiperda. Our comparative genomics shows that a large cluster of CYP9A genes occurs in the two species but with significant differences in its contents, including several species-specific duplicates and substantial sequence divergence, both between orthologs and between duplicates. Bioassays of CRISPR-Cas9 knockouts of the clusters show that, collectively, the CYP9As can detoxify two furanocoumarin plant defense compounds (imperatorin and xanthotoxin) and insecticides representing three different chemotypes (pyrethroids, avermectins, and oxadiazines). However, in vitro metabolic assays of heterologously expressed products of individual genes show several differences between the species in the particular CYP9As with activities against these compounds. We also find that the clusters show tight genetic linkage with high levels of pyrethroid resistance in field strains of the two species. We propose that their divergent amplifications of the CYP9A subfamily have not only contributed to the development of the broad host ranges of these species over long evolutionary timeframes but also supplied them with diverse genetic options for evolving resistance to chemical insecticides in the very recent past.

Evolution of chemosensory and detoxification gene families across herbivorous Drosophilidae

Herbivorous insects are exceptionally diverse, accounting for a quarter of all known eukaryotic species, but the genomic basis of adaptations that enabled this dietary transition remains poorly understood. Many studies have suggested that expansions and contractions of chemosensory and detoxification gene families-genes directly mediating interactions with plant chemical defenses-underlie successful plant colonization. However, this hypothesis has been challenging to test because the origins of herbivory in many insect lineages are ancient (>150 million years ago (mya)), obscuring genomic evolutionary patterns. Here, we characterized chemosensory and detoxification gene family evolution across Scaptomyza, a genus nested within Drosophila that includes a recently derived (<15 mya) herbivore lineage of mustard (Brassicales) specialists and carnation (Caryophyllaceae) specialists, and several nonherbivorous species. Comparative genomic analyses revealed that herbivorous Scaptomyza has among the smallest chemosensory and detoxification gene repertoires across 12 drosophilid species surveyed. Rates of gene turnover averaged across the herbivore clade were significantly higher than background rates in over half of the surveyed gene families. However, gene turnover was more limited along the ancestral herbivore branch, with only gustatory receptors and odorant-binding proteins experiencing strong losses. The genes most significantly impacted by gene loss, duplication, or changes in selective constraint were those involved in detecting compounds associated with feeding on living plants (bitter or electrophilic phytotoxins) or their ancestral diet (fermenting plant volatiles). These results provide insight into the molecular and evolutionary mechanisms of plant-feeding adaptations and highlight gene candidates that have also been linked to other dietary transitions in Drosophila.

Aliens in the CYPome of the black fungus gnat, Bradysia coprophila

The diverse cytochrome P450 enzymes of insects play essential physiological roles and also play important roles in the metabolism of environmental chemicals such as insecticides. We manually curated the complement of P450 (CYP) genes, or CYPome, of the black fungus gnat, Bradysia (Sciara) coprophila (Diptera, Sciaroidea), a species with a variable number of chromosomes. This CYPome carries two types of "alien" P450 genes. The first type of alien P450s was found among the 163 CYP genes of the core genome (autosomes and X). They consist of 28 sequences resulting from horizontal gene transfer, with closest sequences not found in insects, but in other arthropods, often Collembola. These genes are not contaminants, because they are expressed genes with introns, found in synteny with regular dipteran genes, also found in B. odoriphaga and B. hygida. Two such "alien" genes are representatives of CYP clans not otherwise found in insects, a CYP53 sequence related to fungal CYP53 genes, and a CYP19-like sequence similar to some collembolan sequences but of unclear origin. The second type of alien P450s are represented by 99 sequences from germline-restricted chromosomes (GRC). While most are P450 pseudogenes, 33 are apparently intact, with half being more closely related to P450s from Cecidomyiidae than from Sciaridae, thus supporting the hypothesis of a cross-family hybridization origin of the GRC.

Wood feeding and social living: Draft genome of the subterranean termite Reticulitermes lucifugus (Blattodea; Termitoidae)

Termites (Insecta, Blattodea, Termitoidae) are a widespread and diverse group of eusocial insects known for their ability to digest wood matter. Herein, we report the draft genome of the subterranean termite Reticulitermes lucifugus, an economically important species and among the most studied taxa with respect to eusocial organization and mating system. The final assembly (~813 Mb) covered up to 88% of the estimated genome size and, in agreement with the Asexual Queen Succession Mating System, it was found completely homozygous. We predicted 16,349 highly supported gene models and 42% of repetitive DNA content. Transposable elements of R. lucifugus show similar evolutionary dynamics compared to that of other termites, with two main peaks of activity localized at 25% and 8% of Kimura divergence driven by DNA, LINE and SINE elements. Gene family turnover analyses identified multiple instances of gene duplication associated with R. lucifugus diversification, with significant lineage-specific gene family expansions related to development, perception and nutrient metabolism pathways. Finally, we analysed P450 and odourant receptor gene repertoires in detail, highlighting the large diversity and dynamical evolutionary history of these proteins in the R. lucifugus genome. This newly assembled genome will provide a valuable resource for further understanding the molecular basis of termites biology as well as for pest control.

Evolution of chemosensory and detoxification gene families across herbivorous Drosophilidae

Herbivorous insects are exceptionally diverse, accounting for a quarter of all known eukaryotic species, but the genetic basis of adaptations that enabled this dietary transition remains poorly understood. Many studies have suggested that expansions and contractions of chemosensory and detoxification gene families - genes directly mediating interactions with plant chemical defenses - underlie successful plant colonization. However, this hypothesis has been challenging to test because the origins of herbivory in many lineages are ancient (>150 million years ago [mya]), obscuring genomic evolutionary patterns. Here, we characterized chemosensory and detoxification gene family evolution across Scaptomyza, a genus nested within Drosophila that includes a recently derived (<15 mya) herbivore lineage of mustard (Brassicales) specialists and carnation (Caryophyllaceae) specialists, and several non-herbivorous species. Comparative genomic analyses revealed that herbivorous Scaptomyza have among the smallest chemosensory and detoxification gene repertoires across 12 drosophilid species surveyed. Rates of gene turnover averaged across the herbivore clade were significantly higher than background rates in over half of the surveyed gene families. However, gene turnover was more limited along the ancestral herbivore branch, with only gustatory receptors and odorant binding proteins experiencing strong losses. The genes most significantly impacted by gene loss, duplication, or changes in selective constraint were those involved in detecting compounds associated with feeding on plants (bitter or electrophilic phytotoxins) or their ancestral diet (yeast and fruit volatiles). These results provide insight into the molecular and evolutionary mechanisms of plant-feeding adaptations and highlight strong gene candidates that have also been linked to other dietary transitions in Drosophila .

Discovery of malathion resistance QTL in Drosophila melanogaster using a bulked phenotyping approach

Drosophila melanogaster has proved an effective system with which to understand the evolutionary genetics and molecular mechanisms of insecticide resistance. Insecticide use has left signatures of selection in the fly genome, and both functional and quantitative genetic studies in the system have identified genes and variants associated with resistance. Here, we use D. melanogaster and leverage a bulk phenotyping and pooled sequencing "extreme quantitative trait loci" approach to genetically dissect variation in resistance to malathion, an organophosphate insecticide. We resolve 2 quantitative trait loci, one of which implicates allelic variation at the cytochrome P450 gene Cyp6g1, a strong candidate based on previous work. The second shows no overlap with hits from a previous genome-wide association study for malathion resistance, recapitulating other studies showing that different strategies for complex trait dissection in flies can yield apparently different architectures. Notably, we see no genetic signal at the Ace gene. Ace encodes the target of organophosphate insecticide inhibition, and genome-wide association studies have identified strong Ace-linked associations with resistance in flies. The absence of quantitative trait locus implicating Ace here is most likely because our mapping population does not segregate for several of the known functional polymorphisms impacting resistance at Ace, perhaps because our population is derived from flies collected prior to the widespread use of organophosphate insecticides. Our fundamental approach can be an efficient, powerful strategy to dissect genetic variation in resistance traits. Nonetheless, studies seeking to interrogate contemporary insecticide resistance variation may benefit from deriving mapping populations from more recently collected strains.

The cytochrome P450 Cyp6t3 is not required for ecdysone biosynthesis in Drosophila melanogaster

The steroid hormone 20-hydroxyecdysone (20E) is essential for proper development and the timing of intermediary stage transitions in insects. As a result, there is intense interest in identifying and defining the roles of the enzymes and signaling pathways that regulate 20E production in the prothoracic gland (PG), the major endocrine organ of juvenile insect phases. Transcriptomics is one powerful tool that has been used to identify novel genes that are up- or down-regulated in the PG which may contribute to 20E regulation. Additional functional characterization of putative regulatory candidate genes typically involves qRT-PCR and/or RNAi mediated knockdown of the candidate mRNA in the PG to assess whether the gene's expression shows temporal regulation in the PG and whether its expression is essential for proper 20E production and the correct timing of developmental transitions. While these methods have proved fruitful for identifying novel regulators of 20E production, characterizing the null phenotype of putative regulatory genes is the gold standard for assigning gene function since RNAi is known to generate various types of "off target" effects. Here we describe the genetic null mutant phenotype of the Drosophila melanogaster Cyp6t3 gene . Cyp6t3 was originally identified as a differentially regulated gene in a PG microarray screen and assigned a place in the "Black Box" step of the E biosynthetic pathway based on RNAi mediated knockdown phenotypes and rescue experiments involving feeding of various intermediate compounds of the E biosynthetic pathway. In contrast, we find that Crispr generated null mutations in Cyp6t3 are viable and have normal developmental timing. Therefore, we conclude that Cyp6t3 is not required for E production under typical lab growth conditions and therefore is not an obligate enzymatic component of the Black Box.

Involvement of CYP2 and mitochondrial clan P450s of Helicoverpa armigera in xenobiotic metabolism

Insect CYP2 and mitochondrial clan P450s are relatively conserved genes encoding enzymes generally thought to be involved in biosynthesis or metabolism of endobiotics. However, emerging evidence argues they have potential roles in chemical defense as well, but their actual detoxification functions remain largely unknown. Here, we focused on the full complement of 8 CYP2 and 10 mitochondrial P450s in the generalist herbivore, Helicoverpa armigera. Their varied spatiotemporal expression profiles were analyzed and reflected their specific functions. For functional study of the mitochondrial clan P450s, the redox partners, adrenodoxin reductase (AdR) and adrenodoxin (Adx), were identified from genomes of eight insects and an efficient in vitro electron transfer system of mitochondrial P450 was established by co-expression with Adx and AdR of H. armigera. All CYP2 clan P450s and 8 mitochondrial P450s were successfully expressed in Sf9 cells and compared functionally. In vitro metabolism assays showed that two CYP2 clan P450s (CYP305B1 and CYP18A1) and CYP333B3 (mito clan) could epoxidize aldrin to dieldrin, while CYP305B1 and CYP339A1 (mito clan) have limited but significant hydroxylation capacities to esfenvalerate. CYP303A1 of the CYP2 clan exhibits high metabolic efficiency to 2-tridecanone. Screening the xenobiotic metabolism competence of CYP2 and mitochondrial clan P450s not only provides new insights on insect chemical defense but also can give indications on their physiological functions in H. armigera and other insects.

A Population Genomic Assessment of Three Decades of Evolution in a Natural Drosophila Population

Population genetics seeks to illuminate the forces shaping genetic variation, often based on a single snapshot of genomic variation. However, utilizing multiple sampling times to study changes in allele frequencies can help clarify the relative roles of neutral and non-neutral forces on short time scales. This study compares whole-genome sequence variation of recently collected natural population samples of Drosophila melanogaster against a collection made approximately 35 years prior from the same locality—encompassing roughly 500 generations of evolution. The allele frequency changes between these time points would suggest a relatively small local effective population size on the order of 10,000, significantly smaller than the global effective population size of the species. Some loci display stronger allele frequency changes than would be expected anywhere in the genome under neutrality—most notably the tandem paralogs Cyp6a17 and Cyp6a23, which are impacted by structural variation associated with resistance to pyrethroid insecticides. We find a genome-wide excess of outliers for high genetic differentiation between old and new samples, but a larger number of adaptation targets may have affected SNP-level differentiation versus window differentiation. We also find evidence for strengthening latitudinal allele frequency clines: northern-associated alleles have increased in frequency by an average of nearly 2.5% at SNPs previously identified as clinal outliers, but no such pattern is observed at random SNPs. This project underscores the scientific potential of using multiple sampling time points to investigate how evolution operates in natural populations, by quantifying how genetic variation has changed over ecologically relevant timescales.

The Birth-and-Death Evolution of Cytochrome P450 Genes in Bees

The birth-and-death model of multigene family evolution describes how gene families evolve and diversify through duplication and deletion. The cytochrome P450s are one of the most diverse and well-studied multigene families, involved in both physiological and xenobiotic functions. Extensive studies of insect P450 genes have demonstrated their role in insecticide resistance. Bees are thought to experience toxin exposure through their diet of nectar and pollen, as well as the resin-collecting behavior exhibited by some species. Here, we describe the repertoire of P450 genes in the orchid bee Euglossa dilemma. Male orchid bees form perfume bouquets used in courtship displays by collecting volatile compounds, resulting in exposure to compounds known to be toxic. In addition, we conducted phylogenetic and selection analyses across ten bee species encompassing three bee families. We find that social behavior and resin collection are not correlated with the repertoire of P450 present in a bee species. However, our analyses revealed that P450 clades can be classified as stable and unstable, and that genes involved in xenobiotic metabolism are more likely to belong to unstable clades. Furthermore, we find that unstable clades are under more dynamic evolutionary pressures and exhibit signals of adaptive evolution. This work highlights the complexity of multigene family evolution, revealing that multiple factors contribute to the diversification, stability, and dynamics of this gene family. Furthermore, we provide a resource for future detailed studies investigating the function of different P450s in economically important bee species.

De novo transcriptome assemblies of red king crab (Paralithodes camtschaticus) and snow crab (Chionoecetes opilio) molting gland and eyestalk ganglia - Temperature effects on expression of molting and growth regulatory genes in adult red king crab

Red king crab (Paralithodes camtschaticus) and snow crab (Chionoecetes opilio) are deep-sea crustaceans widely distributed in the North Pacific and Northwest Atlantic Oceans. These giant predators have invaded the Barents Sea over the past decades, and climate-driven temperature changes may influence their distribution and abundance in the sub-Arctic region. Molting and growth in crustaceans are strongly affected by temperature, but the underlying molecular mechanisms are little known, particularly in cold-water species. Here, we describe multiple regulatory factors in the two high-latitude crabs by developing de novo transcriptomes from the molting gland (Y-organ or YO) and eye stalk ganglia (ESG), in addition to the hepatopancreas and claw muscle of red king crab. The Halloween genes encoding the ecdysteroidogenic enzymes were expressed in YO, and the ESG contained multiple neuropeptides, including molt-inhibiting hormone (MIH), crustacean hyperglycemic hormone (CHH), and ion-transport peptide (ITP). Both crabs expressed a diversity of growth-related factors, such as mTOR, AKT, Rheb and AMPKα, and stress-responsive factors, including multiple heat shock proteins (HSPs). Temperature effects on the expression of key regulatory genes were quantified by qPCR in adult red king crab males kept at 4 °C or 10 °C for two weeks during intermolt. The Halloween genes tended to be upregulated in YO at high temperature, while the ecdysteroid receptor and several growth regulators showed tissue-specific responses to elevated temperature. Constitutive and heat-inducible HSPs were expressed in an inverse temperature-dependent manner, suggesting that adult red king crabs can acclimate to increased water temperatures.

Identification and characterization of CYPs induced in the Drosophila antenna by exposure to a plant odorant

Members of the cytochrome p450 (CYP) enzyme family are abundantly expressed in insect olfactory tissues, where they are thought to act as Odorant Degrading Enzymes (ODEs). However, their contribution to olfactory signaling in vivo is poorly understood. This is due in part to the challenge of identifying which of the dozens of antennal-expressed CYPs might inactivate a given odorant. Here, we tested a high-throughput deorphanization strategy in Drosophila to identify CYPs that are transcriptionally induced by exposure to odorants. We discovered three CYPs selectively upregulated by geranyl acetate using transcriptional profiling. Although these CYPs are broadly expressed in the antenna in non-neuronal cells, electrophysiological recordings from CYP mutants did not reveal any changes in olfactory neuron responses to this odorant. Neurons were desensitized by pre-exposing flies to the odorant, but this effect was similar in CYP mutants. Together, our data suggest that the induction of a CYP gene by an odorant does not necessarily indicate a role for that CYP in neuronal responses to that odorant. We go on to show that some CYPs have highly restricted expression patterns in the antenna, and suggest that such CYPs may be useful candidates for further studies on olfactory CYP function.

Parallel and Population-specific Gene Regulatory Evolution in Cold-Adapted Fly Populations

Changes in gene regulation at multiple levels may comprise an important share of the molecular changes underlying adaptive evolution in nature. However, few studies have assayed within- and between-population variation in gene regulatory traits at a transcriptomic scale, and therefore inferences about the characteristics of adaptive regulatory changes have been elusive. Here, we assess quantitative trait differentiation in gene expression levels and alternative splicing (intron usage) between three closely related pairs of natural populations of Drosophila melanogaster from contrasting thermal environments that reflect three separate instances of cold tolerance evolution. The cold-adapted populations were known to show population genetic evidence for parallel evolution at the SNP level, and here we find evidence for parallel expression evolution between them, with stronger parallelism at larval and adult stages than for pupae. We also implement a flexible method to estimate cis- vs trans-encoded contributions to expression or splicing differences at the adult stage. The apparent contributions of cis- vs trans-regulation to adaptive evolution vary substantially among population pairs. While two of three population pairs show a greater enrichment of cis-regulatory differences among adaptation candidates, trans-regulatory differences are more likely to be implicated in parallel expression changes between population pairs. Genes with significant cis-effects are enriched for signals of elevated genetic differentiation between cold- and warm-adapted populations, suggesting that they are potential targets of local adaptation. These findings expand our knowledge of adaptive gene regulatory evolution and our ability to make inferences about this important and widespread process.

Functional Diversification, Redundancy, and Epistasis among Paralogs of the Drosophila melanogaster Obp50a-d Gene Cluster

Large multigene families, such as the insect odorant-binding proteins (OBPs), are thought to arise through functional diversification after repeated gene duplications. Whereas many OBPs function in chemoreception, members of this family are also expressed in tissues outside chemosensory organs. Paralogs of the Obp50 gene cluster are expressed in metabolic and male reproductive tissues, but their functions and interrelationships remain unknown. Here, we report the genetic dissection of four members of the Obp50 cluster, which are in close physical proximity without intervening genes. We used CRISPR technology to excise the entire cluster while introducing a PhiC31 reintegration site to reinsert constructs in which different combinations of the constituent Obp genes were either intact or rendered inactive. We performed whole transcriptome sequencing and assessed sexually dimorphic changes in transcript abundances (transcriptional niches) associated with each gene-edited genotype. Using this approach, we were able to estimate redundancy, additivity, diversification, and epistasis among Obp50 paralogs. We analyzed the effects of gene editing of this cluster on organismal phenotypes and found a significant skewing of sex ratios attributable to Obp50a, and sex-specific effects on starvation stress resistance attributable to Obp50d. Thus, there is functional diversification within the Obp50 cluster with Obp50a contributing to development and Obp50d to stress resistance. The deletion-reinsertion approach we applied to the Obp50 cluster provides a general paradigm for the genetic dissection of paralogs of multigene families.

Signatures of selection and evolutionary relevance of cytochrome P450s in plant-insect interactions

Enzymes in the cytochrome P450 (P450) superfamily have important functions ranging from those that are essential for the physiology and development of the individual to those that mediate interactions between individuals and their biotic environment. Until recently the study of P450s had focused on single functions, substrates, or pathways. Recent advances in sequencing, genome assembly, and phylogenetic methods have returned emphasis to the adaptive value of these enzymes in the context of herbivory. Comparisons of whole repertoires of P450s across related species reveal that P450s capable of metabolizing xenobiotics have an increased rate of gains compared to losses after gene duplications. In plants, studies have focused on enzymes and end-functions that have converged to provide increased resistance to herbivory. This review summarizes the latest findings related to the ecological value of P450s in the interactions between phytophagous insects and their host plants.

Diversity and evolution of the P450 family in arthropods

The P450 family (CYP genes) of arthropods encodes diverse enzymes involved in the metabolism of foreign compounds and in essential endocrine or ecophysiological functions. The P450 sequences (CYPome) from 40 arthropod species were manually curated, including 31 complete CYPomes, and a maximum likelihood phylogeny of nearly 3000 sequences is presented. Arthropod CYPomes are assembled from members of six CYP clans of variable size, the CYP2, CYP3, CYP4 and mitochondrial clans, as well as the CYP20 and CYP16 clans that are not found in Neoptera. CYPome sizes vary from two dozen genes in some parasitic species to over 200 in species as diverse as collembolans or ticks. CYPomes are comprised of few CYP families with many genes and many CYP families with few genes, and this distribution is the result of dynamic birth and death processes. Lineage-specific expansions or blooms are found throughout the phylogeny and often result in genomic clusters that appear to form a reservoir of catalytic diversity maintained as heritable units. Among the many P450s with physiological functions, six CYP families are involved in ecdysteroid metabolism. However, five so-called Halloween genes are not universally represented and do not constitute the unique pathway of ecdysteroid biosynthesis. The diversity of arthropod CYPomes has only partially been uncovered to date and many P450s with physiological functions regulating the synthesis and degradation of endogenous signal molecules (including ecdysteroids) and semiochemicals (including pheromones and defense chemicals) remain to be discovered. Sequence diversity of arthropod P450s is extreme, and P450 sequences lacking the universally conserved Cys ligand to the heme have evolved several times. A better understanding of P450 evolution is needed to discern the relative contributions of stochastic processes and adaptive processes in shaping the size and diversity of CYPomes.

Genomic and transcriptomic analyses in Drosophila suggest that the ecdysteroid kinase-like (EcKL) gene family encodes the 'detoxification-by-phosphorylation' enzymes of insects

Phosphorylation is a phase II detoxification reaction that, among animals, occurs near exclusively in insects, but the enzymes responsible have never been cloned or otherwise identified. We propose the hypothesis that members of the arthropod-specific ecdysteroid kinase-like (EcKL) gene family encode detoxicative kinases. To test this hypothesis, we annotated the EcKL gene family in 12 species of Drosophila and explored their evolution within the genus. Many ancestral EcKL clades are evolutionarily unstable and have experienced repeated gene gain and loss events, while others are conserved as single-copy orthologs. Leveraging multiple published gene expression datasets from D. melanogaster, and using the cytochrome P450s-a classical detoxification family-as a test case, we demonstrate relationships between xenobiotic induction, detoxification tissue-enriched expression and evolutionary instability in the EcKLs and the P450s. We devised a systematic method for identifying candidate detoxification genes in large gene families that is concordant with experimentally determined functions of P450 genes in D. melanogaster. Applying this method to the EcKLs suggested a significant proportion of these genes play roles in detoxification, and that the EcKLs may constitute a detoxification gene family in insects. Additionally, we estimate that between 11 and 16 uncharacterised D. melanogaster P450s are strong detoxification candidates. Lastly, we also found previously unreported genomic and transcriptomic variation in a number of EcKLs and P450s associated with toxic stress phenotypes using a targeted phenome-wide association study (PheWAS) approach in D. melanogaster, presenting multiple future avenues of research for detoxification genetics in this species.

Detoxification Genes Differ Between Cactus-, Fruit-, and Flower-Feeding Drosophila

We use annotated genomes of 14 Drosophila species covering diverse host use phenotypes to test whether 4 gene families that often have detoxification functions are associated with host shifts among species. Bark, slime flux, flower, and generalist necrotic fruit-feeding species all have similar numbers of carboxyl/cholinesterase, glutathione S-transferase, cytochrome P450, and UDP-glucuronosyltransferase genes. However, species feeding on toxic Morinda citrifolia fruit and the fresh fruit-feeding Drosophila suzukii have about 30 and 60 more, respectively. ABC transporters show a different pattern, with the flower-feeding D. elegans and the generalist necrotic fruit and cactus feeder D. hydei having about 20 and >100 more than the other species, respectively. Surprisingly, despite the complex secondary chemistry we find that 3 cactophilic specialists in the mojavensis species cluster have variably fewer genes than any of the other species across all 4 families. We also find 82 positive selection events across the 4 families, with the terminal D. suzukii and M. citrifolia-feeding D. sechellia branches again having the highest number of such events in proportion to their respective branch lengths. Many of the genes involved in these host-use-specific gene number differences or positive selection events lie in specific clades of the gene families that have been recurrently associated with detoxification. Several genes are also found to be involved in multiple duplication and/or positive selection events across the species studied regardless of their host use phenotypes; the most frequently involved are the ABC transporter CG1718, which is not in a specific clade associated with detoxification, and the α-esterase gene cluster, which is.

Chemosensation and Evolution of Drosophila Host Plant Selection

The ability to respond to chemosensory cues is critical for survival of most organisms. Among insects, Drosophila melanogaster has the best characterized olfactory system, and the availability of genome sequences of 30 Drosophila species provides an ideal scenario for studies on evolution of chemosensation. Gene duplications of chemoreceptor genes allow for functional diversification of the rapidly evolving chemoreceptor repertoire. Although some species of the genus Drosophila are generalists for host plant selection, rapid evolution of olfactory receptors, gustatory receptors, odorant-binding proteins, and cytochrome P450s has enabled diverse host specializations of different members of the genus. Here, I review diversification of the chemoreceptor repertoire among members of the genus Drosophila along with co-evolution of detoxification mechanisms that may have enabled occupation of diverse host plant ecological niches.

Birth-and-Death Evolution of the Fatty Acyl-CoA Reductase (FAR) Gene Family and Diversification of Cuticular Hydrocarbon Synthesis in Drosophila

The birth-and-death evolutionary model proposes that some members of a multigene family are phylogenetically stable and persist as a single copy over time, whereas other members are phylogenetically unstable and undergo frequent duplication and loss. Functional studies suggest that stable genes are likely to encode essential functions, whereas rapidly evolving genes reflect phenotypic differences in traits that diverge rapidly among species. One such class of rapidly diverging traits are insect cuticular hydrocarbons (CHCs), which play dual roles in chemical communications as short-range recognition pheromones as well as protecting the insect from desiccation. Insect CHCs diverge rapidly between related species leading to ecological adaptation and/or reproductive isolation. Because the CHC and essential fatty acid biosynthetic pathways share common genes, we hypothesized that genes involved in the synthesis of CHCs would be evolutionary unstable, whereas those involved in fatty acid-associated essential functions would be evolutionary stable. To test this hypothesis, we investigated the evolutionary history of the fatty acyl-CoA reductases (FARs) gene family that encodes enzymes in CHC synthesis. We compiled a unique data set of 200 FAR proteins across 12 Drosophila species. We uncovered a broad diversity in FAR content which is generated by gene duplications, subsequent gene losses, and alternative splicing. We also show that FARs expressed in oenocytes and presumably involved in CHC synthesis are more unstable than FARs from other tissues. Taken together, our study provides empirical evidence that a comparative approach investigating the birth-and-death evolution of gene families can identify candidate genes involved in rapidly diverging traits between species.

Transcriptome-Based Identification and Molecular Evolution of the Cytochrome P450 Genes and Expression Profiling under Dimethoate Treatment in Amur Stickleback (Pungitius sinensis)

Cytochrome P450s (CYPs) are a family of membrane-bound mono-oxygenase proteins, which are involved in cell metabolism and detoxification of various xenobiotic substances. In this study, we identified 58 putative CYP genes in Amur stickleback (Pungitius sinensis) based on the transcriptome sequencing. Conserved motif distribution suggested their functional relevance within each group. Some present recombination events have accelerated the evolution of this gene family. Moreover, a few positive selection sites were identified, which may have accelerated the functional divergence of this family of proteins. Expression patterns of these CYP genes were investigated and indicated that most were affected by dimethoate treatment, suggesting that CYPs were involved in the detoxication of dimethoate. This study will provide a foundation for the further functional investigation of CYP genes in fishes.

Beyond drugs: the evolution of genes involved in human response to medications

The genetic variation of our species reflects human demographic history and adaptation to diverse local environments. Part of this genetic variation affects individual responses to exogenous substances, such as food, pollutants and drugs, and plays an important role in drug efficacy and safety. This review provides a synthesis of the evolution of loci implicated in human pharmacological response and metabolism, interpreted within the theoretical framework of population genetics and molecular evolution. In particular, I review and discuss key evolutionary aspects of different pharmacogenes in humans and other species, such as the relationship between the type of substrates and rate of evolution; the selective pressure exerted by landscape variables or dietary habits; expected and observed patterns of rare genetic variation. Finally, I discuss how this knowledge can be translated directly or after the implementation of specific studies, into practical guidelines.

Cytochrome P450s Cyp4p1 and Cyp4p2 associated with the DDT tolerance in the Drosophila melanogaster strain 91-R

Cytochrome P450s are part of a super-gene family that has undergone gene duplication, divergence, over-expression and, in some cases, loss of function. One such case is the 91-R and 91-C strains of common origin, in Drosophila melanogaster, whereby 91-R (DDT resistant strain) overexpresses Cyp4p1 and Cyp4p2 and both genes are lost in 91-C (DDT susceptible strain). In this study, we used a comparative approach to demonstrate that transcription of Cyp4p1 and Cyp4p2 were constitutively up-regulated in the Drosophila melanogaster strain 91-R as compared to another DDT susceptible strain Canton-S which does not have a loss of function of these genes. Furthermore, significantly increased expression of Cyp4p1 and Cyp4p2 was induced in 91-R in response to sublethal DDT exposure, however, such induction did not occur in the DDT treated Canton-S. Additionally, fixed nucleotide variation within putative transcription factor binding sites of Cyp4p1 and Cyp4p2 promoters were observed between 91-R and Canton-S, however, their impact on transcription remains to be determined. Two GAL4/UAS transgenic strains with integrated heat shock-inducible Cyp4p1- or Cyp4p2-RNAi constructs within wild-type genetic backgrounds were developed. Following heat shock induction of Cyp4p1 and Cyp4p2 knockdown, these transgenic lines showed increased DDT mortality as compared to their corresponding non-heat shock controls. These results provide a functional link of Cyp4p1 and Cyp4p2 in conferring tolerance to DDT exposure.

Candidate genes involved in spiroacetal biosynthesis in the oriental fruit fly, Bactrocera dorsalis

Spiroacetals are widespread in nature as components of volatile semiochemical secretions from many insect species. The general pathway for spiroacetal biosynthesis in Bactrocera sp. is preliminarily established, but many genes involved in this pathway remain to be characterized. By analyzing transcriptomes of the rectal glands (RGs) from immature and mature females of the oriental fruit fly, Bactrocera dorsalis, a set of genes encoding two acetyl-CoA carboxylases (ACCs), two fatty acid synthases (FASs), eight desaturases (DESs), twelve fatty acyl-CoA reductases (FARs), seventy-two cytochrome P450 enzymes (CYPs), and twenty-three odorant binding proteins (OBPs) were identified. We investigated the expression of candidate genes in immature and mature stages based on the RNA-seq data and Real-time quantitative PCR. Expression profiling revealed that some of these genes were primarily expressed in female rectal glands among different tissues, and were up-regulated in mature females. Semi-quantitative RT-PCR assays were also adapted to examine tissue-specific expression of selected candidate genes. Additionally, their putative functions in spiroacetal synthesis and transportation are proposed. Our study provided large-scale sequence information for further functional studies on spiroacetal biosynthetic pathways.

Cis- and trans-acting variants contribute to survivorship in a naïve Drosophila melanogaster population exposed to ryanoid insecticides

Around the world insecticides are being deregistered and banned, as their environmental costs are deemed too great or their efficacy against pest insects is reduced through the evolution of insecticide resistance. With the introduction of replacement insecticides comes the responsibility to assess the way new insecticides perturb various levels of biological systems, from insect physiology to ecosystems. We used a systems genetics approach to identify genetic variants affecting survivorship of Drosophila melanogaster exposed to chlorantraniliprole. The study population was completely naïve to this insecticide chemistry and yet we find associations with variants in neuromuscular genes and coregulated detoxification genes. We predict that these variants will increase in populations of this “sentinel species” as these insecticides are applied in the environment.

Insecticide resistance is a paradigm of microevolution, and insecticides are responsible for the strongest cases of recent selection in the genome of Drosophila melanogaster. Here we use a naïve population and a novel insecticide class to examine the ab initio genetic architecture of a potential selective response. Genome-wide association studies (GWAS) of chlorantraniliprole susceptibility reveal variation in a gene of major effect, Stretchin Myosin light chain kinase (Strn-Mlck), which we validate with linkage mapping and transgenic manipulation of gene expression. We propose that allelic variation in Strn-Mlck alters sensitivity to the calcium depletion attributable to chlorantraniliprole’s mode of action. GWAS also reveal a network of genes involved in neuromuscular biology. In contrast, phenotype to transcriptome associations identify differences in constitutive levels of multiple transcripts regulated by cnc, the homolog of mammalian Nrf2. This suggests that genetic variation acts in trans to regulate multiple metabolic enzymes in this pathway. The most outstanding association is with the transcription level of Cyp12d1 which is also affected in cis by copy number variation. Transgenic overexpression of Cyp12d1 reduces susceptibility to both chlorantraniliprole and the closely related insecticide cyantraniliprole. This systems genetics study reveals multiple allelic variants segregating at intermediate frequency in a population that is completely naïve to this new insecticide chemistry and it foreshadows a selective response among natural populations to these chemicals.

Host use and host shifts in Drosophila

Over a thousand Drosophila species have radiated onto a wide range of feeding and breeding sites. These radiations involve adaptations for locating, accepting, and growing in hosts with highly differing characteristics. In a number of species, owing to the availability of sequenced genomes, particular steps in host specialization and genes that control them, are being identified. Many cases of specialization involve the ability to detoxify some component of the host. Examples include Drosophila sechellia and the octanoic acid in Morinda citrifolia, alpha-amanitin in mycophagous drosophilids, and the alkaloids in cactophilic species. Owing to the known ecologies of many species for which genomes exist, the Drosophila model system provides an unprecedented opportunity to simultaneously examine the genes underlying HOST LOCATION, HOST ACCEPTANCE and HOST USE, the types of selection acting upon them and any coevolutionary interactions among the genes underlying these steps.

The genome of the water strider Gerris buenoi reveals expansions of gene repertoires associated with adaptations to life on the water

BACKGROUND

Having conquered water surfaces worldwide, the semi-aquatic bugs occupy ponds, streams, lakes, mangroves, and even open oceans. The diversity of this group has inspired a range of scientific studies from ecology and evolution to developmental genetics and hydrodynamics of fluid locomotion. However, the lack of a representative water strider genome hinders our ability to more thoroughly investigate the molecular mechanisms underlying the processes of adaptation and diversification within this group.

RESULTS

Here we report the sequencing and manual annotation of the Gerris buenoi (G. buenoi) genome; the first water strider genome to be sequenced thus far. The size of the G. buenoi genome is approximately 1,000 Mb, and this sequencing effort has recovered 20,949 predicted protein-coding genes. Manual annotation uncovered a number of local (tandem and proximal) gene duplications and expansions of gene families known for their importance in a variety of processes associated with morphological and physiological adaptations to a water surface lifestyle. These expansions may affect key processes associated with growth, vision, desiccation resistance, detoxification, olfaction and epigenetic regulation. Strikingly, the G. buenoi genome contains three insulin receptors, suggesting key changes in the rewiring and function of the insulin pathway. Other genomic changes affecting with opsin genes may be associated with wavelength sensitivity shifts in opsins, which is likely to be key in facilitating specific adaptations in vision for diverse water habitats.

CONCLUSIONS

Our findings suggest that local gene duplications might have played an important role during the evolution of water striders. Along with these findings, the sequencing of the G. buenoi genome now provides us the opportunity to pursue exciting research opportunities to further understand the genomic underpinnings of traits associated with the extreme body plan and life history of water striders.

Structural Variants and Selective Sweep Foci Contribute to Insecticide Resistance in the Drosophila Genetic Reference Panel

Patterns of nucleotide polymorphism within populations of Drosophila melanogaster suggest that insecticides have been the selective agents driving the strongest recent bouts of positive selection. However, there is a need to explicitly link selective sweeps to the particular insecticide phenotypes that could plausibly account for the drastic selective responses that are observed in these non-target insects. Here, we screen the Drosophila Genetic Reference Panel with two common insecticides; malathion (an organophosphate) and permethrin (a pyrethroid). Genome-wide association studies map survival on malathion to two of the largest sweeps in the D. melanogaster genome; Ace and Cyp6g1. Malathion survivorship also correlates with lines which have high levels of Cyp12d1, Jheh1 and Jheh2 transcript abundance. Permethrin phenotypes map to the largest cluster of P450 genes in the Drosophila genome, however in contrast to a selective sweep driven by insecticide use, the derived allele seems to be associated with susceptibility. These results underscore previous findings that highlight the importance of structural variation to insecticide phenotypes: Cyp6g1 exhibits copy number variation and transposable element insertions, Cyp12d1 is tandemly duplicated, the Jheh loci are associated with a Bari1 transposable element insertion, and a Cyp6a17 deletion is associated with susceptibility.

Signatures of Insecticide Selection in the Genome of Drosophila melanogaster

Resistance to insecticides has evolved in multiple insect species, leading to increased application rates and even control failures. Understanding the genetic basis of insecticide resistance is fundamental for mitigating its impact on crop production and disease control. We performed a GWAS approach with the Drosophila Genetic Reference Panel (DGRP) to identify the mutations involved in resistance to two widely used classes of insecticides: organophosphates (OPs, parathion) and pyrethroids (deltamethrin). Most variation in parathion resistance was associated with mutations in the target gene Ace, while most variation in deltamethrin resistance was associated with mutations in Cyp6a23, a gene encoding a detoxification enzyme never previously associated with resistance. A "nested GWAS" further revealed the contribution of other loci: Dscam1 and trpl were implicated in resistance to parathion, but only in lines lacking WolbachiaCyp6a17, the paralogous gene of Cyp6a23, and CG7627, an ATP-binding cassette transporter, were implicated in deltamethrin resistance. We observed signatures of recent selective sweeps at all of these resistance loci and confirmed that the soft sweep at Ace is indeed driven by the identified resistance mutations. Analysis of allele frequencies in additional population samples revealed that most resistance mutations are segregating across the globe, but that frequencies can vary substantially among populations. Altogether, our data reveal that the widely used OP and pyrethroid insecticides imposed a strong selection pressure on natural insect populations. However, it remains unclear why, in Drosophila, resistance evolved due to changes in the target site for OPs, but due to a detoxification enzyme for pyrethroids.

Complex gene expansion of the CYP2D gene subfamily

Cytochrome P450 (CYP) superfamily genes encode enzymes that play a role in metabolizing endogenous compounds and in detoxifying exogenous chemicals. The CYP2D subfamily is a member of the CYP2 family, and its gene expansion in herbivores is presumably linked with the need to detoxify abundant plant toxins in the diet, which indicates that CYP2D gene expansion is associated with dietary preferences. To test this hypothesis, the dietary information and CYP2D gene number for 73 vertebrates from different taxonomic groups including 22 mammals, 49 birds, 1 reptile, and 1 amphibian were collected, and correlation analysis and ANOVA were conducted. The results showed that most species (45/73) had only one CYP2D gene, despite their different diets, and dietary preferences were not correlated with CYP2D gene numbers. Specifically, the majority of birds and 7 mammals had only 1 CYP2D gene, and the CYP2D gene number of mammals ranged from 1 to 11, irrespective of their feeding habits. Species with a CYP2D gene number ≥5 included carnivores, herbivores, and omnivores. Furthermore, statistical analyses revealed that no significant correlation existed between dietary preferences and CYP2D gene number, and there was no significant CYP2D gene number variation among species with different dietary preferences, regardless of whether all vertebrates or specific lineages were considered. Furthermore, gene dynamics which indicated by gene duplication events and loss events showed that CYP2D gene number variation had no relationship with diet, suggesting that diet was not a driving force of CYP2D gene expansion and that CYP2D gene expansion was more complex than previously recognized.

Comparative CYP-omic analysis between the DDT-susceptible and -resistant Drosophila melanogaster strains 91-C and 91-R

BACKGROUND

Cytochrome P450 monooxygenases (P450s) are involved in the biosynthesis of endogenous intracellular compounds and the metabolism of xenobiotics, including chemical insecticides. We investigated the structural and expression level variance across all P450 genes with respect to the evolution of insecticide resistance under multigenerational dichlorodiphenyltrichloroethane (DDT) selection.

RESULTS

RNA-sequencing (RNA-seq) and reverse transcriptase-quantitative polymerase chain reaction (RT-qPCR) indicated that the transcript levels of seven P450 genes were significantly up-regulated and three P450 genes were down-regulated in the DDT-resistant strain 91-R, as compared to the control strain 91-C. The overexpression of Cyp6g1 was associated with the presence of an Accord and an HMS-Beagle element insertion in the 5' upstream region in conjunction with copy number variation in the 91-R strain, but not in the 91-C strain. A total of 122 (50.2%) fixed nonsynonymous (amino acid-changing) mutations were found between 91-C and 91-R, and 20 (8.2%) resulted in amino acid changes within functional domains. Three P450 proteins were truncated as a result of premature stop codons and fixed between strains.

CONCLUSION

Our results demonstrate that a combination of changes in P450 protein-coding regions and transcript levels are possibly associated with DDT resistance, and thereby suggest that selection for variant function may occur within this gene family in response to chronic DDT exposure. © 2018 Society of Chemical Industry.

Functional characterization of CYP4G11-a highly conserved enzyme in the western honey bee Apis mellifera

Determining the functionality of CYP4G11, the only CYP4G in the genome of the western honey bee Apis mellifera, can provide insight into its reduced CYP4 inventory. Toward this objective, CYP4G11 transcripts were quantified, and CYP4G11 was expressed as a fusion protein with housefly CPR in Sf9 cells. Transcript levels varied with age, task, and tissue type in a manner consistent with the need for cuticular hydrocarbon production to prevent desiccation or with comb wax production. Young larvae, with minimal need for desiccation protection, expressed CYP4G11 at very low levels. Higher levels were observed in nurses, and even higher levels in wax producers and foragers, the latter of which risk desiccation upon leaving the hive. Recombinant CYP4G11 readily converted octadecanal to n-heptadecane in a time-dependent manner, demonstrating its functions as an oxidative decarbonylase. CYP4G11 expression levels are high in antennae; heterologously expressed CYP4G11 converted tetradecanal to n-tridecane, demonstrating that it metabolizes shorter-chain aldehydes. Together, these findings confirm the involvement of CYP4G11 in cuticular hydrocarbon production and suggest a possible role in clearing pheromonal and phytochemical compounds from antennae. This possible dual functionality of CYP4G11, i.e., cuticular hydrocarbon and comb wax production and antennal odorant clearance, may explain how honey bees function with a reduced CYP4G inventory.

Rapid evolution in insect pests: the importance of space and time in population genomics studies

Pest species in agroecosystems often exhibit patterns of rapid evolution to environmental and human-imposed selection pressures. Although the role of adaptive processes is well accepted, few insect pests have been studied in detail and most research has focused on selection at insecticide resistance candidate genes. Emerging genomic datasets provide opportunities to detect and quantify selection in insect pest populations, and address long-standing questions about mechanisms underlying rapid evolutionary change. We examine the strengths of recent studies that stratify population samples both in space (along environmental gradients and comparing ancestral vs. derived populations) and in time (using chronological sampling, museum specimens and comparative phylogenomics), resulting in critical insights on evolutionary processes, and providing new directions for studying pests in agroecosystems.

Evolution and Expansion of the Prokaryote-Like Lipoxygenase Family in the Brown Alga Saccharina japonica

Lipoxygenase (LOX) plays important roles in fatty acid oxidation and lipid mediator biosynthesis. In this study, we give first insights into brown algal LOX evolution. Whole genome searches revealed four, three, and eleven LOXs in Ectocarpus siliculosus, Cladosiphon okamuranus, and Saccharina japonica, respectively. In phylogenetic analyses, LOXs from brown algae form a robust clade with those from prokaryotes, suggesting an ancestral origin and slow evolution. Brown algal LOXs were divided into two clades, C1 and C2 in a phylogenetic tree. Compared to the two species of Ectocarpales, LOX gene expansion occurred in the kelp S. japonica through tandem duplication and segmental duplication. Selection pressure analysis showed that LOX genes in brown algae have undergone strong purifying selection, while the selective constraint in the C2 clade was more relaxed than that in the C1 clade. Furthermore, within each clade, LOXs of S. japonica evolved under more relaxed selection constraints than E. siliculosus and C. okamuranus. Structural modeling showed that unlike LOXs of plants and animals, which contain a β barrel in the N-terminal part of the protein, LOXs in brown algae fold into a single domain. Analysis of previously published transcriptomic data showed that LOXs in E. siliculosus are responsive to hyposaline, hypersaline, oxidative, and copper stresses. Moreover, clear divergence of expression patterns was observed among different life stages, as well as between duplicate gene pairs. In E. siliculosus, all four LOXs are male-biased in immature gametophytes, and mature gametophytes showed significantly higher LOX mRNA levels than immature gametophytes and sporophytes. In S. japonica, however, our RNA-Seq data showed that most LOXs are highly expressed in sporophytes. Even the most recently duplicated gene pairs showed divergent expression patterns, suggesting that functional divergence has likely occurred since LOX genes duplicated, which potentially contributes to the production of various oxylipins in brown algae.

Duplications and Positive Selection Drive the Evolution of Parasitism-Associated Gene Families in the Nematode Strongyloides papillosus

Gene duplication is a major mechanism playing a role in the evolution of phenotypic complexity and in the generation of novel traits. By comparing parasitic and nonparasitic nematodes, a recent study found that the evolution of parasitism in Strongyloididae is associated with a large expansion in the Astacin and CAP gene families.To gain novel insights into the developmental processes in the sheep parasite Strongyloides papillosus, we sequenced transcriptomes of different developmental stages and sexes. Overall, we found that the majority of genes are developmentally regulated and have one-to-one orthologs in the diverged S. ratti genome. Together with the finding of similar expression profiles between S. papillosus and S. ratti, these results indicate a strong evolutionary constraint acting against change at sequence and expression levels. However, the comparison between parasitic and free-living females demonstrates a quite divergent pattern that is mostly due to the previously mentioned expansion in the Astacin and CAP gene families. More detailed phylogenetic analysis of both gene families shows that most members date back to single expansion events early in the Strongyloides lineage and have undergone subfunctionalization resulting in clusters that are highly expressed either in infective larvae or in parasitic females. Finally, we found increased evidence for positive selection in both gene families relative to the genome-wide expectation.In summary, our study reveals first insights into the developmental transcriptomes of S. papillosus and provides a detailed analysis of sequence and expression evolution in parasitism-associated gene families.

Cytochrome P450 diversification and hostplant utilization patterns in specialist and generalist moths: Birth, death and adaptation

Across insect genomes, the size of the cytochrome P450 monooxygenase (CYP) gene superfamily varies widely. CYPome size variation has been attributed to reciprocal adaptive radiations in insect detoxification genes in response to plant biosynthetic gene radiations driven by co-evolution between herbivores and their chemically defended hostplants. Alternatively, variation in CYPome size may be due to random "birth-and-death" processes, whereby exponential increase via gene duplications is limited by random decay via gene death or transition via divergence. We examined CYPome diversification in the genomes of seven Lepidoptera species varying in host breadth from monophagous (Bombyx mori) to highly polyphagous (Amyelois transitella). CYPome size largely reflects the size of Clan 3, the clan associated with xenobiotic detoxification, and to some extent phylogenetic age. Consistently across genomes, families CYP6, CYP9 and CYP321 are most diverse and CYP6AB, CYP6AE, CYP6B, CYP9A and CYP9G are most diverse among subfamilies. Higher gene number in subfamilies is due to duplications occurring primarily after speciation and specialization ("P450 blooms"), and the genes are arranged in clusters, indicative of active duplicating loci. In the parsnip webworm, Depressaria pastinacella, gene expression levels in large subfamilies are high relative to smaller subfamilies. Functional and phylogenetic data suggest a correlation between highly dynamic loci (reflective of extensive gene duplication, functionalization and in some cases loss) and the ability of enzymes encoded by these genes to metabolize hostplant defences, consistent with an adaptive, nonrandom process driven by ecological interactions.

Partitioning the roles of CYP6G1 and gut microbes in the metabolism of the insecticide imidacloprid in Drosophila melanogaster

Resistance to insecticides through enhanced metabolism is a worldwide problem. The Cyp6g1 gene of the vinegar fly, Drosophila melanogaster, is a paradigm for the study of metabolic resistance. Constitutive overexpression of this gene confers resistance to several classes of insecticides, including the neonicotinoid imidacloprid (IMI). The metabolism of IMI in this species has been previously shown to yield oxidative and nitro-reduced metabolites. While levels of the oxidative metabolites are correlated with CYP6G1 expression, nitro-reduced metabolites are not, raising the question of how these metabolites are produced. Some IMI metabolites are known to be toxic, making their fate within the insect a second question of interest. These questions have been addressed by coupling the genetic tools of gene overexpression and CRISPR gene knock-out with the mass spectrometric technique, the Twin-Ion Method (TIM). Analysing axenic larvae indicated that microbes living within D. melanogaster are largely responsible for the production of the nitro-reduced metabolites. Knock-out of Cyp6g1 revealed functional redundancy, with some metabolites produced by CYP6G1 still detected. IMI metabolism was shown to produce toxic products that are not further metabolized but readily excreted, even when produced in the Central Nervous System (CNS), highlighting the significance of transport and excretion in metabolic resistance.

Naturally Segregating Variation at Ugt86Dd Contributes to Nicotine Resistance in Drosophila melanogaster

Identifying the sequence polymorphisms underlying complex trait variation is a key goal of genetics research, since knowing the precise causative molecular events allows insight into the pathways governing trait variation. Genetic analysis of complex traits in model systems regularly starts by constructing QTL maps, but generally fails to identify causative sequence polymorphisms. Previously we mapped a series of QTL contributing to resistance to nicotine in a Drosophila melanogaster multiparental mapping resource and here use a battery of functional tests to resolve QTL to the molecular level. One large-effect QTL resided over a cluster of UDP-glucuronosyltransferases, and quantitative complementation tests using deficiencies eliminating subsets of these detoxification genes revealed allelic variation impacting resistance. RNAseq showed that Ugt86Dd had significantly higher expression in genotypes that are more resistant to nicotine, and anterior midgut-specific RNA interference (RNAi) of this gene reduced resistance. We discovered a segregating 22-bp frameshift deletion in Ugt86Dd, and accounting for the InDel during mapping largely eliminates the QTL, implying the event explains the bulk of the effect of the mapped locus. CRISPR/Cas9 editing of a relatively resistant genotype to generate lesions in Ugt86Dd that recapitulate the naturally occurring putative loss-of-function allele, leads to a large reduction in resistance. Despite this major effect of the deletion, the allele appears to be very rare in wild-caught populations and likely explains only a small fraction of the natural variation for the trait. Nonetheless, this putatively causative coding InDel can be a launchpad for future mechanistic exploration of xenobiotic detoxification.

Resistance evolution in Drosophila: the case of CYP6G1

The massive use of DDT as an insecticide between 1940 and 1970 has resulted in the emergence of a resistant population of insects. One of the main metabolic mechanisms developed by resistant insects involves detoxification enzymes such as cytochrome P450s. These enzymes can metabolise the insecticide to render it less toxic and facilitate its elimination from the organism. The P450 Cyp6g1 was identified as the major factor responsible for DDT resistance in Drosophila melanogaster field populations. In this article, we review the data available for this gene since it was associated with resistance in 2002. The knowledge gained on Cyp6g1 allows a better understanding of the evolution of insecticide resistance mechanisms and highlights the major role of transposable elements in evolutionary processes. © 2016 Society of Chemical Industry.

Transcriptome Analysis of Drosophila melanogaster Third Instar Larval Ring Glands Points to Novel Functions and Uncovers a Cytochrome p450 Required for Development

In Drosophila melanogaster larvae, the ring gland (RG) is a control center that orchestrates major developmental transitions. It is a composite organ, consisting of the prothoracic gland, the corpus allatum, and the corpora cardiaca, each of which synthesizes and secretes a different hormone. Until now, the RG’s broader developmental roles beyond endocrine secretion have not been explored. RNA sequencing and analysis of a new transcriptome resource from D. melanogaster wandering third instar larval RGs has provided a fascinating insight into the diversity of developmental signaling in this organ. We have found strong enrichment of expression of two gene pathways not previously associated with the RG: immune response and fatty acid metabolism. We have also uncovered strong expression for many uncharacterized genes. Additionally, RNA interference against RG-enriched cytochrome p450s Cyp6u1 and Cyp6g2 produced a lethal ecdysone deficiency and a juvenile hormone deficiency, respectively, flagging a critical role for these genes in hormone synthesis. This transcriptome provides a valuable new resource for investigation of roles played by the RG in governing insect development.

Evolution: Sex, Diet and Red Ketocarotenoids

In many birds, bright red plumage in males helps to attract females, a classic example of sexual selection. Studies in zebra finches and canaries have now identified the gene responsible for red coloration.

RNA interference of cytochrome P450 CYP6F subfamily genes affects susceptibility to different insecticides in Locusta migratoria

BACKGROUND

Many insect cytochrome P450s (CYPs) play critical roles in detoxification of insecticides. The CYP6 family is unique to the class Insecta, and its biochemical function has essentially been associated with the metabolism of xenobiotics. In this study, we sequenced and characterised the full-length cDNAs of five CYP genes from Locusta migratoria, a highly destructive agricultural pest worldwide.

RESULTS

The five genes were predominantly expressed in brain, guts, fat bodies or Malpighian tubules. CYP6FE1, CYP6FF1 and CYP6FG1 were expressed at higher levels in fourth-instar nymphs than in other developmental stages. CYPFD2 is specifically expressed in adults, whereas CYP6FD1, CYP6FD2 and CYP6FE1 showed significantly lower expression in eggs than in other developmental stages. Deltamethrin suppressed CYP6FD1 expression in third-instar nymphs and upregulated the expression level of CYP6FD2, CYP6FF1 and CYP6FG1 at the dose of LD₁₀ . Efficient RNA interference-mediated gene silencing was established for four of the five CYP genes. Silencing of CYP6FF1 increased the nymphal mortality from 23 to 50% in response to deltamethrin. Silencing of CYP6FD2 and CYP6FE1 increased the nymphal mortality from 32 to 72 and 66%, respectively, to carbaryl.

CONCLUSION

Three of the four CYP6F subfamily genes in L. migratoria were associated with the detoxification of deltamethrin or carbaryl. The role of CYPs in insecticide detoxification appears to be both gene and insecticide specific. © 2016 Society of Chemical Industry.

Computational Identification of the Paralogs and Orthologs of Human Cytochrome P450 Superfamily and the Implication in Drug Discovery

The human cytochrome P450 (CYP) superfamily consisting of 57 functional genes is the most important group of Phase I drug metabolizing enzymes that oxidize a large number of xenobiotics and endogenous compounds, including therapeutic drugs and environmental toxicants. The CYP superfamily has been shown to expand itself through gene duplication, and some of them become pseudogenes due to gene mutations. Orthologs and paralogs are homologous genes resulting from speciation or duplication, respectively. To explore the evolutionary and functional relationships of human CYPs, we conducted this bioinformatic study to identify their corresponding paralogs, homologs, and orthologs. The functional implications and implications in drug discovery and evolutionary biology were then discussed. GeneCards and Ensembl were used to identify the paralogs of human CYPs. We have used a panel of online databases to identify the orthologs of human CYP genes: NCBI, Ensembl Compara, GeneCards, OMA ("Orthologous MAtrix") Browser, PATHER, TreeFam, EggNOG, and Roundup. The results show that each human CYP has various numbers of paralogs and orthologs using GeneCards and Ensembl. For example, the paralogs of CYP2A6 include CYP2A7, 2A13, 2B6, 2C8, 2C9, 2C18, 2C19, 2D6, 2E1, 2F1, 2J2, 2R1, 2S1, 2U1, and 2W1; CYP11A1 has 6 paralogs including CYP11B1, 11B2, 24A1, 27A1, 27B1, and 27C1; CYP51A1 has only three paralogs: CYP26A1, 26B1, and 26C1; while CYP20A1 has no paralog. The majority of human CYPs are well conserved from plants, amphibians, fishes, or mammals to humans due to their important functions in physiology and xenobiotic disposition. The data from different approaches are also cross-validated and validated when experimental data are available. These findings facilitate our understanding of the evolutionary relationships and functional implications of the human CYP superfamily in drug discovery.