Transcriptome characterization of moso bamboo (Phyllostachys edulis) seedlings in response to exogenous gibberellin applications

Integrating Physiological Features and Proteomic Analyses Provides New Insights in Blue/Red Light-Treated Moso Bamboo (Phyllostachys edulis)

Bamboo is one of the most important nontimber forestry products in the world. Light is not only the most critical source of energy for plant photosynthesis but also involved in regulating the biological processes of plants. However, there are few reports on how blue/red light affects Moso bamboo. This study investigated the growth status and physiological responses of Moso bamboo (Phyllostachys edulis) to blue/red light treatments. The growth status of the bamboo plants was evaluated, revealing that both blue- and red-light treatments promoted plant height and overall growth. Gas exchange parameters, chlorophyll fluorescence, and enzyme activity were measured to assess the photosystem response of Moso bamboo to light treatments. Additionally, the blue light treatment led to a higher chlorophyll content and enzyme activities compared to the red light treatment. A tandem mass tag quantitative proteomics approach identified significant changes in protein abundance under different light conditions with specific response proteins associated with distinct pathways, such as photosynthesis and starch metabolism. Overall, this study provides valuable insights into the physiological and proteomic responses of Moso bamboo to blue/red light treatments, highlighting their potential impact on growth and development.

GATA transcription factor in common bean: A comprehensive genome-wide functional characterization, identification, and abiotic stress response evaluation

The GATA transcription factors (TFs) have been extensively studied for its regulatory role in various biological processes in many plant species. The functional and molecular mechanism of GATA TFs in regulating tolerance to abiotic stress has not yet been studied in the common bean. This study analyzed the functional identity of the GATA gene family in the P. vulgaris genome under different abiotic and phytohormonal stress. The GATA gene family was systematically investigated in the P. vulgaris genome, and 31 PvGATA TFs were identified. The study found that 18 out of 31 PvGATA genes had undergone duplication events, emphasizing the role of gene duplication in GATA gene expansion. All the PvGATA genes were classified into four significant subfamilies, with 8, 3, 6, and 13 members in each subfamily (subfamilies I, II, III, and IV), respectively. All PvGATA protein sequences contained a single GATA domain, but subfamily II members had additional domains such as CCT and tify. A total of 799 promoter cis-regulatory elements (CREs) were predicted in the PvGATAs. Additionally, we used qRT-PCR to investigate the expression profiles of five PvGATA genes in the common bean roots under abiotic conditions. The results suggest that PvGATA01/10/25/28 may play crucial roles in regulating plant resistance against salt and drought stress and may be involved in phytohormone-mediated stress signaling pathways. PvGATA28 was selected for overexpression and cloned into N. benthamiana using Agrobacterium-mediated transformation. Transgenic lines were subjected to abiotic stress, and results showed a significant tolerance of transgenic lines to stress conditions compared to wild-type counterparts. The seed germination assay suggested an extended dormancy of transgenic lines compared to wild-type lines. This study provides a comprehensive analysis of the PvGATA gene family, which can serve as a foundation for future research on the function of GATA TFs in abiotic stress tolerance in common bean plants.

Effects of exogenous GA3 on stem secondary growth of Pinus massoniana seedlings

Gibberellins (GAs) play a crucial role in regulating secondary growth in angiosperms, but their effects on the secondary growth of gymnosperms are rarely reported. In this study, we administered exogenous GA3 to two-year-old P. massoniana seedlings, and examined its effects on anatomical structure, physiological and biochemical changes, and gene expression in stems. The results showed that exogenous GA3 could enhance xylem development in P. massoniana by promoting cell division. The content of endogenous hormone (including auxins, brassinosteroids, and gibberellins) were changed and the genes related to phytohormone biosynthesis and signaling pathway, such as GID1, DELLA, TIR1, ARF, SAUR, CPD, BR6ox1, and CYCD3, were differentially expressed under GA3 treatment. Furthermore, GA3 and BR (brassinosteroid) might act synergistically in promoting secondary growth in P. massoniana. Additionally, lignin content was significantly increased after GA3 treatment accompanied by the express of lignin biosynthesis related genes. PmCAD (TRINITY_DN142116_c0_g1), a crucial gene involved in the lignin biosynthesis, was cloned and overexpressed in Nicotiana benthamiana, significantly promoting the xylem development and enhancing stem lignification. It was regarded as a key candidate gene for improving stem growth of P. massoniana. The findings of this study have demonstrated the impact of GA3 treatment on secondary growth of stems in P. massoniana, providing a foundation for understanding the molecular regulatory mechanism of stem secondary growth in Pinaceae seedlings and offering theoretical guidance for cultivating new germplasm with enhanced growth and yield.

Insight into gene expression associated with DNA methylation and small RNA in the rhizome-root system of Moso bamboo

Moso bamboo (Phyllostachys edulis), typically a monopodial scattering bamboo, is famous for its rapid growth. The rhizome-root system of Moso bamboo plays a crucial role in its clonal growth and spatial distribution. However, few studies have focused on rhizome-root systems. Here we collected LBs, RTs, and RGFNSs, the most important parts of the rhizome-root system, to study the molecular basis of the rapid growth of Moso bamboo due to epigenetic changes, such as DNA modifications and small RNAs. The angle of the shoot apical meristem of LB gradually decreased with increasing distance from the mother plant, and the methylation levels of LB were much higher than those of RT and RGFNS. 24 nt small RNAs and mCHH exhibited similar distribution patterns in transposable elements, suggesting a potential association between these components. The miRNA abundance of LB gradually increased with increasing distance from the mother plant, and a negative correlation was observed between gene expression levels and mCG and mCHG levels in the gene body. This study paves the way for further exploring the effects of epigenetic factors on the physiology of Moso bamboo.

GA20ox Family Genes Mediate Gibberellin and Auxin Crosstalk in Moso bamboo (Phyllostachys edulis)

