Transcriptome profiling reveals the crucial biological pathways involved in cold response in Moso bamboo (Phyllostachys edulis)

CaAOS as a hub gene based on physiological and transcriptomic analyses of cold-resistant and cold-sensitive pepper cultivars

The yield and quality of pepper are considerably influenced by the cold conditions. Herein, we performed morphological, physiological and transcriptomic analyses by using two pepper seedlings, '2379' (cold-resistant) and '2380' (cold-sensitive). Briefly, 60 samples from each cultivar were analyzed at four distinct time points (0, 6, 24 and 48 h) at 5 °C in darkness. The physiological indices and activities of enzymes exhibited marked differences between the two cultivars. Transcriptomic analysis indicated that, compared to the control group, 11,415 DEGs were identified in '2379' and '2380' at 24 h. In the early stage, the number of DEGs in '2379' was 5.68 times higher than that in '2380', potentially explaining the observed differences in tolerance to colds. Processes such as protein targeting to membranes, jasmonic acid (JA)-mediated signalling, cold response and abscisic acid-activated signalling were involved. Subsequently, we identified a hub gene, CaAOS, that is involved in JA biosynthesis, positively influences cold tolerance and is a target of CaMYC2. Variations in the GC-motif of the CaAOS's promoter may influence the expression levels of CaAOS under cold treatment. The result of this study may lead to the development of more effective strategies for enhancing cold tolerance, potentially benefitting pepper breeding in cold regions.

Transcriptome profiling of Bergenia purpurascens under cold stress

Bergenia purpurascens is an important medicinal, edible and ornamental plant. It generally grows in high-altitude areas with complex climates. There have been no reports about how B. purpurascens survives under cold stress. Here, the B. purpurascens under low temperature were subjected to transcriptomics analysis to explore the candidate genes and pathways that involved in the cold tolerance of B. purpurascens. Compared with the control treatment, we found 9,600 up-regulated differentially expressed genes (DEGs) and 7,055 down-regulated DEGs. A significant number of DEGs were involved in the Ca2+ signaling pathway, mitogen-activated protein kinase (MAPK) cascade, plant hormone signaling pathway, and lipid metabolism. A total of 400 transcription factors were found to respond to cold stress, most of which belonged to the MYB and AP2/ERF families. Five novel genes were found to be potential candidate genes involved in the cold tolerance of B. purpurascens. The study provide insights into further investigation of the molecular mechanism of how B. purpurascens survives under cold stress.

Comprehensive genome-wide analysis of the DREB gene family in Moso bamboo (Phyllostachys edulis): evidence for the role of PeDREB28 in plant abiotic stress response

Dehydration response element binding (DREB) proteins are vital for plant abiotic stress responses, but the understanding of DREBs in bamboo, an important sustainable non-timber forest product, is limited. Here we conducted a comprehensive genome-wide analysis of the DREB gene family in Moso bamboo, representing the most important running bamboo species in Asia. In total, 44 PeDREBs were identified, and information on their gene structures, protein motifs, phylogenetic relationships, and stress-related cis-regulatory elements (CREs) was provided. Based on the bioinformatical analysis, we further analyzed PeDREBs from the A5 group and found that four of five PeDREB transcripts were induced by salt, drought, and cold stresses, and their proteins could bind to stress-related CREs. Among these, PeDREB28 was selected as a promising candidate for further functional characterization. PeDREB28 is localized in nucleus, has transcriptional activation activity, and could bind to the DRE- and coupling element 1- (CE1) CREs. Overexpression of PeDREB28 in Arabidopsis and bamboo improved plant abiotic stress tolerance. Transcriptomic analysis showed that broad changes due to the overexpression of PeDREB28. Furthermore, 628 genes that may act as the direct PeDREB28 downstream genes were identified by combining DAP-seq and RNA-seq analysis. Moreover, we confirmed that PeDREB28 could bind to the promoter of pyrabactin-resistance-like gene (DlaPYL3), which is a homolog of abscisic acid receptor in Arabidopsis, and activates its expression. In summary, our study provides important insights into the DREB gene family in Moso bamboo, and contributes to their functional verification and genetic engineering applications in the future.

