Global Pharmacopoeia Genome Database is an integrated and mineable genomic database for traditional medicines derived from eight international pharmacopoeias

Genome-wide identification of YABBY gene family and its expression pattern analysis in Astragalus mongholicus

During plant growth and development, the YABBY gene plays a crucial role in the morphological structure, hormone signaling, stress resistance, crop breeding, and agricultural production of plant lateral organs, leaves, flowers, and fruits. Astragalus mongholicus is a perennial herbaceous plant in the legume family, widely used worldwide due to its high medicinal and edible value. However, there have been no reports of the YABBY gene family in A. mongholicus. This study used bioinformatics methods, combined with databases and analysis websites, to systematically analyze the AmYABBY gene family in the entire genome of A. mongholicus and verified its expression patterns in different tissues of A. mongholicus through transcriptome data and qRT-PCR experiments. A total of seven AmYABBY genes were identified, which can be divided into five subfamilies and distributed on three chromosomes. Two pairs of AmYABBY genes may be involved in fragment duplication on three chromosomes. All AmYABBY proteins have a zinc finger YABBY domain, and members of the same group have similar motif composition and intron - exon structure. In the promoter region of the genes, light-responsive and MeJa-response cis-elements are dominant. AmYABBY is highly expressed in stems and leaves, especially AmYABBY1, AmYABBY2, and AmYABBY3, which play important roles in the growth and development of stems and leaves. The AmYABBY gene family regulates the growth and development of A. mongholicus. In summary, this study provides a theoretical basis for in-depth research on the function of the AmYABBY gene and new insights into the molecular response mechanism of the growth and development of the traditional Chinese medicine A. mongholicus.

Genome-wide identification and expression pattern analysis of the GRF transcription factor family in Astragalus mongholicus

BACKGROUND

Astragalus membranaceus is a plant of the Astragalus genus, which is used as a traditional Chinese herbal medicine with extremely high medicinal and edible value. Astragalus mongholicus, as one of the representative medicinal materials with the same origin of medicine and food, has a rising market demand for its raw materials, but the quality is different in different production areas. Growth-regulating factors (GRF) are transcription factors unique to plants that play important roles in plant growth and development. Up to now, there is no report about GRF in A. mongholicus.

METHODS AND RESULTS

This study conducted a genome-wide analysis of the AmGRF gene family, identifying a total of nine AmGRF genes that were classified into subfamily V based on phylogenetic relationships. In the promoter region of the AmGRF gene, we successfully predicted cis-elements that respond to abiotic stress, growth, development, and hormone production in plants. Based on transcriptomic data and real-time quantitative polymerase chain reaction (qPCR) validation, the results showed that AmGRFs were expressed in the roots, stems, and leaves, with overall higher expression in leaves, higher expression of AmGRF1 and AmGRF8 in roots, and high expression levels of AmGRF1 and AmGRF9 in stems.

CONCLUSIONS

The results of this study provide a theoretical basis for the further exploration of the functions of AmGRFs in plant growth and development.

Pogostemon cablin essential oil affects anxiety- and depressive-like behaviors and the gut microbiota in chronic unpredictable mild stress model rats

The gut microbiota is thought to be an important factor that influences brain processes and behaviors through the gut-brain axis. Pogostemon cablin is used in traditional Chinese medicine (TCM) to treat gastrointestinal symptoms. Patchouli essential oil (PCO), the main active agent in P. cablin, is used in aromatherapy for stress relief. The aim of our study was to investigate the effects of orally administered PCO on anxiety- and depressive-like behaviors and the gut microbiota. We constructed a rat model of chronic unpredictable mild stress (CUMS) and explored the anxiolytic- and antidepressant-like effects of PCO using the open field test (OFT) and forced swim test (FST). Changes in the abundance of the gut microbiota, short-chain fatty acids (SCFAs), and other related molecules were assessed to determine the role of the gut microbiota. Our results showed that CUMS induced an anxiety-like phenotype in the OFT, which was reversed by PCO, and that PCO also significantly mitigated the depression-like behaviors caused by CUMS in the FST. Furthermore, we found that PCO increased the relative abundances of several probiotics, including Bacteroides and Blautia, and decreased the relative abundances of Ruminococcus_1 and Ruminococcus_2, which were increased by CUMS. Regarding SCFAs, the metabolites of the gut microbiota, PCO increased the concentration of propionic acid and decreased that of caproic acid. Finally, PCO restored the serotonin (5-hydroxytryptamine, 5-HT) level in the hippocampus, which had been decreased by CUMS. The results of this study suggested that PCO can improve stress-related anxiety- and depression-like behaviors and might exert its effects on the central nervous system through interactions with the gut microbiota.

Whole-genome sequencing in medicinal plants: current progress and prospect