Moso bamboo (Phyllostachys edulis) is one of the fastest growing plants. Gibberellin (GA) is a key phytohormone regulating growth, but there are few studies on the growth of Moso bamboo regulated by GA. The gibberellin 20 oxidase (GA20ox) gene family was targeted in this study. Chromosomal distribution and collinearity analysis identified 10 GA20ox genes evenly distributed on chromosomes, and the family genes were relatively conservative in evolution. The genetic relationship of GA20ox genes had been confirmed to be closest in different genera of plants in a phylogenetic and selective pressure analysis between Moso bamboo and rice. About 1/3 GA20ox genes experienced positive selective pressure with segmental duplication being the main driver of gene family expansion. Analysis of expression patterns revealed that only six PheGA20ox genes were expressed in different organs of shoot development and flowers, that there was redundancy in gene function. Underground organs were not the main site of GA synthesis in Moso bamboo, and floral organs are involved in the GA biosynthesis process. The auxin signaling factor PheARF47 was located upstream of PheGA20ox3 and PheGA20ox6 genes, where PheARF47 regulated PheGA20ox3 through cis-P box elements and cis-AuxRR elements, based on the result that promoter analysis combined with yeast one-hybrid and dual luciferase detection analysis identified. Overall, we identified the evolutionary pattern of PheGA20ox genes in Moso bamboo and the possible major synthesis sites of GA, screened for key genes in the crosstalk between auxin and GA, and laid the foundation for further exploration of the synergistic regulation of growth by GA and auxin in Moso bamboo.

Single-cell transcriptome atlas reveals spatiotemporal developmental trajectories in the basal roots of moso bamboo (Phyllostachys edulis)

Roots are essential for plant growth and development. Bamboo is a large Poaceae perennial with 1642 species worldwide. However, little is known about the transcriptional atlas that underpins root cell-type differentiation. Here, we set up a modified protocol for protoplast preparation and report single-cell transcriptomes of 14 279 filtered single cells derived from the basal root tips of moso bamboo. We identified four cell types and defined new cell-type-specific marker genes for the basal root. We reconstructed the developmental trajectories of the root cap, epidermis, and ground tissues and elucidated critical factors regulating cell fate determination. According to in situ hybridization and pseudotime trajectory analysis, the root cap and epidermis originated from a common initial cell lineage, revealing the particularity of bamboo basal root development. We further identified key regulatory factors for the differentiation of these cells and indicated divergent root developmental pathways between moso bamboo and rice. Additionally, PheWOX13a and PheWOX13b ectopically expressed in Arabidopsis inhibited primary root and lateral root growth and regulated the growth and development of the root cap, which was different from WOX13 orthologs in Arabidopsis. Taken together, our results offer an important resource for investigating the mechanism of root cell differentiation and root system architecture in perennial woody species of Bambusoideae.

Multi-omics of Circular RNAs and Their Responses to Hormones in Moso Bamboo (Phyllostachys edulis)

Circular RNAs (circRNAs) are endogenous non-coding RNAs with covalently closed structures, which have important functions in plants. However, their biogenesis, degradation, and function upon treatment with gibberellins (GAs) and auxins (1-naphthaleneacetic acid, NAA) remain unknown. Here, we systematically identified and characterized the expression patterns, evolutionary conservation, genomic features, and internal structures of circRNAs using RNase R-treated libraries from moso bamboo (Phyllostachys edulis) seedlings. Moreover, we investigated the biogenesis of circRNAs dependent on both cis- and trans-regulation. We explored the function of circRNAs, including their roles in regulating microRNA (miRNA)-related genes and modulating the alternative splicing of their linear counterparts. Importantly, we developed a customized degradome sequencing approach to detect miRNA-mediated cleavage of circRNAs. Finally, we presented a comprehensive view of the participation of circRNAs in the regulation of hormone metabolism upon treatment of bamboo seedlings with GA and NAA. Collectively, our study provides insights into the biogenesis, function, and miRNA-mediated degradation of circRNAs in moso bamboo.

Phytohormone Crosstalk of Cytokinin Biosynthesis and Signaling Family Genes in Moso Bamboo (Phyllostachys edulis)

Cytokinin is widely involved in the regulation of plant growth, but its pathway-related genes have not been reported in Moso bamboo. In this study, a total of 129 candidate sequences were identified by bioinformatic methods. These included 15 IPT family genes, 19 LOG family genes, 22 HK family genes, 11 HP family genes and 62 RR family genes. Phylogenetic analysis revealed that the cytokinin pathway was closely related to rice, and evolutionary pattern analysis found that most of the genes have syntenic relationship with rice-related genes. The Moso bamboo cytokinin pathway was evolutionarily conservative and mainly underwent purifying selection, and that gene family expansion was mainly due to whole-gene duplication events. Analysis of transcriptome data revealed a tissue-specific expression pattern of Moso bamboo cytokinin family genes, with auxin and gibberellin response patterns. Analysis of co-expression patterns at the developmental stages of Moso bamboo shoots revealed the existence of a phytohormone co-expression pattern centered on cytokinin signaling genes. The auxin signaling factor PheARF52 was identified by yeast one-hybrid assay as regulating the PheRR3 gene through a P-box element in the PheRR3 promoter region. Auxin and cytokinin signaling crosstalk to regulate Moso bamboo growth. Overall, we systematically identified and analyzed key gene families of the cytokinin pathway in Moso bamboo and obtained key factors for auxin and cytokinin crosstalk, laying the foundation for the study of hormone regulation in Moso bamboo.

Plasticity in the Morphology of Growing Bamboo: A Bayesian Analysis of Exogenous Treatment Effects on Plant Height, Internode Length, and Internode Numbers

Sucrose (Suc) and gibberellin (GA) can promote the elongation of certain internodes in bamboo. However, there is a lack of field studies to support these findings and no evidence concerning how Suc and GA promote the plant height of bamboo by regulating the internode elongation and number. We investigated the plant height, the length of each internode, and the total number of internodes of Moso bamboo (Phyllostachys edulis) under exogenous Suc, GA, and control group (CTRL) treatments in the field and analyzed how Suc and GA affected the height of Moso bamboo by promoting the internode length and number. The lengths of the 10th-50th internodes were significantly increased under the exogenous Suc and GA treatments, and the number of internodes was significantly increased by the exogenous Suc treatment. The increased effect of Suc and GA exogenous treatment on the proportion of longer internodes showed a weakening trend near the plant height of 15-16 m compared with the CTRL, suggesting that these exogenous treatments may be more effective in regions where bamboo growth is suboptimal. This study demonstrated that both the exogenous Suc and GA treatments could promote internode elongation of Moso bamboo in the field. The exogenous GA treatment had a stronger effect on internode elongation, and the exogenous Suc treatment had a stronger effect on increasing the internode numbers. The increase in plant height by the exogenous Suc and GA treatments was promoted by the co-elongation of most internodes or the increase in the proportion of longer internodes.

