Origins of P450 diversity

Structure-guided discovery of novel AflG inhibitors for aflatoxin contamination control in aspergillus flavus

Aflatoxins (AFs) are highly carcinogenic metabolites produced by Aspergillus species that can contaminate critical food staples, leading to significant health and economic risks. The cytochrome P450 monooxygenase AflG catalyzes an early step in AF biosynthesis, resulting in the conversion of averantin (AVN) to 5'-hydroxy-averantin. However, the molecular mechanism underlying the AflG-AVN interaction remains unclear. Here, we sought to understand the structural features of AflG in complex with AVN to enable the identification of inhibitors targeting the AflG binding pocket. To achieve this goal, we employed a comprehensive approach combining computational and experimental methods. Structural modeling and microsecond-scale molecular dynamics (MD) simulations yielded new insights into AflG architecture and unveiled unique ligand binding conformations of the AflG-AVN complex. High-throughput virtual screening of more than 1.3 million compounds pinpointed specific subsets with favorable predicted docking scores. The resulting compounds were ranked based on binding free energy calculations and evaluated with MD simulations and in vitro experiments with Aspergillus flavus. Our results revealed two compounds significantly inhibited AF biosynthesis. Comprehensive structural analysis elucidated the binding sites of competitive inhibitors and demonstrated their regulation of AflG dynamics. This structure-guided pipeline successfully enabled the identification of novel AflG inhibitors and provided novel molecular insights that will guide future efforts to develop effective therapeutics that prevent AF contamination.

Sequence similarity networks bear out hierarchical relationships of green cytochrome P450

Land plants have diversified enzyme families. One of the most prominent is the cytochrome P450 (CYP or CYP450) family. With over 443,000 CYP proteins sequenced across the tree of life, CYPs are ubiquitous in archaea, bacteria, and eukaryotes. Here, we focused on land plants and algae to study the role of CYP diversification. CYPs, acting as monooxygenases, catalyze hydroxylation reactions crucial for specialized plant metabolic pathways, including detoxification and phytohormone production; the CYPome consists of one enormous superfamily that is divided into clans and families. Their evolutionary history speaks of high substrate promiscuity; radiation and functional diversification have yielded numerous CYP families. To understand the evolutionary relationships within the CYPs, we employed sequence similarity network analyses. We recovered distinct clusters representing different CYP families, reflecting their diversified sequences that we link to the prediction of functionalities. Hierarchical clustering and phylogenetic analysis further elucidated relationships between CYP clans, uncovering their shared deep evolutionary history. We explored the distribution and diversification of CYP subfamilies across plant and algal lineages, uncovering novel candidates and providing insights into the evolution of these enzyme families. This identified unexpected relationships between CYP families, such as the link between CYP82 and CYP74, shedding light on their roles in plant defense signaling pathways. Our approach provides a methodology that brings insights into the emergence of new functions within the CYP450 family, contributing to the evolutionary history of plants and algae. These insights can be further validated and implemented via experimental setups under various external conditions.

Sex-biased gene expression across mammalian organ development and evolution

Sexually dimorphic traits are common among mammals and are specified during development through the deployment of sex-specific genetic programs. Because little is known about these programs, we investigated them using a resource of gene expression profiles in males and females throughout the development of five organs in five mammals (human, mouse, rat, rabbit, and opossum) and a bird (chicken). We found that sex-biased gene expression varied considerably across organs and species and was often cell-type specific. Sex differences increased abruptly around sexual maturity instead of increasing gradually during organ development. Finally, sex-biased gene expression evolved rapidly at the gene level, with differences between organs in the evolutionary mechanisms used, but more slowly at the cellular level, with the same cell types being sexually dimorphic across species.

Evolution and functional role prediction of the CYP6DE and CYP6DJ subfamilies in Dendroctonus (Curculionidae: Scolytinae) bark beetles

Dendroctonus-bark beetles are natural components and key ecological agents of coniferous forests. They spend most of their lives under the bark, where they are exposed to highly toxic terpenes present in the oleoresin. Cytochrome P450 (CYP) is a multigene family involved in the detoxification of these compounds. It has been demonstrated that CYP6DE and CYP6DJ subfamilies hydroxylate monoterpenes, whose derivatives can act as pheromone synergist compounds or be pheromones themselves in these insects. Given the diversity and functional role of CYPs, we investigated whether these cytochromes have retained their function throughout the evolution of these insects. To test this hypothesis, we performed a Bayesian phylogenetic analysis to determine phylogenetic subgroups of cytochromes in these subfamilies. Subgroups were mapped and reconciled with the Dendroctonus phylogeny. Molecular docking analyses were performed with the cytochromes of each subgroup and enantiomers of α-pinene and β-pinene, (+)-3-carene, β-myrcene and R-(+)-limonene. In addition, functional divergence analysis was performed to identify critical amino acid sites that influence changes in catalytic site conformation and/or protein folding. Three and two phylogenetic subgroups were recovered for the CYP6DE and CYP6DJ subfamilies, respectively. Mapping and reconciliation analysis showed different gain and loss patterns for cytochromes of each subgroup. Functional predictions indicated that the cytochromes analyzed are able to hydroxylate all monoterpenes; however, they showed preferential affinities to different monoterpenes. Functional divergence analyses indicated that the CYP6DE subfamily has experimented type I and II divergence, whereas the CYP6DJ subfamily has evolved under strong functional constraints. Results suggest cytochromes of the CYP6DE subfamily evolve to reinforce their detoxifying capacity hydroxylating mainly α- and β-pinene to (+) and (-)-trans-verbenol, being the negative enantiomer used as a pheromone by several Dendroctonus species; whereas cytochromes of the CYP6DJ subfamily appear to retain their original function related to the detoxification of these compounds.

CYP19A1 is downregulated by BRD4 and suppresses castration-resistant prostate cancer cell invasion and proliferation by decreasing AR expression

