Origins of P450 diversity

Cytochrome P450s from the fungal symbiont of Sirex noctilio, Amylostereum areolatum: characterization, sequence analysis and their response to host terpenoids

Amylostereum areolatum is the main symbiotic fungus of woodwasp (Sirex noctilio) and is an active participant in the large-scale death of Pinus sylvestris. Woodwasps and associated fungi resort to active detoxification mechanisms to overcome the toxicity of host defence chemicals, which consist of a multitude of monoterpenes and diterpenes. Since cytochrome P450 (CYP) is considered to have extensive detoxification capabilities in fungi, 35 CYP genes were identified from A. areolatum by Iso-Seq. Moreover, we analysed the transcriptional levels of 11 chosen CYPs in the mycelia of fungi grown on different carbon sources or sprayed with different terpene mixtures or extracts to explore the relationship between CYPs and utilization of terpenoids or the detoxification capabilities for terpenoids. Molecular modeling and docking were also employed to predict the interaction between the P450 protein structure and substrate. The results showed that the 35 CYPs of A. areolatum belong to 22 families and 12 clans, and all had the typical P450 conserved domains. RT-qPCR revealed that most CYPs were down-regulated with monoterpenes as the sole carbon source and up-regulated with diterpenes. CYPs were induced after spraying with xylem extract, and venom can increase the expression levels of CYPs. Molecular docking predicted that three P450 proteins (CYP5037BM3, CYP5144KC1, and CYP5152T1) bind tightly to diterpenes through hydrogen bonds and π-alkyl interactions, they have a higher binding affinity for diterpenes than for monoterpenes, suggesting that they may preferentially metabolize diterpenes. This indicates that A. areolatum regulates the expression levels of key P450 enzymes through a special response pattern to reduce the toxicity of the chemical defence of the host on itself and its symbiotic insect, Sirex noctilio.

Genome-wide identification and expression analyses of CYP450 genes in Chrysanthemum indicum

BACKGROUND

The cytochrome P450 superfamily comprises a large group of enzymes crucial for the biosynthesis and metabolism of diverse endogenous and exogenous secondary metabolites in plants. Chrysanthemum, an ornamental genus with considerable medicinal value, is one of the most economically important floricultural crops in the world. The characteristics and functions of CYP450 genes in Chrysanthemum species, however, remain largely unknown.

RESULTS

In this study, we identified 371 CYP450 genes in the Chrysanthemum indicum genome, and categorized them into 8 clans and 44 families through phylogenetic analysis. Gene duplication analysis revealed 111 genes in 47 tandem duplicated clusters and 28 genes in 15 syntenic blocks, suggesting that extensive duplication events may account for the rapid expansion of CiCYP450 superfamily. Additionally, extensive variations in gene structure, motif composition, and cis-regulatory element likely enhance the functional diversity of CiCYP450 proteins. Volatile metabolomic analysis detected a total of 53 distinct volatile organic compounds across the leaves, stems, and roots of C. indicum, with 19 and 16 compounds being exclusive to leaves and stems, respectively. Transcriptomic analysis identified 248 expressed CiCYP450 genes, with 31, 40, and 88 specifically or preferentially expressed in leaves, stems, and roots, respectively. Further correlation analyses between gene expression levels and compound contents highlighted 36 candidate CiCYP450 genes potentially responsible for the biosynthesis of 47 volatile organic compounds.

CONCLUSIONS

The genome-wide analyses of cytochrome P450 superfamily offers essential genomic resources for functional studies of CiCYP450 genes, and is significant for the molecular breeding of Chrysanthemum.

Non-synonymous ERG11 mutations in M. restricta and M. arunalokei: impact on azole susceptibility

