The Cytochrome P450 Engineering Database: a navigation and prediction tool for the cytochrome P450 protein family

Genome analysis of the esca-associated Basidiomycetes Fomitiporia mediterranea, Fomitiporia polymorpha, Inonotus vitis, and Tropicoporus texanus reveals virulence factor repertoires characteristic of white-rot fungi

Some Basidiomycete fungi are important plant pathogens, and certain species have been associated with the grapevine trunk disease esca. We present the genomes of 4 species associated with esca: Fomitiporia mediterranea, Fomitiporia polymorpha, Tropicoporus texanus, and Inonotus vitis. We generated high-quality phased genome assemblies using long-read sequencing. The genomic and functional comparisons identified potential virulence factors, suggesting their roles in disease development. Similar to other white-rot fungi known for their ability to degrade lignocellulosic substrates, these 4 genomes encoded a variety of lignin peroxidases and carbohydrate-active enzymes (CAZymes) such as CBM1, AA9, and AA2. The analysis of gene family expansion and contraction revealed dynamic evolutionary patterns, particularly in genes related to secondary metabolite production, plant cell wall decomposition, and xenobiotic degradation. The availability of these genomes will serve as a reference for further studies of diversity and evolution of virulence factors and their roles in esca symptoms and host resistance.

P450Rdb: A manually curated database of reactions catalyzed by cytochrome P450 enzymes

INTRODUCTION

Cytochrome P450 enzymes (P450s) are recognized as the most versatile catalysts worldwide, playing vital roles in numerous biological metabolism and biosynthesis processes across all kingdoms of life. Despite the vast number of P450 genes available in databases (over 300,000), only a small fraction of them (less than 0.2 %) have undergone functional characterization.

OBJECTIVES

To provide a convenient platform with abundant information on P450s and their corresponding reactions, we introduce the P450Rdb database, a manually curated resource compiles literature-supported reactions catalyzed by P450s.

METHODS

All the P450s and Reactions were manually curated from the literature and known databases. Subsequently, the P450 reactions organized and categorized according to their chemical reaction type and site. The website was developed using HTML and PHP languages, with the MySQL server utilized for data storage.

RESULTS

The current version of P450Rdb catalogs over 1,600 reactions, involving more than 590 P450s across a diverse range of over 200 species. Additionally, it offers a user-friendly interface with comprehensive information, enabling easy querying, browsing, and analysis of P450s and their corresponding reactions. P450Rdb is free available at http://www.cellknowledge.com.cn/p450rdb/.

CONCLUSIONS

We believe that this database will significantly promote structural and functional research on P450s, thereby fostering advancements in the fields of natural product synthesis, pharmaceutical engineering, biotechnological applications, agricultural and crop improvement, and the chemical industry.

Genome sequencing of Elaeocarpus spp. stem blight pathogen Pseudocryphonectria elaeocarpicola reveals potential adaptations to colonize woody bark

BACKGROUND

Elaeocarpus spp. stem blight, caused by Pseudocryphonectria elaeocarpicola, is a destructive disease, which will significantly reduce the productivity and longevity of Elaeocarpus spp. plants, especially in the Guangdong Province of China. However, few information is available for P. elaeocarpicola. To unravel the potential adaptation mechanism of stem adaptation, the whole genome of P. elaeocarpicola was sequenced by using the DNBSEQ and PacBio platforms.

RESULTS

P. elaeocarpicola harbors 44.49 Mb genome with 10,894 predicted coding genes. Genome analysis revealed that the P. elaeocarpicola genome encodes a plethora of pathogenicity-related genes. Analysis of carbohydrate-active enzymes (CAZymes) revealed a rich variety of enzymes participated in plant cell wall degradation, which could effectively degrade cellulose, hemicellulose and xyloglucans in the plant cell wall and promote the invasion of the host plant. There are 213 CAZyme families found in P. elaeocarpicola, among which glycoside hydrolase (GH) family has the largest number, far exceeding other tested fungi by 53%. Besides, P. elaeocarpicola has twice as many genes encoding chitin and cellulose degradation as Cryphonectria parasitica, which belong to the same family. The predicted typical secreted proteins of P. elaeocarpicola are numerous and functional, including many known virulence effector factors, indicating that P. elaeocarpicola has great potential to secrete virulence effectors to promote pathogenicity on host plants. AntiSMASH revealed that the genome encoded 61 secondary metabolic gene clusters including 86 secondary metabolic core genes which was much higher than C. parasitica (49). Among them, two gene cluster of P. elaeocarpicola, cluster12 and cluster52 showed 100% similarity with the mycotoxins synthesis clusters from Aspergillus steynii and Alternaria alternata, respectively. In addition, we annotated cytochrome P450 related enzymes, transporters, and transcription factors in P. elaeocarpicola, which are important virulence determinants of pathogenic fungi.

CONCLUSIONS

Taken together, our study represents the first genome assembly for P. elaeocarpicola and reveals the key virulence factors in the pathogenic process of P. elaeocarpicola, which will promote our understanding of its pathogenic mechanism. The acquired knowledge lays a foundation for further exploration of molecular interactions with the host and provide target for management strategies in future research.

A fungal P450 enzyme from Fusarium equiseti HG18 with 7β-hydroxylase activity in biosynthesis of ursodeoxycholic acid

Cytochrome P450 enzyme with 7β-hydroxylation capacity has attracted widespread attentions due to the vital roles in the biosynthesis of ursodeoxycholic acid (UDCA), a naturally active molecule for the treatment of liver and gallbladder diseases. In this study, a novel P450 hydroxylase (P450FE) was screen out from Fusarium equiseti HG18 and identified by a combination of genome and transcriptome sequencing, as well as heterologous expression in Pichia pastoris. The biotransformation of lithocholic acid (LCA) by whole cells of recombinant Pichia pastoris further confirmed the C7β-hydroxylation with 5.2% UDCA yield. It was firstly identified a fungal P450 enzyme from Fusarium equiseti HG18 with the capacity to catalyze the LCA oxidation producing UDCA. The integration of homology modeling and molecular docking discovered the substrate binding to active pockets, and the key amino acids in active center were validated by site-directed mutagenesis, and revealed that Q112, V362 and L363 were the pivotal residues of P450FE in regulating the activity and selectivity of 7β-hydroxylation. Specifically, V362I mutation exhibited 2.6-fold higher levels of UDCA and higher stereospecificity than wild-type P450FE. This advance provided guidance for improving the catalytic efficiency and selectivity of P450FE in LCA hydroxylation, indicative of the great potential in green synthesis of UDCA from biologically toxic LCA.

A polylactic acid degrading lipase from Bacillus safensis: Characterization and structural analysis

A polylactic acid degrading triacylglycerol lipase (TGL) was identified from Bacillus safensis based on genome annotation and validated by real-time quantitative PCR. TGL displayed optimal activity at pH 9.0 and 55 °C. It maintained stability at pH 9.0 and temperatures 45 °C. The activity of TGL was found to benefit from the presence of potassium sodium ions, and low concentrations of Triton X-100. The TGL could erode the surface of polylactic acid films and increase its hydrophilicity. The hydrolysis products of polylactic acid by TGL were lactate monomer and dimer. TGL contains a classical catalytic triad structure of lipase (Ser77, Asp133, and His156) and an Ala-X-Ser-X-Gly sequence. Compared with some lipases produced by the same genus Bacillus, TGL is highly conserved in its amino acid sequence, mainly reflected in the amino acid residues that exercise the enzyme activity, including the catalytic activity center and the substrate binding sites.

Preliminary exploratory research on the application value of oral and intestinal meta-genomics in predicting subjects' occupations-A case study of the distinction between students and migrant workers

Background

In the field of forensic science, accurately determining occupation of an individual can greatly assist in resolving cases such as criminal investigations or disaster victim identifications. However, estimating occupation can be challenging due to the intricate relationship between occupation and various factors, including gender, age, living environment, health status, medication use, and lifestyle habits such as alcohol consumption and smoking. All of these factors can impact the composition of oral or gut microbial community of an individual.

Methods and results