Advancements in genomics have dramatically accelerated the research on medicinal plants, and the development of herbgenomics has promoted the "Project of 1K Medicinal Plant Genome" to decipher their genetic code. However, it is difficult to obtain their high-quality whole genomes because of the prevalence of polyploidy and/or high genomic heterozygosity. Whole genomes of 123 medicinal plants were published until September 2022. These published genome sequences were investigated in this review, covering their classification, research teams, ploidy, medicinal functions, and sequencing strategies. More than 1,000 institutes or universities around the world and 50 countries are conducting research on medicinal plant genomes. Diploid species account for a majority of sequenced medicinal plants. The whole genomes of plants in the Poaceae family are the most studied. Almost 40% of the published papers studied species with tonifying, replenishing, and heat-cleaning medicinal effects. Medicinal plants are still in the process of domestication as compared with crops, thereby resulting in unclear genetic backgrounds and the lack of pure lines, thus making their genomes more difficult to complete. In addition, there is still no clear routine framework for a medicinal plant to obtain a high-quality whole genome. Herein, a clear and complete strategy has been originally proposed for creating a high-quality whole genome of medicinal plants. Moreover, whole genome-based biological studies of medicinal plants, including breeding and biosynthesis, were reviewed. We also advocate that a research platform of model medicinal plants should be established to promote the genomics research of medicinal plants.

Genome-wide characterization of regulator of chromosome condensation 1 (RCC1) gene family in Artemisia annua L. revealed a conservation evolutionary pattern

BACKGROUND

Artemisia annua is the major source for artemisinin production. The artemisinin content in A. annua is affected by different types of light especially the UV light. UVR8, a member of RCC1 gene family was found to be the UV-B receptor in plants. The gene structures, evolutionary history and expression profile of UVR8 or RCC1 genes remain undiscovered in A. annua.

RESULTS

Twenty-two RCC1 genes (AaRCC1) were identified in each haplotype genome of two diploid strains of A. annua, LQ-9 and HAN1. Varied gene structures and sequences among paralogs were observed. The divergence of most RCC1 genes occurred at 46.7 - 51 MYA which overlapped with species divergence of core Asteraceae during the Eocene, while no recent novel RCC1 members were found in A. annua genome. The number of RCC1 genes remained stable among eudicots and RCC1 genes underwent purifying selection. The expression profile of AaRCC1 is analogous to that of Arabidopsis thaliana (AtRCC1) when responding to environmental stress.

CONCLUSIONS

This study provided a comprehensive characterization of the AaRCC1 gene family and suggested that RCC1 genes were conserved in gene number, structures, constitution of amino acids and expression profiles among eudicots.

Herbgenomics meets Papaveraceae: a promising -omics perspective on medicinal plant research

Herbal medicines were widely used in ancient and modern societies as remedies for human ailments. Notably, the Papaveraceae family includes well-known species, such as Papaver somniferum and Chelidonium majus, which possess medicinal properties due to their latex content. Latex-bearing plants are a rich source of diverse bioactive compounds, with applications ranging from narcotics to analgesics and relaxants. With the advent of high-throughput technologies and advancements in sequencing tools, an opportunity exists to bridge the knowledge gap between the genetic information of herbs and the regulatory networks underlying their medicinal activities. This emerging discipline, known as herbgenomics, combines genomic information with other -omics studies to unravel the genetic foundations, including essential gene functions and secondary metabolite biosynthesis pathways. Furthermore, exploring the genomes of various medicinal plants enables the utilization of modern genetic manipulation techniques, such as Clustered Regularly-Interspaced Short Palindromic Repeats (CRISPR/Cas9) or RNA interference. This technological revolution has facilitated systematic studies of model herbs, targeted breeding of medicinal plants, the establishment of gene banks and the adoption of synthetic biology approaches. In this article, we provide a comprehensive overview of the recent advances in genomic, transcriptomic, proteomic and metabolomic research on species within the Papaveraceae family. Additionally, it briefly explores the potential applications and key opportunities offered by the -omics perspective in the pharmaceutical industry and the agrobiotechnology field.

Comparative chloroplast genomes provided insights into the evolution and species identification on the Datureae plants