Castration-resistant prostate cancer (CRPC) is the final stage of prostate cancer (PCa). As the main androgen in males, testosterone, and its androgen receptor (AR) play an important role in CRPC. The enzyme that catalyzes testosterone, aromatase, can be influenced by CYP19A1 activity - thus possibly affecting both AR expression and CRPC. However, the function of CYP19A1 in CRPC remains unclear. Using data derived from public databases and clinical samples, we analyzed the expression of CYP19A1 in PCa and CRPC specimens. The effect of CYP19A1 on cell invasion and proliferation was investigated in vitro and in vivo; while its function in metabolizing testosterone was detected in vitro. The effect of BRD4 on CYP19A1 and AR was investigated by qRT-PCR and western blot; whereas the effect of JQ1 on cells was assessed based on the IC50 value. We found that CYP19A1 was downregulated in CRPC samples and cells which correlated with a decrease in CRPC cell invasion and proliferation, and an increase in AR expression. Inversely, CYP19A1 affected CRPC cell invasion and proliferation by suppressing the expression of AR which may be attributed to the metabolism of testosterone by CYP19A1. Moreover, the BRD4 inhibitor JQ1 induced the CYP19A1 expression and suppressed the AR expression. Following BRD4 knockdown, CYP19A1 showed higher expression while AR expression was decreased. Our findings demonstrated that CYP19A1 could reduce CRPC cell invasion and proliferation by targeting AR, and this process could be regulated by BRD4. CYP19A1 may be a potential therapeutic target and enhance BRD4 inhibition in treating CRPC.

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.

Identifying biomarkers associated with the diagnosis of ulcerative colitis via bioinformatics and machine learning

BACKGROUND

Ulcerative colitis (UC) is an idiopathic inflammatory disease with an increasing incidence. This study aimed to identify potential UC biomarkers and associated immune infiltration characteristics.

METHODS

Two datasets (GSE87473 and GSE92415) were merged to obtain 193 UC samples and 42 normal samples. Using R, differentially expressed genes (DEGs) between UC and normal samples were filtered out, and their biological functions were investigated using Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses. Promising biomarkers were identified using least absolute shrinkage selector operator regression and support vector machine recursive feature elimination, and their diagnostic efficacy was evaluated through receiver operating characteristic (ROC) curves. Finally, CIBERSORT was used to investigate the immune infiltration characteristics in UC, and the relationship between the identified biomarkers and various immune cells was examined.

RESULTS

We found 102 DEGs, of which 64 were significantly upregulated, and 38 were significantly downregulated. The DEGs were enriched in pathways associated with interleukin-17, cytokine-cytokine receptor interaction and viral protein interactions with cytokines and cytokine receptors, among others. Using machine learning methods and ROC tests, we confirmed DUOX2, DMBT1, CYP2B7P, PITX2 and DEFB1 to be essential diagnostic genes for UC. Immune cell infiltration analysis revealed that all five diagnostic genes were correlated with regulatory T cells, CD8 T cells, activated and resting memory CD4 T cells, activated natural killer cells, neutrophils, activated and resting mast cells, activated and resting dendritic cells and M0, M1 and M2 macrophages.

CONCLUSIONS

DUOX2, DMBT1, CYP2B7P, PITX2 and DEFB1 were identified as prospective biomarkers for UC. A new perspective on understanding the progression of UC may be provided by these biomarkers and their relationship with immune cell infiltration.

Insights into the Effects of Insecticides on Aphids (Hemiptera: Aphididae): Resistance Mechanisms and Molecular Basis

With the passage of time and indiscreet usage of insecticides on crops, aphids are becoming resistant to their effect. The different classes of insecticides, including organophosphates, carbamates, pyrethroids and neonicotinoids, have varied effects on insects. Furthermore, the molecular effects of these insecticides in aphids, including effects on the enzymatic machinery and gene mutation, are resulting in aphid resistance to the insecticides. In this review, we will discuss how aphids are affected by the overuse of pesticides, how resistance appears, and which mechanisms participate in the resistance mechanisms in various aphid species as significant crop pests. Gene expression studies were analyzed using the RNA-Seq technique. The stress-responsive genes were analyzed, and their expression in response to insecticide administration was determined. Putative insecticide resistance-related genes, cytochrome P450, glutathione S-transferase, carboxylesterase CarEs, ABC transporters, cuticle protein genes, and trypsin-related genes were studied. The review concluded that if insecticide-susceptible aphids interact with ample dosages of insecticides with sublethal effects, this will result in the upregulation of genes whose primary role is to detoxify insecticides. In the past decade, certain advancements have been observed regarding insecticide resistance on a molecular basis. Even so, not much is known about how aphids detoxify the insecticides at molecular level. Thus, to attain equilibrium, it is important to observe the manipulation of pest and insect species with the aim of restoring susceptibility to insecticides. For this purpose, this review has included critical insights into insecticide resistance in aphids.

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.

Promiscuity, a Driver of Plant Cytochrome P450 Evolution?

Plant cytochrome P450 monooxygenases were long considered to be highly substrate-specific, regioselective and stereoselective enzymes, in this respect differing from their animal counterparts. The functional data that have recently accumulated clearly counter this initial dogma. Highly promiscuous P450 enzymes have now been reported, mainly in terpenoid pathways with functions in plant adaptation, but also some very versatile xenobiotic/herbicide metabolizers. An overlap and predictable interference between endogenous and herbicide metabolism are starting to emerge. Both substrate preference and permissiveness vary between plant P450 families, with high promiscuity seemingly favoring retention of gene duplicates and evolutionary blooms. Yet significant promiscuity can also be observed in the families under high negative selection and with essential functions, usually enhanced after gene duplication. The strategies so far implemented, to systematically explore P450 catalytic capacity, are described and discussed.

An efficient genotyper and star-allele caller for pharmacogenomics