In this study, we collected saliva and feces samples from individuals representing different occupational sectors, specifically students and manual laborers. We then performed metagenomic sequencing on the DNA extracted from these samples to obtain data that could be analyzed for taxonomic and functional annotations in five different databases. The correlation between occupation with microbial information was assisted from the perspective of α and β diversity, showing that individuals belonging to the two occupations hold significantly different oral and gut microbial communities, and that this correlation is basically not affected by gender, drinking, and smoking in our datasets. Finally, random forest (RF) models were built with recursive feature elimination (RFE) processes. Models with 100% accuracy in both training and testing sets were constructed based on three species in saliva samples or on a single pathway annotated by the KEGG database in fecal samples, namely, "ko04145" or Phagosome.

Conclusion

Although this study may have limited representativeness due to its small sample size, it provides preliminary evidence of the potential of using microbiome information for occupational inference.

Comparative Pangenomic Insights into the Distinct Evolution of Virulence Factors Among Grapevine Trunk Pathogens

The permanent organs of grapevines (Vitis vinifera L.), like those of other woody perennials, are colonized by various unrelated pathogenic ascomycete fungi secreting cell wall-degrading enzymes and phytotoxic secondary metabolites that contribute to host damage and disease symptoms. Trunk pathogens differ in the symptoms they induce and the extent and speed of damage. Isolates of the same species often display a wide virulence range, even within the same vineyard. This study focuses on Eutypa lata, Neofusicoccum parvum, and Phaeoacremonium minimum, causal agents of Eutypa dieback, Botryosphaeria dieback, and Esca, respectively. We sequenced 50 isolates from viticulture regions worldwide and built nucleotide-level, reference-free pangenomes for each species. Through examination of genomic diversity and pangenome structure, we analyzed intraspecific conservation and variability of putative virulence factors, focusing on functions under positive selection and recent gene family dynamics of contraction and expansion. Our findings reveal contrasting distributions of putative virulence factors in the core, dispensable, and private genomes of each pangenome. For example, carbohydrate active enzymes (CAZymes) were prevalent in the core genomes of each pangenome, whereas biosynthetic gene clusters were prevalent in the dispensable genomes of E. lata and P. minimum. The dispensable fractions were also enriched in Gypsy transposable elements and virulence factors under positive selection (polyketide synthase genes in E. lata and P. minimum, glycosyltransferases in N. parvum). Our findings underscore the complexity of the genomic architecture in each species and provide insights into their adaptive strategies, enhancing our understanding of the underlying mechanisms of virulence. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Integrated Genome Sequencing and Transcriptome Analysis Identifies Candidate Pathogenicity Genes from Ustilago crameri

Ustilago crameri is a pathogenic basidiomycete fungus that causes foxtail millet kernel smut (FMKS), a devastating grain disease in most foxtail-millet-growing regions of the world. Here, we report an assembled high-quality genome sequence of U. crameri strain SCZ-6 isolated from the diseased grains of foxtail millet in Changzhi, Shanxi Province, China. The genome size is 19.55 Mb, consisting of 73 contigs (N50 = 840,209 bp) with a G + C content of 54.09%, and encoding 6576 predicted genes and 6486 genes supported by RNA-seq. Evolutionarily, U. crameri lies close to the barley smut U. hordei, and an obvious co-linearity was observed between these two smut fungi. We annotated the genome of U. crameri strain SCZ-6 using databases, identifying 1827 pathogen-host interaction (PHI)-associated genes, 1324 genes encoding fungal virulence factors, 259 CAZy-related genes, 80 genes encoding transporters, and 206 putative cytochrome P450 genes; their expression profiles at different inoculation time points were also detected. Additionally, 70 candidate pathogen effectors were identified according to their expression patterns and predicted functions. In summary, our results provide important insights into the pathogenic mechanisms of the pathogenesis-related genes of U. crameri and a robust foundation for further investigation.

Rumen microbiome-driven insight into bile acid metabolism and host metabolic regulation

Gut microbes play a crucial role in transforming primary bile acids (BAs) into secondary forms, which influence systemic metabolic processes. The rumen, a distinctive and critical microbial habitat in ruminants, boasts a diverse array of microbial species with multifaceted metabolic capabilities. There remains a gap in our understanding of BA metabolism within this ecosystem. Herein, through the analysis of 9371 metagenome-assembled genomes and 329 cultured organisms from the rumen, we identified two enzymes integral to BA metabolism: 3-dehydro-bile acid delta4,6-reductase (baiN) and the bile acid:Na + symporter family (BASS). Both in vitro and in vivo experiments were employed by introducing exogenous BAs. We revealed a transformation of BAs in rumen and found an enzyme cluster, including L-ribulose-5-phosphate 3-epimerase and dihydroorotate dehydrogenase. This cluster, distinct from the previously known BA-inducible operon responsible for 7α-dehydroxylation, suggests a previously unrecognized pathway potentially converting primary BAs into secondary BAs. Moreover, our in vivo experiments indicated that microbial BA administration in the rumen can modulate amino acid and lipid metabolism, with systemic impacts underscored by core secondary BAs and their metabolites. Our study provides insights into the rumen microbiome's role in BA metabolism, revealing a complex microbial pathway for BA biotransformation and its subsequent effect on host metabolic pathways, including those for glucose, amino acids, and lipids. This research not only advances our understanding of microbial BA metabolism but also underscores its wider implications for metabolic regulation, offering opportunities for improving animal and potentially human health.

Chromosome-level assembly of Dictyophora rubrovolvata genome using third-generation DNA sequencing and Hi-C analysis

Dictyophora rubrovolvata, a rare edible mushroom with both nutritional and medicinal values, was regarded as the "queen of the mushroom" for its attractive appearance. Dictyophora rubrovolvata has been widely cultivated in China in recent years, and many researchers were focusing on its nutrition, culture condition, and artificial cultivation. Due to a lack of genomic information, research on bioactive substances, cross breeding, lignocellulose degradation, and molecular biology is limited. In this study, we report a chromosome-level reference genome of D. rubrovolvata using the PacBio single-molecule real-time-sequencing technique and high-throughput chromosome conformation capture (Hi-C) technologies. A total of 1.83 Gb circular consensus sequencing reads representing ∼983.34 coverage of the D. rubrovolvata genome were generated. The final genome was assembled into 136 contigs with a total length of 32.89 Mb. The scaffold and contig N50 length were 2.71 and 2.48 Mb, respectively. After chromosome-level scaffolding, 11 chromosomes with a total length of 28.24 Mb were constructed. Genome annotation further revealed that 9.86% of the genome was composed of repetitive sequences, and a total of 508 noncoding RNA (rRNA: 329, tRNA: 150, ncRNA: 29) were annotated. In addition, 9,725 protein-coding genes were predicted, among which 8,830 (90.79%) genes were predicted using homology or RNA-seq. Benchmarking Universal Single-Copy Orthologs results further revealed that there were 80.34% complete single-copy fungal orthologs. In this study, a total of 360 genes were annotated as belonging to the carbohydrate-active enzymes family. Further analysis also predicted 425 cytochromes P450 genes, which can be classified into 41 families. This highly accurate, chromosome-level reference genome of D. rubrovolvata will provide essential genomic information for understanding the molecular mechanism in its fruiting body formation during morphological development and facilitate the exploitation of medicinal compounds produced by this mushroom.

Genomic insights into Aspergillus sydowii 29R-4-F02: unraveling adaptive mechanisms in subseafloor coal-bearing sediment environments

Introduction

Aspergillussydowii is an important filamentous fungus that inhabits diverse environments. However, investigations on the biology and genetics of A. sydowii in subseafloor sediments remain limited.

Methods

Here, we performed de novo sequencing and assembly of the A. sydowii 29R-4-F02 genome, an isolate obtained from approximately 2.4 km deep, 20-million-year-old coal-bearing sediments beneath the seafloor by employing the Nanopore sequencing platform.

Results and Discussion

