Plant genomes remain highly fragmented and are often characterized by hundreds to thousands of assembly gaps. Here, we report chromosome-level reference and phased genome assembly ofOphiorrhiza pumila, a camptothecin-producing medicinal plant, through an ordered multi-scaffolding and experimental validation approach. With 21 assembly gaps and a contig N50 of 18.49 Mb,Ophiorrhizagenome is one of the most complete plant genomes assembled to date. We also report 273 nitrogen-containing metabolites, including diverse monoterpene indole alkaloids (MIAs). A comparative genomics approach identifies strictosidine biogenesis as the origin of MIA evolution. The emergence of strictosidine biosynthesis-catalyzing enzymes precede downstream enzymes’ evolution post γ whole-genome triplication, which occurred approximately 110 Mya inO. pumila, and before the whole-genome duplication inCamptotheca acuminataidentified here. Combining comparative genome analysis, multi-omics analysis, and metabolic gene-cluster analysis, we propose a working model for MIA evolution, and a pangenome for MIA biosynthesis, which will help in establishing a sustainable supply of camptothecin.
Background: Transcription elongation is frequently interrupted by pausing signals in DNA, with downstream effects on gene expression. Transcription errors also induce prolonged pausing, which can lead to a destabilized genome by interfering with DNA replication. Mechanisms of pausing associated with translocation blocks and misincorporation have been characterized in vitro, but not in vivo. Results: We investigate the pausing pattern of RNA polymerase (RNAP) in Escherichia coli by a novel approach, combining native elongating transcript sequencing (NET-seq) with RNase footprinting of the transcripts (RNET-seq). We reveal that the G-dC base pair at the 5′ end of the RNA-DNA hybrid interferes with RNAP translocation. The distance between the 5′ G-dC base pair and the 3′ end of RNA fluctuates over a three-nucleotide width. Thus, the G-dC base pair can induce pausing in post-translocated, pre-translocated, and backtracked states of RNAP. Additionally, a CpG sequence of the template DNA strand spanning the active site of RNAP inhibits elongation and induces G-to-A errors, which leads to backtracking of RNAP. Gre factors efficiently proofread the errors and rescue the backtracked complexes. We also find that pausing events are enriched in the 5′ untranslated region and antisense transcription of mRNA genes and are reduced in rRNA genes. Conclusions: In E. coli, robust transcriptional pausing involves RNAP interaction with G-dC at the upstream end of the RNA-DNA hybrid, which interferes with translocation. CpG DNA sequences induce transcriptional pausing and G-to-A errors.
Aconitum carmichaelii is an important medicinal herb used widely in China, Japan, India, Korea, and other Asian countries. While extensive research on the characterization of metabolic extracts of A. carmichaelii has shown accumulation of numerous bioactive metabolites including aconitine and aconitine-type diterpene alkaloids, its biosynthetic pathway remains largely unknown. Biosynthesis of these secondary metabolites is tightly controlled and mostly occurs in a tissue-specific manner; therefore, transcriptome analysis across multiple tissues is an attractive method to identify the molecular components involved for further functional characterization. In order to understand the biosynthesis of secondary metabolites, Illumina-based deep transcriptome profiling and analysis was performed for four tissues (flower, bud, leaf, and root) of A. carmichaelii, resulting in 5.5 Gbps clean RNA-seq reads assembled into 128,183 unigenes. Unigenes annotated as possible rate-determining steps of an aconitine-type biosynthetic pathway were highly expressed in the root, in accordance with previous reports describing the root as the accumulation site for these metabolites. We also identified 21 unigenes annotated as cytochrome P450s and highly expressed in roots, which represent candidate unigenes involved in the diversification of secondary metabolites. Comparative transcriptome analysis of A. carmichaelii with A. heterophyllum identified 20,232 orthogroups, representing 30,633 unigenes of A. carmichaelii, gene ontology enrichment analysis of which revealed essential biological process together with a secondary metabolic process to be highly enriched. Unigenes identified in this study are strong candidates for aconitine-type diterpene alkaloid biosynthesis, and will serve as useful resources for further validation studies.
Aspergillus fumigatus is a human fungal pathogen that causes aspergillosis in immunocompromised hosts. A. fumigatus is believed to be exposed to diverse environmental stresses in the host cells. The adaptation mechanisms are critical for infections in human bodies. Transcriptional networks in response to diverse environmental challenges remain to be elucidated. To gain insights into the adaptation to environmental stresses in A. fumigatus mycelia, we conducted time series transcriptome analyses. With the aid of RNA-seq, we explored the global gene expression profiles of mycelia in A. fumigatus upon exposure to diverse environmental changes, including heat, superoxide, and osmotic stresses. From the perspective of global transcriptomes, transient responses to superoxide and osmotic stresses were observed while responses to heat stresses were gradual. We identified the stress-responsive genes for particular stresses, and the 266 genes whose expression levels drastically fluctuated upon exposure to all tested stresses. Among these, the 77 environmental stress response genes are conserved in S. cerevisiae, suggesting that these genes might be more general prerequisites for adaptation to environmental stresses. Finally, we revealed the strong correlations among expression profiles of genes related to ‘rRNA processing’. The time series transcriptome analysis revealed the stress-responsive genes underlying the adaptation mechanisms in A. fumigatus mycelia. These results will shed light on the regulatory networks underpinning the adaptation of the filamentous fungi.
Aspergillus fumigatus is the Aspergillus species most commonly associated with aspergillosis. Of the various presentations of aspergillosis, one of the most frequently observed in cases involving A. fumigatus pulmonary infections is aspergilloma (PA). In such infections one finds a fungus ball composed of fungal hyphae, inflammatory cells, fibrin, mucus, and tissue debris. Chronic necrotizing pulmonary aspergillosis (CNPA), also known as semi-invasive or invasive aspergillosis, is locally invasive and predominantly seen in patients with mild immunodeficiency or with a chronic lung disease. In the present study, with the aid of a next-generation sequencer, we conducted whole genome sequence (WGS) analyses of 17 strains isolated from patients in Japan with PA and CNPA. A total of 99,088 SNPs were identified by mapping the reads to A. fumigatus genome reference strain Af293, and according to genome-wide phylogenetic analysis, there were no correlations between the whole genome sequence typing results and pathologic conditions of patients. Here, we conducted the first multi-genome WGS study to focus on the A. fumigatus strains isolated from patients with PA and CNPA, and comprehensively characterized genetic variations of strains. WGS approach will help in better understanding of molecular mechanisms of aspergillosis cases caused by A. fumigatus.