Study of the whole genome, methylome and transcriptome of Cordyceps militaris

Genomic and Transcriptome Analysis Reveals the Biosynthesis Network of Cordycepin in Cordyceps militaris

Cordycepin is the primary active compound of Cordyceps militaris. However, the definitive genetic mechanism governing cordycepin synthesis in fruiting body growth and development remains elusive, necessitating further investigation. This study consists of 64 C. militaris strains collected from northeast China. The high-yielding cordycepin strain CMS19 was selected for the analysis of cordycepin production and the genetic basis of cordycepin anabolism. First, the whole-genome sequencing of CMS19 yielded a final size of 30.96 Mb with 8 contigs and 9781 protein-coding genes. The genome component revealed the presence of four additional secondary metabolite gene clusters compared with other published genomes, suggesting the potential for the production of new natural products. The analyses of evolutionary and genetic differentiation revealed a close relationship between C. militaris and Beauveria bassiana. The population of strains distributed in northeast China exhibited the significant genetic variation. Finally, functional genes associated with cordycepin synthesis were identified using a combination of genomic and transcriptomic analyses. A large number of functional genes associated with energy and purine metabolism were significantly enriched, facilitating the reconstruction of a hypothetical cordycepin metabolic pathway. Therefore, our speculation of the cordycepin metabolism pathway involved 24 genes initiating from the glycolysis and pentose phosphate pathways, progressing through purine metabolism, and culminating in the core region of cordycepin synthesis. These findings could offer fundamental support for scientific utilizations of C. militaris germplasm resources and standardized cultivation for cordycepin production.

The Sordariomycetes: an expanding resource with Big Data for mining in evolutionary genomics and transcriptomics

Advances in genomics and transcriptomics accompanying the rapid accumulation of omics data have provided new tools that have transformed and expanded the traditional concepts of model fungi. Evolutionary genomics and transcriptomics have flourished with the use of classical and newer fungal models that facilitate the study of diverse topics encompassing fungal biology and development. Technological advances have also created the opportunity to obtain and mine large datasets. One such continuously growing dataset is that of the Sordariomycetes, which exhibit a richness of species, ecological diversity, economic importance, and a profound research history on amenable models. Currently, 3,574 species of this class have been sequenced, comprising nearly one-third of the available ascomycete genomes. Among these genomes, multiple representatives of the model genera Fusarium, Neurospora, and Trichoderma are present. In this review, we examine recently published studies and data on the Sordariomycetes that have contributed novel insights to the field of fungal evolution via integrative analyses of the genetic, pathogenic, and other biological characteristics of the fungi. Some of these studies applied ancestral state analysis of gene expression among divergent lineages to infer regulatory network models, identify key genetic elements in fungal sexual development, and investigate the regulation of conidial germination and secondary metabolism. Such multispecies investigations address challenges in the study of fungal evolutionary genomics derived from studies that are often based on limited model genomes and that primarily focus on the aspects of biology driven by knowledge drawn from a few model species. Rapidly accumulating information and expanding capabilities for systems biological analysis of Big Data are setting the stage for the expansion of the concept of model systems from unitary taxonomic species/genera to inclusive clusters of well-studied models that can facilitate both the in-depth study of specific lineages and also investigation of trait diversity across lineages. The Sordariomycetes class, in particular, offers abundant omics data and a large and active global research community. As such, the Sordariomycetes can form a core omics clade, providing a blueprint for the expansion of our knowledge of evolution at the genomic scale in the exciting era of Big Data and artificial intelligence, and serving as a reference for the future analysis of different taxonomic levels within the fungal kingdom.

Biological roles of adenine methylation in RNA

N6-Methyladenosine (m6A) is one of the most abundant modifications of the epitranscriptome and is found in cellular RNAs across all kingdoms of life. Advances in detection and mapping methods have improved our understanding of the effects of m6A on mRNA fate and ribosomal RNA function, and have uncovered novel functional roles in virtually every species of RNA. In this Review, we explore the latest studies revealing roles for m6A-modified RNAs in chromatin architecture, transcriptional regulation and genome stability. We also summarize m6A functions in biological processes such as stem-cell renewal and differentiation, brain function, immunity and cancer progression.

Construction of Cordycepin High-Production Strain and Optimization of Culture Conditions