Phylogeny, transcriptional profile, and auxin-induced phosphorylation modification characteristics of conserved PIN proteins in Moso bamboo (Phyllostachys edulis)

Auxin polar transport is an important way for auxin to exercise its function, and auxin plays an irreplaceable role in the rapid growth of Moso bamboo. We identified and performed the structural analysis of PIN-FORMED auxin efflux carriers in Moso bamboo and obtained a total of 23 PhePIN genes from five gene subfamilies. We also performed chromosome localization and intra- and inter-species synthesis analysis. Phylogenetic analyses of 216 PIN genes showed that PIN genes are relatively conserved in the evolution of the Bambusoideae and have undergone intra-family segment replication in Moso bamboo. The PIN genes' transcriptional patterns showed that the PIN1 subfamily plays a major regulatory role. PIN genes and auxin biosynthesis maintain a high degree of consistency in spatial and temporal distribution. Phosphoproteomics analysis identified many phosphorylated protein kinases that respond to auxin regulation through autophosphorylation and phosphorylation of PIN proteins. The protein interaction network showed that there is a plant hormone interaction regulatory network with PIN protein as the core. We provide a comprehensive PIN protein analysis that complements the auxin regulatory pathway in Moso bamboo and paves the way for further auxin regulatory studies in bamboo.

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.

Rapid growth of Moso bamboo (Phyllostachys edulis): Cellular roadmaps, transcriptome dynamics, and environmental factors

Moso bamboo (Phyllostachys edulis) shows remarkably rapid growth (114.5 cm/day), but the underlying biological mechanisms remain unclear. After examining more than 12,750 internodes from more than 510 culms from 17 Moso populations, we identified internode 18 as a representative internode for rapid growth. This internode includes a 2-cm cell division zone (DZ), a cell elongation zone up to 12 cm, and a secondary cell wall (SCW) thickening zone. These zones elongated 11.8 cm, produced approximately 570,000,000 cells, and deposited ∼28 mg g-1 dry weight (DW) lignin and ∼44 mg g-1 DW cellulose daily, far exceeding vegetative growth observed in other plants. We used anatomical, mathematical, physiological, and genomic data to characterize development and transcriptional networks during rapid growth in internode 18. Our results suggest that (1) gibberellin may directly trigger the rapid growth of Moso shoots, (2) decreased cytokinin and increased auxin accumulation may trigger cell DZ elongation, and (3) abscisic acid and mechanical pressure may stimulate rapid SCW thickening via MYB83L. We conclude that internode length involves a possible tradeoff mediated by mechanical pressure caused by rapid growth, possibly influenced by environmental temperature and regulated by genes related to cell division and elongation. Our results provide insight into the rapid growth of Moso bamboo.

Genome-wide identification and characterization of GATA family genes in wheat

BACKGROUND

Transcription factors GATAs were a member of zinc finger protein, which could bind DNA regulatory regions to control expression of target genes, thus influencing plant growth and development either in normal condition or environmental stresses. Recently, GATA genes have been found and functionally characterized in a number of plant species. However, little information of GATA genes were annotated in wheat.

RESULTS

In the current study, 79 GATA genes were identified in wheat, which were unevenly located on 21 chromosomes. According to the analysis of phylogenetic tree and functional domain structures, TaGATAs were classified into four subfamilies (I, II, III, and IV), consist of 35, 21, 12, and 11 genes, respectively. Meanwhile, the amino acids of 79 TaGATAs exhibited apparent difference in four subfamilies according to GATA domains comparison, gene structures and conserved motif analysis. We then analyze the gene duplication and synteny between the genomes of wheat and Arabidopsis, rice and barley, which provided insights into evolutionary characteristics. In addition, expression patterns of TaGATAs were analyzed, and they showed obvious difference in diverse tissues and abiotic stresses.

CONCLUSION

In general, these results provide useful information for future TaGATA gene function analysis, and it helps to better understand molecular breeding and stress response in wheat.

Genome-Wide Identification of NAP1 and Function Analysis in Moso Bamboo (Phyllostachys edulis)

The nucleosome assembly protein 1 (NAP1) family is the main histone chaperone of histone H2A–H2B. To explore the function of NAP1 family genes in moso bamboo (Phyllostachys edulis), characterized by extremely rapid growth and a long flowering cycle, we originally conducted a genome-wide analysis of the PheNAP1 gene. The phylogenetic relationship, gene expression pattern, DNA methylation, and histone modification were analyzed. Eventually, 12 PheNAP1 genes were recognized from the Phyllostachys edulis genome, divided into two sorts: the NRP subfamily (four members) and the NAP subfamily (eight members). Highly conserved motifs exist in each subfamily, which are distinct between subfamilies. PheNAP1 was distributed homogeneously on 10 out of 24 chromosomes, and gene duplication contributed significantly to the enhancement of the PheNAP1 gene in the genome. Cis-acting element analysis showed that PheNAP1 family genes are involved in light, hormone, and abiotic stress responses and may play an important role in the rapid growth and flowering. PheNAP1 exhibited the highest expression level in fast-growing shoots, indicating it is closely associated with the rapid growth of moso bamboo. Besides, PheNAP1 can rescue the early-flowering phenotype of nrp1-1 nrp2-2, and it affected the expression of genes related to the flowering pathway, like BSU1, suggesting the vital role that PheNAP1 may take in the flowering process of moso bamboo. In addition, histone modification results showed that PheNAP1 could bind to phosphorylation-, acetylation-, and methylation-modified histones to further regulate gene expression. A sketch appears: that PheNAP1 can accompany histones to regulate fast-growth- and flowering-related genes in moso bamboo. The consequences of this study enrich the understanding of the epigenetic regulation mechanism of bamboo plants and lays a foundation for further studies on the role of the NAP1 gene in Phyllostachys edulis and the function of chromatin regulation in forest growth and development.