The generated genome was 37.19 Mb with GC content of 50.05%. The final assembly consisted of 11 contigs with N50 of 4.6 Mb, encoding 12,488 putative genes. Notably, the subseafloor strain 29R-4-F02 showed a higher number of carbohydrate-active enzymes (CAZymes) and distinct genes related to vesicular fusion and autophagy compared to the terrestrial strain CBS593.65. Furthermore, 257 positively selected genes, including those involved in DNA repair and CAZymes were identified in subseafloor strain 29R-4-F02. These findings suggest that A. sydowii possesses a unique genetic repertoire enabling its survival in the extreme subseafloor environments over tens of millions of years.

Differential RNA-Seq Analysis Predicts Genes Related to Terpene Tailoring in Caryopteris × clandonensis

Enzymatic terpene functionalization is an essential part of plant secondary metabolite diversity. Within this, multiple terpene-modifying enzymes are required to enable the chemical diversity of volatile compounds essential in plant communication and defense. This work sheds light on the differentially transcribed genes within Caryopteris × clandonensis that are capable of functionalizing cyclic terpene scaffolds, which are the product of terpene cyclase action. The available genomic reference was subjected to further improvements to provide a comprehensive basis, where the number of contigs was minimized. RNA-Seq data of six cultivars, Dark Knight, Grand Bleu, Good as Gold, Hint of Gold, Pink Perfection, and Sunny Blue, were mapped on the reference, and their distinct transcription profile investigated. Within this data resource, we detected interesting variations and additionally genes with high and low transcript abundancies in leaves of Caryopteris × clandonensis related to terpene functionalization. As previously described, different cultivars vary in their modification of monoterpenes, especially limonene, resulting in different limonene-derived molecules. This study focuses on predicting the cytochrome p450 enzymes underlying this varied transcription pattern between investigated samples. Thus, making them a reasonable explanation for terpenoid differences between these plants. Furthermore, these data provide the basis for functional assays and the verification of putative enzyme activities.

Cytochrome P450 superfamily in spotted sea bass: Genome-wide identification and expression profiles under trichlorfon and environmental stresses

Cytochrome P450s (CYPs), as one of the most diverse enzyme superfamilies in nature, play critical functions in antioxidant reactions against endogenous and exogenous compounds. In this study, we performed genome-wide characterization of CYP superfamily members and analyzed their expression patterns under several abiotic stresses in spotted sea bass, which is known as an economically important fish species in the Chinese aquaculture industry. A total of 55 CYP genes were identified and divided into 17 families within 10 clans. The analysis of phylogeny, gene structure, and syntenic relationships provided evidence for the evolution of CYP genes and confirmed their annotation and orthology. The expression of CYP genes was examined in the liver during trichlorfon stress using quantitative real-time PCR. The results showed that 20 tested CYP genes displayed significant mRNA expression changes, indicating that they may play crucial roles in the metabolism of trichlorfon and can be potential biomarkers for trichlorfon pollution. Moreover, by screening transcriptomic databases, 10, 3 and 19 CYP genes exhibited differential expression patterns in response to hypoxia, alkalinity and heat stress, respectively. Taken together, this study provided insights into the regulation of CYP genes by toxicological and environmental stresses, laid basis for extensive functional studies of the CYP superfamily in spotted sea bass and other teleost species.

Artificial Intelligence in Drug Metabolism and Excretion Prediction: Recent Advances, Challenges, and Future Perspectives

Drug metabolism and excretion play crucial roles in determining the efficacy and safety of drug candidates, and predicting these processes is an essential part of drug discovery and development. In recent years, artificial intelligence (AI) has emerged as a powerful tool for predicting drug metabolism and excretion, offering the potential to speed up drug development and improve clinical success rates. This review highlights recent advances in AI-based drug metabolism and excretion prediction, including deep learning and machine learning algorithms. We provide a list of public data sources and free prediction tools for the research community. We also discuss the challenges associated with the development of AI models for drug metabolism and excretion prediction and explore future perspectives in the field. We hope this will be a helpful resource for anyone who is researching in silico drug metabolism, excretion, and pharmacokinetic properties.

Intraspecific Comparative Analysis Reveals Genomic Variation of Didymella arachidicola and Pathogenicity Factors Potentially Related to Lesion Phenotype

Didymella arachidicola is one of the most important fungal pathogens, causing foliar disease and leading to severe yield losses of peanuts (Arachis hypogaea L.) in China. Two main lesion phenotypes of peanut web blotch have been identified as reticulation type (R type) and blotch type (B type). As no satisfactory reference genome is available, the genomic variations and pathogenicity factors of D. arachidicola remain to be revealed. In the present study, we collected 41 D. arachidicola isolates from 26 geographic locations across China (33 for R type and 8 for B type). The chromosome-scale genome of the most virulent isolate (YY187) was assembled as a reference using PacBio and Hi-C technologies. In addition, we re-sequenced 40 isolates from different sampling sites. Genome-wide alignments showed high similarity among the genomic sequences from the 40 isolates, with an average mapping rate of 97.38%. An average of 3242 SNPs and 315 InDels were identified in the genomic variation analysis, which revealed an intraspecific polymorphism in D. arachidicola. The comparative analysis of the most and least virulent isolates generated an integrated gene set containing 512 differential genes. Moreover, 225 genes individually or simultaneously harbored hits in CAZy-base, PHI-base, DFVF, etc. Compared with the R type reference, the differential gene sets from all B type isolates identified 13 shared genes potentially related to lesion phenotype. Our results reveal the intraspecific genomic variation of D. arachidicola isolates and pathogenicity factors potentially related to different lesion phenotypes. This work sets a genomic foundation for understanding the mechanisms behind genomic diversity driving different pathogenic phenotypes of D. arachidicola.

The genome and transcriptome of Sarocladium terricola provide insight into ergosterol biosynthesis

Sarocladium terricola is a species of ascomycete fungus that has been recognized as a biocontrol agent for managing animal and plant pathogens, and exhibits significant potential as a feed additive. In this study, we utilized a combination of short-read Illumina sequencing and long-read PacBio sequencing to sequence, assemble, and analyze the genome of S. terricola. The resulting genome consisted of 11 scaffolds encompassing 30.27 Mb, with a GC content of 54.07%, and 10,326 predicted protein coding gene models. We utilized 268 single-copy ortholog genes to reconstruct the phylogenomic relationships among 26 ascomycetes, and found that S. terricola was closely related to two Acremonium species. We also determined that the ergosterol content of S. terricola was synthesized to nearly double levels when cultured in potato dextrose media compared to bean media (4509 mg/kg vs. 2382 mg/kg). Furthermore, transcriptome analyses of differentially expressed genes suggested that the ergosterol synthesis genes ERG3, ERG5, and ERG25 were significantly up-regulated in potato dextrose media. These results will help us to recognize metabolic pathway of ergosterol biosynthesis of S. terricloa comprehensivelly.

Molecular Characterization of an Isoflavone 2'-Hydroxylase Gene Revealed Positive Insights into Flavonoid Accumulation and Abiotic Stress Tolerance in Safflower

Flavonoids with significant therapeutic properties play an essential role in plant growth, development, and adaptation to various environments. The biosynthetic pathway of flavonoids has long been studied in plants; however, its regulatory mechanism in safflower largely remains unclear. Here, we carried out comprehensive genome-wide identification and functional characterization of a putative cytochrome P45081E8 gene encoding an isoflavone 2'-hydroxylase from safflower. A total of 15 CtCYP81E genes were identified from the safflower genome. Phylogenetic classification and conserved topology of CtCYP81E gene structures, protein motifs, and cis-elements elucidated crucial insights into plant growth, development, and stress responses. The diverse expression pattern of CtCYP81E genes in four different flowering stages suggested important clues into the regulation of secondary metabolites. Similarly, the variable expression of CtCYP81E8 during multiple flowering stages further highlighted a strong relationship with metabolite accumulation. Furthermore, the orchestrated link between transcriptional regulation of CtCYP81E8 and flavonoid accumulation was further validated in the yellow- and red-type safflower. The spatiotemporal expression of CtCYP81E8 under methyl jasmonate, polyethylene glycol, light, and dark conditions further highlighted its likely significance in abiotic stress adaption. Moreover, the over-expressed transgenic Arabidopsis lines showed enhanced transcript abundance in OE-13 line with approximately eight-fold increased expression. The upregulation of AtCHS, AtF3'H, and AtDFR genes and the detection of several types of flavonoids in the OE-13 transgenic line also provides crucial insights into the potential role of CtCYP81E8 during flavonoid accumulation. Together, our findings shed light on the fundamental role of CtCYP81E8 encoding a putative isoflavone 2'-hydroxylase via constitutive expression during flavonoid biosynthesis.

