Haimei Chen

Genome sequence of the model medicinal mushroom Ganoderma lucidum

Ganoderma lucidum is a widely used medicinal macrofungus in traditional Chinese medicine that creates a diverse set of bioactive compounds. Here we report its 43.3-Mb genome, encoding 16,113 predicted genes, obtained using next-generation sequencing and optical mapping approaches. The sequence analysis reveals an impressive array of genes encoding cytochrome P450s (CYPs), transporters and regulatory proteins that cooperate in secondary metabolism. The genome also encodes one of the richest sets of wood degradation enzymes among all of the sequenced basidiomycetes. In all, 24 physical CYP gene clusters are identified. Moreover, 78 CYP genes are coexpressed with lanosterol synthase, and 16 of these show high similarity to fungal CYPs that specifically hydroxylate testosterone, suggesting their possible roles in triterpenoid biosynthesis. The elucidation of the G. lucidum genome makes this organism a potential model system for the study of secondary metabolic pathways and their regulation in medicinal fungi.

Complete Mitochondrial Genome of the Medicinal Mushroom Ganoderma lucidum

Ganoderma lucidum is one of the well-known medicinal basidiomycetes worldwide. The mitochondrion, referred to as the second genome, is an organelle found in most eukaryotic cells and participates in critical cellular functions. Elucidating the structure and function of this genome is important to understand completely the genetic contents of G. lucidum. In this study, we assembled the mitochondrial genome of G. lucidum and analyzed the differential expressions of its encoded genes across three developmental stages. The mitochondrial genome is a typical circular DNA molecule of 60,630 bp with a GC content of 26.67%. Genome annotation identified genes that encode 15 conserved proteins, 27 tRNAs, small and large rRNAs, four homing endonucleases, and two hypothetical proteins. Except for genes encoding trnW and two hypothetical proteins, all genes were located on the positive strand. For the repeat structure analysis, eight forward, two inverted, and three tandem repeats were detected. A pair of fragments with a total length around 5.5 kb was found in both the nuclear and mitochondrial genomes, which suggests the possible transfer of DNA sequences between two genomes. RNA-Seq data for samples derived from three stages, namely, mycelia, primordia, and fruiting bodies, were mapped to the mitochondrial genome and qualified. The protein-coding genes were expressed higher in mycelia or primordial stages compared with those in the fruiting bodies. The rRNA abundances were significantly higher in all three stages. Two regions were transcribed but did not contain any identified protein or tRNA genes. Furthermore, three RNA-editing sites were detected. Genome synteny analysis showed that significant genome rearrangements occurred in the mitochondrial genomes. This study provides valuable information on the gene contents of the mitochondrial genome and their differential expressions at various developmental stages of G. lucidum. The results contribute to the understanding of the functions and evolution of fungal mitochondrial DNA.

Computational Identification and Systematic Classification of Novel Cytochrome P450 Genes in Salvia miltiorrhiza