Generally, chloroplast genomes of angiosperms are always highly conserved but carry a certain number of variation among species. In this study, chloroplast genomes of 13 species from Datureae tribe that are of importance both in ornamental gardening and medicinal usage were studied. In addition, seven chloroplast genomes from Datureae together with two from Solanaceae species retrieved from the National Center for Biotechnology Information (NCBI) were integrated into this study. The chloroplast genomes ranged in size from 154,686 to 155,979 and from 155,497 to 155,919 bp for species of Datura and Brugmansia, respectively. As to Datura and Brugmansia, a total of 128 and 132 genes were identified, in which 83 and 87 protein coding genes were identified, respectively; Furthermore, 37 tRNA genes and 8 rRNA genes were both identified in Datura and Brugmansia. Repeats analysis indicated that the number and type varied among species for Simple sequence repeat (SSR), long repeats, and tandem repeats ranged in number from 53 to 59, 98 to 99, and 22 to 30, respectively. Phylogenetic analysis based on the plastid genomes supported the monophyletic relationship among Datura and Brugmansia and Trompettia, and a refined phylogenic relationships among each individual was resolved. In addition, a species-specific marker was designed based on variation spot that resulted from a comparative analysis of chloroplast genomes and verified as effective maker for identification of D. stramonium and D. stramonium var. inermis. Interestingly, we found that 31 genes were likely to be under positive selection, including genes encoding ATP protein subunits, photosystem protein subunit, ribosome protein subunits, NAD(P)H dehydrogenase complex subunits, and clpP, petB, rbcL, rpoCl, ycf4, and cemA genes. These genes may function as key roles in the adaption to diverse environment during evolution. The diversification of Datureae members was dated back to the late Oligocene periods. These chloroplast genomes are useful genetic resources for taxonomy, phylogeny, and evolution for Datureae.

Manual correction of genome annotation improved alternative splicing identification of Artemisia annua

Gene annotation is essential for genome-based studies. However, algorithm-based genome annotation is difficult to fully and correctly reveal genomic information, especially for species with complex genomes. Artemisia annua L. is the only commercial resource of artemisinin production though the content of artemisinin is still to be improved. Genome-based genetic modification and breeding are useful strategies to boost artemisinin content and therefore, ensure the supply of artemisinin and reduce costs, but better gene annotation is urgently needed. In this study, we manually corrected the newly released genome annotation of A. annua using second- and third-generation transcriptome data. We found that incorrect gene information may lead to differences in structural, functional, and expression levels compared to the original expectations. We also identified alternative splicing events and found that genome annotation information impacted identifying alternative splicing genes. We further demonstrated that genome annotation information and alternative splicing could affect gene expression estimation and gene function prediction. Finally, we provided a valuable version of A. annua genome annotation and demonstrated the importance of gene annotation in future research.

Chromatin architecture of two different strains of Artemisia annua reveals the alterations in interaction and gene expression

MAIN CONCLUSION

The hierarchical architecture of chromatins affects the gene expression level of glandular secreting trichomes and the artemisinin biosynthetic pathway-related genes, consequently bringing on huge differences in the content of artemisinin and its derivatives of A. annua. The plant of traditional Chinese medicine "Qinghao" is called Artemisia annua L. in Chinese Pharmacopoeia. High content and the total amount of artemisinin is the main goal of A. annua breeding, nevertheless, the change of chromatin organization during the artemisinin synthesis process has not been discovered yet. This study intended to find the roles of chromatin structure in the production of artemisinin through bioinformatics and experimental validation. Chromosome conformation capture analysis was used to scrutinize the interactions among chromosomes and categorize various scales of chromatin during artemisinin synthesis in A. annua. To confirm the effect of the changes in chromatin structure, Hi-C and RNA-sequencing were performed on two different strains to find the correlation between chromatin structure and gene expression levels on artemisinin synthesis progress and regulation. Our results revealed that the frequency of intra-chromosomal interactions was higher in the inter-chromosomal interactions between the root and leaves on a high artemisinin production strain (HAP) compared to a low artemisinin production strain (LAP). We found that compartmental transition was connected with interactions among different chromatins. Interestingly, glandular secreting trichomes (GSTs) and the artemisinin biosynthetic pathway (ABP) related genes were enriched in the areas which have the compartmental transition, reflecting the regulation of artemisinin synthesis. Topologically associated domain boundaries were associated with various distributions of genes and expression levels. Genes associated with ABP and GST in the adjacent loop were highly expressed, suggesting that epigenetic regulation plays an important role during artemisinin synthesis and glandular secreting trichomes production process. Chromatin structure could show an important status in the mechanisms of artemisinin synthesis process in A. annua.

Functional characterization of nine critical genes encoding rate-limiting enzymes in the flavonoid biosynthesis of the medicinal herb Grona styracifolia

Grona styracifolia is a photophilous legume that contains abundant flavonoids with multiple pharmacological activities, which is used to cure urethral and biliary calculus in China for thousands of years. The authentication of the rate-limiting enzymes involved in the flavonoids biosynthesis pathway enabled a better understanding of the molecular aspect of quality formation and modulation of this medicinal herb. In this study, the chemical distribution characteristics and content of flavonoids in different tissues of Grona styracifolia were analyzed using ultraperormance liquid chromatography coupled with Q-TOF mass spectrometry and showed that active flavonoids were primarily synthesized and stored in the leaves. Subsequently, RNA sequencing (RNA-seq)-based transcriptome profiling of the different tissues revealed that the flavonoids biosynthesis in the leaves was the most active. Meanwhile, 27 full-length transcripts inferred encoding vital enzymes involved in the flavonoids biosynthesis were preliminarily excavated. Finally, four CHSs, four CHIs, and one FNSII were successfully characterized by heterologous expression, which involved in three rate-limiting steps of the flavonoid biosynthetic pathway. In conclusion, these results laid a foundation for further investigation of the molecular mechanism of the biosynthesis and modulation of active flavonoids in Grona styracifolia.

