Full Transcriptome Analysis of Callus Suspension Culture System of Bletilla striata

Cytochemical localization and synthesis mechanism of the glucomannan in pseudobulbs of Bletilla striata Reichb. f

The dried pseudobulbs of Bletilla striata, an important traditional Chinese medicine named BaiJi, have an extraordinary polysaccharide content and excellent prospects for medicinal effects. However, the distribution and molecular mechanism underlying biosynthesis are poorly understood. In this study, chemical and immunologic analyses were performed in representative tissues of B. striata, and the results showed that what are conventionally termed Bletilla striata polysaccharides (BSPs) are water-soluble polysaccharides deposited only in pseudobulbs. The structural component of BSPs is glucomannan, with a mannose:glucose mass ratio of ~3:2. BSPs are present in the parenchyma of the pseudobulbs in cells known as glucomannan idioblasts and distributed in the cytoplasm within cellular membranes, but are not contained in the vacuole. Comparative transcriptomics and bioinformatics analyses mapped the pathway from sucrose to BSP and identified BsGPI, BsmanA, and BsCSLAs as the key genes of BSP biosynthesis, suggesting that the functional differentiation of the cellulose synthase-like family A (CSLA) may be critical for the flow of glucomannan to the BSP or cell wall. Subsequently, virus-mediated gene silencing showed that silencing of two CSLAs (Bs03G11846 and Bs03G11849) led to a decrease in BSP content, and yeast two-hybrid and luciferase complementation experiments confirmed that four CSLAs (Bs03G11846, Bs03G11847, Bs03G11848, and Bs03G11849) can form homo- or heterodimers, suggesting that multiple CSLAs may form a large complex that functions in BSP synthesis. Our results provide cytological evidence of BSP and describe the isolation and characterization of candidate genes involved in BSP synthesis, laying a solid foundation for further research on its regulation mechanisms and the genetic engineering breeding of B. striata.

Transcriptional Changes in Damask Rose Suspension Cell Culture Revealed by RNA Sequencing

Damask roses (Rosa x damascena) are widely used in cosmetics and pharmaceutics. Here, we established an in vitro suspension cell culture for calli derived from damask rose petals. We analyzed rose suspension cell transcriptomes obtained at two different time points by RNA sequencing to reveal transcriptional changes during rose suspension cell culture. Of the 580 coding RNAs (1.3%) highly expressed in the suspension rose cells, 68 encoded cell wall-associated proteins. However, most RNAs encoded by the chloroplasts and mitochondria are not expressed. Many highly expressed coding RNAs are involved in translation, catalyzing peptide synthesis in ribosomes. Moreover, the amide metabolic process producing naturally occurring alkaloids was the most abundant metabolic process during the propagation of rose suspension cells. During rose cell propagation, coding RNAs involved in the stress response were upregulated at an early stage, while coding RNAs associated with detoxification and transmembrane transport were upregulated at the late stage. We used transcriptome analyses to reveal important biological processes and molecular mechanisms during rose suspension cell culture. Most non-coding (nc) RNAs were not expressed in the rose suspension cells, but a few ncRNAs with unknown functions were highly expressed. The expression of ncRNAs and their target coding RNAs was highly correlated. Taken together, we revealed significant biological processes and molecular mechanisms occurring during rose suspension cell culture using transcriptome analyses.

De novo full length transcriptome analysis and gene expression profiling to identify genes involved in phenylethanol glycosides biosynthesis in Cistanche tubulosa

Background

The dried stem of Cistanche, is a famous Chinese traditional medicine. The main active pharmacodynamic components are phenylethanol glycosides (PhGs). Cistanche tubulosa produces higher level of PhGs in its stems than that of Cistanche deserticola. However, the key genes in the PhGs biosynthesis pathway is not clear in C. tubulosa.

Results

In this study, we performed the full-length transcriptome sequencing and gene expression profiling of C. tubulosa using PacBio combined with BGISEQ-500 RNA-seq technology. Totally, 237,772 unique transcripts were obtained, ranging from 199 bp to 31,857 bp. Among the unique transcripts, 188,135 (79.12%) transcripts were annotated. Interestingly, 1080 transcripts were annotated as 22 enzymes related to PhGs biosynthesis. We measured the content of echinacoside, acteoside and total PhGs at two development stages, and found that the content of PhGs was 46.74% of dry matter in young fleshy stem (YS1) and then decreased to 31.22% at the harvest stage (HS2). To compare with YS1, 13,631 genes were up-regulated, and 15,521 genes were down regulated in HS2. Many differentially expressed genes (DEGs) were identified to be involved in phenylpropanoid biosynthesis pathway, phenylalanine metabolism pathway, and tyrosine metabolism pathway.

Conclusions