Multi-omics analysis identifies a CYP9K1 haplotype conferring pyrethroid resistance in the malaria vector Anopheles funestus in East Africa

Metabolic resistance to pyrethroids is a menace to the continued effectiveness of malaria vector controls. Its molecular basis is complex and varies geographically across Africa. Here, we used a multi‐omics approach, followed‐up with functional validation to show that a directionally selected haplotype of a cytochrome P450, CYP9K1 is a major driver of resistance in Anopheles funestus. A PoolSeq GWAS using mosquitoes alive and dead after permethrin exposure, from Malawi and Cameroon, detected candidate genomic regions, but lacked consistency across replicates. Targeted sequencing of candidate resistance genes detected several SNPs associated with known pyrethroid resistance QTLs. The most significant SNPs were in the cytochrome P450 CYP304B1 (Cameroon), CYP315A1 (Uganda) and the ABC transporter gene ABCG4 (Malawi). However, when comparing field resistant mosquitoes to laboratory susceptible, the pyrethroid resistance locus rp1 and SNPs around the ABC transporter ABCG4 were consistently significant, except for Uganda where SNPs in the P450 CYP9K1 was markedly significant. In vitro heterologous metabolism assays with recombinant CYP9K1 revealed that it metabolises type II pyrethroid (deltamethrin; 64% depletion) but not type I (permethrin; 0%), while moderately metabolising DDT (17%). CYP9K1 exhibited reduced genetic diversity in Uganda underlying an extensive selective sweep. Furthermore, a glycine to alanine (G454A) amino acid change in CYP9K1 was fixed in Ugandan mosquitoes but not in other An. funestus populations. This study sheds further light on the evolution of metabolic resistance in a major malaria vector by implicating more genes and variants that can be used to design field‐applicable markers to better track resistance Africa‐wide.

Whole genome sequencing of an edible and medicinal mushroom, Russula griseocarnosa, and its association with mycorrhizal characteristics

Russula griseocarnosa is a well-known ectomycorrhizal mushroom, which is mainly distributed in the Southern China. Although several scholars have attempted to isolate and cultivate fungal strains, no accurate method for culture of artificial fruiting bodies has been presented owing to difficulties associated with mycelium growth on artificial media. Herein, we sequenced R. griseocarnosa genome using the second- and third-generation sequencing technologies, followed by de novo assembly of high-throughput sequencing reads, and GeneMark-ES, BLAST, CAZy, and other databases were utilized for functional gene annotation. We also constructed a phylogenetic tree using different species of fungi, and also conducted comparative genomics analysis of R. griseocarnosa against its four representative species. In addition, we evaluated the accuracy of one already sequenced genome of R. griseocarnosa based on the internal transcribed spacer (ITS) sequencing of that type of species. The assembly process resulted in identification of 230 scaffolds with a total genome size of 50.67 Mbp. The gene prediction showed that R. griseocarnosa genome included 14,229 coding sequences (CDs). In addition, 470 RNAs were predicted with 155 transfer RNAs (tRNAs), 49 ribosomal RNAs (rRNAs), 41 small noncoding RNAs (sRNAs), 42 small nuclear RNAs (snRNAs), and 183 microRNAs (miRNAs). The predicted protein sequences of R. griseocarnosa were analyzed to indicate the existence of carbohydrate-active enzymes (CAZymes), and the results revealed that 153 genes encoded CAZymes, which were distributed in 58 CAZyme families. These enzymes included 78 glycoside hydrolases (GHs), 34 glycosyl transferases (GTs), 30 auxiliary activities (AAs), 2 carbohydrate esterases (CEs), 8 carbohydrate-binding modules (CBMs), and only one polysaccharide lyase (PL). Compared with other fungi, R. griseocarnosa had fewer CAZymes, and the number and distribution of CAZymes were similar to other mycorrhizal fungi, such as Tricholoma matsutake and Suillus luteus. Well-defined effector proteins that were associated with mycorrhiza-induced small-secreted proteins (MiSSPs) were not found in R. griseocarnosa, which indicated that there may be some special effector proteins to interact with host plants in R. griseocarnosa. The genome of R. griseocarnosa may provide new insights into the energy metabolism of ectomycorrhizal (ECM) fungi, a reference to study ecosystem and evolutionary diversification of R. griseocarnosa, as well as promoting the study of artificial domestication.

Whole genome sequences of two Trichophyton indotineae clinical isolates from India emerging as threats during therapeutic treatment of dermatophytosis

In the current study, we report the genome sequence of two different clinical isolates from India, Trichophyton indotineae UCMS-IGIB-CI12 and Trichophyton indotineae UCMS-IGIB-CI14. The resulting genome assembly achieved a 143-fold coverage in 824 contigs for T. indotineae UCMS-IGIB-CI12 and a 136-fold coverage in 904 contigs for T. indotineae UCMS-IGIB-CI14. Both the clinical isolates contain a c.1342G>A mutation corresponding to Ala448Thr amino acid substitution in erg1 and exhibit an intermittent drug response to terbinafine. Comparative genomics analysis with available genomes of Trichophyton interdigitale/Trichophyton mentagrophytes species complex revealed a similar genome architecture and identified large number of genes associated with virulence and pathogenicity, namely, lipases, proteases, LysM domain-containing factors, carbon metabolism enzymes and cytochrome P450 enzymes, in all the genomes. An analysis of single amino acid polymorphisms (SAPs) in the protein sequences of subtilisin and lipase enzyme families identified a higher frequency of SAPs in functionally important proteins, Sub3 and Sub6 and their possible use in multilocus phylogenetic analysis of T. interdigitale/T. mentagrophytes species complex. The whole genome sequences of T. indotineae clinical isolates provided in this report will, hence, serve as a key reference point for investigation of clinical strains and emerging drug resistance among dermatophytes originating from different parts of the world.

Extended many-item similarity indices for sets of nucleotide and protein sequences

Quantification of similarities between protein sequences or DNA/RNA strands is a (sub-)task that is ubiquitously present in bioinformatics workflows, and is usually accomplished by pairwise comparisons of sequences, utilizing simple (e.g. percent identity) or more intricate concepts (e.g. substitution scoring matrices). Complex tasks (such as clustering) rely on a large number of pairwise comparisons under the hood, instead of a direct quantification of set similarities. Based on our recently introduced framework that enables multiple comparisons of binary molecular fingerprints (i.e., direct calculation of the similarity of fingerprint sets), here we introduce novel symmetric similarity indices for analogous calculations on sets of character sequences with more than two (t) possible items (e.g. DNA/RNA sequences with t = 4, or protein sequences with t = 20). The features of these new indices are studied in detail with analysis of variance (ANOVA), and demonstrated with three case studies of protein/DNA sequences with varying degrees of similarity (or evolutionary proximity). The Python code for the extended many-item similarity indices is publicly available at: https://github.com/ramirandaq/tn_Comparisons.

Assembly and annotation of whole-genome sequence of Fusarium equiseti

Fusarium equiseti is a plant pathogen with a wide range of hosts and diverse effects, including probiotic effects. However, the molecular mechanisms underlying these effects remain unclear, hindering its effective utilization. The final assembly included 16 scaffolds of contiguous sequence without gaps. The total sequence length was 40,776,005 bp, and the GC content of 48.01%. In total, we annotated the putative function of 13,134 genes, accounting for 94.97% of the candidate genes. We identified two and 23 candidate genes that are likely involved in the production of mycotoxins zearalenone and trichothecene, respectively. A comparative genomic analysis supported the high quality of the F. equiseti assembly. Our comprehensive analysis of whole-genome sequence will serve as a valuable resource for future studies of expression, regulation, function and evolution of the genes of F. equiseti as well as studies into disease prevention and control.