A reference-grade genome assembly for Astragalus mongholicus and insights into the biosynthesis and high accumulation of triterpenoids and flavonoids in its roots

Astragalus membranaceus var. mongholicus (AMM), a member of the Leguminosae, is one of the most important medicinal plants worldwide. The dried roots of AMM have a wide range of pharmacological effects and are a traditional Chinese medicine. Here, we report the first chromosome-level reference genome of AMM, comprising nine pseudochromosomes with a total size of 1.47 Gb and 27 868 protein-encoding genes. Comparative genomic analysis reveals that AMM has not experienced an independent whole-genome duplication (WGD) event after the WGD event shared by the Papilionoideae species. Analysis of long terminal repeat retrotransposons suggests a recent burst of these elements at approximately 0.13 million years ago, which may explain the large size of the AMM genome. Multiple gene families involved in the biosynthesis of triterpenoids and flavonoids were expanded, and our data indicate that tandem duplication has been the main driver for expansion of these families. Among the expanded families, the phenylalanine ammonia-lyase gene family was primarily expressed in the roots of AMM, suggesting their roles in the biosynthesis of phenylpropanoid compounds. The functional versatility of 2,3-oxidosqualene cyclase genes in cluster III may play a critical role in the diversification of triterpenoids in AMM. Our findings provide novel insights into triterpenoid and flavonoid biosynthesis and can facilitate future research on the genetics and medical applications of AMM.

Orally administered BZL-sRNA-20 oligonucleotide targeting TLR4 effectively ameliorates acute lung injury in mice

The global COVID-19 pandemic emerged at the end of December 2019. Acute respiratory distress syndrome (ARDS) and acute lung injury (ALI) are common lethal outcomes of bacterial lipopolysaccharide (LPS), avian influenza virus, and SARS-CoV-2. Toll-like receptor 4 (TLR4) is a key target in the pathological pathway of ARDS and ALI. Previous studies have reported that herbal small RNAs (sRNAs) are a functional medical component. BZL-sRNA-20 (Accession number: B59471456; Family ID: F2201.Q001979.B11) is a potent inhibitor of Toll-like receptor 4 (TLR4) and pro-inflammatory cytokines. Furthermore, BZL-sRNA-20 reduces intracellular levels of cytokines induced by lipoteichoic acid (LTA) and polyinosinic-polycytidylic acid (poly (I:C)). We found that BZL-sRNA-20 rescued the viability of cells infected with avian influenza H5N1, SARS-CoV-2, and several of its variants of concern (VOCs). Acute lung injury induced by LPS and SARS-CoV-2 in mice was significantly ameliorated by the oral medical decoctosome mimic (bencaosome; sphinganine (d22:0)+BZL-sRNA-20). Our findings suggest that BZL-sRNA-20 could be a pan-anti-ARDS ALI drug.

Integration of high-throughput omics technologies in medicinal plant research: The new era of natural drug discovery

Medicinal plants are natural sources to unravel novel bioactive compounds to satisfy human pharmacological potentials. The world's demand for herbal medicines is increasing year by year; however, large-scale production of medicinal plants and their derivatives is still limited. The rapid development of modern technology has stimulated multi-omics research in medicinal plants, leading to a series of breakthroughs on key genes, metabolites, enzymes involved in biosynthesis and regulation of active compounds. Here, we summarize the latest research progress on the molecular intricacy of medicinal plants, including the comparison of genomics to demonstrate variation and evolution among species, the application of transcriptomics, proteomics and metabolomics to explore dynamic changes of molecular compounds, and the utilization of potential resources for natural drug discovery. These multi-omics research provide the theoretical basis for environmental adaptation of medicinal plants and allow us to understand the chemical diversity and composition of bioactive compounds. Many medicinal herbs' phytochemical constituents and their potential health benefits are not fully explored. Given their large diversity and global distribution as well as the impacts of growth duration and environmental factors on bioactive phytochemicals in medicinal plants, it is crucial to emphasize the research needs of using multi-omics technologies to address basic and applied problems in medicinal plants to aid in developing new and improved medicinal plant resources and discovering novel medicinal ingredients.

Essential Oils and Their Compounds as Potential Anti-Influenza Agents