Gibberellin inhibition of taproot formation by modulation of DELLA-NAC complex activity in turnip (Brassica rapa var. rapa)

Turnip is a member of the Brassica rapa species and is characterized by a swollen taproot that develops from the hypocotyl and part of the root. Gibberellins (GAs) are plant growth regulators involved in promoting cell elongation and play important roles in many aspects of plant growth and development. Interestingly, exogenous application of GA₃ was found to significantly inhibit taproot formation in turnip. Moreover, endogenous GA contents decreased during the early developmental stages of taproot formation, suggesting that GA plays a negative role in taproot formation. We examined the anatomical structure of the taproot and found that lignification of the xylem cell wall was enhanced after treatment with GA₃. Yeast two-hybrid assays suggested the occurrence of protein interactions between DELLAs and NACs in turnip. We also found that the expression of NAC-targeted genes involved in lignification of the secondary cell wall was significantly upregulated upon GA₃ treatment. Taken together, these results supported the hypothesis that GA induced DELLA proteins degradation to release NAC proteins and induced xylem lignification, therefore inhibiting taproot formation, providing new insight into the molecular mechanism underlying turnip taproot formation.

MicroRNAs play important roles in regulating the rapid growth of the Phyllostachys edulis culm internode

Moso bamboo (Phyllostachys edulis) is a fast-growing species with uneven growth and lignification from lower to upper segments within one internode. MicroRNAs (miRNAs) play a vital role in post-transcriptional regulation in plants. However, how miRNAs regulate fast growth in bamboo internodes is poorly understood. In this study, one moso bamboo internode was divided during early rapid growth into four segments called F4 (bottom) to F1 (upper) and these were then analysed for transcriptomes, miRNAs and degradomes. The F4 segment had a higher number of actively dividing cells as well as a higher content of auxin (IAA), cytokinin (CK) and gibberellin (GA) compared with the F1 segment. RNA-seq analysis showed DNA replication and cell division-associated genes highly expressed in F4 rather than in F1. In total, 63 miRNAs (DEMs) were identified as differentially expressed between F4 and F1. The degradome and the transcriptome indicated that many downstream transcription factors and hormonal responses genes were modulated by DEMs. Several miR-target interactions were further validated by tobacco co-infiltration. Our findings give new insights into miRNA-mediated regulatory pathways in bamboo, and will contribute to a comprehensive understanding of the molecular mechanisms governing rapid growth.

Study on the Regulatory Effects of GA3 on Soybean Internode Elongation

Excessive plant height is an important factor that can lead to lodging, which is closely related to soybean yield. Gibberellins are widely used as plant growth regulators in agricultural production. Gibberellic acid (GA₃), one of the most effective active gibberellins, has been used to regulate plant height and increase yields. The mechanism through which GA₃ regulates internode elongation has been extensively investigated. In 2019 and 2020, we applied GA₃ to the stems, leaves, and roots of two soybean cultivars, Heinong 48 (a high-stalk cultivar) and Henong 60 (a dwarf cultivar), and GA₃ was also applied to plants whose apical meristem was removed or to girded plants to compare the internode length and stem GA₃ content of soybean plants under different treatments. These results suggested that the application of GA₃ to the stems, leaves, and roots of soybean increased the internode length and GA₃ content in the stems. Application of GA₃ decreased the proportion of the pith in the soybean stems and primary xylem while increasing the proportion of secondary xylem. The apical meristem is an important site of GA₃ synthesis in soybean stems and is involved in the regulation of stem elongation. GA₃ was shown to be transported acropetally through the xylem and laterally between the xylem and phloem in soybean stems. We conclude that the GA₃ level in stems is an important factor affecting internode elongation.

Production of purple Ma bamboo (Dendrocalamus latiflorus Munro) with enhanced drought and cold stress tolerance by engineering anthocyanin biosynthesis

Overexpression of the leaf color (Lc) gene in Ma bamboo substantially increased the accumulation level of anthocyanin, and improved plant tolerance to cold and drought stresses, probably due to the increased antioxidant capacity. Most bamboos, including Ma bamboo (Dendrocalamus latiflorus Munro), are naturally evergreen and sensitive to cold and drought stresses, while it's nearly impossible to make improvements through conventual breeding due to their long and irregular flowering habit. Moreover, few studies have reported bamboo germplasm innovation through genetic engineering as bamboo genetic transformation remains difficult. In this study, we have upregulated anthocyanin biosynthesis in Ma bamboo, to generate non-green Ma bamboo with increased abiotic stress tolerance. By overexpressing the maize Lc gene, a bHLH transcription activator involved in the anthocyanin biosynthesis in Ma bamboo, we generated purple bamboos with increased anthocyanin levels including cyanidin-3-O-rutinoside, peonidin 3-O-rutinoside, and an unknown cyanidin pentaglycoside derivative. The expression levels of 9 anthocyanin biosynthesis genes were up-regulated. Overexpression of the Lc gene improved the plant tolerance to cold and drought stress, probably due to increased antioxidant capacity. The levels of the cold- and drought-related phytohormone jasmonic acid in the transgenic plants were also enhanced, which may also contribute to the plant stress-tolerant phenotypes. High anthocyanin accumulation level did not affect plant growth. Transcriptomic analysis showed higher expressions of genes involved in the flavonoid pathway in Lc transgenic bamboos compared with those in wild-type ones. The anthocyanin-rich bamboos generated here provide an example of ornamental and multiple agronomic trait improvements by genetic engineering in this important grass species.

The BBX gene family in Moso bamboo (Phyllostachys edulis): identification, characterization and expression profiles

BACKGROUND

The BBX (B-box) family are zinc finger protein (ZFP) transcription factors that play an essential role in plant growth, development and response to abiotic stresses. Although BBX genes have been characterized in many model organisms, genome-wide identification of the BBX family genes have not yet been reported in Moso bamboo (Phyllostachys edulis), and the biological functions of this family remain unknown.

RESULT