Overexpression of Efflux Pumps Mediate Pan Resistance of Klebsiella pneumoniae Sequence Type 11

A clinically isolated pan-resistant Klebsiella pneumoniae strain (ST11), KPN142 was subjected to whole-genome sequencing. Genomic sequence of KPN142 showed that limited antibiotic resistances (β-lactams [bla_shv-11], sulfonamides [sul1 and dfrA22], bacitracin [bacA], tetracycline [tet34], aminoglycosides [ksgA, kdpE, aph(3)Ia, aac(3)III, and ant(3)Ia], and chloramphenicol [catA1]) were mediated by enzymes, and efflux pumps contributed most to pan resistance. Five types of multidrug resistance efflux pump families were identified, including the resistance nodulation division superfamily (AcrAB-TolC, AcrD, MdtABC, and KexD), the ATP-binding cascade superfamily (MacAB), the small multidrug resistance family (KpnEF), the multidrug and toxic compound extrusion family (KdeA), and the major facilitator superfamily (EmrAB). There was an AcrAB-TolC efflux pump system, and inhibitory regulatory gene acrR and ramR of system carried deletion mutation, which lead to overexpression of AcrAB-TolC efflux pump, and in turn plays key role in the pan resistance of KPN142. Moreover, we did not find mgrb, a suppressor in the expression of phoPQ, overexpression of which may confer the resistance of KPN142 to colistin B. In addition, K. pneumoniae KPN142 carries IS1, IS3, and IntI1, which means that KPN142 is able to transfer drug-resistance genes. Of note, we detected the overexpression of acrB, ramA, phoP, and phoQ by real-time quantitative reverse transcription-polymerase chain reaction, and carbonyl cyanide chlorophenylhydrazone was able to reverse the resistance patterns of K. pneumoniae KPN142. In conclusion, we consider that the overexpression of AcrAB-TolC efflux pump mediates the resistance to most common clinical antimicrobial agents, and the overexpression of phoPQ mediates the resistance to colistin B in K. pneumoniae KPN142.

Phylogenomics of Plant-Associated Botryosphaeriaceae Species

The Botryosphaeriaceae is a fungal family that includes many destructive vascular pathogens of woody plants (e.g., Botryosphaeria dieback of grape, Panicle blight of pistachio). Species in the genera Botryosphaeria, Diplodia, Dothiorella, Lasiodiplodia, Neofusicoccum, and Neoscytalidium attack a range of horticultural crops, but they vary in virulence and their abilities to infect their hosts via different infection courts (flowers, green shoots, woody twigs). Isolates of seventeen species, originating from symptomatic apricot, grape, pistachio, and walnut were tested for pathogenicity on grapevine wood after 4 months of incubation in potted plants in the greenhouse. Results revealed significant variation in virulence in terms of the length of the internal wood lesions caused by these seventeen species. Phylogenomic comparisons of the seventeen species of wood-colonizing fungi revealed clade-specific expansion of gene families representing putative virulence factors involved in toxin production and mobilization, wood degradation, and nutrient uptake. Statistical analyses of the evolution of the size of gene families revealed expansions of secondary metabolism and transporter gene families in Lasiodiplodia and of secreted cell wall degrading enzymes (CAZymes) in Botryosphaeria and Neofusicoccum genomes. In contrast, Diplodia, Dothiorella, and Neoscytalidium generally showed a contraction in the number of members of these gene families. Overall, species with expansions of gene families, such as secreted CAZymes, secondary metabolism, and transporters, were the most virulent (i.e., were associated with the largest lesions), based on our pathogenicity tests and published reports. This study represents the first comparative phylogenomic investigation into the evolution of possible virulence factors from diverse, cosmopolitan members of the Botryosphaeriaceae.

Pathogenic Adaptations Revealed by Comparative Genome Analyses of Two Colletotrichum spp., the Causal Agent of Anthracnose in Rubber Tree

Colletotrichum siamense and Colletotrichum australisinense cause Colletotrichum leaf disease that differ in their symptoms in rubber tree (Hevea brasiliensis), and pathogenicity of these two fungal species is also not identical on different cultivars of rubber tree. This divergence is often attributed to pathogen virulence factors, namely carbohydrate-active enzymes (CAZymes), secondary metabolites (SM), and small-secreted protein (SSP) effectors. The draft genome assembly and functional annotation of potential pathogenicity genes of both species obtained here provide an important and timely genomic resource for better understanding the biology and lifestyle of Colletotrichum spp. This should pave the way for designing more efficient disease control strategies in plantations of rubber tree. In this study, the genes associated with these categories were manually annotated in the genomes of C. australisinense GX1655 and C. siamense HBCG01. Comparative genomic analyses were performed to address the evolutionary relationships among these gene families in the two species. First, the size of genome assembly, number of predicted genes, and some of the functional categories differed significantly between the two congeners. Second, from the comparative genomic analyses, we identified some specific genes, certain higher abundance of gene families associated with CAZymes, CYP450, and SM in the genome of C. siamense, and Nep1-like proteins (NLP) in the genome of C. australisinense.

Discovery and modification of cytochrome P450 for plant natural products biosynthesis

Cytochrome P450s are widespread in nature and play key roles in the diversification and functional modification of plant natural products. Over the last few years, there has been remarkable progress in plant P450s identification with the rapid development of sequencing technology, "omics" analysis and synthetic biology. However, challenges still persist in respect of crystal structure, heterologous expression and enzyme engineering. Here, we reviewed several research hotspots of P450 enzymes involved in the biosynthesis of plant natural products, including P450 databases, gene mining, heterologous expression and protein engineering.

Genome Sequencing and Analysis of the Hypocrellin-Producing Fungus Shiraia bambusicola S4201

Shiraia bambusicola has long been used as a traditional Chinese medicine and its major medicinal active metabolite is hypocrellin, which exhibits outstanding antiviral and antitumor properties. Here we report the 32 Mb draft genome sequence of S. bambusicola S4201, encoding 11,332 predicted genes. The genome of S. bambusicola is enriched in carbohydrate-active enzymes (CAZy) and pathogenesis-related genes. The phylogenetic tree of S. bambusicola S4201 and nine other sequenced species was constructed and its taxonomic status was supported (Pleosporales, Dothideomycetes). The genome contains a rich set of secondary metabolite biosynthetic gene clusters, suggesting that strain S4201 has a remarkable capacity to produce secondary metabolites. Overexpression of the zinc finger transcription factor zftf, which is involved in hypocrellin A (HA) biosynthesis, increases HA production when compared with wild type. In addition, a new putative HA biosynthetic pathway is proposed. These results provide a framework to study the mechanisms of infection in bamboo and to understand the phylogenetic relationships of S. bambusicola S4201. At the same time, knowledge of the genome sequence may potentially solve the puzzle of HA biosynthesis and lead to the discovery of novel genes and secondary metabolites of importance in medicine and agriculture.

Comprehensive analysis of full genome sequence and Bd-milRNA/target mRNAs to discover the mechanism of hypovirulence in Botryosphaeria dothidea strains on pear infection with BdCV1 and BdPV1