Orthogroup and phylotranscriptomic analyses identify transcription factors involved in the plant cold response: A case study of Arabidopsis BBX29

C-repeat binding factors (CBFs) are well-known transcription factors (TFs) that regulate plant cold acclimation. RNA sequencing (RNA-seq) data from diverse plant species provide opportunities to identify other TFs involved in the cold response. However, this task is challenging because gene gain and loss has led to an intertwined community of co-orthologs and in-paralogs between and within species. Using orthogroup (closely related homologs) analysis, we identified 10,549 orthogroups in five representative eudicots. A phylotranscriptomic analysis of cold-treated seedlings from eudicots identified 35 high-confidence conserved cold-responsive transcription factor orthogroups (CoCoFos). These 35 CoCoFos included the well-known cold-responsive regulators CBFs, HSFC1, ZAT6/10, and CZF1 among others. We used Arabidopsis BBX29 for experimental validation. Expression and genetic analyses showed that cold-induction of BBX29 is CBF- and abscisic acid-independent, and BBX29 is a negative regulator of cold tolerance. Integrative RNA-seq and Cleavage Under Targets and Tagmentation followed by sequencing analyses revealed that BBX29 represses a set of cold-induced TFs (ZAT12, PRR9, RVE1, MYB96, etc.). Altogether, our analysis yielded a library of eudicot CoCoFos and demonstrated that BBX29 is a negative regulator of cold tolerance in Arabidopsis.

A phylotranscriptomic dataset of angiosperm species under cold stress

Angiosperms are one of the most diverse and abundant plant groups that are widely distributed on Earth, from tropical to temperate and polar zones. The wide distribution of angiosperms may be attributed to the evolution of sophisticated mechanisms of environmental adaptability, including cold tolerance. Since the development of high-throughput sequencing, transcriptome has been widely utilized to gain insights into the molecular mechanisms of plants in response to cold stress. However, previous studies generally focused on single or two species, and comparative transcriptome analyses for multispecies responding to cold stress were limited. In this study, we selected 11 representative angiosperm species, performed phylotranscriptome experiments at four time points before and after cold stress, and presented a profile of cold-induced transcriptome changes in angiosperms. Our multispecies cold-responsive RNA-seq datasets provide valuable references for exploring conserved and evolutionary mechanisms of angiosperms in adaptation to cold stress.

Data-Independent Acquisition Proteomics Reveals the Effects of Red and Blue Light on the Growth and Development of Moso Bamboo (Phyllostachys edulis) Seedlings

Moso bamboo is a rapidly growing species with significant economic, social, and cultural value. Transplanting moso bamboo container seedlings for afforestation has become a cost-effective method. The growth and development of the seedlings is greatly affected by the quality of light, including light morphogenesis, photosynthesis, and secondary metabolite production. Therefore, studies on the effects of specific light wavelengths on the physiology and proteome of moso bamboo seedlings are crucial. In this study, moso bamboo seedlings were germinated in darkness and then exposed to blue and red light conditions for 14 days. The effects of these light treatments on seedling growth and development were observed and compared through proteomics analysis. Results showed that moso bamboo has higher chlorophyll content and photosynthetic efficiency under blue light, while it displays longer internode and root length, more dry weight, and higher cellulose content under red light. Proteomics analysis reveals that these changes under red light are likely caused by the increased content of cellulase CSEA, specifically expressed cell wall synthetic proteins, and up-regulated auxin transporter ABCB19 in red light. Additionally, blue light is found to promote the expression of proteins constituting photosystem II, such as PsbP and PsbQ, more than red light. These findings provide new insights into the growth and development of moso bamboo seedlings regulated by different light qualities.

Identification and expression profiles of xylogen-like arabinogalactan protein (XYLP) gene family in Phyllostachys edulis in different developmental tissues and under various abiotic stresses