Essential oils (EOs) are chemical substances, mostly produced by aromatic plants in response to stress, that have a history of medicinal use for many diseases. In the last few decades, EOs have continued to gain more attention because of their proven therapeutic applications against the flu and other infectious diseases. Influenza (flu) is an infectious zoonotic disease that affects the lungs and their associated organs. It is a public health problem with a huge health burden, causing a seasonal outbreak every year. Occasionally, it comes as a disease pandemic with unprecedentedly high hospitalization and mortality. Currently, influenza is managed by vaccination and antiviral drugs such as Amantadine, Rimantadine, Oseltamivir, Peramivir, Zanamivir, and Baloxavir. However, the adverse side effects of these drugs, the rapid and unlimited variabilities of influenza viruses, and the emerging resistance of new virus strains to the currently used vaccines and drugs have necessitated the need to obtain more effective anti-influenza agents. In this review, essential oils are discussed in terms of their chemistry, ethnomedicinal values against flu-related illnesses, biological potential as anti-influenza agents, and mechanisms of action. In addition, the structure-activity relationships of lead anti-influenza EO compounds are also examined. This is all to identify leading agents that can be optimized as drug candidates for the management of influenza. Eucalyptol, germacrone, caryophyllene derivatives, eugenol, terpin-4-ol, bisabolene derivatives, and camphecene are among the promising EO compounds identified, based on their reported anti-influenza activities and plausible molecular actions, while nanotechnology may be a new strategy to achieve the efficient delivery of these therapeutically active EOs to the active virus site.

Comparative chloroplast genome analyses of Amomum: insights into evolutionary history and species identification

BACKGROUND

Species in genus Amomum always have important medicinal and economic values. Classification of Amomum using morphological characters has long been a challenge because they exhibit high similarity. The main goals of this study were to mine genetic markers from cp genomes for Amomum species identification and discover their evolutionary history through comparative analysis.

RESULTS

Three species Amomum villosum, Amomum maximum and Amomum longipetiolatum were sequenced and annotated for the complete chloroplast (cp) genomes, and the cp genomes of A. longipetiolatum and A. maximum were the first reported. Three cp genomes exhibited typical quadripartite structures with 163,269-163,591 bp in length. Each genome encodes 130 functional genes including 79 protein-coding, 26 tRNAs and 3 rRNAs genes. 113-152 SSRs and 99 long repeats were identified in the three cp genomes. By designing specific primers, we amplified the highly variable loci and the mined genetic marker ccsA exhibited a relatively high species identification resolution in Amomum. The nonsynonymous and synonymous substitution ratios (Ka/Ks) in Amomum and Alpinia showed that most genes were subjected to a purifying selection. Phylogenetic analysis revealed the evolutionary relationships of Amomum and Alpinia species and proved that Amomum is paraphyletic. In addition, the sequenced sample of A. villosum was found to be a hybrid, becoming the first report of natural hybridization of this genus. Meanwhile, the high-throughput sequencing-based ITS2 analysis was proved to be an efficient tool for interspecific hybrid identification and with the help of the chloroplast genome, the hybrid parents can be also be determined.

CONCLUSION

The comparative analysis and mined genetic markers of cp genomes were conducive to species identification and evolutionary relationships of Amomum.

Natural products of medicinal plants: biosynthesis and bioengineering in post-genomic era

Globally, medicinal plant natural products (PNPs) are a major source of substances used in traditional and modern medicine. As we human race face the tremendous public health challenge posed by emerging infectious diseases, antibiotic resistance and surging drug prices etc., harnessing the healing power of medicinal plants gifted from mother nature is more urgent than ever in helping us survive future challenge in a sustainable way. PNP research efforts in the pre-genomic era focus on discovering bioactive molecules with pharmaceutical activities, and identifying individual genes responsible for biosynthesis. Critically, systemic biological, multi- and inter-disciplinary approaches integrating and interrogating all accessible data from genomics, metabolomics, structural biology, and chemical informatics are necessary to accelerate the full characterization of biosynthetic and regulatory circuitry for producing PNPs in medicinal plants. In this review, we attempt to provide a brief update on the current research of PNPs in medicinal plants by focusing on how different state-of-the-art biotechnologies facilitate their discovery, the molecular basis of their biosynthesis, as well as synthetic biology. Finally, we humbly provide a foresight of the research trend for understanding the biology of medicinal plants in the coming decades.

Allele-aware chromosome-level genome assembly of Artemisia annua reveals the correlation between ADS expansion and artemisinin yield