High-throughput sequencing provides sufficient means for determining genotypes of clinically important pharmacogenes that can be used to tailor medical decisions to individual patients. However, pharmacogene genotyping, also known as star-allele calling, is a challenging problem that requires accurate copy number calling, structural variation identification, variant calling, and phasing within each pharmacogene copy present in the sample. Here we introduce Aldy 4, a fast and efficient tool for genotyping pharmacogenes that uses combinatorial optimization for accurate star-allele calling across different sequencing technologies. Aldy 4 adds support for long reads and uses a novel phasing model and improved copy number and variant calling models. We compare Aldy 4 against the current state-of-the-art star-allele callers on a large and diverse set of samples and genes sequenced by various sequencing technologies, such as whole-genome and targeted Illumina sequencing, barcoded 10x Genomics, and Pacific Biosciences (PacBio) HiFi. We show that Aldy 4 is the most accurate star-allele caller with near-perfect accuracy in all evaluated contexts, and hope that Aldy remains an invaluable tool in the clinical toolbox even with the advent of long-read sequencing technologies.

Analysis of Cyp51 protein sequences shows 4 major Cyp51 gene family groups across fungi

Azole drugs target fungal sterol biosynthesis and are used to treat millions of human fungal infections each year. Resistance to azole drugs has emerged in multiple fungal pathogens including Candida albicans, Cryptococcus neoformans, Histoplasma capsulatum, and Aspergillus fumigatus. The most well-studied resistance mechanism in A. fumigatus arises from missense mutations in the coding sequence combined with a tandem repeat in the promoter of cyp51A, which encodes a cytochrome P450 enzyme in the fungal sterol biosynthesis pathway. Filamentous members of Ascomycota such as A. fumigatus have either 1 or 2 of 3 Cyp51 paralogs (Cyp51A, Cyp51B, and Cyp51C). Most previous research in A. fumigatus has focused on Cyp51A due to its role in azole resistance. We used the A. fumigatus Cyp51A protein sequence as the query in database searches to identify Cyp51 proteins across fungi. We found 435 Cyp51 proteins in 295 species spanning from early-diverging fungi (Blastocladiomycota, Chytridiomycota, Zoopagomycota, and Mucormycota) to late-diverging fungi (Ascomycota and Basidiomycota). We found these sequences formed 4 major Cyp51 groups: Cyp51, Cyp51A, Cyp51B, and Cyp51C. Surprisingly, we found all filamentous Ascomycota had a Cyp51B paralog, while only 50% had a Cyp51A paralog. We created maximum likelihood trees to investigate the evolution of Cyp51 in fungi. Our results suggest Cyp51 is present in all fungi with 3 paralogs emerging in Pezizomycotina, including Cyp51C which appears to have diverged from the progenitor of the Cyp51A and Cyp51B groups.

Downregulation of NADPH-cytochrome P450 reductase via RNA interference increases the susceptibility of Acyrthosiphon pisum to desiccation and insecticides

Nicotinamide adenine dinucleotide phosphate (NADPH)-cytochrome P450 reductase (CPR) is involved in the metabolism of endogenous and exogenous substances, and detoxification of insecticides. RNA interference (RNAi) of CPR in certain insects causes developmental defects and enhanced susceptibility to insecticides. However, the CPR of Acyrthosiphon pisum has not been characterized, and its function is still not understood. In this study, we investigated the biochemical functions of A. pisum CPR (ApCPR). ApCPR was found to be transcribed in all developmental stages and was abundant in the embryo stage, and in the gut, head, and abdominal cuticle. After optimizing the dose and silencing duration of RNAi for downregulating ApCPR, we found that ApCPR suppression resulted in a significant decrease in the production of cuticular and internal hydrocarbon contents, and of cuticular waxy coatings. Deficiency in cuticular hydrocarbons (CHCs) decreased the survival rate of A. pisum under desiccation stress and increased its susceptibility to contact insecticides. Moreover, desiccation stress induced a significant increase in ApCPR mRNA levels. We further confirmed that ApCPR participates in CHC production. These results indicate that ApCPR modulates CHC production, desiccation tolerance, and insecticide susceptibility in A. pisum, and presents a novel target for pest control.

Single amino acid variations drive functional divergence of cytochrome P450s in Helicoverpa species

Divergence of gene function is a hallmark of evolution, but assessing such divergence in one species or between species requires information on functional alterations of the alleles and homologs. Here, we explore the functional divergence of two paralogs, CYP6AE19 and CYP6AE20, from Helicoverpa armigera, and two close orthologs, CYP6B8 and CYP6B7, from two related species (Helicoverpa zea and H. armigera); although there is high sequence identity within each pair of enzymes, the latter P450 of each pair has lost metabolic competence towards the plant allelochemical xanthotoxin. Multiple chimeric and single/double site mutants were created by exchanging the diverse substrate recognition sites (SRSs) and amino acids within each pair of P450s. Heterologous expression in Sf9 cells and in vitro metabolism studies showed that the exchange of SRS4 swapped the activity of CYP6AE19 and CYP6AE20, and subsequent site-directed mutagenesis demonstrated that the CYP6AE20 V318M substitution causes a gain-of-function towards xanthotoxin. Meanwhile, a single amino acid substitution (L489P) in SRS6 was found to swap activity between the CYP6B orthologs. Sequence alignments of CYP6AE paralogs and all reported insect xanthotoxin-metabolizing P450s suggest M318 and P489 are essential for the catalytic activities of CYP6AE paralogs and CYP6B orthologs, respectively, but P450s in different subfamilies may have different mechanisms towards the same substrate. Our findings demonstrate that a single amino acid substitution can suffice to alter substrate metabolism and this functional divergence resulting from natural mutations will help to further our understanding of the process of natural selection of P450 genes and their role in insect-host plant interactions.

Evaluation of biotransformation capacity of transplastomic plants and hairy roots of Nicotiana tabacum expressing human cytochrome P450 2D6