Pear ring rot disease, mainly caused by Botryosphaeria dothidea, is widespread in most pear and apple-growing regions. Mycoviruses are used for biocontrol, especially in fruit tree disease. BdCV1 (Botryosphaeria dothidea chrysovirus 1) and BdPV1 (Botryosphaeria dothidea partitivirus 1) influence the biological characteristics of B. dothidea strains. BdCV1 is a potential candidate for the control of fungal disease. Therefore, it is vital to explore interactions between B. dothidea and mycovirus to clarify the pathogenic mechanisms of B. dothidea and hypovirulence of B. dothidea in pear. A high-quality full-length genome sequence of the B. dothidea LW-Hubei isolate was obtained using Single Molecule Real-Time sequencing. It has high repeat sequence with 9.3% and DNA methylation existence in the genome. The 46.34 Mb genomes contained 14,091 predicted genes, which of 13,135 were annotated. B. dothidea was predicted to express 3833 secreted proteins. In bioinformatics analysis, 351 CAZy members, 552 transporters, 128 kinases, and 1096 proteins associated with plant-host interaction (PHI) were identified. RNA-silencing components including two endoribonuclease Dicer, four argonaute (Ago) and three RNA-dependent RNA polymerase (RdRp) molecules were identified and expressed in response to mycovirus infection. Horizontal transfer of the LW-C and LW-P strains indicated that BdCV1 induced host gene silencing in LW-C to suppress BdPV1 transmission. To investigate the role of RNA-silencing in B. dothidea defense, we constructed four small RNA libraries and sequenced B. dothidea micro-like RNAs (Bd-milRNAs) produced in response to BdCV1 and BdPV1 infection. Among these, 167 conserved and 68 candidate novel Bd-milRNAs were identified, of which 161 conserved and 20 novel Bd-milRNA were differentially expressed. WEGO analysis revealed involvement of the differentially expressed Bd-milRNA-targeted genes in metabolic process, catalytic activity, cell process and response to stress or stimulus. BdCV1 had a greater effect on the phenotype, virulence, conidiomata, vertical and horizontal transmission ability, and mycelia cellular structure biological characteristics of B. dothidea strains than BdPV1 and virus-free strains. The results obtained in this study indicate that mycovirus regulates biological processes in B. dothidea through the combined interaction of antiviral defense mediated by RNA-silencing and milRNA-mediated regulation of target gene mRNA expression.

Whole genome sequencing and genome annotation of the wild edible mushroom, Russula griseocarnosa

Russula griseocarnosa is a species of edible ectomycorrhizal fungi with medicinal properties that grows in southern China. Total DNA was isolated from a fresh fruiting body of R. griseocarnosa and subjected to sequencing using Illumina Hiseq with the PacBio RS sequencing platform. Here, we present the 64.81 Mb draft genome map of R. griseocarnosa based on 471 scaffolds and 16,128 coding protein genes. The gene annotation of protein coding genes was used to obtain corresponding annotations by blastp. Phylogenetic analysis revealed a close evolutionary relationship of R. griseocarnosa to Heterobasidion irregulare and Stereum hirsutum in the core Russulales clade. The R. griseocarnosa genome encodes a repertoire of enzymes engaged in carbohydrate and polysaccharide metabolism, along with cytochrome P450s and secondary metabolite biosynthesis. The genome content of R. griseocarnosa provides insights into the genetic basis of its reported medicinal properties and serves as a reference for comparative genomics of fungi.

Phylogenomic analysis of cytochrome P450 multigene family and their differential expression analysis in Solanum lycopersicum L. suggested tissue specific promoters

Background

Cytochrome P450 (P450) is a functionally diverse and multifamily class of enzymes which catalyses vast variety of biochemical reactions. P450 genes play regulatory role in growth, development and secondary metabolite biosynthesis. Solanum lycopersicum L. (Tomato) is an economically important crop plant and model system for various studies with massive genomic data. The comprehensive identification and characterization of P450 genes was lacking. Probing tomato genome for P450 identification would provide valuable information about the functions and evolution of the P450 gene family.

Results

In the present study, we have identified 233 P450 genes from tomato genome along with conserved motifs. Through the phylogenetic analysis of Solanum lycopersicum P450 (SlP450) protein sequences, they were classified into two major clades and nine clans further divided into 42 families. RT-qPCR analysis of selected six candidate genes were corroborated with digital expression profile. Out of 233 SlP450 genes, 73 showed expression evidence in 19 tissues of tomato. Out of 22 intron gain/loss positions, two positions were conserved in tomato P450 genes supporting intron late theory of intron evolution in SlP450 families. The comparison between tomato and other related plant P450s families showed that CYP728, CYP733, CYP80, CYP92, CYP736 and CYP749 families have been evolved in tomato and few higher plants whereas lost from Arabidopsis. The global promoter analysis of SlP450 against all the protein coding genes, coupled with expression data, revealed statistical overrepresentation of few promoter motifs in SlP450 genes which were highly expressed in specific tissue of tomato. Hence, these identified promoter motifs can be pursued further as tissue specific promoter that are driving expression of respective SlP450.

Conclusions

The phylogenetic analysis and expression profiles of tomato P450 gene family offers essential genomic resource for their functional characterization. This study allows comparison of SlP450 gene family with other Solanaceae members which are also economically important and attempt to classify functionally important SlP450 genes into groups and families. This report would enable researchers working on Tomato P450 to select appropriate candidate genes from huge repertoire of P450 genes depending on their phylogenetic class, tissue specific expression and promoter prevalence.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5483-x) contains supplementary material, which is available to authorized users.

Enzymatic Synthesis of ω-Hydroxydodecanoic Acid By Employing a Cytochrome P450 from Limnobacter sp. 105 MED

ω-Hydroxylated fatty acids are valuable and versatile building blocks for the production of various adhesives, lubricants, cosmetic intermediates, etc. The biosynthesis of ω-hydroxydodecanoic acid from vegetable oils is one of the important green pathways for their chemical-based synthesis. In the present study, the novel monooxygenase CYP153AL.m from Limnobacter sp. 105 MED was used for the whole-cell biotransformations. We constructed three-component system that was comprised of CYP153AL.m, putidaredoxin and putidaredoxin reductase from Pseudomonas putida. This in vivo study demonstrated that CYP153AL.m is a powerful catalyst for the biosynthesis of ω-hydroxydodecanoic acid. Under optimized conditions, the application of a solid-state powdered substrate rather than a substrate dissolved in DMSO significantly enhanced the overall reaction titer of the process. By employing this efficient system, 2 g/L of 12-hydroxydodecanoic acid (12-OHDDA) was produced from 4 g/L of its corresponding fatty acid, which was namely dodecanoic acid. Furthermore, the system was extended to produce 3.28 g/L of 12-OHDDA using 4 g/L of substrate by introducing native redox partners. These results demonstrate the utility of CYP153AL.m-catalyzed biotransformations in the industrial production of 12-OHDDA and other valuable building blocks.

Kauniolide synthase is a P450 with unusual hydroxylation and cyclization-elimination activity

Guaianolides are an important class of sesquiterpene lactones with unique biological and pharmaceutical properties. They have been postulated to be derived from germacranolides, but for years no progress has been made in the elucidation of their biosynthesis that requires an unknown cyclization mechanism. Here we demonstrate the isolation and characterization of a cytochrome P450 from feverfew (Tanacetum parthenium), kauniolide synthase. Kauniolide synthase catalyses the formation of the guaianolide kauniolide from the germacranolide substrate costunolide. Unlike most cytochrome P450s, kauniolide synthase combines stereoselective hydroxylation of costunolide at the C3 position, with water elimination, cyclization and regioselective deprotonation. This unique mechanism of action is supported by in silico modelling and docking experiments. The full kauniolide biosynthesis pathway is reconstructed in the heterologous hosts Nicotiana benthamiana and yeast, paving the way for biotechnological production of guaianolide-type sesquiterpene lactones.

The crystal structure of P450-TT heme-domain provides the first structural insights into the versatile class VII P450s

The first crystal structure of a class VII P450, CYP116B46 from Tepidiphilus thermophilus, has been solved at 1.9 Å resolution. The structure reveals overall conservation of the P450-fold and a water conduit around the I-helix. Active site residues have been identified and sequence comparisons have been made with other class VII enzymes. A structure similarity search demonstrated that the P450-TT structure is similar to enzymes capable of oxy-functionalization of fatty acids, terpenes, macrolides, steroids and statins. The insight gained from solving this structure will provide a guideline for future engineering and modelling studies on this catalytically promiscuous class of enzymes.

Extreme genome diversity in the hyper-prevalent parasitic eukaryote Blastocystis