Xylogen-like arabinogalactan protein (XYLP) is an atypical lipid transport protein. In this study, 23 Phyllostachys edulis XYLPs were identified, and their proteins contain characteristic structures of AGP and nsLTP domain. All PeXYLPs can be divided into four clades, and their genes were unevenly distributed on 11 chromosome scaffolds. Collinear analysis revealed that segmental duplication was the main driver for PeXYLP family expansion. The cis-acting elements presented in the promoter are involved in various regulations of PeXYLPs expression. G.O. annotation revealed that PeXYLPs are mainly interested in lipid transport and synthesis and primarily function at the plasma membrane. Transcriptome analysis revealed that PeXYLPs were spatiotemporally expressed and displayed significant variability during various tissue development. Besides that, some PeXYLPs also respond to multiple phytohormones and abiotic stresses. By semi-quantitative RT-PCR, the response of some PeXYLPs to MeJA was confirmed, and the proteins were shown to localize to the plasma membrane mainly. WGCNA in defined regions of fast-growing bamboo shoots revealed that 5 PeXYLPs in 4 gene co-expression modules showed a positive module-trait relationship with three fast-growing regions. This systematic analysis of the PeXYLP family will provide a foundation for further insight into the functions of individual PeXYLP in a specific tissue or organ development, phytohormone perception, and stress responses in the future.

Transcriptome-wide analysis reveals core transcriptional regulators associated with culm development and variation in Dendrocalamus sinicus, the strongest woody bamboo in the world

Transcription factors (TFs) play indispensable roles in plant development and stress responses. As the largest woody bamboo species in the world, Dendrocalamus sinicus is endemic to Yunnan Province, China, and possesses two natural variants characterized by culm shape, namely straight or bent culms. Understanding the transcriptional regulation network of D. sinicus provides a unique opportunity to clarify the growth and development characteristics of woody bamboos. In this study, 10,236 TF transcripts belonging to 57 families were identified from transcriptome data of two variants at different developmental stages, from which we constructed a transcriptional regulatory network and unigene-coding protein-TFs interactive network of culm development for this attractive species. Gene function enrichment analysis revealed that hormone signaling and MAPK signaling pathways were two most enriched pathways in TF-regulated network. Based on PPI analysis, 50 genes interacting with nine TFs were screened as the core regulation components related to culm development. Of them, 18 synergistic genes of seven TFs, including nuclear cap-binding protein subunit 1, transcription factor GTE9-like, and ATP-dependent DNA helicase DDX11 isoform X1, involved in culm-shape variation. Most of these genes would interact with MYB, C3H, and ARF transcription factors. Six members with two each from ARF, C3H, and MYB transcription factor families and six key interacting genes (IAA3, IAA19, leucine-tRNA ligase, nuclear cap-binding protein subunit 1, elongation factor 2, and coiled-coil domain-containing protein 94) cooperate with these transcription factors were differentially expressed at development stage of young culms, and were validated by quantitative PCR. Our results represent a crucial step towards understanding the regulatory mechanisms of TFs involved in culm development and variation of D. sinicus.

Histological, metabolomic and transcriptomic analyses reveal mechanisms of cold acclimation of the Moso bamboo (Phyllostachys edulis) leaf

The Moso bamboo (Phyllostachys edulis) leaf copes well with cold winters in southeastern China. However, until now, there has been almost no research on its adaptation mechanisms to cold weather. Herein, we found that the Moso bamboo leaf has evolved several anatomical structures that may play a role in enhancing its cold tolerance. These structures include fewer fusiform cells, smaller bulliform cells, lower stomata density and many more trichomes, as well as lower relative water content than in the leaf of a cold-sensitive bamboo species, Bambusa ventricosa. Untargeted metabolomic analysis revealed that the winter leaf of Moso bamboo had 10- to 1000-fold higher stress-resistant metabolites such as glutathione, trehalose and ascorbic acid than the leaf of B. ventricosa on both warm and cold days. In contrast to the leaves that grew on a warm day, some metabolites such as glutathione and trehalose increased dramatically in the leaves of Moso bamboo that grew on a cold day. However, they unexpectedly decreased in the leaf of B. ventricosa growing at cold temperatures. Transcriptome analysis revealed a cold stress response network that includes trehalose, glutathione, flavonoid metabolism, DNA repair, reactive oxygen species degradation, stress-associated genes and abiotic stress-related plant hormones such as jasmonic acid, abscisic acid and ethylene. The potential mediator transcription factors, such as EREBP, HSF, MYB, NAC and WRYK, were also significantly upregulated in Moso bamboo leaves growing at cold temperatures. Interestingly, many newly identified genes were involved in the transcriptome of the winter leaf of the Moso bamboo. Most of these new genes have not even been annotated yet. The above results indicate that the Moso bamboo leaf has evolved special histological structures, metabolic pathways and a cold stress-tolerant transcriptome to adapt to the cold weather in its distribution areas.

