Building an octaploid genome and transcriptome of the medicinal plant Pogostemon cablin from Lamiales

Allopolyploid origin and diversification of the Hawaiian endemic mints

Island systems provide important contexts for studying processes underlying lineage migration, species diversification, and organismal extinction. The Hawaiian endemic mints (Lamiaceae family) are the second largest plant radiation on the isolated Hawaiian Islands. We generated a chromosome-scale reference genome for one Hawaiian species, Stenogyne calaminthoides, and resequenced 45 relatives, representing 34 species, to uncover the continental origins of this group and their subsequent diversification. We further resequenced 109 individuals of two Stenogyne species, and their purported hybrids, found high on the Mauna Kea volcano on the island of Hawai'i. The three distinct Hawaiian genera, Haplostachys, Phyllostegia, and Stenogyne, are nested inside a fourth genus, Stachys. We uncovered four independent polyploidy events within Stachys, including one allopolyploidy event underlying the Hawaiian mints and their direct western North American ancestors. While the Hawaiian taxa may have principally diversified by parapatry and drift in small and fragmented populations, localized admixture may have played an important role early in lineage diversification. Our genomic analyses provide a view into how organisms may have radiated on isolated island chains, settings that provided one of the principal natural laboratories for Darwin's thinking about the evolutionary process.

Polyploidization of Indotyphlops braminus: evidence from isoform-sequencing

BACKGROUND

Indotyphlops braminus, the only known triploid parthenogenetic snake, is a compelling species for revealing the mechanism of polyploid emergence in vertebrates.

METHODS

In this study, we applied PacBio isoform sequencing technology to generate the first full-length transcriptome of I. braminus, aiming to improve the understanding of the molecular characteristics of this species.

RESULTS

A total of 51,849 nonredundant full-length transcript assemblies (with an N50 length of 2980 bp) from I. braminus were generated and fully annotated using various gene function databases. Our analysis provides preliminary evidence supporting a recent genome duplication event in I. braminus. Phylogenetic analysis indicated that the divergence of I. braminus subgenomes occurred approximately 11.5 ~ 15 million years ago (Mya). The full-length transcript resource generated as part of this research will facilitate transcriptome analysis and genomic evolution studies in the future.

A Genome-Wide Identification and Expression Analysis of the Casparian Strip Membrane Domain Protein-like Gene Family in Pogostemon cablin in Response to p-HBA-Induced Continuous Cropping Obstacles

Casparian strip membrane domain protein-like (CASPL) genes are key genes for the formation and regulation of the Casparian strip and play an important role in plant abiotic stress. However, little research has focused on the members, characteristics, and biological functions of the patchouli PatCASPL gene family. In this study, 156 PatCASPL genes were identified at the whole-genome level. Subcellular localization predicted that 75.6% of PatCASPL proteins reside on the cell membrane. A phylogenetic analysis categorized PatCASPL genes into five subclusters alongside Arabidopsis CASPL genes. In a cis-acting element analysis, a total of 16 different cis-elements were identified, among which the photo-responsive element was the most common in the CASPL gene family. A transcriptome analysis showed that p-hydroxybenzoic acid, an allelopathic autotoxic substance, affected the expression pattern of PatCASPLs, including a total of 27 upregulated genes and 30 down-regulated genes, suggesting that these PatCASPLs may play an important role in the regulation of patchouli continuous cropping obstacles by affecting the formation and integrity of Casparian strip bands. These results provided a theoretical basis for exploring and verifying the function of the patchouli PatCASPL gene family and its role in continuous cropping obstacles.

An integrated transcriptome and metabolome analysis reveals the gene network regulating flower development in Pogostemon cablin