Malassezia are commensal lipid-dependent yeasts and opportunistic pathogens that cause superficial mycoses and systemic infection. Azole antifungals target cell wall ergosterol synthesis and are the first line of antifungal treatment. ERG11 gene mutations and overexpression are major mechanisms conferring azole resistance and resulting in antifungal therapy failure. Malassezia restricta is found ubiquitously on healthy and diseased skin, with azole-resistant isolates described. Malassezia arunalokei is a relatively new, closely related common skin species. Ketoconazole and itraconazole were the most effective at inhibiting both species. Isolates of M. restricta and M. arunalokei from healthy skin of Singapore subjects were cultured, evaluated, and generally susceptible to common over-the-counter azoles, including clotrimazole, except for select less-susceptible strains. Some less-susceptible strains have novel or reported non-synonymous mutations in the ERG11 gene, such as R88C. The QK178RQ ERG11 sequence variation was observed to be associated with differences in M. restricta and M. arunalokei as independent species. In the absence of identified ERG11 mutations, strains with elevated MICs were observed to have elevated ERG11 expression and drug efflux pump expression/activity. We conclude that antifungal susceptibility is determined by a combination of intrinsic (e.g., mutations, gene expression, efflux pump activity) and extrinsic (e.g., skin condition, prior antifungal exposure) factors and that the skin microbiome serves as a reference for the emergence of new mutations and strain phenotypes.

IMPORTANCE

Malassezia over colonization is associated with conditions such as dandruff and seborrheic dermatitis, which give rise to unpleasant itching and swelling on the skin. Azole antifungals such as ketoconazole, clotrimazole, and miconazole are the primary treatments of choice available as over-the-counter creams or shampoos. However, the emergence of antifungal resistance leads to a loss of treatment efficacy and persistent fungal infection. To understand the mechanisms underlying antifungal resistance, we profiled the susceptibility profiles of commensal Malassezia isolates from the skin and identified novel ERG11 mutations. Our results indicate that antifungal susceptibility is determined by a combination of factors (mutations, efflux pump activity, gene expression, copy number) and suggest that the healthy skin microbiome serves as a reference for the emergence of new mutations and strain phenotypes.

Properties of "Stable" Mosquito Cytochrome P450 Enzymes

The use of insecticides is widespread in the control of debilitating mosquito-borne diseases. P450 enzymes (CYPs) play essential roles in mosquito physiological function but also in the enzymatic detoxification of xenobiotics. Broadly speaking, CYPs can be classified as "stable", meaning those that have no or very few paralogs, and "labile", constituting gene families with many paralogous members. The evolutionary dichotomy between "stable" and "labile" P450 genes is fuzzy and there is not a clear phylogenetic demarcation between P450s involved in detoxification and P450s involved in essential metabolic processes. In this study, bioinformatic methods were used to explore differences in the sequences of "stable" and "labile" P450s that may facilitate their functional classification. Genomic and sequence data of Anopheles gambiae (Agam), Aedes aegypti (Aaeg), and Culex quinquefasciatus (Cqui) CYPs were obtained from public databases. The results of this study show that "stable" CYPs are encoded by longer genes, have longer introns and more exons, and contain a higher proportion of hydrophobic amino acids than "labile" CYPs. Compared to "labile" CYPs, a significantly higher proportion of "stable" CYPs are associated with biosynthetic and developmental processes.

High Antennal Expression of CYP6K1 and CYP4V2 Participate in the Recognition of Alarm Pheromones by Solenopsis invicta Buren

Insects have highly developed olfactory systems in which cytochrome P450s (CYPs) were involved as odor-degrading enzymes throughout the olfactory recognition of odor compounds by insects to avoid continuous stimulation of signaling molecules and thus damage to the olfactory nervous. To understand whether the highly expressed CYPs in the antennae play an olfactory function in Solenopsis invicta worker, in this study, we find six highly expressed antennal CYPs from the transcriptome of S. invicta. Multiple sequence alignment and phylogenetic analysis divided them into two families: the CYP3 family (SinvCYP6K1, SinvCYP6K1-1) and the CYP4 family (SinvCYP4C1, SinvCYP4C1-1, SinvCYP4C1-2, SinvCYP4V2). The expression patterns of these six CYPs were analyzed by RT-qPCR, which revealed that SinvCYP6K1 and SinvCYP4V2 were only highly expressed in the antennae of adult workers. The expression of SinvCYP6K1 and SinvCYP4V2 in workers was markedly diminished after feeding with dsRNA. The electroantennography (EAG) assay demonstrated that the silencing of either SinvCYP6K1 or SinvCYP4V2 resulted in a notable reduction in the EAG response of workers to 2-ethyl-3,6(5)-dimethylpyrazine (EDMP). Furthermore, the trajectory behavior assay showed that the worker's range and speed of movement in response to EDMP significant decreased after the silencing of SinvCYP6K1 and SinvCYP4V2. The findings indicated that both SinvCYP6K1 and SinvCYP4V2 were implicated in the recognition of EDMP by S. invicta.