Artemisia annua is the major natural source of artemisinin, an anti-malarial medicine commonly used worldwide. Here, we present chromosome-level haploid maps for two A. annua strains with different artemisinin contents to explore the relationships between genomic organization and artemisinin production. High-fidelity sequencing, optical mapping, and chromatin conformation capture sequencing were used to assemble the heterogeneous and repetitive genome and resolve the haplotypes of A. annua. Approximately 50,000 genes were annotated for each haplotype genome, and a triplication event that occurred approximately 58.12 million years ago was examined for the first time in this species. A total of 3,903,467-5,193,414 variants (SNPs, indels, and structural variants) were identified in the 1.5-Gb genome during pairwise comparison between haplotypes, consistent with the high heterozygosity of this species. Genomic analyses revealed a correlation between artemisinin concents and the copy number of amorpha-4,11-diene synthase genes. This correlation was further confirmed by resequencing of 36 A. annua samples with varied artemisinin contents. Circular consensus sequencing of transcripts facilitated the detection of paralog expression. Collectively, our study provides chromosome-level allele-aware genome assemblies for two A. annua strains and new insights into the biosynthesis of artemisinin and its regulation, which will contribute to conquering malaria worldwide.

Gllac7 Is Induced by Agricultural and Forestry Residues and Exhibits Allelic Expression Bias in Ganoderma lucidum

Ganoderma lucidum has a wide carbon spectrum, while the expression profile of key genes relevant to carbon metabolism on different carbon sources has been seldom studied. Here, the transcriptomes of G. lucidum mycelia cultured on each of 19 carbon sources were conducted. In comparison with glucose, 16 to 1,006 genes were upregulated and 7 to 1,865 genes were downregulated. Significant gene expression dynamics and induced activity were observed in laccase genes when using agricultural and forestry residues (AFRs) as solo carbon sources. Furthermore, study of laccase gene family in two haploids of G. lucidum GL0102 was conducted. Totally, 15 and 16 laccase genes were identified in GL0102_53 and GL0102_8, respectively, among which 15 pairs were allelic genes. Gene structures were conserved between allelic laccase genes, while sequence variations (most were SNPs) existed. Nine laccase genes rarely expressed on all the tested carbon sources, while the other seven genes showed high expression level on AFRs, especially Gllac2 and Gllac7, which showed 5- to 1,149-fold and 4- to 94-fold upregulation in mycelia cultured for 5 days, respectively. The expression of H53lac7 was consistently higher than that of H8lac7_1 on all the carbon sources except XM, exhibiting a case of allelic expression bias. A total of 47 SNPs and 3 insertions/deletions were observed between promoters of H53lac7 and H8lac7_1, which lead to differences in predicted binding sites of zinc fingers. These results provide scientific data for understanding the gene expression profile and regulatory role on different carbon sources and may support further functional research of laccase.

Strategies for molecular authentication of herbal products: from experimental design to data analysis

Molecular herbal authentication has gained worldwide popularity in the past decade. DNA-based methods, including DNA barcoding and species-specific amplification, have been adopted for herbal identification by various pharmacopoeias. Development of next-generating sequencing (NGS) drastically increased the throughput of sequencing process and has sped up sequence collection and assembly of organelle genomes, making more and more reference sequences/genomes available. NGS allows simultaneous sequencing of multiple reads, opening up the opportunity of identifying multiple species from one sample in one go. Two major experimental approaches have been applied in recent publications of identification of herbal products by NGS, the PCR-dependent DNA metabarcoding and PCR-free genome skimming/shotgun metagenomics. This review provides a brief introduction of the use of DNA metabarcoding and genome skimming/shotgun metagenomics in authentication of herbal products and discusses some important considerations in experimental design for botanical identification by NGS, with a specific focus on quality control, reference sequence database and different taxon assignment programs. The potential of quantification or abundance estimation by NGS is discussed and new scientific findings that could potentially interfere with accurate taxon assignment and/or quantification is presented.

Endophytic bacterial and fungal community compositions in different organs of ginseng (Panax ginseng)

Panax ginseng (Panax ginseng C. A. Mey.) is a perennial herb of the genus ginseng, which is used as medicine with dried roots and rhizomes. With the deepening of research on ginseng, the chemical components and pharmacological effects of ginseng have gradually been discovered. Endophytes are beneficial to host plants. However, the composition of endophytes in different organs from ginseng is poorly elucidated. The report of ginsenoside production by endophytic microbes isolated from Panax sp., motivated us to explore the endophytic microbial diversity related to the roots, stems, and leaves. In this study, the V5-V7 variable region of endophytic bacteria 16S rRNA gene and V1 variable region of endophytic fungi ITS gene in different organs were analyzed by high-throughput sequencing. The diversity and abundance of endophytic microbes in the three organs are different and are affected by the organs. For example, the most abundant endophytic bacterial genus in roots was Mycobacterium, while, the stems and leaves were Ochrobactrum. Similarly, the fungal endophytes, Coniothyrium and Cladosporium, were also found in high abundance in stems, in comparison to roots and leaves. The Shannon index shows that the diversity of endophytic bacteria in roots is the highest, and the richness of endophytic bacterial was root > stem (p < 0.05). Principal coordinate analysis showed that there were obvious microbial differences among the three groups, and the endophytic bacterial composition of the leaves was closer to that of the roots. This study provides an important reference for the study of endophytic microorganisms in ginseng.