This study aimed to increase cordycepin production by over-expressing bio-synthetic enzyme genes, including the adenylosuccinate synthase, adenylosuccinate lyase, and 5'-nucleotidase genes. Research data showed that the extracellular and intracellular cordycepin concent of 24 recombinant strains were higher than those of C. militaris WT, indicating that over-expression of key enzyme genes increased cordycepin production. Among them, the CM-adss-5 strain had highest cordycepin production, and the extracellular and intracellular cordycepin concent were 1119.75 ± 1.61 and 65.56 ± 0.97 mg/L, which were 1.26 and 2.61 times that of C. militaris WT. This study also optimized the culture conditions of CM-adss-5 strain through single factor experiments to obtain the best culture conditions. The best culture condition was 25 °C constant temperature, 180-rpm shaking culture, fermentation period 12 days, inoculate amount 5%, initial pH 6, seed age 108 h, and liquid volume 110/250 mL. Then, the extracellular and intracellular cordycepin content of CM-adss-5 strain reached 2581.96 ± 21.07 and 164.08 ± 1.44 mg/L, which were higher by 130.6% and 150.3%, respectively. Therefore, our research provides a way to efficiently produce cordycepin for the development of cordycepin and its downstream products.

iLoc-miRNA: extracellular/intracellular miRNA prediction using deep BiLSTM with attention mechanism