The impact of genetic modified Ma bamboo on soil microbiome

Evaluating the potential alteration of microbial communities is a vital step for biosafety of genetic modified plants. Recently, we have produced genetic modified Ma bamboo with increased cold and drought tolerance by anthocyanin accumulation. In this work, we aim to study the potential effects on microbial communities in rhizosphere soils during the cultivation of genetic modified bamboo. Rhizosphere and surrounding soil were collected at 3-month post-transplant. The amplicon (16S rDNA and ITS1) were sequenced for analysis of bacterial and fungal communities. Multiple software and database (Picrust2, FAPROTAX and FUNGulid) were applied to predict and compare the microbial functions involving basic metabolisms, nitrogen usage and presence of plant pathogens. There were no substantial change of the structure and abundance of rhizosphere soil microbial communities between genetic modified and wild type bamboo. For the surrounding soil, the bacterial biota α-diversity increased (chao1: 1,001 ± 80-1,276 ± 84, observed species: 787 ± 52-1,194 ± 137, PD whole tree: 75 ± 4-117 ± 18) and fungal biota α-diversity decreased (chao1: 187 ± 18-145 ± 10) in samples of genetic modified bamboo compared to those of wild type bamboo. The microbiota predicted functions did not change or had no negative alteration between genetic modified and wild type bamboo, in both rhizosphere and surrounding soils. As a conclusion, the growth of genetic modified bamboo had no substantial change on rhizosphere soil microbial communities, while minor alteration on bamboo surrounding soil microbial communities with no harmful effects. Moreover, the genetic modified bamboo had no negative effect on the predicted functions of microbiota in soil.

Transcriptome Analysis of Moso Bamboo (Phyllostachys edulis) Reveals Candidate Genes Involved in Response to Dehydration and Cold Stresses

Bamboo (Bambusoideae) belongs to the grass family (Poaceae) and has been utilized as one of the most important nontimber forest resources in the world. Moso bamboo (Phyllostachys edulis) is a large woody bamboo with high ecological and economic values. Global climate change brings potential challenges to the normal growth of moso bamboo, and hence its production. Despite the release of moso bamboo genome sequence, the knowledge on genome-wide responses to abiotic stress is still limited. In this study, we generated a transcriptome data set with respect to dehydration and cold responses of moso bamboo using RNA-seq technology. The differentially expressed genes (DEGs) under treatments of dehydration and cold stresses were identified. By combining comprehensive gene ontology (GO) analysis, time-series analysis, and co-expression analysis, candidate genes involved in dehydration and cold responses were identified, which encode abscisic acid (ABA)/water deficit stress (WDS)-induced protein, late embryogenesis abundant (LEA) protein, 9-cis-epoxycarotenoid dioxygenase (NCED), anti-oxidation enzymes, transcription factors, etc. Additionally, we used PeLEA14, a dehydration-induced gene encoding an "atypical" LEA protein, as an example to validate the function of the identified stress-related gene in tolerance to abiotic stresses, such as drought and salt. In this study, we provided a valuable genomic resource for future excavation of key genes involved in abiotic stress responses and genetic improvement of moso bamboo to meet the requirement for environmental resilience and sustainable production.

Transcriptome analysis of the winter wheat Dn1 in response to cold stress

BACKGROUND