Cytochrome P450 monooxygenases (CYPs) are important tools for regio- and stereoselective oxidation of target molecules or engineering of metabolic pathways. Functional heterologous expression of eukaryotic CYPs is often problematic due to their dependency on the specific redox partner and the necessity of correct association with the membranes for displaying enzymatic activity. Plant hosts offer advantages of accessibility of reducing partners and a choice of membranes to insert heterologous CYPs. For the evaluation of plant systems for efficient CYP expression, we established transplastomic plants and hairy root cultures of Nicotiana tabacum carrying the gene encoding human CYP2D6 with broad substrate specificity. The levels of CYP2D6 transcript accumulation and enzymatic activity were estimated and compared with the data of CYP2D6 transient expression in N. benthamiana. The relative level of CYP2D6 transcripts in transplastomic plants was 2-3 orders of magnitude higher of that observed after constitutive or transient expression from the nucleus. CYP2D6 expressed in chloroplasts converted exogenous synthetic substrate loratadine without the need for co-expression of the cognate CYP reductase. The loratadine conversion rate in transplastomic plants was comparable to that in N. benthamiana plants transiently expressing a chloroplast targeted CYP2D6 from the nucleus, but was lower than the value reported for transiently expressed CYP2D6 with the native endoplasmic reticulum signal-anchor sequence. Hairy roots showed the lowest substrate conversion rate, but demonstrated the ability to release the product into the culture medium. The obtained results illustrate the potential of plant-based expression systems for exploiting the enzymatic activities of eukaryotic CYPs with broad substrate specificities.

Simple and Robust Detection of CYP2D6 Gene Deletions and Duplications Using CYP2D8P as Reference

Genotyping of the CYP2D6 gene is the most commonly applied pharmacogenetic test globally. Significant economic interests have led to the development of a plurality of assays, available for almost any genotyping platform or DNA detection chemistry. Of all the genetic variants, copy number variations are particular difficult to detect by polymerase chain reaction. Here, we present two simple novel approaches for the identification of samples carrying either deletions or duplications of the CYP2D6 gene; by relative quantification using a singleplex 5′nuclease real-time PCR assay, and by high-resolution melting of PCR products. These methods make use of universal primers, targeting both the CYP2D6 and the reference gene CYP2D8P, which is necessary for the analysis. The assays were validated against a reference method using a large set of samples. The singleplex nature of the 5′nuclease real-time PCR ensures that the primers anneal with equal affinity to both the sequence of the CYP2D6 and the reference gene. This facilitates robust identification of gene deletions and duplications based on the cycle threshold value. In contrast, the high-resolution melting assay is an end-point PCR, where the identification relies on variations between the amount of product generated from each of the two genes.

The Role of Cytochrome P450s in Insect Toxicology and Resistance

Insect cytochrome P450 monooxygenases (P450s) perform a variety of important physiological functions, but it is their role in the detoxification of xenobiotics, such as natural and synthetic insecticides, that is the topic of this review. Recent advances in insect genomics and postgenomic functional approaches have provided an unprecedented opportunity to understand the evolution of insect P450s and their role in insect toxicology. These approaches have also been harnessed to provide new insights into the genomic alterations that lead to insecticide resistance, the mechanisms by which P450s are regulated, and the functional determinants of P450-mediated insecticide resistance. In parallel, an emerging body of work on the role of P450s in defining the sensitivity of beneficial insects to insecticides has been developed. The knowledge gained from these studies has applications for the management of P450-mediated resistance in insect pests and can be leveraged to safeguard the health of important beneficial insects.

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.

Assembling the P450 puzzle: on the sources of nonadditivity in drug metabolism

There is an increasing number of indications of an oversimplification in the premise that the cumulative properties of the human drug-metabolizing ensemble represent a simple aggregate of the properties of the constituting enzymes. Recent studies of the functional effects of hetero-association of multiple cytochrome P450 species and their interactions with metabolically related enzymes revealed a tight integration in the drug-metabolizing ensemble. In our opinion, the sources of interindividual variability in drug metabolism can be elucidated only when considering this ensemble as a multienzyme system, the functional parameters of which are determined by interactions between its constituents. In this article, we present a conceptual model providing a mechanistic explanation for the functional effects of the interactions between multiple P450 species and propose a clue to understanding the nonadditive behavior of the drug-metabolizing ensemble.

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.

Comparative analysis of the detoxification gene inventory of four major Spodoptera pest species in response to xenobiotics

The genus Spodoptera (Lepidoptera: Noctuidae) comprises some of the most polyphagous and destructive agricultural pests worldwide. The success of many species of this genus is due to their striking abilities to adapt to a broad range of host plants. Superfamilies of detoxification genes play a crucial role in the adaption to overcome plant defense mechanisms mediated by numerous secondary metabolites and toxins. Over the past decade, a substantial amount of expression data in Spodoptera larvae was produced for those genes in response to xenobiotics such as plant secondary metabolites, but also insecticide exposure. However, this information is scattered throughout the literature and in most cases does not allow to clearly identify candidate genes involved in host-plant adaptation and insecticide resistance. In the present review, we analyzed and compiled information on close to 600 pairs of inducers (xenobiotics) and induced genes from four main Spodoptera species: S. exigua, S. frugiperda, S. littoralis and S. litura. The cytochrome P450 monooxygenases (P450s; encoded by CYP genes) were the most upregulated detoxification genes across the literature for all four species. Most of the data was provided from studies on S. litura, followed by S. exigua, S. frugiperda and S. littoralis. We examined whether these detoxification genes were reported for larval survival under xenobiotic challenge in forward and reverse genetic studies. We further analyzed whether biochemical assays were carried out showing the ability of corresponding enzymes and transporters to breakdown and excrete xenobiotics, respectively. This revealed a clear disparity between species and the lack of genetic and biochemical information in S. frugiperda. Finally, we discussed the biological importance of detoxification genes for this genus and propose a workflow to study the involvement of these enzymes in an ecological and agricultural context.

Pyrethroid metabolism by eleven Helicoverpa armigera P450s from the CYP6B and CYP9A subfamilies