The location of microRNAs (miRNAs) in cells determines their function in regulation activity. Studies have shown that miRNAs are stable in the extracellular environment that mediates cell-to-cell communication and are located in the intracellular region that responds to cellular stress and environmental stimuli. Though in situ detection techniques of miRNAs have made great contributions to the study of the localization and distribution of miRNAs, miRNA subcellular localization and their role are still in progress. Recently, some machine learning-based algorithms have been designed for miRNA subcellular location prediction, but their performance is still far from satisfactory. Here, we present a new data partitioning strategy that categorizes functionally similar locations for the precise and instructive prediction of miRNA subcellular location in Homo sapiens. To characterize the localization signals, we adopted one-hot encoding with post padding to represent the whole miRNA sequences, and proposed a deep bidirectional long short-term memory with the multi-head self-attention algorithm to model. The algorithm showed high selectivity in distinguishing extracellular miRNAs from intracellular miRNAs. Moreover, a series of motif analyses were performed to explore the mechanism of miRNA subcellular localization. To improve the convenience of the model, a user-friendly web server named iLoc-miRNA was established (http://iLoc-miRNA.lin-group.cn/).

Enhancing Genome-Scale Model by Integrative Exometabolome and Transcriptome: Unveiling Carbon Assimilation towards Sphingolipid Biosynthetic Capability of Cordyceps militaris

Cordyceps militaris is an industrially important fungus, which is often used in Asia as traditional medicine. There has been a published genome-scale metabolic model (GSMM) of C. militaris useful for predicting its growth behaviors; however, lipid metabolism, which plays a vital role in cellular functions, remains incomplete in the GSMM of C. militaris. A comprehensive study on C. militaris was thus performed by enhancing GSMM through integrative analysis of metabolic footprint and transcriptome data. Through the enhanced GSMM of C. militaris (called iPC1469), it contained 1469 genes, 1904 metabolic reactions and 1229 metabolites. After model evaluation, in silico growth simulation results agreed well with the experimental data of the fungal growths on different carbon sources. Beyond the model-driven integrative data analysis, interestingly, we found key metabolic responses in alteration of lipid metabolism in C. militaris upon different carbon sources. The sphingoid bases (e.g., sphinganine, sphingosine, and phytosphingosine) and ceramide were statistically significant accumulated in the xylose culture when compared with other cultures; this study suggests that the sphingolipid biosynthetic capability in C. militaris was dependent on the carbon source assimilated for cell growth; this finding provides a comprehensive basis for the sphingolipid biosynthesis in C. militaris that can help to further redesign its metabolic control for medicinal and functional food applications.

Application of High-Throughput Sequencing on the Chinese Herbal Medicine for the Data-Mining of the Bioactive Compounds

The Chinese Herbal Medicine (CHM) has been used worldwide in clinic to treat the vast majority of human diseases, and the healing effect is remarkable. However, the functional components and the corresponding pharmacological mechanism of the herbs are unclear. As one of the main means, the high-throughput sequencing (HTS) technologies have been employed to discover and parse the active ingredients of CHM. Moreover, a tremendous amount of effort is made to uncover the pharmacodynamic genes associated with the synthesis of active substances. Here, based on the genome-assembly and the downstream bioinformatics analysis, we present a comprehensive summary of the application of HTS on CHM for the synthesis pathways of active ingredients from two aspects: active ingredient properties and disease classification, which are important for pharmacological, herb molecular breeding, and synthetic biology studies.

Epitranscriptomics of Ischemic Heart Disease-The IHD-EPITRAN Study Design and Objectives

Epitranscriptomic modifications in RNA can dramatically alter the way our genetic code is deciphered. Cells utilize these modifications not only to maintain physiological processes, but also to respond to extracellular cues and various stressors. Most often, adenosine residues in RNA are targeted, and result in modifications including methylation and deamination. Such modified residues as N-6-methyl-adenosine (m6A) and inosine, respectively, have been associated with cardiovascular diseases, and contribute to disease pathologies. The Ischemic Heart Disease Epitranscriptomics and Biomarkers (IHD-EPITRAN) study aims to provide a more comprehensive understanding to their nature and role in cardiovascular pathology. The study hypothesis is that pathological features of IHD are mirrored in the blood epitranscriptome. The IHD-EPITRAN study focuses on m6A and A-to-I modifications of RNA. Patients are recruited from four cohorts: (I) patients with IHD and myocardial infarction undergoing urgent revascularization; (II) patients with stable IHD undergoing coronary artery bypass grafting; (III) controls without coronary obstructions undergoing valve replacement due to aortic stenosis and (IV) controls with healthy coronaries verified by computed tomography. The abundance and distribution of m6A and A-to-I modifications in blood RNA are charted by quantitative and qualitative methods. Selected other modified nucleosides as well as IHD candidate protein and metabolic biomarkers are measured for reference. The results of the IHD-EPITRAN study can be expected to enable identification of epitranscriptomic IHD biomarker candidates and potential drug targets.

DNA Methylation Regulates Transcription Factor-Specific Neurodevelopmental but Not Sexually Dimorphic Gene Expression Dynamics in Zebra Finch Telencephalon

Song learning in zebra finches (Taeniopygia guttata) is a prototypical example of a complex learned behavior, yet knowledge of the underlying molecular processes is limited. Therefore, we characterized transcriptomic (RNA-sequencing) and epigenomic (RRBS, reduced representation bisulfite sequencing; immunofluorescence) dynamics in matched zebra finch telencephalon samples of both sexes from 1 day post hatching (1 dph) to adulthood, spanning the critical period for song learning (20 and 65 dph). We identified extensive transcriptional neurodevelopmental changes during postnatal telencephalon development. DNA methylation was very low, yet increased over time, particularly in song control nuclei. Only a small fraction of the massive differential expression in the developing zebra finch telencephalon could be explained by differential CpG and CpH DNA methylation. However, a strong association between DNA methylation and age-dependent gene expression was found for various transcription factors (i.e., OTX2, AR, and FOS) involved in neurodevelopment. Incomplete dosage compensation, independent of DNA methylation, was found to be largely responsible for sexually dimorphic gene expression, with dosage compensation increasing throughout life. In conclusion, our results indicate that DNA methylation regulates neurodevelopmental gene expression dynamics through steering transcription factor activity, but does not explain sexually dimorphic gene expression patterns in zebra finch telencephalon.

Experimental approaches to tracking mobile genetic elements in microbial communities

Horizontal gene transfer is an important mechanism of microbial evolution and is often driven by the movement of mobile genetic elements between cells. Due to the fact that microbes live within communities, various mechanisms of horizontal gene transfer and types of mobile elements can co-occur. However, the ways in which horizontal gene transfer impacts and is impacted by communities containing diverse mobile elements has been challenging to address. Thus, the field would benefit from incorporating community-level information and novel approaches alongside existing methods. Emerging technologies for tracking mobile elements and assigning them to host organisms provide promise for understanding the web of potential DNA transfers in diverse microbial communities more comprehensively. Compared to existing experimental approaches, chromosome conformation capture and methylome analyses have the potential to simultaneously study various types of mobile elements and their associated hosts. We also briefly discuss how fermented food microbiomes, given their experimental tractability and moderate species complexity, make ideal models to which to apply the techniques discussed herein and how they can be used to address outstanding questions in the field of horizontal gene transfer in microbial communities.

Optimizing cultivation of Cordyceps militaris for fast growth and cordycepin overproduction using rational design of synthetic media

Cordyceps militaris is an entomopathogenic fungus which is often used in Asia as a traditional medicine developed from age-old wisdom. Presently, cordycepin from C. militaris is a great interest in medicinal applications. However, cellular growth of C. militaris and the association with cordycepin production remain poorly understood. To explore the metabolism of C. militaris as potential cell factories in medical and biotechnology applications, this study developed a high-quality genome-scale metabolic model of C. militaris, iNR1329, based on its genomic content and physiological data. The model included a total of 1329 genes, 1821 biochemical reactions, and 1171 metabolites among 4 different cellular compartments. Its in silico growth simulation results agreed well with experimental data on different carbon sources. iNR1329 was further used for optimizing the growth and cordycepin overproduction using a novel approach, POPCORN, for rational design of synthetic media. In addition to the high-quality GEM iNR1329, the presented POPCORN approach was successfully used to rationally design an optimal synthetic medium with C:N ratio of 8:1 for enhancing 3.5-fold increase in cordycepin production. This study thus provides a novel insight into C. militaris physiology and highlights a potential GEM-driven method for synthetic media design and metabolic engineering application. The iNR1329 and the POPCORN approach are available at the GitHub repository: https://github.com/sysbiomics/Cordyceps_militaris-GEM.