Pogostemon cablin is a well-known protected species widely used in medicine and spices, however the underlying molecular mechanisms and metabolite dynamics of P. cablin flower development remain unclear due to the difficulty in achieving flowering in this species. A comparison of the transcriptome and widely targeted metabolome during P. cablin flower development was first performed in this study. Results showed that a total of 13,469 differentially expressed unigenes (DEGs) and 371 differentially accumulated metabolites (DAMs) were identified. Transcriptomic analysis revealed that the DEGs were associated with starch and sucrose metabolism, terpenoid biosynthesis and phenylpropanoid biosynthesis. Among these DEGs, 75 MIKC-MADS unigenes were associated with the development of floral organs. Gibberellins (GAs), auxin, and aging signaling might form a cross-regulatory network to regulate flower development in P. cablin. According to the metabolic profile, the predominant DAMs were amino acids, flavonoids, terpenes, phenols, and their derivatives. The accumulation patterns of these predominant DAMs were closely associated with the flower developmental stage. The integration analysis of DEGs and DAMs indicated that phenylpropanoids, flavonoids, and amino acids might be accumulated due to the activation of starch and sucrose metabolism. Our results provide some important insights for elucidating the reproductive process, floral organ, and color formation of P. cablin flowers at the molecular level. These results will improve our understanding of the molecular and genetic mechanisms involved in the floral development of P. cablin.

Uncovering a miltiradiene biosynthetic gene cluster in the Lamiaceae reveals a dynamic evolutionary trajectory

The spatial organization of genes within plant genomes can drive evolution of specialized metabolic pathways. Terpenoids are important specialized metabolites in plants with diverse adaptive functions that enable environmental interactions. Here, we report the genome assemblies of Prunella vulgaris, Plectranthus barbatus, and Leonotis leonurus. We investigate the origin and subsequent evolution of a diterpenoid biosynthetic gene cluster (BGC) together with other seven species within the Lamiaceae (mint) family. Based on core genes found in the BGCs of all species examined across the Lamiaceae, we predict a simplified version of this cluster evolved in an early Lamiaceae ancestor. The current composition of the extant BGCs highlights the dynamic nature of its evolution. We elucidate the terpene backbones generated by the Callicarpa americana BGC enzymes, including miltiradiene and the terpene (+)-kaurene, and show oxidization activities of BGC cytochrome P450s. Our work reveals the fluid nature of BGC assembly and the importance of genome structure in contributing to the origin of metabolites.

The Current Developments in Medicinal Plant Genomics Enabled the Diversification of Secondary Metabolites' Biosynthesis

Medicinal plants produce important substrates for their adaptation and defenses against environmental factors and, at the same time, are used for traditional medicine and industrial additives. Plants have relatively little in the way of secondary metabolites via biosynthesis. Recently, the whole-genome sequencing of medicinal plants and the identification of secondary metabolite production were revolutionized by the rapid development and cheap cost of sequencing technology. Advances in functional genomics, such as transcriptomics, proteomics, and metabolomics, pave the way for discoveries in secondary metabolites and related key genes. The multi-omics approaches can offer tremendous insight into the variety, distribution, and development of biosynthetic gene clusters (BGCs). Although many reviews have reported on the plant and medicinal plant genome, chemistry, and pharmacology, there is no review giving a comprehensive report about the medicinal plant genome and multi-omics approaches to study the biosynthesis pathway of secondary metabolites. Here, we introduce the medicinal plant genome and the application of multi-omics tools for identifying genes related to the biosynthesis pathway of secondary metabolites. Moreover, we explore comparative genomics and polyploidy for gene family analysis in medicinal plants. This study promotes medicinal plant genomics, which contributes to the biosynthesis and screening of plant substrates and plant-based drugs and prompts the research efficiency of traditional medicine.

Integrated Analysis of Physiological, mRNA Sequencing, and miRNA Sequencing Data Reveals a Specific Mechanism for the Response to Continuous Cropping Obstacles in Pogostemon cablin Roots