Lepidopteran P450s of the CYP6B and CYP9A subfamilies are thought to play important roles in host plant adaptation and insecticide resistance. An increasing number of paralogs and orthologs with high levels of sequence identity have been found in these subfamilies by mining recent genome projects. However, the biochemical function of most of them remains unknown. A better understanding of the evolution of P450 genes and of the catalytic competence of the enzymes they encode is needed to facilitate studies of host plant use and insecticide resistance. Here, we focused on the full complement of CYP6B (4 genes) and CYP9A (7 genes) in the generalist herbivore, Helicoverpa armigera. These P450s were heterologously expressed in Sf9 cells and compared functionally. In vitro assays showed that all CYP6B and CYP9A P450s can metabolize esfenvalerate efficiently, except for the evolutionarily divergent CYP6B43. A new 2'-hydroxy-metabolite of esfenvalerate was identified and found to be the main metabolite produced by CYP9A12. All tested P450s showed only low induction responses to esfenvalerate. To put these results from H. armigera P450s in perspective, 158 complete CYP6B and 100 complete CYP9A genes from 34 ditrysian species were manually curated. The CYP9A subfamily was more widespread than the CYP6B subfamily and the latter showed dramatic gains and losses, with ten species lacking CYP6B genes. Two adjacent CYP6B loci were found on chromosome 21, with different fates during the evolution of Lepidoptera. The diversity and functional redundancy of CYP6B and CYP9A genes challenge resistance management and pest control strategies as many P450s are available to insects to cope with chemical stresses they encounter.

Concerning P450 Evolution: Structural Analyses Support Bacterial Origin of Sterol 14α-Demethylases

Sterol biosynthesis, primarily associated with eukaryotic kingdoms of life, occurs as an abbreviated pathway in the bacterium Methylococcus capsulatus. Sterol 14α-demethylation is an essential step in this pathway and is catalyzed by cytochrome P450 51 (CYP51). In M. capsulatus, the enzyme consists of the P450 domain naturally fused to a ferredoxin domain at the C-terminus (CYP51fx). The structure of M. capsulatus CYP51fx was solved to 2.7 Å resolution and is the first structure of a bacterial sterol biosynthetic enzyme. The structure contained one P450 molecule per asymmetric unit with no electron density seen for ferredoxin. We connect this with the requirement of P450 substrate binding in order to activate productive ferredoxin binding. Further, the structure of the P450 domain with bound detergent (which replaced the substrate upon crystallization) was solved to 2.4 Å resolution. Comparison of these two structures to the CYP51s from human, fungi, and protozoa reveals strict conservation of the overall protein architecture. However, the structure of an "orphan" P450 from nonsterol-producing Mycobacterium tuberculosis that also has CYP51 activity reveals marked differences, suggesting that loss of function in vivo might have led to alterations in the structural constraints. Our results are consistent with the idea that eukaryotic and bacterial CYP51s evolved from a common cenancestor and that early eukaryotes may have recruited CYP51 from a bacterial source. The idea is supported by bioinformatic analysis, revealing the presence of CYP51 genes in >1,000 bacteria from nine different phyla, >50 of them being natural CYP51fx fusion proteins.

Function and Evolution of Nuclear Receptors in Environmental-Dependent Postembryonic Development

Nuclear receptors (NRs) fulfill key roles in the coordination of postembryonal developmental transitions in animal species. They control the metamorphosis and sexual maturation in virtually all animals and by that the two main environmental-dependent developmental decision points. Sexual maturation and metamorphosis are controlled by steroid receptors and thyroid receptors, respectively in vertebrates, while both processes are orchestrated by the ecdysone receptor (EcR) in insects. The regulation of these processes depends on environmental factors like nutrition, temperature, or photoperiods and by that NRs form evolutionary conserved mediators of phenotypic plasticity. While the mechanism of action for metamorphosis and sexual maturation are well studied in model organisms, the evolution of these systems is not entirely understood and requires further investigation. We here review the current knowledge of NR involvement in metamorphosis and sexual maturation across the animal tree of life with special attention to environmental integration and evolution of the signaling mechanism. Furthermore, we compare commonalities and differences of the different signaling systems. Finally, we identify key gaps in our knowledge of NR evolution, which, if sufficiently investigated, would lead to an importantly improved understanding of the evolution of complex signaling systems, the evolution of life history decision points, and, ultimately, speciation events in the metazoan kingdom.

Transcriptome analysis unravels RNAi pathways genes and putative expansion of CYP450 gene family in cotton leafhopper Amrasca biguttula (Ishida)

Cotton Leafhopper, Amrasca biguttula is an important pest of cotton and okra in the Indian subcontinent. Presently limited genomic/transcriptomic information is available for this insect in any of open source databases. The present study reports the first assembled and annotated de novo transcriptome of cotton leafhopper. Out of 75,551 transcripts, 39,613 CDS (Coding Sequence) were predicted with 35,282 showing positive blast hits with NCBI nr database. The Gene ontology (GO) analysis annotated 7431 CDS  with KEGG pathway categorizing these CDS into 22 different functional groups. The majority of CDS were annotated in signal transduction and transport catabolism pathways. The sequence data was screened for RNAi pathway genes and presence of 37 transcripts associated with this process confirmed the existence of robust RNAi machinery. The role of core RNAi machinery genes (Dicer-2, Ago-2, Piwi and Staufen) has been validated through dsRNA feeding studies. The data resource has also been used to identify potential RNAi targets and genes associated with insecticide detoxification specifically CYP 450 family. The current study provides a useful sequence resource which can be used to initiate molecular studies in this insect with emphasis on insecticide resistance, RNAi and functional genomics.

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.

Characterization of Shed genes encoding ecdysone 20-monooxygenase (CYP314A1) in the Y-organ of the blackback land crab, Gecarcinus lateralis