In the present study, we identified 27 BBX genes in the genome of Moso bamboo, and analysis of their conserved motifs and multiple sequence alignments revealed that they all shared highly similar structures. Additionally, phylogenetic and homology analyses indicated that PeBBX genes were divided into three clusters, with whole-genome duplication (WGD) events having facilitated the expansion of this gene family. Light-responsive and stress-related cis-elements were identified by analyzing cis-elements in the promoters of all PeBBX genes. Short time-series expression miner (STEM) analysis revealed that the PeBBX genes had spatiotemporal-specific expression patterns and were likely involved in the growth and development of bamboo shoots. We further explored the downstream target genes of PeBBXs, and GO/KEGG enrichment analysis predicted multiple functions of BBX target genes, including those encoding enzymes involved in plant photosynthesis, pyruvate metabolism and glycolysis/gluconeogenesis.

CONCLUSIONS

In conclusion, we analyzed the PeBBX genes at multiple different levels, which will contribute to further studies of the BBX family and provide valuable information for the functional validation of this family.

Transcriptome analysis reveals key genes regulating signaling and metabolic pathways during the growth of moso bamboo (Phyllostachys edulis) shoots

Moso bamboo (Phyllostachys edulis), a high-value bamboo used to produce food (young shoots), building, and industrial goods. To explore key candidate genes regulating signal transduction and metabolic processes during the initiation of stem elongation in moso bamboo, a transcriptome analysis of the shoots during three successive early elongation stages was performed. From cluster and differential expression analyses, 2984 differentially expressed genes (DEGs) were selected for an enrichment analysis. The DEGs were significantly enriched in the plant hormone signal transduction, sugar and starch metabolism, and energy metabolism pathways. Consequently, the DEG expression patterns of these pathways were analyzed, and the plant endogenous hormone and carbon metabolite (including sucrose, total soluble sugar, and starch) contents for each growth stage, of the shoot, were determined. The cytokinin-signaling pathway was continuously active in the three successive elongation stages, in which several cytokinin-signaling genes played indispensable roles. Additionally, many key DEGs regulating sugar, starch metabolism, and energy conversion, which are actively involved in energy production and substrate synthesis during the continuous growth of the shoots, were found. In summary, our study lays a foundation for understanding the mechanisms of moso bamboo growth and provides useful gene resources for breeding through genetic engineering.

Genome-wide identification, evolution and expression analysis of the aspartic protease gene family during rapid growth of moso bamboo (Phyllostachys edulis) shoots

BACKGROUND

Aspartic proteases (APs) are a class of aspartic peptidases belonging to nine proteolytic enzyme families whose members are widely distributed in biological organisms. APs play essential functions during plant development and environmental adaptation. However, there are few reports about APs in fast-growing moso bamboo.

RESULT

In this study, we identified a total of 129 AP proteins (PhAPs) encoded by the moso bamboo genome. Phylogenetic and gene structure analyses showed that these 129 PhAPs could be divided into three categories (categories A, B and C). The PhAP gene family in moso bamboo may have undergone gene expansion, especially the members of categories A and B, although homologs of some members in category C have been lost. The chromosomal location of PhAPs suggested that segmental and tandem duplication events were critical for PhAP gene expansion. Promoter analysis revealed that PhAPs in moso bamboo may be involved in plant development and responses to environmental stress. Furthermore, PhAPs showed tissue-specific expression patterns and may play important roles in rapid growth, including programmed cell death, cell division and elongation, by integrating environmental signals such as light and gibberellin signals.

CONCLUSION

Comprehensive analysis of the AP gene family in moso bamboo suggests that PhAPs have experienced gene expansion that is distinct from that in rice and may play an important role in moso bamboo organ development and rapid growth. Our results provide a direction and lay a foundation for further analysis of plant AP genes to clarify their function during rapid growth.

Profiling of circular RNA N6 -methyladenosine in moso bamboo (Phyllostachys edulis) using nanopore-based direct RNA sequencing

N⁶ -methyladenosine (m⁶ A) is a prevalent modification in messenger RNAs and circular RNAs that play important roles in regulating various aspects of RNA metabolism. However, the occurrence of the m⁶ A modification in plant circular RNAs has not been reported. A widely used method to identify m⁶ A modifications relies on m⁶ A-specific antibodies followed by next-generation sequencing of precipitated RNAs (MeRIP-Seq). However, one limitation of MeRIP-Seq is that it does not provide the precise location of m⁶ A at single-nucleotide resolution. Although more recent sequencing techniques such as Nanopore-based direct RNA sequencing (DRS) can overcome such limitations, the technology does not allow sequencing of circular RNAs, as these molecules lack a poly(A) tail. Here, we developed a novel method to detect the precise location of m⁶ A modifications in circular RNAs using Nanopore DRS. We first enriched our samples for circular RNAs, which we then fragmented and sequenced on the Nanopore platform with a customized protocol. Using this method, we identified 470 unique circular RNAs from DRS reads based on the back-spliced junction region. Among exonic circular RNAs, about 10% contained m⁶ A sites, which mainly occurred around acceptor and donor splice sites. This study demonstrates the utility of our antibody-independent method in identifying total and methylated circular RNAs using Nanopore DRS. This method has the additional advantage of providing the exact location of m⁶ A sites at single-base resolution in circular RNAs or linear transcripts from non-coding RNA without poly(A) tails.

Enhanced Activity of Genes Associated With Photosynthesis, Phytohormone Metabolism and Cell Wall Synthesis Is Involved in Gibberellin-Mediated Sugarcane Internode Growth

Internode elongation is an important trait in sugarcane as it affects the sugarcane yield. Gibberellin (GA) is a key modulator of internode elongation in sugarcane. Understanding the gene expression features of GA-mediated internode elongation has both scientific and practical significance. This study aimed to examine the transcriptomic changes in the internode elongation of sugarcane following GA treatment. Eighteen cDNA libraries from the internode tissues on days of 0, 3, and 6 in control and GA treatment groups were sequenced and their gene expression were studied. RNA-seq analysis revealed 1,338,723,248 reads and 70,821 unigenes from elongating internodes of sugarcane. Comparative studies discovered a large number of transcripts that were differentially expressed in GA-treated samples compared to the control. Further analysis revealed that the differentially expressed genes were enriched in the metabolic process, one-carbon compound transport, and single-organism process. Kyoto Encyclopedia of Genes and Genomes pathway annotation showed significant enrichment in photosynthesis and plant hormone signal transduction, indicating its involvement in internode elongation. The function analysis suggested that metabolic pathways and biosynthesis of secondary metabolites, plant hormones, and cell wall components were enriched in the internodes of the GA-treated plants. The hub genes were identified, with the function of cellulose synthesis. The results of this study provide a global view of mRNA changes during sugarcane internode elongation and extend our knowledge of the GA-mediated cellular processes involved in sugarcane stem growth.