Blastocystis is the most prevalent eukaryotic microbe colonizing the human gut, infecting approximately 1 billion individuals worldwide. Although Blastocystis has been linked to intestinal disorders, its pathogenicity remains controversial because most carriers are asymptomatic. Here, the genome sequence of Blastocystis subtype (ST) 1 is presented and compared to previously published sequences for ST4 and ST7. Despite a conserved core of genes, there is unexpected diversity between these STs in terms of their genome sizes, guanine-cytosine (GC) content, intron numbers, and gene content. ST1 has 6,544 protein-coding genes, which is several hundred more than reported for ST4 and ST7. The percentage of proteins unique to each ST ranges from 6.2% to 20.5%, greatly exceeding the differences observed within parasite genera. Orthologous proteins also display extreme divergence in amino acid sequence identity between STs (i.e., 59%-61% median identity), on par with observations of the most distantly related species pairs of parasite genera. The STs also display substantial variation in gene family distributions and sizes, especially for protein kinase and protease gene families, which could reflect differences in virulence. It remains to be seen to what extent these inter-ST differences persist at the intra-ST level. A full 26% of genes in ST1 have stop codons that are created on the mRNA level by a novel polyadenylation mechanism found only in Blastocystis. Reconstructions of pathways and organellar systems revealed that ST1 has a relatively complete membrane-trafficking system and a near-complete meiotic toolkit, possibly indicating a sexual cycle. Unlike some intestinal protistan parasites, Blastocystis ST1 has near-complete de novo pyrimidine, purine, and thiamine biosynthesis pathways and is unique amongst studied stramenopiles in being able to metabolize α-glucans rather than β-glucans. It lacks all genes encoding heme-containing cytochrome P450 proteins. Predictions of the mitochondrion-related organelle (MRO) proteome reveal an expanded repertoire of functions, including lipid, cofactor, and vitamin biosynthesis, as well as proteins that may be involved in regulating mitochondrial morphology and MRO/endoplasmic reticulum (ER) interactions. In sharp contrast, genes for peroxisome-associated functions are absent, suggesting Blastocystis STs lack this organelle. Overall, this study provides an important window into the biology of Blastocystis, showcasing significant differences between STs that can guide future experimental investigations into differences in their virulence and clarifying the roles of these organisms in gut health and disease.

Gene-enriched draft genome of the cattle tick Rhipicephalus microplus: assembly by the hybrid Pacific Biosciences/Illumina approach enabled analysis of the highly repetitive genome

The genome of the cattle tick Rhipicephalus microplus, an ectoparasite with global distribution, is estimated to be 7.1Gbp in length and consists of approximately 70% repetitive DNA. We report the draft assembly of a tick genome that utilized a hybrid sequencing and assembly approach to capture the repetitive fractions of the genome. Our hybrid approach produced an assembly consisting of 2.0Gbp represented in 195,170 scaffolds with a N50 of 60,284bp. The Rmi v2.0 assembly is 51.46% repetitive with a large fraction of unclassified repeats, short interspersed elements, long interspersed elements and long terminal repeats. We identified 38,827 putative R. microplus gene loci, of which 24,758 were protein coding genes (≥100 amino acids). OrthoMCL comparative analysis against 11 selected species including insects and vertebrates identified 10,835 and 3,423 protein coding gene loci that are unique to R. microplus or common to both R. microplus and Ixodes scapularis ticks, respectively. We identified 191 microRNA loci, of which 168 have similarity to known miRNAs and 23 represent novel miRNA families. We identified the genomic loci of several highly divergent R. microplus esterases with sequence similarity to acetylcholinesterase. Additionally we report the finding of a novel cytochrome P450 CYP41 homolog that shows similar protein folding structures to known CYP41 proteins known to be involved in acaricide resistance.

The Genome of Undifilum oxytropis Provides Insights into Swainsonine Biosynthesis and Locoism

Undifilum oxytropis is a fungal endophyte of locoweeds. It produces swainsonine, which is the principal toxic ingredient of locoweeds. However, the genes, pathways and mechanisms of swainsonine biosynthesis are not known. In this study, the genome of U. oxytropis was firstly sequenced and assembled into a 70.05 megabases (Mb) draft genome, which encoded 11,057 protein-coding genes, and 54% of them were similar to current publicly available sequences. U. oxytropis genes were annotated and 164 putative genes were annotated into enzymes, such as Saccharopine dehydrogenase, Saccharopine oxidase, and Pyrroline-5-carboxylate reductase, hypothesized to be involved in the biosynthesis pathway of swainsonine. The genome sequence and gene annotation of U. oxytropis will provide new insights into functional analyses. The characterization of genes in swainsonine biosynthesis will greatly facilitate locoweed poisoning research and help direct locoism management.

Mannitol Stress Directs Flavonoid Metabolism toward Synthesis of Flavones via Differential Regulation of Two Cytochrome P450 Monooxygenases in Coleus forskohlii

Cytochrome P450 monooxygenases (CYP450s) are known to play important roles in biosynthesis of all secondary metabolites, including flavonoids. Despite this, few CYP450s have been functionally characterized in model plants and roles of fewer CYP450s are known in non-model, medicinal, and aromatic plants. Our study in Coleus forskohlii indicates that flavone synthase (CYP93B) and flavonoid 3' monooxygenase (CYP706C) are key enzymes positioned at a metabolic junction, to execute the biosynthesis of different sub-classes of flavonoids (flavones, flavonol, anthocynanin, isoflavones etc.) from a common precursor. Such branch points are favored targets for artificially modulating the metabolic flux toward specific metabolites, through genetic manipulation or use of elicitors that differentially impact the expression of branch point genes. Genkwanin, the only flavone reported from C. forskohlii, is known to possess anti-inflammatory activity. It is biosynthesized from the general flavonoid precursor: naringenin. Two differentially expressed cytochrome P450 genes (CfCYP93B, CfCYP706C), exhibiting maximum expression in leaf tissues, were isolated from C. forskohlii. Mannitol treatment resulted in increased expression of CfCYP93B and decrease in expression of CfCYP706C. Metabolite quantification data showed that genkwanin content increased and anthocyanin levels decreased in response to mannitol treatment. Alignment, phylogenetic analysis, modeling, and molecular docking analysis of protein sequences suggested that CfCYP93B may be involved in conversion of naringenin to flavones (possibly genkwanin via apigenin), while CfCYP706C may act on common precursors of flavonoid metabolism and channel the substrate toward production of flavonols or anthocynanins. Decrease in expression of CfCYP706C and increase in accumulation of genkwanin suggested that mannitol treatment may possibly lead to accumulation of genkwanin via suppression of a competitive branch of flavonoids in C. forskohlii.

Identification of imine reductase-specific sequence motifs

Chiral amines are valuable building blocks for the production of a variety of pharmaceuticals, agrochemicals and other specialty chemicals. Only recently, imine reductases (IREDs) were discovered which catalyze the stereoselective reduction of imines to chiral amines. Although several IREDs were biochemically characterized in the last few years, knowledge of the reaction mechanism and the molecular basis of substrate specificity and stereoselectivity is limited. To gain further insights into the sequence-function relationships, the Imine Reductase Engineering Database (www.IRED.BioCatNet.de) was established and a systematic analysis of 530 putative IREDs was performed. A standard numbering scheme based on R-IRED-Sk was introduced to facilitate the identification and communication of structurally equivalent positions in different proteins. A conservation analysis revealed a highly conserved cofactor binding region and a predominantly hydrophobic substrate binding cleft. Two IRED-specific motifs were identified, the cofactor binding motif GLGxMGx(5 )[ATS]x(4) Gx(4) [VIL]WNR[TS]x(2) [KR] and the active site motif Gx[DE]x[GDA]x[APS]x(3){K}x[ASL]x[LMVIAG]. Our results indicate a preference toward NADPH for all IREDs and explain why, despite their sequence similarity to β-hydroxyacid dehydrogenases (β-HADs), no conversion of β-hydroxyacids has been observed. Superfamily-specific conservations were investigated to explore the molecular basis of their stereopreference. Based on our analysis and previous experimental results on IRED mutants, an exclusive role of standard position 187 for stereoselectivity is excluded. Alternatively, two standard positions 139 and 194 were identified which are superfamily-specifically conserved and differ in R- and S-selective enzymes.

Generate a bioactive natural product library by mining bacterial cytochrome P450 patterns