Pogostemon cablin (patchouli) is a commercially important medicinal and industrial crop grown worldwide for its medicinal and aromatic properties. Patchoulol and pogostone, derived from the essential oil of patchouli, are considered valuable components in the cosmetic and pharmaceutical industries. Due to its high application value in the clinic and industry, the demand for patchouli is constantly growing. Unfortunately, patchouli cultivation has suffered due to severe continuous cropping obstacles, resulting in a significant decline in yield and quality. Moreover, the physiological and transcriptional changes in patchouli in response to continuous cropping obstacles remain unclear. This has greatly restricted the development of the patchouli industry. To explore the mechanism underlying the rapid response of patchouli roots to continuous cropping stress, integrated analysis of the transcriptome and miRNA profiles of patchouli roots under continuous and noncontinuous cropping conditions in different growth periods was conducted using RNA sequencing (RNA-seq) and miRNA-seq and complemented with physiological data. The physiological and biochemical results showed that continuous cropping significantly inhibited root growth, decreased root activity, and increased the activity of antioxidant enzymes (superoxide dismutase, peroxidase, and catalase) and the levels of osmoregulators (malondialdehyde, soluble protein, soluble sugar, and proline). Subsequently, we found 4,238, 3,494, and 7,290 upregulated and 4,176, 3,202, and 8,599 downregulated differentially expressed genes (DEGs) in the three growth periods of continuously cropped patchouli, many of which were associated with primary carbon and nitrogen metabolism, defense responses, secondary metabolite biosynthesis, and transcription factors. Based on miRNA-seq, 927 known miRNAs and 130 novel miRNAs were identified, among which 67 differentially expressed miRNAs (DEMIs) belonging to 24 miRNA families were induced or repressed by continuous cropping. By combining transcriptome and miRNA profiling, we obtained 47 miRNA-target gene pairs, consisting of 18 DEMIs and 43 DEGs, that likely play important roles in the continuous cropping response of patchouli. The information provided in this study will contribute to clarifying the intricate mechanism underlying the patchouli response to continuous cropping obstacles. In addition, the candidate miRNAs and genes can provide a new strategy for breeding continuous cropping-tolerant patchouli.

Transcriptome analysis of two Pogostemon cablin chemotypes reveals genes related to patchouli alcohol biosynthesis

Pogostemon cablin, a medicinally and economically important perennial herb, is cultivated around the world due to its medicinal and aromatic properties. Different P. cablin cultivars exhibit different morphological traits and patchouli oil components and contents (especially patchouli alcohol (PA) and pogostone (PO)). According to the signature constituent of the leaf, P. cablin was classified into two different chemotypes, including PA-type and PO-type. To better understand the molecular mechanisms of PA biosynthesis, the transcriptomes of Chinese-cultivated P. cablin cv. PA-type “Nanxiang” (NX) and PO-type “Paixiang” (PX) were analyzed and compared with ribonucleic acid sequencing (RNA-Seq) technology. We obtained a total of 36.83 G clean bases from the two chemotypes, compared them with seven databases and revealed 45,394 annotated unigenes. Thirty-six candidate unigenes participating in the biosynthesis of PA were found in the P. cablin transcriptomes. Overall, 8,390 differentially expressed unigenes were identified between the chemotypes, including 2,467 upregulated and 5,923 downregulated unigenes. Furthermore, six and nine differentially expressed genes (DEGs) were mapped to the terpenoid backbone biosynthetic and sesquiterpenoid and triterpenoid biosynthetic pathways, respectively. One key sesquiterpene synthase gene involved in the sesquiterpenoid and triterpenoid biosynthetic pathways, encoding patchoulol synthase variant 1, was significantly upregulated in NX. Additionally, GC-MS analysis of the two chemotypes in this study showed that the content of PA in NX was significantly higher than that of PX, while the content of PO showed the opposite phenotype. Quantitative real-time polymerase chain reaction (qRT-PCR) analysis showed that the DEG expression tendency was consistent with the transcriptome sequencing results. Overall, 23 AP2/ERF, 13 bHLH, 11 MYB, 11 NAC, three Trihelix, 10 WRKY and three bZIP genes that were differentially expressed may act as regulators of terpenoid biosynthesis. Altogether, 8,314 SSRs were recognized within 6,825 unigenes, with a distribution frequency of 18.32%, among which 1,202 unigenes contained more than one SSR. The transcriptomic characteristics of the two P. cablin chemotypes are comprehensively reported in this study, and these results will contribute to a better understanding of the molecular mechanism of PA biosynthesis. Our transcriptome data also provide a valuable genetic resource for further studies on P. cablin.

Genome-wide analyses of tandem repeats and transposable elements in patchouli