Molting in decapod crustaceans is controlled by ecdysteroid hormones synthesized and secreted by the molting gland, or Y-organ (YO). Halloween genes encode cytochrome P450 enzymes in the ecdysteroid synthetic pathway. The current paradigm is that YOs secrete an inactive precursor (e.g., ecdysone or E), which is hydroxylated at the #20 carbon to form an active hormone (20-hydroxyecdysone or 20E) by a mitochonrial 20-monooxygenase (CYP314A1) in peripheral tissues. 20-Monooxygenase is encoded by Shed in decapods and Shade in insects. We used eastern spiny lobster Shed sequences to extract six orthologs in the G. lateralis YO transcriptome. Phylogenetic analysis of the deduced amino acid sequences from six decapod species organized the Sheds into seven classes (Sheds 1-7), resulting in the assignment of the G. lateralis Sheds to Gl-Shed1, 2, 4A, 4B, 5A, and 5B. The mRNA levels of the six Gl-Sheds in the YO of intermolt animals were comparable to those in nine other tissues that included hepatopancreas and muscle. qPCR was used to compare the effects of molt induction by multiple leg autotomy (MLA) and eyestalk ablation (ESA) on Gl-Shed mRNA levels in the YO. Molt stage had little effect on Gl-Shed1 and Gl-Shed5B expression in the YO of MLA animals. Gl-Shed5A was expressed at the highest mRNA levels in the YO and was significantly increased during early and mid premolt stages. By contrast, ESA ± SB431542 had no effect on Gl-Shed expression at 1, 3, 5, and 7 days post-ESA. SB431542, which inhibits Transforming Growth Factor-β/activin signaling and blocks YO commitment, decreased Gl-Shed2 and Gl-Shed4A mRNA levels at 14 days post-ESA. A targeted metabolomic analysis showed that YOs cultured in vitro secreted E and 20E to the medium. These data suggest that the YO expresses 20-monooygenases that can convert E to 20E, which may contribute to the increase in active hormone in the hemolymph during premolt.

Cytochrome P450 expression and regulation in the brain

The regulation of brain cytochrome P450 enzymes (CYPs) is different compared with respective hepatic enzymes. This may result from anatomical bases and physiological functions of the two organs. The brain is composed of a variety of functional structures built of different interconnected cell types endowed with specific receptors that receive various neuronal signals from other brain regions. Those signals activate transcription factors or alter functioning of enzyme proteins. Moreover, the blood-brain barrier (BBB) does not allow free penetration of all substances from the periphery into the brain. Differences in neurotransmitter signaling, availability to endogenous and exogenous active substances, and levels of transcription factors between neuronal and hepatic cells lead to differentiated expression and susceptibility to the regulation of CYP genes in the brain and liver. Herein, we briefly describe the CYP enzymes of CYP1-3 families, their distribution in the brain, and discuss brain-specific regulation of CYP genes. In parallel, a comparison to liver CYP regulation is presented. CYP enzymes play an essential role in maintaining the levels of bioactive molecules within normal ranges. These enzymes modulate the metabolism of endogenous neurochemicals, such as neurosteroids, dopamine, serotonin, melatonin, anandamide, and exogenous substances, including psychotropics, drugs of abuse, neurotoxins, and carcinogens. The role of these enzymes is not restricted to xenobiotic-induced neurotoxicity, but they are also involved in brain physiology. Therefore, it is crucial to recognize the function and regulation of CYP enzymes in the brain to build a foundation for future medicine and neuroprotection and for personalized treatment of brain diseases.

How Can Drug Metabolism and Transporter Genetics Inform Psychotropic Prescribing?

Many genetic variants in drug metabolizing enzymes and transporters have been shown to be relevant for treating psychiatric disorders. Associations are strong enough to feature on drug labels and for prescribing guidelines based on such data. A range of commercial tests are available; however, there is variability in included genetic variants, methodology, and interpretation. We herein provide relevant background for understanding clinical associations with specific variants, other factors that are relevant to consider when interpreting such data (such as age, gender, drug-drug interactions), and summarize the data relevant to clinical utility of pharmacogenetic testing in psychiatry and the available prescribing guidelines. We also highlight areas for future research focus in this field.

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.

Research Progress of the Biosynthesis of Natural Bio-Antibacterial Agent Pulcherriminic Acid in Bacillus

Pulcherriminic acid is a cyclic dipeptide found mainly in Bacillus and yeast. Due to the ability of pulcherriminic acid to chelate Fe3+ to produce reddish brown pulcherrimin, microorganisms capable of synthesizing pulcherriminic acid compete with other microorganisms for environmental iron ions to achieve bacteriostatic effects. Therefore, studying the biosynthetic pathway and their enzymatic catalysis, gene regulation in the process of synthesis of pulcherriminic acid in Bacillus can facilitate the industrial production, and promote the wide application in food, agriculture and medicine industries. After initially discussing, this review summarizes current research on the synthesis of pulcherriminic acid by Bacillus, which includes the crystallization of key enzymes, molecular catalytic mechanisms, regulation of synthetic pathways, and methods to improve efficiency in synthesizing pulcherriminic acid and its precursors. Finally, possible applications of pulcherriminic acid in the fermented food, such as Chinese Baijiu, applying combinatorial biosynthesis will be summarized.

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.

Phylogeny, evolution, and potential ecological relationship of cytochrome CYP52 enzymes in Saccharomycetales yeasts

Cytochrome P450s from the CYP52 family participate in the assimilation of alkanes and fatty acids in fungi. In this work, the evolutionary history of a set of orthologous and paralogous CYP52 proteins from Saccharomycetales yeasts was inferred. Further, the phenotypic assimilation profiles were related with the distribution of cytochrome CYP52 members among species. The maximum likelihood phylogeny of CYP52 inferred proteins reveled a frequent ancient and modern duplication and loss events that generated orthologous and paralogous groups. Phylogeny and assimilation profiles of alkanes and fatty acids showed a family expansion in yeast isolated from hydrophobic-rich environments. Docking analysis of deduced ancient CYP52 proteins suggests that the most ancient function was the oxidation of C4-C11 alkanes, while the oxidation of >10 carbon alkanes and fatty acids is a derived character. The ancient CYP52 paralogs displayed partial specialization and promiscuous interaction with hydrophobic substrates. Additionally, functional optimization was not evident. Changes in the interaction of ancient CYP52 with different alkanes and fatty acids could be associated with modifications in spatial orientations of the amino acid residues that comprise the active site. The extended family of CYP52 proteins is likely evolving toward functional specialization, and certain redundancy for substrates is being maintained.