This is the first report of transcriptome study of C. tubulosa which provided the foundation for understanding of PhGs biosynthesis. Based on these results, we proposed a potential model for PhGs biosynthesis in C. tubulosa.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08921-x.

Haplotype-resolved genome assembly of Bletilla striata (Thunb.) Reichb.f. to elucidate medicinal value

Bletilla striata, commonly known as baiji, is a species used in traditional Chinese medicine; it is highly regarded for its medicinal applications and therefore has high economic value. Here, we report a high-quality haplotype-resolved genome of B. striata, haplotype A (2.37 Gb, with a scaffold N50 of 146.39 Mb and a contig N50 of 1.65 Mb) and haplotype B (2.43 Gb, with a scaffold N50 of 150.22 Mb and a contig N50 of 1.66 Mb), assembled from high-fidelity (HiFi) reads and chromosome conformation capture (Hi-C) reads. We find evidence that B. striata has undergone two whole-genome duplication (WGD) events: an ancient WGD event shared by most monocots and a recent WGD event unique to all orchids. We also reconstructed the ancestral orchid karyotype (AOK) of 18 ancient chromosomes and the evolutionary trajectories of 16 modern B. striata chromosomes. Comparative genomic analysis suggests that the expanded gene families of B. striata might play important roles in secondary metabolite biosynthesis and environmental adaptation. By combining genomic and transcriptomic data, we identified the 10 core members from nine gene families that were probably involved in B. striata polysaccharide (BSP) biosynthesis. Based on virus-induced gene silencing (VIGS) and yeast two-hybrid experiments, we present an MYB transcription factor (TF), BsMYB2, that can regulate BSP biosynthesis by directly interacting with eight key BSP-related genes: sacA1, HK1, scrK1, scrK2, GPI1, manA1, GMPP1 and UGP2_1. Our study will enhance the understanding of orchid evolution and accelerate the molecular-assisted breeding of B. striata for improving traits of medicinal value.

Genome-Wide Identification of Genes Related to Biosynthesis of Phenolic Acid Derivatives in Bletilla striata at Different Suspension Culture Stages

To screen the genes regulating the biosynthesis of phenolic acid derivatives from the genome of Bletilla striata, we designed a suspension culture system to sample the cells for the following experiments. The contents of four phenolic acid derivatives were determined by high-performance liquid chromatography, and several full-length transcriptome sequencings of RNA samples at 10 time points were performed for bioinformatics analysis. The correlation analysis was used to identify and verify the key DEGs involved in the biosynthesis of the four phenolic acid derivatives. The results showed that the contents of p-hydroxybenzylalcohol (HBA), Dactylorhin A, Militarine, and Coelonin peaked at 33 days postinoculation (Dpi), 18 Dpi, 39 Dpi, and 39 Dpi of the culture system, respectively. Based on transcriptome data, 80 DEGs involved in the biosynthesis of phenolic acid derivatives were obtained. The KEGG pathway enrichment analysis classified them mostly into five metabolic pathways: phenylpropane biosynthesis, starch and sucrose metabolic, cyanoamino acid metabolism, gluconeogenesis and glycolysis, and phenylalanine metabolism. qPCR analysis revealed that the relative gene expression levels were consistent with the overall trend of transcriptome sequencing results. Among them, 14, 18, 23, and 41 unigenes were found to be involved in the synthesis of HBA, Dactylorhin A, Coelonin, and Militarine, respectively. These unigenes laid a solid foundation for elucidating the biosynthesis mechanism of phenolic acid derivatives in suspension cells of B. striata.

An organ-specific transcriptomic atlas of the medicinal plant Bletilla striata: Protein-coding genes, microRNAs, and regulatory networks

As one of the important species belonging to the Bletilla genus of Orchidaceae, Bletilla striata (Thunb.) Rchb. f., possesses both ornamental and medicinal values. Its dried tubers are used as a traditional Chinese medicine, and several secondary metabolites have been indicated to be the active ingredients. However, the molecular mechanisms related to the regulation of secondary metabolism have not been characterized in B. striata. In this study, integrated analysis of RNA sequencing (RNA-seq), small RNA sequencing (sRNA-seq), and degradome sequencing (degradome-seq) data from three organs (leaf, root, and tuber) of B. striata provided us with a comprehensive view of the microRNA (miRNA)-mediated regulatory network. Firstly, based on the RNA-seq data, the organ-specific expression patterns of the protein-coding genes, especially for those related to secondary metabolism, were investigated. Secondly, 342 conserved miRNA candidates were identified from B. striata. These miRNAs were assigned to 88 families, some of which were selected for expression pattern analysis. Additionally, 31 hairpin-structured precursors encoding 23 novel miRNAs were uncovered from the transcriptome assembly. Thirdly, based on the degradome signatures, 1,142 validated miRNA-target pairs (involving 167 conserved miRNAs and six novel miRNAs and 51 target genes) were included in the regulatory network. Organ-specific expression level comparison between the miRNAs and their targets revealed some interesting miRNA-target pairs. Fourthly, some valuable subnetworks were extracted for further functional studies. Additionally, some regulatory pathways were indicated to be monocot specific. Summarily, our results lay a solid basis for in-depth studies on the regulatory mechanisms underlying the production of the medicinal ingredients in B. striata.