Genome-Wide Identification, Characterization and Expression Analysis of Lipoxygenase Gene Family in Artemisia annua L

Lipoxygenase (LOX) is a ubiquitous oxygenase found in animals and plants and plays a pivotal role in diverse biological processes, including defense and development. Artemisinin, which can only be obtained from Artemisia annua L., is the most effective therapeutic drug for malaria without serious side effects. This study identified and analyzed LOX gene family members in the A. annua genome at the chromosomal level. Twenty LOX genes with various molecular weights, isoelectric points, and amino acid numbers were identified and named AaLOX, which were located in the cytoplasm or chloroplast. The average protein length of all AaLOX was 850 aa. Phylogenetic tree analysis revealed that the AaLOX was divided into two major groups, 9-LOX and 13-LOX. The exon numbers ranged from 1 to 12, indicating that different AaLOX genes have different functions. The secondary structure was mainly composed of alpha helix and random coil, and the tertiary structure was similar for most AaLOX. Upstream promoter region analysis revealed that a large number of cis-acting elements were closely related to plant growth and development, light response, hormone, and other stress responses. Transcriptome data analysis of different tissues suggested that the gene family was differently expressed in the roots, stems, leaves, and flowers of two A. annua strains HAN1 and LQ9. qRT-PCR confirmed that AaLOX5 and AaLOX17 had the highest expression in flowers and leaves. This study provides a theoretical basis for the further functional analysis of the AaLOX gene family.

TCMPG: an integrative database for traditional Chinese medicine plant genomes