Comparative Sequence Analysis of TRI1 of Fusarium

Trichothecene mycotoxins are a class of secondary metabolites produced by multiple genera of fungi, including certain plant pathogenic Fusarium species. Functional variation in the TRI1 gene produces a novel Type A trichothecene called NX-2 in strains of F. graminearum. Using a bioinformatics approach, a systematic analysis of 52 translated TRI1 sequences of Fusarium species, including five F. graminearum NX-2 producers and four F. graminearum non-NX-2 producers, was conducted to explain the functional difference of TRI1p of FGNX-2. An assessment of several signature motifs of fungal P450s revealed amino acid substitutions in addition to the post-translational N-X-S/T sequons motif, which is indicative of N-linked glycosylation of this TRI1-encoded protein characteristic of NX-2 producers. There was evidence of selection bias, where TRI1 gene sequences were found to be under positive selection and, therefore, under functional constraints. The cumulative amino acid changes in the TRI1p sequences were reflected in the phylogenetic analyses which revealed species-specific clustering with a distinct separation of FGNX-2 from FG-non-NX-2 producers with high bootstrap support. Together, our findings provide insight into the amino acid sequence features responsible for the functional diversification of this TRI1p.

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.

A cytochrome P450 from the mustard leaf beetles hydroxylates geraniol, a key step in iridoid biosynthesis

Larvae of the leaf beetle Phaedon cochleariae synthesize the iridoid chysomelidial via the mevalonate pathway to repel predators. The normal terpenoid biosynthesis is integrated into the dedicated defensive pathway by the ω-hydroxylation of geraniol to (2E,6E)-2,6-dimethylocta-2,6-diene-1,8-diol (ω-OH-geraniol). Here we identify and characterize the P450 monooxygenase CYP6BH5 as the geraniol hydroxylase using integrated transcriptomics, proteomics and RNA interference (RNAi). In the fat body, 73 cytochrome P450s were identified, and CYP6BH5 was among those that were expressed specifically in fat body. Double stranded RNA mediated knockdown of CYP6BH5 led to a significant reduction of ω-hydroxygeraniol glucoside in the hemolymph and, later, of the chrysomelidial in the defensive secretion. Heterologously expressed CYP6BH5 converted geraniol to ω-OH-geraniol. In addition to geraniol, CYP6BH5 also catalyzes hydroxylation of other monoterpenols, such as nerol and citronellol to the corresponding α,ω-dihydroxy compounds.

Knockdown of cytochrome P450 CYP6 family genes increases susceptibility to carbamates and pyrethroids in the migratory locust, Locusta migratoria

Insect cytochrome P450 monooxygenase (CYP) plays a key role in the detoxification of insecticides. In this study, four cDNA sequences of CYP6 genes were identified and characterized. Transcription levels of LmCYP6HC1 and LmCYP6HCL1 were high in first- and fourth-instar nymph stages, respectively. LmCYP6HN1 was primarily expressed in the egg to third-instar nymph stages, while LmCYP6HQ1 was predominantly expressed in the stages from fourth-instar nymph to the adult. The four CYP6 genes were predominantly distributed in the antenna, brain, fat body, integument, and hemolymph. Piperonyl butoxide exposure inhibited total CYP activity and synergized the toxicity of carbamates and pyrethroids. Knockdown of LmCYP6HL1, LmCYP6HN1, and LmCYP6HQ1 increased nymph mortality following exposure to carbaryl, and silencing of LmCYP6HC1, LmCYP6HL1, LmCYP6HN1, and LmCYP6HQ1 comprehensively raised nymph mortality following exposure to fluvalinate. Knockdown of LmCYP6HL1 or LmCYP6HN1 significantly increased nymph mortality following exposure to cypermethrin or fenvalerate, respectively. These results suggest that the CYP6 family plays a key role in determining the susceptibility of Locusta migratoria to both carbamates and pyrethroids.

Mining insect genomes for functionally affiliated genes

Several hundred insect genome assemblies are already publicly available, and this total grows on a weekly basis. A major challenge now confronting insect science is how best to use genomic data to improve our understanding of insect biology. We consider a framework for genome analysis based on functional affiliation, that is, groups of genes involved in the same biological process or pathway, and explore how such an approach furthers our understanding of several aspects of insect phenotype. We anticipate that this approach will prove useful for future research across the breadth of insect studies, whatever organism or trait it involves.

Combined Toxicity of Insecticides and Fungicides Applied to California Almond Orchards to Honey Bee Larvae and Adults

Beekeepers providing pollination services for California almond orchards have reported observing dead or malformed brood during and immediately after almond bloom—effects that they attribute to pesticide exposure. The objective of this study was to test commonly used insecticides and fungicides during almond bloom on honey bee larval development in a laboratory bioassay. In vitro rearing of worker honey bee larvae was performed to test the effect of three insecticides (chlorantraniliprole, diflubenzuron, and methoxyfenozide) and three fungicides (propiconazole, iprodione, and a mixture of boscalid-pyraclostrobin), applied alone or in insecticide-fungicide combinations, on larval development. Young worker larvae were fed diets contaminated with active ingredients at concentration ratios simulating a tank-mix at the maximum label rate. Overall, larvae receiving insecticide and insecticide-fungicide combinations were less likely to survive to adulthood when compared to the control or fungicide-only treatments. The insecticide chlorantraniliprole increased larval mortality when combined with the fungicides propiconazole or iprodione, but not alone; the chlorantraniliprole-propiconazole combination was also found to be highly toxic to adult workers treated topically. Diflubenzuron generally increased larval mortality, but no synergistic effect was observed when combined with fungicides. Neither methoxyfenozide nor any methoxyfenozide-fungicide combination increased mortality. Exposure to insecticides applied during almond bloom has the potential to harm honey bees and this effect may, in certain instances, be more damaging when insecticides are applied in combination with fungicides.

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.