Phenotype correlation analysis and excellent germplasm screening of herb Bletilla Rchb. f. based on comprehensive evaluation from thirty-three geographic populations

BACKGROUND

The Bletilla genus of Orchidaceae includes plants with great economic value, among which B. striata is the main traditional medicinal plant, and its pseudobulb, known as BaiJi, was first recorded in Shennong's Classic of Materia Medica. However, there has been little systemic evaluation of the germplasm quality of Bletilla plants in China. In order to comprehensive evaluate the Bletilla resources in China and screen out the candidate phenotypic traits determining yield and/or quality of Bletilla, the variation of phenotypic indicators (pseudobulb, leaf, stem, inflorescence, flower) and active ingredients contents (polysaccharide, total phenolics and militarine) in different populations of B. striata and B. ochracea were investigated through 4 years' common-garden experiment.

RESULTS

There were abundant phenotypic variations and significant differences among different populations in the morphological phenotypes, pseudobulb weight and main active ingredient contents. AHBZ, HBLT and HBSN populations showed good prospects for industrial development, presenting higher quality in terms of yield and main active ingredient content. Pseudobulb yield, polysaccharide and total phenol content are positively correlated with phenotypic traits. Militarine content is negatively correlated with almost all indexes. Plant height, leaf width and stem diameter may be important indicators of potential excellent germplasms.

CONCLUSIONS

Bletilla is not strictly geoauthentic medicinal plants. B. ochracea could be accepted as an alternative resource to B. striata. The best harvest period of Bletilla is the third year after cultivation. Plant height, leaf width and stem diameter may be important indicators of potential excellent germplasms. These results provide important information required for the efficient screening and utilization of Bletilla germplasm resources.

Full-length transcriptome sequencing reveals the molecular mechanism of potato seedlings responding to low-temperature

BACKGROUND

Potato (Solanum tuberosum L.) is one of the world's most important crops, the cultivated potato is frost-sensitive, and low-temperature severely influences potato production. However, the mechanism by which potato responds to low-temperature stress is unclear. In this research, we apply a combination of second-generation sequencing and third-generation sequencing technologies to sequence full-length transcriptomes in low-temperature-sensitive cultivars to identify the important genes and main pathways related to low-temperature resistance.

RESULTS

In this study, we obtained 41,016 high-quality transcripts, which included 15,189 putative new transcripts. Amongst them, we identified 11,665 open reading frames, 6085 simple sequence repeats out of the potato dataset. We used public available genomic contigs to analyze the gene features, simple sequence repeat, and alternative splicing event of 24,658 non-redundant transcript sequences, predicted the coding sequence and identified the alternative polyadenylation. We performed cluster analysis, GO, and KEGG functional analysis of 4518 genes that were differentially expressed between the different low-temperature treatments. We examined 36 transcription factor families and identified 542 transcription factors in the differentially expressed genes, and 64 transcription factors were found in the AP2 transcription factor family which was the most. We measured the malondialdehyde, soluble sugar, and proline contents and the expression genes changed associated with low temperature resistance in the low-temperature treated leaves. We also tentatively speculate that StLPIN10369.5 and StCDPK16 may play a central coordinating role in the response of potatoes to low temperature stress.

CONCLUSIONS

Overall, this study provided the first large-scale full-length transcriptome sequencing of potato and will facilitate structure-function genetic and comparative genomics studies of this important crop.

De novo full length transcriptome analysis of Arachis glabrata provides insights into gene expression dynamics in response to biotic and abiotic stresses

The perennial ornamental peanut Arachis glabrata represents one of the most adaptable wild Arachis species. This study used PacBio combined with BGISEQ-500 RNA-seq technology to study the transcriptome and gene expression dynamics of A. glabrata. Of the total 109,747 unique transcripts obtained, >90,566 transcripts showed significant homology to known proteins and contained the complete coding sequence (CDS). RNA-seq revealed that 1229, 1039, 1671, 3923, 1521 and 1799 transcripts expressed specifically in the root, stem, leaf, flower, peg and pod, respectively. We also identified thousands of differentially expressed transcripts in response to drought, salt, cold and leaf spot disease. Furthermore, we identified 30 polyphenol oxidase encoding genes associated with the quality of forage, making A. glabrata suitable as a forage crop. Our findings presented the first transcriptome study of A. glabrata which will facilitate genetic and genomics studies and lays the groundwork for a deeper understanding of the A. glabrata genome.