Cytochrome P450 monooxygenase systems: Diversity and plasticity for adaptive stress response

Superfamily of cytochromes P450 (CYPs) is composed of heme-thiolate-containing monooxygenase enzymes, which play crucial roles in the biosynthesis, bioactivation, and detoxification of a variety of organic compounds, both endogenic and exogenic. Majority of CYP monooxygenase systems are multi-component and contain various redox partners, cofactors and auxiliary proteins, which contribute to their diversity in both prokaryotes and eukaryotes. Recent progress in bioinformatics and computational biology approaches make it possible to undertake whole-genome and phylogenetic analyses of CYPomes of a variety of organisms. Considerable variations in sequences within and between CYP families and high similarity in secondary and tertiary structures between all CYPs along with dramatic conformational changes in secondary structure elements of a substrate binding site during catalysis have been reported. This provides structural plasticity and substrate promiscuity, which underlie functional diversity of CYPs. Gene duplication and mutation events underlie CYP evolutionary diversity and emergence of novel selectable functions, which provide the involvement of CYPs in high adaptability to changing environmental conditions and dietary restrictions. In our review, we discuss the recent advancements and challenges in the elucidating the evolutionary origin and mechanisms underlying the CYP monooxygenase system diversity and plasticity. Our review is in the view of hypothesis that diversity of CYP monooxygenase systems is translated into the broad metabolic profiles, and this has been acquired during the long evolutionary time to provide structural plasticity leading to high adaptative capabilities to environmental stress conditions.

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.

Functional expression and regulation of eukaryotic cytochrome P450 enzymes in surrogate microbial cell factories

Cytochrome P450 (CYP) enzymes play crucial roles during the evolution and diversification of ancestral monocellular eukaryotes into multicellular eukaryotic organisms due to their essential functionalities including catalysis of housekeeping biochemical reactions, synthesis of diverse metabolites, detoxification of xenobiotics, and contribution to environmental adaptation. Eukaryotic CYPs with versatile functionalities are undeniably regarded as promising biocatalysts with great potential for biotechnological, pharmaceutical and chemical industry applications. Nevertheless, the modes of action and the challenges associated with these membrane-bound proteins have hampered the effective utilization of eukaryotic CYPs in a broader range. This review is focused on comprehensive and consolidated approaches to address the core challenges in heterologous expression of membrane-bound eukaryotic CYPs in different surrogate microbial cell factories, aiming to provide key insights for better studies and applications of diverse eukaryotic CYPs in the future. We also highlight the functional significance of the previously underrated cytochrome P450 reductases (CPRs) and provide a rational justification on the progression of CPR from auxiliary redox partner to function modulator in CYP catalysis.

The P450 genes of the cat flea, Ctenocephalides felis: a CYPome in flux

The genome of the cat flea, an ectoparasite of major veterinary importance and the first representative of the Siphonaptera, is highly unusual among arthropod genomes in showing a variable size and a very large number of gene duplications (Driscoll et al., 2020). The cat flea is the target of several classes of insecticides, justifying the description of its CYPome, the complement of P450s that are an important family of detoxification enzymes. 103 P450 genes were annotated on the nine chromosomes, with an additional 12 genes on small, extrachromosomal scaffolds. Only 34 genes were found as single sequences, with 47 duplicated two to four-fold. This included duplication of genes that are mostly single copy P450 genes in other arthropods. Large clusters of mitochondrial clan P450s were observed, resulting in a CYP12 bloom within this clan to 34 genes, a number of mitochondrial P450s not seen in other animals so far. The variable geometry of the cat flea CYPome poses a challenge to the study of P450 function in this species, and raises the question of the underlying causes of single copy control versus multicopy licence of P450 genes.

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.