Full Transcriptome Analysis of Callus Suspension Culture System of Bletilla striata

Background

Bletilla striata has been widely used in the pharmacology industry. To effectively produce the secondary metabolites through suspension cultured cells of B. striata, it is important to exploring the full-length transcriptome data and the genes related to cell growth and chemical producing of all culture stages. We applied a combination of Real-Time Sequencing of Single Molecule (SMRT) and second-generation sequencing (SGS) to generate the complete and full-length transcriptome of B. striata suspension cultured cells.

Methods

The B. striata transcriptome was formed in de novo way by using PacBio isoform sequencing (Iso-Seq) on a pooled RNA sample derived from 23 samples of 10 culture stages, to explore the potential for capturing full-length transcript isoforms. All unigenes were obtained after splicing, assembling, and clustering, and corrected by the SGS results. The obtained unigenes were compared with the databases, and the functions were annotated and classified.

Results and conclusions

A total of 100,276 high-quality full-length transcripts were obtained, with an average length of 2530 bp and an N50 of 3302 bp. About 52% of total sequences were annotated against the Gene Ontology, 53,316 unigenes were hit by KOG annotations and divided into 26 functional categories, 80,020 unigenes were mapped by KEGG annotations and clustered into 363 pathways. Furthermore, 15,133 long-chain non-coding RNAs (lncRNAs) were detected. And 68,996 coding sequences were identified based on SSR analysis, among which 31 pairs of primers selected at random were amplified and obtained stable bands. In conclusion, our results provide new full-length transcriptome data and genetic resources for identifying growth and metabolism-related genes, which provide a solid foundation for further research on its growth regulation mechanisms and genetic engineering breeding mechanisms of B. striata.

Physiological and transcriptomic analyses of brassinosteroid function in moso bamboo (Phyllostachys edulis) seedlings

This study demonstrates that brassinosteroid is essential for seedling and shoot growth in moso bamboo. The shoot of moso bamboo is known to grow extremely fast. The roles of phytohormones in such fast growth of bamboo shoot remain unclear. Here we reported that endogenous brassinosteroid (BR) is a major factor promoting bamboo shoot internode elongation. Reducing endogenous brassinosteroid level by its biosynthesis inhibitor propiconazole stunted shoot growth in seedling stage, whereas exogenous BR application promoted scale leaf elongation and the inclination of lamina joint of leaves and scale leaves. Genome-wide transcriptome analysis identified hundreds of genes whose expression levels are altered by BR and propiconazole in shoots and roots of bamboo seedling. The data show that BR regulates cell wall-related genes, hydrogen peroxide catabolic genes, and auxin-related genes. Our study demonstrates an essential role of BR in fast growth bamboo shoots and identifies a large number of BR-responsive genes in bamboo seedlings.

Transcriptome analysis of lateral buds from Phyllostachys edulis rhizome during germination and early shoot stages

BACKGROUND

The vegetative growth is an important stage for plants when they conduct photosynthesis, accumulate and collect all resources needed and prepare for reproduction stage. Bamboo is one of the fastest growing plant species. The rapid growth of Phyllostachys edulis results from the expansion of intercalary meristem at the basal part of nodes, which are differentiated from the apical meristem of rhizome lateral buds. However, little is known about the major signaling pathways and players involved during this rapid development stage of bamboo. To study this question, we adopted the high-throughput sequencing technology and compared the transcriptomes of Moso bamboo rhizome buds in germination stage and late development stage.

RESULTS

We found that the development of Moso bamboo rhizome lateral buds was coordinated by multiple pathways, including meristem development, sugar metabolism and phytohormone signaling. Phytohormones have fundamental impacts on the plant development. We found the evidence of several major hormones participating in the development of Moso bamboo rhizome lateral bud. Furthermore, we showed direct evidence that Gibberellic Acids (GA) signaling participated in the Moso bamboo stem elongation.

CONCLUSION

Significant changes occur in various signaling pathways during the development of rhizome lateral buds. It is crucial to understand how these changes are translated to Phyllostachys edulis fast growth. These results expand our knowledge on the Moso bamboo internodes fast growth and provide research basis for further study.

PSBR1, encoding a mitochondrial protein, is regulated by brassinosteroid in moso bamboo (Phyllostachys edulis)

PSBR1 is a moso bamboo gene negatively regulated by brassinosteroid, which encodes a mitochondrial localized protein. Overexpression of PSBR1 leads to growth inhibition in various growth progresses in Arabidopsis. The young shoot of moso bamboo (Phyllostachys edulis) is known as one of the fastest growing plant organs. The roles of phytohormones in the fast-growth of bamboo shoot are not fully understood. Brassinosteroids (BRs) are a group of growth-promoting steroid hormones that play important roles in cell elongation and division. While BR related genes are highly enriched in fast-growing internodes in moso bamboo, the functions of BR in the fast-growth process is not understood at the molecular level. Here, we identified a poaceae specific gene, PSBR1 (Poaceae specific and BR responsive gene 1) from the moso bamboo genome. PSBR1 was highly expressed in the stem and leaves of bamboo seedling, and the elongating nodes of fast-growing bamboo shoot. PSBR1's expression is increased by BR biosynthesis inhibitor propiconazole but decreased by BR treatment. PSBR1 encodes a novel protein that is localized to the mitochondria in tobacco and bamboo protoplast. The Arabidopsis transgenic plants overexpressing PSBR1 show growth inhibition in both vegetative and reproductive stages. This study suggests that PSBR1 is a BR regulated mitochondrial protein in bamboo, which inhibits plant growth when overexpressed in Arabidopsis.

Genome-Wide Characterization and Gene Expression Analyses of GATA Transcription Factors in Moso Bamboo (Phyllostachys edulis)