Salvia miltiorrhiza is one of the most economically important medicinal plants. Cytochrome P450 (CYP450) genes have been implicated in the biosynthesis of its active components. However, only a dozen full-length CYP450 genes have been described, and there is no systematic classification of CYP450 genes in S. miltiorrhiza. We obtained 77,549 unigenes from three tissue types of S. miltiorrhiza using RNA-Seq technology. Combining our data with previously identified CYP450 sequences and scanning with the CYP450 model from Pfam resulted in the identification of 116 full-length and 135 partial-length CYP450 genes. The 116 genes were classified into 9 clans and 38 families using standard criteria. The RNA-Seq results showed that 35 CYP450 genes were co-expressed with CYP76AH1, a marker gene for tanshinone biosynthesis, using r≥0.9 as a cutoff. The expression profiles for 16 of 19 randomly selected CYP450 obtained from RNA-Seq were validated by qRT-PCR. Comparing against the KEGG database, 10 CYP450 genes were found to be associated with diterpenoid biosynthesis. Considering all the evidence, 3 CYP450 genes were identified to be potentially involved in terpenoid biosynthesis. Moreover, we found that 15 CYP450 genes were possibly regulated by antisense transcripts (r≥0.9 or r≤–0.9). Lastly, a web resource (SMCYP450, http://www.herbalgenomics.org/samicyp450) was set up, which allows users to browse, search, retrieve and compare CYP450 genes and can serve as a centralized resource.

Intraspecific and heteroplasmic variations, gene losses and inversions in the chloroplast genome of Astragalus membranaceus

Astragalus membranaceus is an important medicinal plant in Asia. Several of its varieties have been used interchangeably as raw materials for commercial production. High resolution genetic markers are in urgent need to distinguish these varieties. Here, we sequenced and analyzed the chloroplast genome of A. membranaceus (Fisch.) Bunge var. mongholicus (Bunge) P.K. Hsiao using the next generation DNA sequencing technology. The genome was assembled using Abyss and then subjected to gene prediction using CPGAVAS and repeat analysis using MISA, Tandem Repeats Finder, and REPuter. Finally, the genome was subjected phylogenetic and comparative genomic analyses. The complete genome is 123,582u2009bp long, containing only one copy of the inverted repeat. Gene prediction revealed 110u2009genes encoding 76 proteins, 30u2009tRNAs, and four rRNAs. Five intra-specific hypermutation loci were identified, three of which are heteroplasmic. Furthermore, three gene losses and two large inversions were identified. Comparative genomic analyses demonstrated the dynamic nature of the Papilionoideae chloroplast genomes, which showed occurrence of numerous hypermutation loci, frequent gene losses, and fragment inversions. Results obtained herein elucidate the complex evolutionary history of chloroplast genomes and have laid the foundation for the identification of genetic markers to distinguish A. membranaceus varieties.

Complex interplay among DNA modification, noncoding RNA expression and protein-coding RNA expression in Salvia miltiorrhiza chloroplast genome.

Salvia miltiorrhiza is one of the most widely used medicinal plants. As a first step to develop a chloroplast-based genetic engineering method for the over-production of active components from S. miltiorrhiza, we have analyzed the genome, transcriptome, and base modifications of the S. miltiorrhiza chloroplast. Total genomic DNA and RNA were extracted from fresh leaves and then subjected to strand-specific RNA-Seq and Single-Molecule Real-Time (SMRT) sequencing analyses. Mapping the RNA-Seq reads to the genome assembly allowed us to determine the relative expression levels of 80 protein-coding genes. In addition, we identified 19 polycistronic transcription units and 136 putative antisense and intergenic noncoding RNA (ncRNA) genes. Comparison of the abundance of protein-coding transcripts (cRNA) with and without overlapping antisense ncRNAs (asRNA) suggest that the presence of asRNA is associated with increased cRNA abundance (p<0.05). Using the SMRT Portal software (v1.3.2), 2687 potential DNA modification sites and two potential DNA modification motifs were predicted. The two motifs include a TATA box–like motif (CPGDMM1, “TATANNNATNA”), and an unknown motif (CPGDMM2 “WNYANTGAW”). Specifically, 35 of the 97 CPGDMM1 motifs (36.1%) and 91 of the 369 CPGDMM2 motifs (24.7%) were found to be significantly modified (p<0.01). Analysis of genes downstream of the CPGDMM1 motif revealed the significantly increased abundance of ncRNA genes that are less than 400 bp away from the significantly modified CPGDMM1motif (p<0.01). Taking together, the present study revealed a complex interplay among DNA modifications, ncRNA and cRNA expression in chloroplast genome.

The complete mitochondrial genome of the white-rot fungus Ganoderma meredithiae (Polyporales, Basidiomycota)

Abstract Complete nucleotide sequence of the 78,447u2009bp mitochondrial genome of the white-rotting fungus Ganoderma meredithiae Adask. & Gilb. has been determined by next-generation sequencing technology. The circular molecule encodes a set of mitochondrial protein and RNA genes, including 15 conserved proteins, 29 tRNAs, large and small ribosomal RNAs, and 18 homing endonucleases, with a GC content of 26.14%. All structural genes are located on the same strand except trnW-CCA. Compared with previously sequenced mtDNAs of G. lucidum and G. sinense, the gene order of protein and rRNA genes among the three mitogenomes is highly conserved; however, the tRNA composition is slightly different. The mitochondrial genome of G. meredithiae will contribute to understanding the phylogeny and evolution of Ganoderma and Ganodermataceae, the group containing many species with high medicinal values.

Genome-wide Identification and Characterization of Natural Antisense Transcripts by Strand-specific RNA Sequencing in Ganoderma lucidum

Ganoderma lucidum is a white-rot fungus best-known for its medicinal and ligninolytic activities. To discover the underlying genes responsible for these activities, we identified and characterized the natural antisense transcripts (NATs) using strand-specific (ss) RNA-seq data obtained from the mycelia, primordia and fruiting bodies. NATs were identified using a custom pipeline and then subjected to functional enrichment and differential expression analyses. A total of 1613 cis- and 244 trans- sense and antisense transcripts were identified. Mapping to GO terms and KEGG pathways revealed that NATs were frequently associated with genes of particular functional categories in particular stages. ssRT-qPCR experiments showed that the expression profiles of 30 of 50 (60%) transcripts were highly correlated with those of the RNA-seq results (ru2009≥u20090.9). Expression profiles of 22 of 25 (88%) pairs of NATs and STs were highly correlated (pu2009≤u20090.01), with 15 having ru2009≥u20090.8 and 4 having ru2009≤u2009-0.8. Six lignin-modifying genes and their NATs were analyzed in detail. Diverse patterns of differential expression among different stages and positive and negative correlations were observed. These results suggested that NATs were implicated in gene expression regulation in a function-group and developmental-stage specific manner through complex mechanisms.

Angiotensin II promotes the anticoagulant effects of rivaroxaban via angiotensin type 2 receptor signaling in mice

Rivaroxaban is an oral direct factor Xa inhibitor approved for the treatment of stroke and systemic thromboembolism in patients with non-valvular atrial fibrillation. Despite its efficacy, rivaroxaban therapy results in adverse effects and complications, such as bleeding. Angiotensin II (AngII) is implicated in many cardiovascular conditions, such as hypertension and heart failure. In this study, we investigate whether AngII influences anticoagulant effects of rivaroxaban by using an experimental mouse model with type 2 diabetes mellitus and advanced glycation end product (AGE)-exposed human umbilical vein endothelial cells (HUVECs). We found that AngII promoted the anticoagulant effects of rivaroxaban in KKAy mice. The combination of rivaroxaban and AngII enhanced in vivo tissue factor pathway inhibitor (TFPI) activity and induced TFPI expression and activity in AGE-exposed HUVECs. Angiotensin type 2 receptor (AT2R) and Mas antagonists attenuated the AngII-enhanced anticoagulant action of rivaroxaban in vivo, and abolished the increased endothelial TFPI expression and activity. However, angiotensin type 1 receptor (AT1R) antagonist exerted no effects. Additionally, combination of rivaroxaban and AngII induced aortic AT2R and Mas expression. Our data suggest that the anticoagulant effects of rivaroxaban are promoted by AngII via AT2R and Mas signaling. These findings are significant for the clinical administration of rivaroxaban.

Diversity of mitochondrial plastid DNAs (MTPTs) in seed plants

Abstract Mitochondrial plastid DNAs (MTPTs) refer to plastid-derived DNA fragments in mitochondrial genomes. While the MTPTs have been described for numerous species, its overall patterns have not been examined in details. Here, we carried out a systematic analysis of MTPTs among 73 plant species, including 28 algae, 1 liverwort, 2 moss, 1 lycophyte, 1 gymnosperm, 1 magnoliid, 12 monocots, 26 eudicots and 1 relic angiosperm Amborella trichopoda. A total of 300 MTPT gene clusters were found in 39 seed plants, which represented 144 MTPT gene cluster types. The detected MTPT gene clusters were evaluated in seven aspects, and they were found to be enriched particularly in monocots and asterids of eudicots. Some MTPT gene clusters were found to be shared by closely related species. All chloroplast genes were found in MTPTs, suggesting that there is no functional relevancy for genes that were transferred. However, after calculation of the frequency of the 115 chloroplast genes, five hot spots and three cold spots were discovered in chloroplast genome. In summary, this study demonstrated the high degree of diversity in MTPTs. The discovered MTPTs would facilitate the accurate assembly of chloroplast and mitochondrial genomes as well as the understanding of organelle genome evolution.

Macrophage CGI-58 Attenuates Inflammatory Responsiveness via Promotion of PPARγ Signaling.

Background/Aims: Comparative gene identification-58 (CGI-58), an adipose triglyceride lipase (ATGL) coactivator, strongly promotes ATGL-mediated triglyceride (TG) catabolism. Beyond its function in promoting lipolysis, other features of CGI-58 have been proposed. Here, we investigated the role of CGI-58 in the regulation of inflammatory responsiveness in macrophages. Methods: Macrophage-specific GCI-58 transgenic mice (TG) and wild type mice (WT) were fed a high fat diet (HFD), and RAW264.7 cells were treated with lipopolysaccharide (LPS). The peroxisome proliferator-activated receptor (PPAR) signaling was detected. The inflammatory responsiveness and mitochondrial function were examined. Results: TG mice showed lower serum levels of proinflammatory cytokines and better mitochondrial function in macrophages compared with WT control. Knockdown of CGI-58 in RAW264.7 cells aggravated LPS-induced inflammation and mitochondrial dysfunction. CGI-58 overexpression and silencing in macrophages induced and inhibited PPARγ expression and activity, respectively. Most importantly, the PPARγ-specific agonist rosiglitazone significantly suppressed inflammation and mitochondrial dysfunction induced by CGI-58 deficiency. Furthermore, knockdown of PPARγ in macrophages significantly dampened the role of CGI-58 in suppression of inflammation and mitochondrial dysfunction. Interestingly, CGI-58 inhibited histone deacetylation and the recruitment of histone deacetylase (HDAC) to the PPARγ promoter. Finally, ATGL deficiency did not affect inflammatory responsiveness and PPARγ signaling in macrophages. Conclusion: These results demonstrate that macrophage CGI-58 enhances PPARγ signaling and thus suppresses inflammatory responsiveness and mitochondrial dysfunction.