The increased number of annotated bacterial genomes provides a vast resource for genome mining. Several bacterial natural products with epoxide groups have been identified as pre-mRNA spliceosome inhibitors and antitumor compounds through genome mining. These epoxide-containing natural products feature a common biosynthetic characteristic that cytochrome P450s (CYPs) and its patterns such as epoxidases are employed in the tailoring reactions. The tailoring enzyme patterns are essential to both biological activities and structural diversity of natural products, and can be used for enzyme pattern-based genome mining. Recent development of direct cloning, heterologous expression, manipulation of the biosynthetic pathways and the CRISPR-CAS9 system have provided molecular biology tools to turn on or pull out nascent biosynthetic gene clusters to generate a microbial natural product library. This review focuses on a library of epoxide-containing natural products and their associated CYPs, with the intention to provide strategies on diversifying the structures of CYP-catalyzed bioactive natural products. It is conceivable that a library of diversified bioactive natural products will be created by pattern-based genome mining, direct cloning and heterologous expression as well as the genomic manipulation.

Genomes to natural products PRediction Informatics for Secondary Metabolomes (PRISM)

Microbial natural products are an invaluable source of evolved bioactive small molecules and pharmaceutical agents. Next-generation and metagenomic sequencing indicates untapped genomic potential, yet high rediscovery rates of known metabolites increasingly frustrate conventional natural product screening programs. New methods to connect biosynthetic gene clusters to novel chemical scaffolds are therefore critical to enable the targeted discovery of genetically encoded natural products. Here, we present PRISM, a computational resource for the identification of biosynthetic gene clusters, prediction of genetically encoded nonribosomal peptides and type I and II polyketides, and bio- and cheminformatic dereplication of known natural products. PRISM implements novel algorithms which render it uniquely capable of predicting type II polyketides, deoxygenated sugars, and starter units, making it a comprehensive genome-guided chemical structure prediction engine. A library of 57 tailoring reactions is leveraged for combinatorial scaffold library generation when multiple potential substrates are consistent with biosynthetic logic. We compare the accuracy of PRISM to existing genomic analysis platforms. PRISM is an open-source, user-friendly web application available at http://magarveylab.ca/prism/.

Identification and functional characterization of the CYP51 gene from the yeast Xanthophyllomyces dendrorhous that is involved in ergosterol biosynthesis

BACKGROUND

Xanthophyllomyces dendrorhous is a basidiomycetous yeast that synthesizes astaxanthin, a carotenoid with great biotechnological impact. The ergosterol and carotenoid synthetic pathways derive from the mevalonate pathway and involve cytochrome P450 enzymes. Among these enzymes, the CYP51 family, which is involved in ergosterol biosynthesis, is one of the most remarkable that has C14-demethylase activity.

RESULTS

In this study, the CYP51 gene from X. dendrorhous was isolated and its function was analyzed. The gene is composed of ten exons and encodes a predicted 550 amino acid polypeptide that exhibits conserved cytochrome P450 structural characteristics and shares significant identity with the sterol C14-demethylase from other fungi. The functionality of this gene was confirmed by heterologous complementation in S. cerevisiae. Furthermore, a CYP51 gene mutation in X. dendrorhous reduced sterol production by approximately 40% and enhanced total carotenoid production by approximately 90% compared to the wild-type strain after 48 and 120 h of culture, respectively. Additionally, the CYP51 gene mutation in X. dendrorhous increased HMGR (hydroxy-methylglutaryl-CoA reductase, involved in the mevalonate pathway) and crtR (cytochrome P450 reductase) transcript levels, which could be associated with reduced ergosterol production.

CONCLUSIONS

These results suggest that the CYP51 gene identified in X. dendrorhous encodes a functional sterol C14-demethylase that is involved in ergosterol biosynthesis.

The crystal structure of the versatile cytochrome P450 enzyme CYP109B1 from Bacillus subtilis

The crystal structure of the versatile CYP109B1 enzyme from Bacillus subtilis has been solved at 1.8 Å resolution. This is the first structure of an enzyme from this CYP family, whose members are prevalent across diverse species of bacteria. In the crystal structure the enzyme has an open conformation with an access channel leading from the heme to the surface. The substrate-free structure reveals the location of the key residues in the active site that are responsible for binding the substrate in the correct orientation for regioselective oxidation. Importantly, there are significant differences among these residues in members of the CYP109 and closely related CYP106 families and these likely account for the variations in substrate binding and oxidation profiles observed with these enzymes. A whole-cell oxidation biosystem was developed, which contains CYP109B1 and a phthalate family oxygenase reductase (PFOR), from Pseudomonas putida KT24440, as the electron transfer partner. This electron transfer system is able to support CYP109B1 activity resulting in the regioselective hydroxylation of both α- and β-ionone in vivo and in vitro. The PFOR is therefore a versatile electron transfer partner that is able to support the activity of CYP enzymes from other bacterium. The crystal structure of CYP109B1 has a positively charged proximal face and this explains why it can interact with PFOR and adrenodoxin which are predominantly negatively charged around their [2Fe-2S] clusters.

MutaCYP: Classification of missense mutations in human cytochromes P450

BACKGROUND

Cytochrome P450 monooxygenases (CYPs) represent a large and diverse family of enzymes involved in various biological processes in humans. Individual genome sequencing has revealed multiple mutations in human CYPs, and many missense mutations have been associated with variety of diseases. Since 3D structures are not resolved for most human CYPs, there is a need for a reliable sequence-based prediction that discriminates benign and disease causing mutations.

METHODS

A new prediction method (MutaCYP) has been developed for scoring de novo missense mutations to have a deleterious effect. The method utilizes only five features, all of which are sequence-based: predicted relative solvent accessibility (RSA), variance of predicted RSA among the residues in close sequence proximity, Z-score of Shannon entropy for a given position, difference in similarity scores and weighted difference in size between wild type and new amino acids. The method is based on a single neural network.

RESULTS

MutaCYP achieves MCC = 0.70, Q2 = 88.52%, Recall = 93.40% with Precision = 91.09%, and AUC = 0.909. Comparative evaluation with other existing methods indicates that MutaCYP outperforms SIFT and PolyPhen-2. Predictions by MutaCYP appear to be orthogonal to predictions by the evaluated methods. Potential issues on reliability of annotations of mutations in the existing databases are discussed.

CONCLUSIONS

A new accurate method, MutaCYP, for classification of missense mutations in human CYPs is presented. The prediction model consists of only five sequence-based features, including a real-valued predicted relative solvent accessibility. The method is publicly available at http://research.cchmc.org/MutaSense/.

Comparative analysis of P450 signature motifs EXXR and CXG in the large and diverse kingdom of fungi: identification of evolutionarily conserved amino acid patterns characteristic of P450 family

Cytochrome P450 monooxygenases (P450s) are heme-thiolate proteins distributed across the biological kingdoms. P450s are catalytically versatile and play key roles in organisms primary and secondary metabolism. Identification of P450s across the biological kingdoms depends largely on the identification of two P450 signature motifs, EXXR and CXG, in the protein sequence. Once a putative protein has been identified as P450, it will be assigned to a family and subfamily based on the criteria that P450s within a family share more than 40% homology and members of subfamilies share more than 55% homology. However, to date, no evidence has been presented that can distinguish members of a P450 family. Here, for the first time we report the identification of EXXR- and CXG-motifs-based amino acid patterns that are characteristic of the P450 family. Analysis of P450 signature motifs in the under-explored fungal P450s from four different phyla, ascomycota, basidiomycota, zygomycota and chytridiomycota, indicated that the EXXR motif is highly variable and the CXG motif is somewhat variable. The amino acids threonine and leucine are preferred as second and third amino acids in the EXXR motif and proline and glycine are preferred as second and third amino acids in the CXG motif in fungal P450s. Analysis of 67 P450 families from biological kingdoms such as plants, animals, bacteria and fungi showed conservation of a set of amino acid patterns characteristic of a particular P450 family in EXXR and CXG motifs. This suggests that during the divergence of P450 families from a common ancestor these amino acids patterns evolve and are retained in each P450 family as a signature of that family. The role of amino acid patterns characteristic of a P450 family in the structural and/or functional aspects of members of the P450 family is a topic for future research.