Moso bamboo is well-known for its rapid-growth shoots and widespread rhizomes. However, the regulatory genes of these two processes are largely unexplored. GATA transcription factors regulate many developmental processes, but their roles in moso bamboo height control and rhizome development remains unexplored. Here, thirty-one bamboo GATA factors (PeGATAs) were identified, which are evolutionarily closer to rice than Arabidopsis, and their gene expression patterns were analyzed in bamboo development and phytohormone response with bioinformatics and molecular methods. Interestingly, PeGATAs could only be classified into three groups. Phytohormone responsive cis-elements were found in PeGATA promoters and the expression profiles showed that PeGATA genes might respond to gibberellin acid and abscisic acid but not to auxin at the transcriptional level. Furthermore, PeGATA genes have a tissue-specific expression pattern in bamboo rhizomes. Interestingly, most PeGATA genes were down-regulated during the rapid-growth of bamboo shoots. In addition, over-expressing one of the PeGATA genes, PeGATA26, significantly repressed the primary root length and plant height of transgenic Arabidopsis plants, which may be achieved by promoting the gibberellin acid turnover. Overall, our results provide insight into the function of GATA transcription factors in bamboo, and into genetic resources for engineering plant height.

Large Scale Profiling of Protein Isoforms Using Label-Free Quantitative Proteomics Revealed the Regulation of Nonsense-Mediated Decay in Moso Bamboo (Phyllostachys edulis)

Moso bamboo is an important forest species with a variety of ecological, economic, and cultural values. However, the gene annotation information of moso bamboo is only based on the transcriptome sequencing, lacking the evidence of proteome. The lignification and fiber in moso bamboo leads to a difficulty in the extraction of protein using conventional methods, which seriously hinders research on the proteomics of moso bamboo. The purpose of this study is to establish efficient methods for extracting the total proteins from moso bamboo for following mass spectrometry-based quantitative proteome identification. Here, we have successfully established a set of efficient methods for extracting total proteins of moso bamboo followed by mass spectrometry-based label-free quantitative proteome identification, which further improved the protein annotation of moso bamboo genes. In this study, 10,376 predicted coding genes were confirmed by quantitative proteomics, accounting for 35.8% of all annotated protein-coding genes. Proteome analysis also revealed the protein-coding potential of 1015 predicted long noncoding RNA (lncRNA), accounting for 51.03% of annotated lncRNAs. Thus, mass spectrometry-based proteomics provides a reliable method for gene annotation. Especially, quantitative proteomics revealed the translation patterns of proteins in moso bamboo. In addition, the 3284 transcript isoforms from 2663 genes identified by Pacific BioSciences (PacBio) single-molecule real-time long-read isoform sequencing (Iso-Seq) was confirmed on the protein level by mass spectrometry. Furthermore, domain analysis of mass spectrometry-identified proteins encoded in the same genomic locus revealed variations in domain composition pointing towards a functional diversification of protein isoform. Finally, we found that part transcripts targeted by nonsense-mediated mRNA decay (NMD) could also be translated into proteins. In summary, proteomic analysis in this study improves the proteomics-assisted genome annotation of moso bamboo and is valuable to the large-scale research of functional genomics in moso bamboo. In summary, this study provided a theoretical basis and technical support for directional gene function analysis at the proteomics level in moso bamboo.

Selection of reliable reference genes for quantitative RT-PCR in garlic under salt stress

Quantitative real-time reverse-transcriptase PCR (qRT-PCR) has been frequently used for detecting gene expression. To obtain reliable results, selection of suitable reference genes is a fundamental and necessary step. Garlic (Allium sativum), a member from Alliaceae family, has been used both as a food flavoring and as a traditional medicine. In the present study, garlic plants were exposed to salt stress (200 mM NaCl) for 0, 1, 4 and 12 h, and garlic roots, bulbs, and leaves were harvested for subsequent analysis. The expression stability of eight candidate reference genes, eukaryotic translation initiation factor 4α (eIF-4α), actin (ACTIN), tubulin β-7 (TUB7), TAP42-interacting protein of 41 kDa (TIP41), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), SAND family protein (SAND), elongation factor 1 alpha (EF-1α), and protein phosphatase 2A (PP2A) were evaluated by geNorm, NormFinder, and BestKeeper. All genes tested displayed variable expression profiles under salt stress. In the leaf and root group, ACTIN was the best reference gene for normalizing gene expression. In garlic clove, ACTIN and SAND were the least variable, and were suitable for gene expression studies under salt stress; these two genes also performed well in all samples tested. Based on our results, we recommend that it is essential to use specific reference genes in different situations to obtain accurate results. Using a combination of multiple stable reference genes, such as ACTIN and SAND, to normalize gene expression is encouraged. The results from the study will be beneficial for accurate determination of gene expression in garlic and other plants.

GSK3/shaggy-like kinase 1 ubiquitously regulates cell growth from Arabidopsis to Moso bamboo (Phyllostachys edulis)

Moso bamboo (Phyllostachys edulis) is one of the fastest growing species with a maximum growth rate of 1 m/day. However, the regulator genes for this explosive growth phenomenon have not been functionally studied. Here, we found that Moso bamboo GSK3/shaggy-like kinase 1 (PeGSK1) acts as a negative regulator of cell growth. Over-expression of PeGSK1 in Arabidopsis showed significant growth arrest phenotypes, including dwarfism, small leaves, reduced cell length, and disturbed cell elongation of petiole. Furthermore, Overexpression of PeGSK1 fully inhibited the longer hypocotyl phenotype of Arabidopsis atgsk1 mutants. In addition, PeGSK1-overexpressing lines were resistant to exogenous BR treatment and PeGSK1 interacted with the brassinosteroid signal transduction key regulator BZR1. The BZR1-dependent cell growth genes were down-regulated in PeGSK1-overexpressing lines. These results indicated that PeGSK1 is functionally similar to AtGSK1 and inhibited cell growth via the brassinosteroid signaling pathway. Importantly, PeGSK1 also interacted with PeBZR1, and the expression pattern of PeGSK1 was negatively correlated with the internode elongation of bamboo, indicating that PeGSK1 is involved in the cell growth of bamboo. In summary, our results provide insight into the role of brassinosteroids in the rapid-growth of bamboo culms and identifying target genes for the genetic manipulation of plant height.