Heilongjiang Province has a long and cold winter season (the minimum temperature can reach -30 ℃), and few winter wheat varieties can safely overwinter. Dongnongdongmai1 (Dn1) is the first winter wheat variety that can safely overwinter in Heilongjiang Province. This variety fills the gap for winter wheat cultivation in the frigid region of China and greatly increases the land utilization rate. To understand the molecular mechanism of the cold response, we conducted RNA-sequencing analysis of Dn1 under cold stress.

RESULTS

Approximately 120,000 genes were detected in Dn1 under cold stress. The numbers of differentially expressed genes (DEGs) in the six comparison groups (0 ℃ vs. 5 ℃, -5 ℃ vs. 5 ℃, -10 ℃ vs. 5 ℃, -15 ℃ vs. 5 ℃, -20 ℃ vs. 5 ℃ and -25 ℃ vs. 5 ℃) were 11,313, 8313, 15,636, 13,671, 14,294 and 13,979, respectively. Gene Ontology functional annotation suggested that the DEGs under cold stress mainly had "binding", "protein kinase" and "catalytic" activities and were involved in "oxidation-reduction", "protein phosphorylation" and "carbohydrate metabolic" processes. Kyoto Encyclopedia of Genes and Genomes enrichment analysis indicated that the DEGs performed important functions in cold signal transduction and carbohydrate metabolism. In addition, major transcription factors (AP2/ERF, bZIP, NAC, WRKY, bHLH and MYB) participating in the Dn1 cold stress response were activated by low temperature.

CONCLUSION

This is the first study to explore the Dn1 transcriptome under cold stress. Our study comprehensively analysed the key genes involved in cold signal transduction and carbohydrate metabolism in Dn1 under cold stress. The results obtained by transcriptome analysis could help to further explore the cold resistance mechanism of Dn1 and provide basis for breeding of cold-resistant crops.

Comparative Transcriptome Analysis Reveals New Insight of Alfalfa (Medicago sativa L.) Cultivars in Response to Abrupt Freezing Stress

Freezing stress is a major limiting environmental factor that affects the productivity and distribution of alfalfa (Medicago sativa L.). There is growing evidence that enhancing freezing tolerance through resistance-related genes is one of the most efficient methods for solving this problem, whereas little is known about the complex regulatory mechanism of freezing stress. Herein, we performed transcriptome profiling of the leaves from two genotypes of alfalfa, freezing tolerance "Gannong NO.3" and freezing-sensitive "WL326GZ" exposure to -10°C to investigate which resistance-related genes could improve the freezing tolerance. Our results showed that a total of 121,366 genes were identified, and there were 7,245 differentially expressed genes (DEGs) between the control and treated leaves. In particular, the DEGs in "Gannong NO.3" were mainly enriched in the metabolic pathways and biosynthesis of secondary metabolites, and most of the DEGs in "WL326GZ" were enriched in the metabolic pathways, the biosynthesis of secondary metabolites, and plant-pathogen interactions. Moreover, the weighted gene co-expression network analysis (WGCNA) showed that ATP-binding cassette (ABC) C subfamily genes were strongly impacted by freezing stress, indicating that ABCC8 and ABCC3 are critical to develop the freezing tolerance. Moreover, our data revealed that numerous Ca2+ signal transduction and CBF/DREB1 pathway-related genes were severely impacted by the freezing resistance, which is believed to alleviate the damage caused by freezing stress. Altogether, these findings contribute the comprehensive information to understand the molecular mechanism of alfalfa adaptation to freezing stress and further provide functional candidate genes that can adapt to abiotic stress.

Characterization of the early gene expression profile in Populus ussuriensis under cold stress using PacBio SMRT sequencing integrated with RNA-seq reads