Firefly genomes illuminate parallel origins of bioluminescence in beetles

Fireflies and their luminous courtships have inspired centuries of scientific study. Today firefly luciferase is widely used in biotechnology, but the evolutionary origin of bioluminescence within beetles remains unclear. To shed light on this long-standing question, we sequenced the genomes of two firefly species that diverged over 100 million-years-ago: the North American Photinus pyralis and Japanese Aquatica lateralis. To compare bioluminescent origins, we also sequenced the genome of a related click beetle, the Caribbean Ignelater luminosus, with bioluminescent biochemistry near-identical to fireflies, but anatomically unique light organs, suggesting the intriguing hypothesis of parallel gains of bioluminescence. Our analyses support independent gains of bioluminescence in fireflies and click beetles, and provide new insights into the genes, chemical defenses, and symbionts that evolved alongside their luminous lifestyle.

Genome-wide analysis of cytochrome P450s of Trichoderma spp.: annotation and evolutionary relationships

Background

Cytochrome P450s form an important group of enzymes involved in xenobiotics degradation and metabolism, both primary and secondary. These enzymes are also useful in industry as biotechnological tools for bioconversion and a few are reported to be involved in pathogenicity. Trichoderma spp. are widely used in industry and agriculture and are known for their biosynthetic potential of a large number of secondary metabolites. For realising the full biosynthetic potential of an organism, it is important to do a genome-wide annotation and cataloguing of these enzymes.

Results

Here, we have studied the genomes of seven species (T. asperellum, T. atroviride, T. citrinoviride, T. longibrachiatum, T. reesei , T. harzianum and T. virens) and identified a total of 477 cytochrome P450s. We present here the classification, evolution and structure as well as predicted function of these proteins. This study would pave the way for functional characterization of these groups of enzymes and will also help in realization of their full economic potential.

Conclusion

Our CYPome annotation and evolutionary studies of the seven Trichoderma species now provides opportunities for exploration of research-driven strategies to select Trichoderma species for various applications especially in relation to secondary metabolism and degradation of environmental pollutants.

Retained duplications and deletions of CYP2C genes among primates

The human genome encodes about 60 functional enzymes of the cytochrome P450 superfamily, including four functional enzymes of the cytochrome P450 2C (CYP2C) subfamily. These enzymes have been shown to metabolize drugs and xenobiotic toxins, such as those in the diet, and are therefore of great importance for biomedical research and applications. While the pharmacology of P450 enzymes has been studied extensively, our understanding of molecular evolution of this gene family is incomplete, in part because a great variation exists in the number of CYP2C homologs across genomes. In humans, the enzymes encoded by these genes are responsible for the metabolism of more than 20% of clinical drugs, but this is not the naturalistic function of these enzymes, which is the metabolism of xenobiotics such as plant secondary metabolites. In this paper, we sought to correlate evolutionary relationships among primate CYP2C genes with known dietary profiles from these species, testing the hypothesis that these genes have evolved under the pressure of dietary toxins. Aside from a small number of deeply divergent genes, primate CYP2C paralogs form three separate clades: CYP2C18, CYP2C9/CYP2C19, and CYP2C8/CYP2C20. Our results showed that the CYP2C18 gene has been separately lost in Nomascus leucogenys and the Panini genomes, and there is no evidence that this gene has been under any positive selection among primates. While CYP2C20 has been retained in cercopithecoids, orthologous loci were separately lost in platyrrhines and hominoids. Notably, nine codons exhibited signature of positive selection. Finally, the CYP2C19 locus was duplicated in basal catarrhines, resulting in the birth of CYP2C9; but the ancestral locus was only retained in hominoid taxa. Overall, our findings support the hypothesis that primate CYP2C genes have evolved in response to selective pressures provided by dietary toxins, although not all gene clusters have evolved in the same manner. Our results may indicate an evolutionarily deep difference in ecology or physiology among higher-order primate taxa.

Intraspecific comparative genomics of isolates of the Norway spruce pathogen (Heterobasidion parviporum) and identification of its potential virulence factors

BACKGROUND

Heterobasidion parviporum is an economically most important fungal forest pathogen in northern Europe, causing root and butt rot disease of Norway spruce (Picea abies (L.) Karst.). The mechanisms underlying the pathogenesis and virulence of this species remain elusive. No reference genome to facilitate functional analysis is available for this species.

RESULTS

To better understand the virulence factor at both phenotypic and genomic level, we characterized 15 H. parviporum isolates originating from different locations across Finland for virulence, vegetative growth, sporulation and saprotrophic wood decay. Wood decay capability and latitude of fungal origins exerted interactive effects on their virulence and appeared important for H. parviporum virulence. We sequenced the most virulent isolate, the first full genome sequences of H. parviporum as a reference genome, and re-sequenced the remaining 14 H. parviporum isolates. Genome-wide alignments and intrinsic polymorphism analysis showed that these isolates exhibited overall high genomic similarity with an average of at least 96% nucleotide identity when compared to the reference, yet had remarkable intra-specific level of polymorphism with a bias for CpG to TpG mutations. Reads mapping coverage analysis enabled the classification of all predicted genes into five groups and uncovered two genomic regions exclusively present in the reference with putative contribution to its higher virulence. Genes enriched for copy number variations (deletions and duplications) and nucleotide polymorphism were involved in oxidation-reduction processes and encoding domains relevant to transcription factors. Some secreted protein coding genes based on the genome-wide selection pressure, or the presence of variants were proposed as potential virulence candidates.

CONCLUSION

Our study reported on the first reference genome sequence for this Norway spruce pathogen (H. parviporum). Comparative genomics analysis gave insight into the overall genomic variation among this fungal species and also facilitated the identification of several secreted protein coding genes as putative virulence factors for the further functional analysis. We also analyzed and identified phenotypic traits potentially linked to its virulence.