Patchouli, Pogostemon cablin (Blanco) Benth., is a traditional Chinese medicinal plant from the order Lamiales. It is considered a valuable herb due to its essential oil content and range of therapeutic effects. This study aimed to explore the evolutionary history of repetitive sequences in the patchouli genome by analyzing tandem repeats and transposable elements (TEs). We first retrieved genomic data for patchouli and four other Lamiales species from the GenBank database. Next, the content of tandem repeats with different period sizes was identified. Long terminal repeats (LTRs) were then identified with LTR_STRUC. Finally, the evolutionary landscape of TEs was explored using an in-house PERL program. The analysis of repetitive sequences revealed that tandem repeats constitute a higher proportion of the patchouli genome compared to the four other species. Analyses of TE families showed that most of the repetitive sequences in the patchouli genome are TEs, and that recently inserted TEs make up a comparatively larger proportion than older ones. Our analyses of LTR retrotransposons in their host genome indicated the existence of ancient LTR retrotransposon expansion, and the escape of these elements from natural selection revealed their ages. Our identification and analyses of repetitive sequences should provide new insights for further investigation of patchouli evolution.

Review on the Development and Applications of Medicinal Plant Genomes

With the development of sequencing technology, the research on medicinal plants is no longer limited to the aspects of chemistry, pharmacology, and pharmacodynamics, but reveals them from the genetic level. As the price of next-generation sequencing technology becomes affordable, and the long-read sequencing technology is established, the medicinal plant genomes with large sizes have been sequenced and assembled more easily. Although the review of plant genomes has been reported several times, there is no review giving a systematic and comprehensive introduction about the development and application of medicinal plant genomes that have been reported until now. Here, we provide a historical perspective on the current situation of genomes in medicinal plant biology, highlight the use of the rapidly developing sequencing technologies, and conduct a comprehensive summary on how the genomes apply to solve the practical problems in medicinal plants, like genomics-assisted herb breeding, evolution history revelation, herbal synthetic biology study, and geoherbal research, which are important for effective utilization, rational use and sustainable protection of medicinal plants.

Gene and genome duplications in the evolution of chemodiversity: perspectives from studies of Lamiaceae

Plants are reservoirs of extreme chemical diversity, yet biosynthetic pathways remain underexplored in the majority of taxa. Access to improved, inexpensive genomic and computational technologies has recently enhanced our understanding of plant specialized metabolism at the biochemical and evolutionary levels including the elucidation of pathways leading to key metabolites. Furthermore, these approaches have provided insights into the mechanisms of chemical evolution, including neofunctionalization and subfunctionalization, structural variation, and modulation of gene expression. The broader utilization of genomic tools across the plant tree of life, and an expansion of genomic resources from multiple accessions within species or populations, will improve our overall understanding of chemodiversity. These data and knowledge will also lead to greater insight into the selective pressures contributing to and maintaining this diversity, which in turn will enable the development of more accurate predictive models of specialized metabolism in plants.

Global analysis of lysine succinylation in patchouli plant leaves

Lysine succinylation is a novel, naturally occurring posttranslational modification (PTM) in living organisms. Global lysine succinylation identification has been performed at the proteomic level in various species; however, the study of lysine succinylation in plant species is relatively limited. Patchouli plant (P. cablin (Blanco) Benth., Lamiaceae) is a globally important industrial plant and medicinal herb. In the present study, lysine succinylome analysis was carried out in patchouli plants to determine the potential regulatory role of lysine succinylation in patchouli growth, development, and physiology. The global succinylation sites and proteins in patchouli plants were screened with an immunoprecipitation affinity enrichment technique and advanced mass spectrometry-based proteomics. Several bioinformatic analyses, such as function classification and enrichment, subcellular location predication, metabolic pathway enrichment and protein-protein interaction networking, were conducted to characterize the functions of the identified sites and proteins. In total, 1097 succinylation sites in 493 proteins were detected in patchouli plants, among which 466 succinylation sites in 241 proteins were repeatedly identified within three independent experiments. The functional characterization of these proteins indicated that the tricarboxylic acid (TCA) cycle, oxidative phosphorylation, photosynthesis processes, and amino acid biosynthesis may be regulated by lysine succinylation. In addition, these succinylated proteins showed a wide subcellular location distribution, although the chloroplast and cytoplasm were the top two preferred cellular components. Our study suggested the important role of lysine succinylation in patchouli plant physiology and biology and could serve as a useful reference for succinylation studies in other medicinal plants.