Populus ussuriensis is an important and fast-growing afforestation plant species in north-eastern China. The whole-genome sequencing of P. ussuriensis has not been completed. Also, the transcriptional network of P. ussuriensis response to cold stress remains unknown. To unravel the early response of P. ussuriensis to chilling (3 °C) stress and freezing (-3 °C) stresses at the transcriptional level, we performed single-molecule real-time (SMRT) and Illumina RNA sequencing for P. ussuriensis. The SMRT long-read isoform sequencing led to the identification of 29,243,277 subreads and 575,481 circular consensus sequencing reads. Approximately 50,910 high-quality isoforms were generated, and 2272 simple sequence repeats and 8086 long non-coding RNAs were identified. The Ca2+ content and abscisic acid (ABA) content in P. ussuriensis were significantly increased under cold stresses, while the value in the freezing stress treatment group was significantly higher than the chilling stress treatment group. A total of 49 genes that are involved in the signal transduction pathways related to perception and transmission of cold stress signals, such as the Ca2+ signaling pathway, ABA signaling pathway and MAPK signaling cascade, were found to be differentially expressed. In addition, 158 transcription factors from 21 different families, such as MYB, WRKY and AP2/ERF, were differentially expressed during chilling and freezing treatments. Moreover, the measurement of physiological indicators and bioinformatics observations demonstrated the altered expression pattern of genes involved in reactive oxygen species balance and the sugar metabolism pathway during chilling and freezing stresses. This is the first report of the early responses of P. ussuriensis to cold stress, which lays the foundation for future studies on the regulatory mechanisms in cold-stress response. In addition the full-length reference transcriptome of P. ussuriensis deciphered could be used in future studies on P. ussuriensis.

Bamboo Forest Mapping in China Using the Dense Landsat 8 Image Archive and Google Earth Engine

It is of great significance to understand the extent and distribution of bamboo for its valuable ecological services and economic benefits. However, it is challenging to map bamboo using remote sensing images over a large area because of the similarity between bamboo and other vegetation types, the availability of clear optical images, huge workload of image processing, and sample collection. In this study, we use the Landsat 8 times series images archive to map bamboo forests in China via the Google Earth engine. Several spectral indices were calculated and used as classification features, including the normalized difference vegetation index (NDVI), the normalized difference moisture index (NDMI) and textural features of the gray-level co-occurrence matrix (GLCM). We found that the bamboo forest covered an area of 709.92 × 104 hectares, with the provinces of Fujian, Jiangxi, and Zhejiang containing the largest area concentrations. The bamboo forest map was accurate and reliable with an average producer’s accuracy of 89.97%, user’s accuracy of 78.45% and kappa coefficient of 0.7789. In addition, bamboo was mainly distributed in forests with an elevation of 300–1200 m above sea level, average annual precipitation of 1200–1500 mm and average day land surface temperature of 19–25 °C. The NDMI is particularly useful in differentiating bamboo from other vegetation because of the clear difference in canopy moisture content, whilst NDVI and elevation are also helpful to improve the bamboo classification accuracy. The bamboo forest map will be helpful for bamboo forest industry planning and could be used for evaluating the ecological service of the bamboo forest.

Multi-omics analysis of cellular pathways involved in different rapid growth stages of moso bamboo

Moso bamboo (Phyllostachys edulis (Carriere) J. Houzeau) is a rapidly growing grass of industrial and ecological importance. However, the molecular mechanisms of its remarkable growth are not well understood. In this study, we investigated the early-stage growth of moso bamboo shoots and defined three different growth stages based on histological and biochemical analyses, namely, starting of cell division (SD), rapid division (RD) and rapid elongation (RE). Further analyses on potentially relevant cellular pathways in these growth stages using multi-omics approaches such as transcriptomics and proteomics revealed the involvement of multiple cellular pathways, including DNA replication, repair and ribosome biogenesis. A total of 8045 differentially expressed genes (DEGs) and 1053 differentially expressed proteins (DEPs) were identified in our analyses. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses of detected DEGs identified several key biological pathways such as phytohormone metabolism, signal transduction, cell wall development and carbohydrate metabolism. The comparative analysis of proteins displayed that a total of 213 DEPs corresponded with DEGs and 3 significant expression profiles that could be promoting the fast growth of bamboo internodes. Moreover, protein-protein interaction network prediction analysis is suggestive of the involvement of five major proteins of signal transduction, DNA synthesis and RNA transcription, and may act as key elements responsible for the rapid shoot growth. Our work exploits multi-omics and bioinformatic approaches to unfurl the complexity of molecular networks involved in the rapid growth of moso bamboo and opens up questions related to the interactions between the functions played by individual molecular pathway.