Genome-Wide Profiling of Circular RNAs in the Rapidly Growing Shoots of Moso Bamboo (Phyllostachys edulis)

Circular RNAs, including circular exonic RNAs (circRNA), circular intronic RNAs (ciRNA) and exon-intron circRNAs (EIciRNAs), are a new type of noncoding RNAs. Growing shoots of moso bamboo (Phyllostachys edulis) represent an excellent model of fast growth and their circular RNAs have not been studied yet. To understand the potential regulation of circular RNAs, we systematically characterized circular RNAs from eight different developmental stages of rapidly growing shoots. Here, we identified 895 circular RNAs including a subset of mutually inclusive circRNA. These circular RNAs were generated from 759 corresponding parental coding genes involved in cellulose, hemicellulose and lignin biosynthetic process. Gene co-expression analysis revealed that hub genes, such as DEFECTIVE IN RNA-DIRECTED DNA METHYLATION 1 (DRD1), MAINTENANCE OF METHYLATION (MOM), dicer-like 3 (DCL3) and ARGONAUTE 1 (AGO1), were significantly enriched giving rise to circular RNAs. The expression level of these circular RNAs presented correlation with its linear counterpart according to transcriptome sequencing. Further protoplast transformation experiments indicated that overexpressing circ-bHLH93 generating from transcription factor decreased its linear transcript. Finally, the expression profiles suggested that circular RNAs may have interplay with miRNAs to regulate their cognate linear mRNAs, which was further supported by overexpressing miRNA156 decreasing the transcript of circ-TRF-1 and linear transcripts of TRF-1. Taken together, the overall profile of circular RNAs provided new insight into an unexplored category of long noncoding RNA regulation in moso bamboo.

Analysis of Transcriptome and Epitranscriptome in Plants Using PacBio Iso-Seq and Nanopore-Based Direct RNA Sequencing

Nanopore sequencing from Oxford Nanopore Technologies (ONT) and Pacific BioSciences (PacBio) single-molecule real-time (SMRT) long-read isoform sequencing (Iso-Seq) are revolutionizing the way transcriptomes are analyzed. These methods offer many advantages over most widely used high-throughput short-read RNA sequencing (RNA-Seq) approaches and allow a comprehensive analysis of transcriptomes in identifying full-length splice isoforms and several other post-transcriptional events. In addition, direct RNA-Seq provides valuable information about RNA modifications, which are lost during the PCR amplification step in other methods. Here, we present a comprehensive summary of important applications of these technologies in plants, including identification of complex alternative splicing (AS), full-length splice variants, fusion transcripts, and alternative polyadenylation (APA) events. Furthermore, we discuss the impact of the newly developed nanopore direct RNA-Seq in advancing epitranscriptome research in plants. Additionally, we summarize computational tools for identifying and quantifying full-length isoforms and other co/post-transcriptional events and discussed some of the limitations with these methods. Sequencing of transcriptomes using these new single-molecule long-read methods will unravel many aspects of transcriptome complexity in unprecedented ways as compared to previous short-read sequencing approaches. Analysis of plant transcriptomes with these new powerful methods that require minimum sample processing is likely to become the norm and is expected to uncover novel co/post-transcriptional gene regulatory mechanisms that control biological outcomes during plant development and in response to various stresses.

De novo genome assembly of the stress tolerant forest species Casuarina equisetifolia provides insight into secondary growth

Casuarina equisetifolia (C. equisetifolia), a conifer-like angiosperm with resistance to typhoon and stress tolerance, is mainly cultivated in the coastal areas of Australasia. C. equisetifolia, making it a valuable model to study secondary growth associated genes and stress-tolerance traits. However, the genome sequence is unavailable and therefore wood-associated growth rate and stress resistance at the molecular level is largely unexplored. We therefore constructed a high-quality draft genome sequence of C. equisetifolia by a combination of Illumina second-generation sequencing reads and Pacific Biosciences single-molecule real-time (SMRT) long reads to advance the investigation of this species. Here, we report the genome assembly, which contains approximately 300 megabases (Mb) and scaffold size of N50 is 1.06 Mb. Additionally, gene annotation, assisted by a combination of prediction and RNA-seq data, generated 29 827 annotated protein-coding genes and 1983 non-coding genes, respectively. Furthermore, we found that the total number of repetitive sequences account for one-third of the genome assembly. Here we also construct the genome-wide map of DNA modification, such as two novel forms N⁶ -adenine (6mA) and N4-methylcytosine (4mC) at the level of single-nucleotide resolution using single-molecule real-time (SMRT) sequencing. Interestingly, we found that 17% of 6mA modification genes and 15% of 4mC modification genes also included alternative splicing events. Finally, we investigated cellulose, hemicellulose, and lignin-related genes, which were associated with secondary growth and contained different DNA modifications. The high-quality genome sequence and annotation of C. equisetifolia in this study provide a valuable resource to strengthen our understanding of the diverse traits of trees.

Increased Gibberellins and Light Levels Promotes Cell Wall Thickness and Enhance Lignin Deposition in Xylem Fibers

Light intensity and hormones (gibberellins; GAs) alter plant growth and development. A fine regulation triggered by light and GAs induces changes in stem cell walls (CW). Cross-talk between light-stimulated and GAs-induced processes as well as the phenolic compounds metabolism leads to modifications in lignin formation and deposition on cell walls. How these factors (light and GAs) promote changes in lignin content and composition. In addition, structural changes were evaluated in the stem anatomy of tobacco plants. GA₃ was sprayed onto the leaves and paclobutrazol (PAC), a GA biosynthesis inhibitor, via soil, at different irradiance levels. Fluorescence microscopy techniques were applied to detect lignin, and electron microscopy (SEM and TEM) was used to obtain details on cell wall structure. Furthermore, determination of total lignin and monomer contents were analyzed. Both light and GAs induces increased lignin content and CW thickening as well as greater number of fiber-like cells but not tracheary elements. The assays demonstrate that light exerts a role in lignification under GA₃ supplementation. In addition, the existence of an exclusive response mechanism to light was detected, that GAs are not able to replace.