Because of their great therapeutic and economic value, medicinal plants have attracted increasing scientific attention. With the rapid development of high-throughput sequencing technology, the genomes of many medicinal plants have been sequenced. Storing and analyzing the increasing volume of genomic data has become an urgent task. To solve this challenge, we have proposed the Traditional Chinese Medicine Plant Genome database (TCMPG, http://cbcb.cdutcm.edu.cn/TCMPG/), an integrative database for storing the scattered genomes of medicinal plants. TCMPG currently includes 160 medicinal plants, 195 corresponding genomes, and 255 herbal medicines. Detailed information on plant species, genomes, and herbal medicines is also integrated into TCMPG. Popular genomic analysis tools are embedded in TCMPG to facilitate the systematic analysis of medicinal plants. These include BLAST for identifying orthologs from different plants, SSR Finder for identifying simple sequence repeats, JBrowse for browsing genomes, Synteny Viewer for displaying syntenic blocks between two genomes, and HmmSearch for identifying protein domains. TCMPG will be continuously updated by integrating new data and tools for comparative and functional genomic analysis.

Licorice Germplasm Resources Identification Using DNA Barcodes Inner-Variants

Based on the gradual transformation from wild growth to artificial cultivation, the accurate authentication of licorice seeds contributes to the first committed step of its quality control and is pivotal to ensure the clinical efficacy of licorice. However, it is still challenging to obtain genetically stable licorice germplasm resources due to the multi-source, multi-heterozygous, polyploid, and hybrid characteristics of licorice seeds. Here, a new method for determining the heterozygosity of licorice seed mixture, based on the various sites, and finding the composition characteristics of licorice seed is preliminarily designed and proposed. Namely, high-throughput full-length multiple DNA barcodes(HFMD), based on ITS multi-copy variation exist, the full-length amplicons of ITS2, psbA-trnH and ITS are directly sequenced by rDNA through the next-generation sequence(NGS) and single-molecule real-time (SMRT) technologies. By comparing the three sequencing methods, our results proved that SMRT sequencing successfully identified the complete gradients of complex mixed samples with the best performance. Meanwhile, HFMD is a brilliant and feasible method for evaluating the heterozygosity of licorice seeds. It shows a perfect interpretation of DNA barcoding and can be applied in multi-base multi-heterozygous and polyploid species.

Hierarchical chromatin features reveal the toxin production in Bungarus multicinctus

BACKGROUND

Bungarus multicinctus, from which a classical Chinese medicine is produced, is known as the most venomous land snake in the world, but the chromatin organization and transcription factor activity during venom replenishment progress have not been explored yet. This study aimed to determine the roles of chromatin structure in toxin activity via bioinformatics and experimental validation.

METHODS

Chromosome conformation capture (Hi-C) analysis was used to examine interactions among chromosomes and identify different scales of chromatin during envenomation in B. multicinctus. Correlations between epigenetic modifications and chromatin structure were verified through ChIP-seq analysis. RNA-seq was used to validate the influence of variation in chromatin structure and gene expression levels on venom production and regulation.

RESULTS

Our results suggested that intra-chromosomal interactions are more intense than inter-chromosomal interactions among the control group, 3-day group of venom glands and muscles. Through this, we found that compartmental transition was correlated with chromatin interactions. Interestingly, the up-regulated genes in more compartmental switch regions reflect the function of toxin activity. Topologically associated domain (TAD) boundaries enriched with histone modifications are associated with different distributions of genes and the expression levels. Toxin-coding genes in the same loop are highly expressed, implying that the importance of epigenetic regulation during envenomination. On a smaller scale, the epigenetic markers affect transcriptional regulation by controlling the recruitment/inhibition of transcription initiation complexes.

CONCLUSIONS

Chromatin structure and epigenetic modifications could play a vital status role in the mechanisms of venom regulation in B. multicinctus.

Haploid Genome Analysis Reveals a Tandem Cluster of Four HSP20 Genes Involved in the High-Temperature Adaptation of Coriolopsis trogii

Coriolopsis trogii is a typical thermotolerant basidiomycete fungus, but its thermotolerance mechanisms are currently unknown. In this study, two monokaryons of C. trogii strain Ct001 were assembled: Ct001_29 had a genome assembly size of 38.85 Mb and encoded 13,113 genes, while Ct001_31 was 40.19 Mb in length and encoded 13,309 genes. Comparative intra- and interstrain genomic analysis revealed the rich genetic diversity of C. trogii, which included more than 315,194 single-nucleotide polymorphisms (SNPs), 30,387 insertion/deletions (indels), and 1,460 structural variations. Gene family analysis showed that the expanded families of C. trogii were functionally enriched in lignocellulose degradation activities. Furthermore, a total of 14 allelic pairs of heat shock protein 20 (HSP20) genes were identified in the C. trogii genome. The expression profile obtained from RNA sequencing (RNA-Seq) showed that four tandem-duplicated allelic pairs, HSP20.5 to HSP20.8, had more than 5-fold higher expression at 35°C than at 25°C. In particular, HSP20.5 and HSP20.8 were the most highly expressed HSP20 genes. Allelic expression bias was found for HSP20.5 and HSP20.8; the expression of Ct29HSP20.8 was at least 1.34-fold higher than that of Ct31HSP20.8, and that of Ct31HSP20.5 was at least 1.5-fold higher than that of Ct29HSP20.5. The unique structural and expression profiles of the HSP20 genes revealed by these haplotype-resolved genomes provide insight into the molecular mechanisms of high-temperature adaptation in C. trogii. IMPORTANCE Heat stress is one of the most frequently encountered environmental stresses for most mushroom-forming fungi. Currently available fungal genomes are mostly haploid because high heterozygosity hinders diploid genome assembly. Here, two haplotype genomes of C. trogii, a thermotolerant basidiomycete, were assembled separately. A conserved tandem cluster of four HSP20 genes showing allele-specific expression was found to be closely related to high-temperature adaptation in C. trogii. The obtained haploid genomes and their comparison offer a more thorough understanding of the genetic background of C. trogii. In addition, the responses of HSP20 genes at 35°C, which may contribute to the growth and survival of C. trogii at high temperatures, could inform the selection and breeding of elite strains in the future.

Coriolopsis trogii hydrophobin genes favor a clustering distribution and are widely involved in mycelial growth and primordia formation

Hydrophobins are small, secreted proteins with important physiological functions in mycelial growth and fungal development. Here, 1 nucleus-specific and 35 allelic hydrophobin genes were identified in the genome of a white rot fungus, Coriolopsis trogii. Among these, 22 were eight-cysteine class I hydrophobin genes and the other 14 were uncommon six-cysteine hydrophobin genes. The six-cysteine hydrophobins were speculated to have originated from a common ancestor. The hydrophobin genes favored a clustering distribution and two recent duplication pairs were identified. The genes had conserved gene structures with three exons and two introns. Cthyd18, Cthyd19, and Cthyd32 were constitutively highly expressed in all developmental stages. Cthyd20, Cthyd21, Cthyd22, Cthyd28, Cthyd30, Cthyd31, and Cthyd33 were highly expressed in mycelia, and Cthyd12 and Cthyd35 in the reproductive stages. Sixteen hydrophobin genes were regulated differently in the transition from mycelia to primordia; Cthyd35 showed maximal upregulation of 1922-fold, and Cthyd23 showed maximal downregulation of 552-fold. Most (32) hydrophobin genes showed significant differential expression between mycelia cultured in different media (potato dextrose agar or broth). Weighted gene co-expression network analysis and promoter analysis revealed that C2H2 zinc finger proteins may regulate hydrophobin genes. These results may support further research into the function and evolution of hydrophobins.