Characterization of the developmental dynamics of the elongation of a bamboo internode during the fast growth stage

Identification and characterization of FBA genes in moso bamboo reveals PeFBA8 related to photosynthetic carbon metabolism

Fructose 1,6-bisphosphate aldolase (FBA) is a pivotal enzyme, which plays a critical role in fixing CO2 through the process of in the Calvin cycle. In this study, a comprehensive exploration of the FBA family genes in moso bamboo (Phyllostachys edulis) was conducted by the bioinformatics and biological analyses. A total of nine FBA genes (PeFBA1-PeFBA9) were identified in the moso bamboo genome. The expression patterns of PeFBAs across diverse tissues of moso bamboo suggested that they have multifaceted functionality. Notably, PeFBA8 might play an important role in regulating photosynthetic carbon metabolism. Co-expression and cis-element analyses demonstrated that PeFBA8 was regulated by a photosynthetic regulatory transcription factor (PeGLK1), which was confirmed by yeast one-hybrid and dual-luciferase assays. In-planta gene editing analysis revealed that the edited PeFBA8 mutants displayed compromised photosynthetic functionality, characterized by reduced electron transport rate and impaired photosystem I, leading to decreased photosynthesis rate overall, compared to the unedited control. The recombinant protein of PeFBA8 from prokaryotic expression exhibited enzymatic catalytic function. The findings suggest that the expression of PeFBA8 can affect photosynthetic efficiency of moso bamboo leaves, which underlines the potential of leveraging PeFBA8's regulatory mechanism to breed bamboo varieties with enhanced carbon fixation capability.

Unveiling the Biological Function of Phyllostachys edulis FBA6 (PeFBA6) through the Identification of the Fructose-1,6-Bisphosphate Aldolase Gene

Fructose-1,6-bisphosphate aldolase (FBA) is a pivotal enzyme in various metabolic pathways, including glycolysis, gluconeogenesis, and the Calvin cycle. It plays a critical role in CO2 fixation. Building on previous studies on the FBA gene family in Moso bamboo, our study revealed the biological function of PeFBA6. To identify CSN5 candidate genes, this study conducted a yeast two-hybrid library screening experiment. Subsequently, the interaction between CSN5 and PeFBA6 was verified using yeast two-hybrid and LCI experiments. This investigation uncovered evidence that FBA may undergo deubiquitination to maintain glycolytic stability. To further assess the function of PeFBA6, it was overexpressed in rice. Various parameters were determined, including the light response curve, CO2 response curve, and the levels of glucose, fructose, sucrose, and starch in the leaves of overexpressing rice. The results demonstrated that overexpressed rice exhibited a higher saturation light intensity, net photosynthetic rate, maximum carboxylation rate, respiration rate, and increased levels of glucose, fructose, and starch than wild-type rice. These findings indicated that PeFBA6 not only enhanced the photoprotection ability of rice but also improved the photosynthetic carbon metabolism. Overall, this study enhanced our understanding of the function of FBA and revealed the biological function of PeFBA6, thereby providing a foundation for the development of excellent carbon fixation bamboo varieties through breeding.

Genome-wide analysis of the KNOX gene family in Moso bamboo: insights into their role in promoting the rapid shoot growth

BACKGROUND

KNOTTED1-like homeobox (KNOX) genes, plant-specific homologous box transcription factors (TFs), play a central role in regulating plant growth, development, organ formation, and response to biotic and abiotic stresses. However, a comprehensive genome-wide identification of the KNOX genes in Moso bamboo (Phyllostachys edulis), the fastest growing plant, has not yet been conducted, and the specific biological functions of this family remain unknown.

RESULTS

The expression profiles of 24 KNOX genes, divided into two subfamilies, were determined by integrating Moso bamboo genome and its transcriptional data. The KNOX gene promoters were found to contain several light and stress-related cis-acting elements. Synteny analysis revealed stronger similarity with rice KNOX genes than with Arabidopsis KNOX genes. Additionally, several conserved structural domains and motifs were identified in the KNOX proteins. The expansion of the KNOX gene family was primarily regulated by tandem duplications. Furthermore, the KNOX genes were responsive to naphthaleneacetic acid (NAA) and gibberellin (GA) hormones, exhibiting distinct temporal expression patterns in four different organs of Moso bamboo. Short Time-series Expression Miner (STEM) analysis and quantitative real-time PCR (qRT-PCR) assays demonstrated that PeKNOX genes may play a role in promoting rapid shoot growth. Additionally, Gene Ontology (GO) and Protein-Protein Interaction (PPI) network enrichment analyses revealed several functional annotations for PeKNOXs. By regulating downstream target genes, PeKNOXs are involved in the synthesis of AUX /IAA, ultimately affecting cell division and elongation.

CONCLUSIONS

In the present study, we identified and characterized a total of 24 KNOX genes in Moso bamboo and investigated their physiological properties and conserved structural domains. To understand their functional roles, we conducted an analysis of gene expression profiles using STEM and RNA-seq data. This analysis successfully revealed regulatory networks of the KNOX genes, involving both upstream and downstream genes. Furthermore, the KNOX genes are involved in the AUX/IAA metabolic pathway, which accelerates shoot growth by influencing downstream target genes. These results provide a theoretical foundation for studying the molecular mechanisms underlying the rapid growth and establish the groundwork for future research into the functions and transcriptional regulatory networks of the KNOX gene family.

Anatomical Mechanisms of Leaf Blade Morphogenesis in Sasaella kogasensis 'Aureostriatus'

There are limited studies on the cytology of bamboo leaf development from primordium to maturity. This study delves into the leaf morphological characteristics and growth patterns of Sasaella kogasensis 'Aureostriatus' and provides a three-dimensional anatomical analysis of cell division, expansion, and degradation. Leaves on the same branch develop bottom-up, while individual leaves develop the other way around. Like bamboo shoots and culms, the leaves follow a "slow-fast-slow" growth pattern, with longitudinal growth being predominant during their development. The growth zones of individual leaves included division, elongation, and maturation zones based on the distribution of growth space. By measuring 13,303 epidermal long cells and 3293 mesophyll cells in longitudinal sections of rapidly elongating leaves, we observed that in the rapid elongation phase (S4-S5), the division zone was located in the 1-2 cm segment at the bottom of the leaf blade and maintained a constant size, continuously providing new cells for leaf elongation, whereas in the late rapid elongation phase (S6), when the length of the leaf blade was approaching that of a mature leaf, its cells at the bottom of the blade no longer divided and were replaced by the ability to elongate. Furthermore, to gain an insight into the dynamic changes in the growth of the S. kogasensis 'Aureostriatus' leaves in the lateral and periclinal directions, the width and thickness of 1459 epidermal and 2719 mesophyll cells were counted in the mid-cross section of leaves at different developmental stages. The results showed that during the early stages of development (S1-S3), young leaves maintained vigorous division in the lateral direction, while periplasmic division gradually expanded from the bottom to the top of the leaf blade and the number of cell layers stabilized at S4. The meristematic tissues on both sides of the leaf were still able to divide at S4 but the frequency of the division gradually decreased, while cell division and expansion occurred simultaneously between the veins. At S6, the cells at the leaf margins and between the veins were completely differentiated and the width of the leaf blade no longer expanded. These findings revealed changes in cell growth anisotropically during the leaf development of S. kogasensis 'Aureostriatus' and demonstrated that leaf elongation was closely related to the longitudinal expansion of epidermal cells and proliferative growth of mesophyll cells, whereas the cell division of meristematic tissues and expansion of post-divisional cells contributed to the increases in blade width and thickness. The presented framework will facilitate a further exploration of the molecular regulatory mechanisms of leaf development in S. kogasensis 'Aureostriatus' and provide relevant information for developmental and taxonomic studies of bamboo plants.

Chromosomal-level genome and metabolome analyses of highly heterozygous allohexaploid Dendrocalamus brandisii elucidate shoot quality and developmental characteristics

Dendrocalamus brandisii (Munro) Kurz is a sympodial bamboo species with inimitable taste and flavorful shoots. Its rapid growth and use as high-quality material make this bamboo species highly valued for both food processing and wood applications. However, genome information for D. brandisii is lacking, primarily due to its polyploidy and large genome size. Here, we assembled a high-quality genome for hexaploid D. brandisii, which comprises 70 chromosomes with a total size of 2,756 Mb, using long-read HiFi sequencing. Furthermore, we accurately separated the genome into its three constituent subgenomes. We used Oxford Nanopore Technologies long reads to construct a transcriptomic dataset covering 15 tissues for gene annotation to complement our genome assembly, revealing differential gene expression and post-transcriptional regulation. By integrating metabolome analysis, we unveiled that well-balanced lignin formation, as well as abundant flavonoid and fructose contents, contribute to the superior quality of D. brandisii shoots. Integrating genomic, transcriptomic, and metabolomic datasets provided a solid foundation for enhancing bamboo shoot quality and developing efficient gene-editing techniques. This study should facilitate research on D. brandisii and enhance its use as a food source and wood material by providing crucial genomic resources.

Anatomical and Transcriptome Analyses of Moso Bamboo Culm Neck Growth: Unveiling Key Insights

The Moso bamboo culm neck, connected with the rhizome and the shoot bud, is an important hub for connecting and transporting the aboveground and belowground systems of bamboo for the shoot bud development and rapid growth. Our previous study revealed that the culm neck generally undergoes six different developmental stages (CNS1-CNS6), according to the primary thickening growth of the underground shoot bud. However, the molecular mechanism of the culm neck development remains unknown. The present study focused on the developmental process of the CNS3-CNS5 stages, representing the early, middle, and late elongation stages, respectively. These stages are densely packed with vascular tissues and consist of epidermis, hypodermis, cortex, and ground tissue. Unlike the hollow structure of the culms, the culm necks are solid structures. As the culm neck continues to grow, the lignin deposition increases noticeably, contributing to its progressive strengthening. For the transcriptome analysis, a total of 161,160 transcripts with an average length of 2373 were obtained from these stages using both PacBio and Illumina sequencing. A total of 92.2% of the reads mapped to the Moso bamboo reference genome. Further analysis identified a total of 5524 novel genes and revealed a dynamic transcriptome. Secondary-metabolism- and transport-related genes were upregulated particularly with the growth of the culm neck. Further analysis revealed the molecular processes of lignin accumulation in the culm neck, which include differentially expressed genes (DEGs) related to cell wall loosening and remodeling and secondary metabolism. Moreover, the upregulations of transcription factors such as MYBH and RSM in the MYB family play crucial roles during critical transitions in the culm neck development, such as changes in the angle between the rhizome and the culm neck. Our new findings provide essential insights into the cellular roadmaps, transcriptional networks, and key genes involved in the culm neck development.

Rhizosphere Microbe Affects Soil Available Nitrogen and Its Implication for the Ecological Adaptability and Rapid Growth of Dendrocalamus sinicus, the Strongest Bamboo in the World

The interaction between soil microbes and plants has a significant effect on soil microbial structure and function, as well as plant adaptability. However, the effect of soil micro-organisms on ecological adaption and rapid growth of woody bamboos remains unclear. Here, 16S rRNA and ITS rRNA genes of rhizosphere micro-organisms were sequenced, and the soil properties of three different types of Dendrocalamus sinicus were determined at the dormancy and germination stages of rhizome buds. The result showed that each type of D. sinicus preferred to absorb ammonia nitrogen (NH4+-N) rather than nitrate nitrogen (NO3--N) and required more NH4+-N at germination or rapid growth period than during the dormancy period. In total, nitrogen fixation capacity of soil bacteria in the straight type was significantly higher than that in the introduced straight type, while the ureolysis capacity had an opposite trend. Saprophytic fungi were the dominant fungal functional taxa in habitat soils of both straight and introduced straight type. Our findings are of great significance in understanding how soil microbes affect growth and adaptation of woody bamboos, but also for soil management of bamboo forests in red soil.

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.

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.

Transcriptomic Complexity of Culm Growth and Development in Different Types of Moso Bamboo

Moso bamboo is capable of both sexual and asexual reproduction during natural growth, resulting in four distinct types of culms: the bamboo shoot-culm, the seedling stem, the leptomorph rhizome, and a long-ignored culm-the outward-rhizome. Sometimes, when the outward rhizomes break through the soil, they continue to grow longitudinally and develop into a new individual. However, the roles of alternative transcription start sites (aTSS) or termination sites (aTTS) as well as alternative splicing (AS) have not been comprehensively studied for their development. To re-annotate the moso bamboo genome and identify genome-wide aTSS, aTTS, and AS in growing culms, we utilized single-molecule long-read sequencing technology. In total, 169,433 non-redundant isoforms and 14,840 new gene loci were identified. Among 1311 lncRNAs, most of which showed a positive correlation with their target mRNAs, one-third of these IncRNAs were preferentially expressed in winter bamboo shoots. In addition, the predominant AS type observed in moso bamboo was intron retention, while aTSS and aTTS events occurred more frequently than AS. Notably, most genes with AS events were also accompanied by aTSS and aTTS events. Outward rhizome growth in moso bamboo was associated with a significant increase in intron retention, possibly due to changes in the growth environment. As different types of moso bamboo culms grow and develop, a significant number of isoforms undergo changes in their conserved domains due to the regulation of aTSS, aTTS, and AS. As a result, these isoforms may play different roles than their original functions. These isoforms then performed different functions from their original roles, contributing to the transcriptomic complexity of moso bamboo. Overall, this study provided a comprehensive overview of the transcriptomic changes underlying different types of moso bamboo culm growth and development.

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.

Network of miR396-mRNA in Tissue Differentiation in Moso Bamboo (Phyllostachys edulis)

MiR396 plays an essential role in various developmental processes. However, the miR396-mRNA molecular network in bamboo vascular tissue differentiation during primary thickening has not been elucidated. Here, we revealed that three of the five members from the miR396 family were overexpressed in the underground thickening shoots collected from Moso bamboo. Furthermore, the predicted target genes were up/down-regulated in the early (S2), middle (S3) and late (S4) developmental samples. Mechanistically, we found that several of the genes encoding protein kinases (PKs), growth-regulating factors (GRF), transcription factors (TFs), and transcription regulators (TRs) were the potential targets of miR396 members. Moreover, we identified QLQ (Gln, Leu, Gln) and WRC (Trp, Arg, Cys) d omains in five PeGRF homologs and a Lipase_3 domain and a K_trans domain in another two potential targets, where the cleavage targets were identified via degradome sequencing (p < 0.05). The sequence alignment indicated many mutations in the precursor sequence of miR396d between Moso bamboo and rice. Our dual-luciferase assay revealed that ped-miR396d-5p binds to a PeGRF6 homolog. Thus, the miR396-GRF module was associated with Moso bamboo shoot development. Fluorescence in situ hybridization localized miR396 in the vascular tissues of the leaves, stems, and roots of pot Moso bamboo seedlings at the age of two months. Collectively, these experiments revealed that miR396 functions as a regulator of vascular tissue differentiation in Moso bamboo. Additionally, we propose that miR396 members are targets for bamboo improvement and breeding.

Comprehensive profiling of epigenetic modifications in fast-growing Moso bamboo shoots

The fast growth of Moso bamboo (Phyllostachys edulis) shoots is caused by the rapid elongation of each internode. However, the key underlying cellular processes and epigenetic mechanisms remain largely unexplored. We used microscopy and multi-omics approaches to investigate two regions (bottom and middle) of the 18th internode from shoots of two different heights (2 and 4 m). We observed that internode cells become longer, and that lignin biosynthesis and glycosyltransferase family 43 (GT43) genes are substantially upregulated with shoot height. Nanopore direct RNA sequencing (DRS) revealed a higher N6-methyladenine (m6A) modification rate in 2-m shoots than in 4-m shoots. In addition, different specific m6A modification sites were enriched at different growth stages. Global DNA methylation profiling indicated that DNA methylation levels are higher in 4-m shoots than in 2-m shoots. We also detected shorter poly(A) tail lengths (PALs) in 4-m shoots compared with 2-m shoots. Genes showing differential PAL were mainly enriched in the functional terms of protein translation and vesicle fusion. An association analysis between PALs and DNA methylation strongly suggested that gene body CG methylation levels are positively associated with PAL. This study provides valuable information to better understand post-transcriptional regulations responsible for fast-growing shoots in Moso 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.

Partitioning of respired CO2 in newly sprouted Moso bamboo culms

Stem respiration (R _s) plays a vital role in ecosystem carbon cycling. However, the measured efflux on the stem surface (E _s) is not always in situ R _s but only part of it. A previously proposed mass balance framework (MBF) attempted to explore the multiple partitioning pathways of R _s, including sap-flow-transported and internal storage of R _s, in addition to E _s. This study proposed stem photosynthesis as an additional partitioning pathway to the MBF. Correspondingly, a double-chamber apparatus was designed and applied on newly sprouted Moso bamboo (Phyllostachys edulis) in leafless and leaved stages. R _s of newly sprouted bamboo were twice as high in the leafless stage (7.41 ± 2.66 μmol m^-2 s^-1) than in the leaved stage (3.47 ± 2.43 μmol m^-2 s^-1). E _s accounted for ~80% of R _s, while sap flow may take away ~2% of R _s in both leafless and leaved stages. Culm photosynthesis accounted for ~9% and 13% of R _s, respectively. Carbon sequestration from culm photosynthesis accounted for approximately 2% of the aboveground bamboo biomass in the leafless stage. High culm photosynthesis but low sap flow during the leafless stage and vice versa during the leaved stage make bamboo an outstanding choice for exploring the MBF.

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.

A spatiotemporal transcriptomic network dynamically modulates stalk development in maize

Maize (Zea mays) is an important cereal crop with suitable stalk formation which is beneficial for acquiring an ideal agronomic trait to resist lodging and higher planting density. The elongation pattern of stalks arises from the variable growth of individual internodes driven by cell division and cell expansion comprising the maize stalk. However, the spatiotemporal dynamics and regulatory network of the maize stalk development and differentiation process remain unclear. Here, we report spatiotemporally resolved transcriptomes using all internodes of the whole stalks from developing maize at the elongation and maturation stages. We identified four distinct groups corresponding to four developmental zones and nine specific clusters with diverse spatiotemporal expression patterns among individual internodes of the stalk. Through weighted gene coexpression network analysis, we constructed transcriptional regulatory networks at a fine spatiotemporal resolution and uncovered key modules and candidate genes involved in internode maintenance, elongation, and division that determine stalk length and thickness in maize. Further CRISPR/Cas9-mediated knockout validated the function of a cytochrome P450 gene, ZmD1, in the regulation of stalk length and thickness as predicted by the WGCN. Collectively, these results provide insights into the high genetic complexity of stalk development and the potentially valuable resources with ideal stalk lengths and widths for genetic improvements in maize.

Genome-Wide Identification of the Highly Conserved INDETERMINATE DOMAIN (IDD) Zinc Finger Gene Family in Moso Bamboo (Phyllostachys edulis)

INDETERMINATE DOMAIN (IDD) proteins, a family of transcription factors unique to plants, function in multiple developmental processes. Although the IDD gene family has been identified in many plants, little is known about it in moso bamboo. In this present study, we identified 32 PheIDD family genes in moso bamboo and randomly sequenced the full-length open reading frames (ORFs) of ten PheIDDs. All PheIDDs shared a highly conserved IDD domain that contained two canonical C2H2-ZFs, two C2HC-ZFs, and a nuclear localization signal. Collinearity analysis showed that segmental duplication events played an important role in expansion of the PheIDD gene family. Synteny analysis indicated that 30 PheIDD genes were orthologous to those of rice (Oryza sativa). Thirty PheIDDs were expressed at low levels, and most PheIDDs exhibited characteristic organ-specific expression patterns. Despite their diverse expression patterns in response to exogenous plant hormones, 8 and 22 PheIDDs responded rapidly to IAA and 6-BA treatments, respectively. The expression levels of 23 PheIDDs were closely related to the outgrowth of aboveground branches and 20 PheIDDs were closely related to the awakening of underground dormant buds. In addition, we found that the PheIDD21 gene generated two products by alternative splicing. Both isoforms interacted with PheDELLA and PheSCL3. Furthermore, both isoforms could bind to the cis-elements of three genes (PH02Gene17121, PH02Gene35441, PH02Gene11386). Taken together, our work provides valuable information for studying the molecular breeding mechanism of lateral organ development in moso bamboo.

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.

Characterization of anatomical features, developmental roadmaps, and key genes of bamboo leaf epidermis

The exact developmental roadmaps of bamboo leaf epidermis and the regulating genes are largely unknown. In this study, we comprehensively investigated the morphological features of the leaf epidermis of bamboo, Pseudosasa japonica. We also established the developmental roadmaps of the abaxial epidermis along the linearly growing leaf. A variant of P. japonica, P. japonica var. tsutsumiana, with smaller stomata and higher stomata density, was identified. Further analysis revealed that the higher stomata density of the variant was due to the abnormal increase in stomata columns within the single stomata band. This abnormal development of stomata bands was observed as early as the guard mother cell stage in the leaf division zone (DZ). Interestingly, the developmental pattern of the single stomata was similar in P. japonica and the variant. Molecular data showed that PjDLT (Dwarf and Low Tillering) was significantly downregulated in leaves DZ of the variant. Overexpression of PjDLT in Arabidopsis and rice results in smaller plants with lower stomata density, whereas downregulation or mutation of OsDLT results in increased stomata density. Our results highlight the morphological features and developmental schedule of the leaf epidermis of bamboo and provide evidence that DLT plays an important role in regulating stomata in bamboo and rice.

Photosynthesis, Phytohormone Signaling and Sugar Catabolism in the Culm Sheaths of Phyllostachys edulis

Culm sheaths play an important role in supporting and protecting bamboo shoots during the growth and development period. The physiological and molecular functions of bamboo sheaths during the growth of bamboo shoots remain unclear. In this study, we investigated the morphological anatomy of culm sheaths, photosynthesis in sheath blades, storage and distribution of sugars, and the transcriptome of the sheath. Respiration in the base of the culm sheath was higher than that in the sheath blades; chloroplasts matured with the development of the sheath blades, the fluorescence efficiency Fv/Fm value increased from 0.3 to 0.82; and sucrose and hexose accumulated in the sheath blade and the culm sheath. The sucrose, glucose, and fructose contents of the middle sheath blades were 10.66, 5.73, and 8.84 mg/g FW, respectively. Starches accumulated in parenchymal cells close to vascular bundles. Genes related to the plant hormone signaling pathway and sugar catabolism were highly expressed in the culm sheath base. These findings provide a research basis for further understanding the possible role of bamboo sheaths in the growth and development of bamboo shoots.

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.

Identification of the 14-3-3 Gene Family in Bamboo and Characterization of Pe14-3-3b Reveals Its Potential Role in Promoting Growth

The 14-3-3 protein family plays an important role in regulating plant growth and development. The genes of the 14-3-3 family have been reported in multiple species. However, little is known about the 14-3-3 gene family in bamboo. In this study, a total of 58 genes belonging to the 14-3-3 family were identified in three representative bamboo species, i.e., Olyra latifolia, Phyllostachys edulis, and Bonia amplexicaulis, whose encoding proteins were grouped into ε and non-ε groups by phylogeny analysis with 14-3-3 proteins from Arabidopsis thaliana and Oryza sativa. The 14-3-3s had diverse gene structures and motif characteristics among the three bamboo species. Collinearity analysis suggested that the genes of the 14-3-3 family in bamboo had undergone a strong purification selection during evolution. Tissue-specific expression analysis showed the expression of Pe14-3-3s varied in different tissues of P. edulis, suggesting that they had functional diversity during growth and development. Co-expression analysis showed that four Pe14-3-3s co-expressed positively with eight ribosomal genes. Yeast two-hybrid (Y2H) assays showed that Pe14-3-3b/d could interact with Pe_ribosome-1/5/6, and qPCR results demonstrated that Pe14-3-3b/d and Pe_ribosome-1/5/6 had similar expression trends with the increase in shoot height, which further confirmed that they would work together to participate in the shoot growth and development of bamboo. Additionally, the transgenic Arabidopsis plants overexpressing Pe14-3-3b had longer roots, a larger stem diameter, an earlier bolting time and a faster growth rate than wild-type Arabidopsis, indicating that Pe14-3-3b acted as a growth promoter. Our results provide comprehensive information on 14-3-3 genes in bamboo and highlight Pe14-3-3b as a potential target for bamboo improvement.

Cellular differentiation, hormonal gradient, and molecular alternation between the division zone and the elongation zone of bamboo internodes

Bamboo differentiates a cell division zone (DZ) and a cell elongation zone (EZ) to promote internode elongation during rapid growth. However, the biological mechanisms underlying this sectioned growth behavior are still unknown. Using histological, physiological, and genomic data, we found that the cell wall and other subcellular organelles such as chloroplasts are more developed in the EZ. Abundant hydrogen peroxide accumulated in the pith cells of the EZ, and stomata formed completely in the EZ. In contrast, most cells in the DZ were in an undifferentiated state with wrinkled cell walls and dense cytoplasm. Hormone detection revealed that the levels of gibberellin, auxin, cytokinin, and brassinosteroid were higher in the DZ than in the EZ. However, the levels of salicylic acid and jasmonic acid were higher in the EZ than in the DZ. Transcriptome analysis with qRT-PCR quantification revealed that the transcripts for cell division and primary metabolism had higher expression in the DZ, whereas the genes for photosynthesis, cell wall growth, and secondary metabolism were dramatically upregulated in the EZ. Overexpression of a MYB transcription factor, BmMYB83, promotes cell wall lignification in transgenic plants. BmMYB83 is specifically expressed in cells that may have lignin deposits, such as protoxylem vessels and fiber cells. Our results indicate that hormone gradient and transcriptome reprogramming, as well as specific expression of key genes such as BmMYB83, may lead to differentiation of cell growth in the bamboo internode.

Overexpression of DsEXLA2 gene from Dendrocalamus sinicus accelerates the plant growth rate of Arabidopsis

Expansins play crucial roles in cell wall loosening and a range of life activities involving cell wall modification. Nevertheless, the biological functions of expansin genes during fast growth of bamboo remain unclear. In this study, Dendrocalamus sinicus, the largest and fastest growing bamboo species in the world, was used as the research material, and the full length of DsEXLA2 was cloned. Bioinformatics analysis revealed that DsEXLA2 contained expansin family typical domains (DPBB_1 and Pollen_allerg_1, CDRC motif) and amino acid sequence was highly conserved among different species. The expression level of DsEXLA2 increased from top section to basal section in different internodes. Subcellular localization verified that DsEXLA2 protein was located in the cell wall. Further genetic transformation studies in Arabidopsis indicated that compared with the wild type, DsEXLA2 overexpressed transgenic plants exhibited higher plant height, thicker stem, larger leaf, and less epidermal hair number and smaller stomatal aperture in the prophase and metaphase of growth. In addition, the cellulose content in the stem of transgenic plants was increased, and cell wall was thickened significantly. Moreover, a total of 1656 differentially expressed genes (DEGs) were identified by RNA-seq. The upregulated genes were predominantly enriched in the plant-pathogen interaction, MAPK signaling pathway-plant, plant hormone signal transduction, lipid metabolism and amino acid metabolism, while the downregulated genes were mainly enriched in energy metabolism, carbohydrate metabolism, plant hormone signal transduction and ribosome. These data implied that overexpression of DsEXLA2 gene accelerates the plant growth rate of Arabidopsis. This study is helpful to reveal the molecular mechanism of DsEXLA2 in culm growth and development of D. sinicus, and to understand the rapid growth of bamboos.

New Insights Into the Local Auxin Biosynthesis and Its Effects on the Rapid Growth of Moso Bamboo (Phyllostachys edulis)

Auxin plays a crucial regulatory role in higher plants, but systematic studies on the location of auxin local biosynthesis are rare in bamboo and other graminaceous plants. We studied moso bamboo (Phyllostachys edulis), which can grow up to 1 m/day and serves as a reference species for bamboo and other fast-growing species. We selected young tissues such as root tips, shoot tips, young culm sheaths, sheath blades, and internode divisions for local auxin biosynthesis site analysis. IAA immunofluorescence localization revealed that auxin was similarly distributed in different stages of 50-cm and 300-cm bamboo shoots. Shoot tips had the highest auxin content, and it may be the main site of auxin biosynthesis in the early stage of rapid growth. A total of 22 key genes in the YUCCA family for auxin biosynthesis were identified by genome-wide identification, and these had obvious tissue-specific and spatio-temporal expression patterns. In situ hybridization analysis revealed that the localization of YUCCA genes was highly consistent with the distribution of auxin. Six major auxin synthesis genes, PheYUC3-1, PheYUC6-1, PheYUC6-3, PheYUC9-1, PheYUC9-2, and PheYUC7-3, were obtained that may have regulatory roles in auxin accumulation during moso bamboo growth. Culm sheaths were found to serve as the main local sites of auxin biosynthesis and the auxin required for internode elongation may be achieved mainly by auxin transport.

Genome-Wide Identification, Expansion, and Evolution Analysis of Homeobox Gene Family Reveals TALE Genes Important for Secondary Cell Wall Biosynthesis in Moso Bamboo (Phyllostachys edulis)

The TALE gene family is a subfamily of the homeobox gene family and has been implicated in regulating plant secondary growth. However, reports about the evolutionary history and function of the TALE gene family in bamboo are limited. Here, the homeobox gene families of moso bamboo Olyra latifolia and Bonia amplexicaulis were identified and compared. Many duplication events and obvious expansions were found in the TALE family of woody bamboo. PhTALEs were found to have high syntenies with TALE genes in rice. Through gene co-expression analysis and quantitative real-time PCR analysis, the candidate PhTALEs were thought to be involved in regulating secondary cell wall development of moso bamboo during the fast-growing stage. Among these candidate PhTALEs, orthologs of OsKNAT7, OSH15, and SH5 in moso bamboo may regulate xylan synthesis by regulating the expression of IRX-like genes. These results suggested that PhTALEs may participate in the secondary cell wall deposition in internodes during the fast-growing stage of moso bamboo. The expansion of the TALE gene family may be implicated in the increased lignification of woody bamboo when divergent from herbaceous bamboos.

Total and Mitochondrial Transcriptomic and Proteomic Insights into Regulation of Bioenergetic Processes for Shoot Fast-Growth Initiation in Moso Bamboo

As a fast-growing, woody grass plant, Moso bamboo (Phyllostachys edulis) can supply edible shoots, building materials, fibrous raw material, raw materials for crafts and furniture and so on within a relatively short time. Rapid growth of Moso bamboo occurs after the young bamboo shoots are covered with a shell and emerge from the ground. However, the molecular reactions of bioenergetic processes essential for fast growth remain undefined. Herein, total and mitochondrial transcriptomes and proteomes were compared between spring and winter shoots. Numerous key genes and proteins responsible for energy metabolism were significantly upregulated in spring shoots, including those involved in starch and sucrose catabolism, glycolysis, the pentose phosphate pathway, the tricarboxylic acid cycle and oxidative phosphorylation. Accordingly, significant decreases in starch and soluble sugar, higher ATP content and higher rates of respiration and glycolysis were identified in spring shoots. Further, the upregulated genes and proteins related to mitochondrial fission significantly increased the number of mitochondria, indirectly promoting intracellular energy metabolism. Moreover, enhanced alternate-oxidase and uncoupled-protein pathways in winter shoots showed that an efficient energy-dissipating system was important for winter shoots to adapt to the low-temperature environment. Heterologous expression of PeAOX1b in Arabidopsis significantly affected seedling growth and enhanced cold-stress tolerance. Overall, this study highlights the power of comparing total and mitochondrial omics and integrating physiochemical data to understand how bamboo initiates fast growth through modulating bioenergetic processes.

Cellular and molecular characterizations of the irregular internode division zone formation of a slow-growing bamboo variant

The key molecular mechanisms underlying the sectionalized growth within bamboo or other grass internodes remain largely unknown. Here, we genetically and morphologically compared the culm and rhizome internode division zones (DZs) of a slow-growing bamboo variant (sgv) having dwarf internodes, with those of the corresponding wild type (WT). Histological analysis discovers that the sgv has an irregular internode DZ. However, the shoot apical meristems in height, width, outside shape, cell number and cell width of the sgv and the WT were all similar. The DZ irregularities first appeared post apical meristem development, in 1-mm sgv rhizome internodes. Thus, the sgv is a DZ irregularity bamboo variant, which has been first reported in bamboo according to our investigation. Transcriptome sequencing analysis finds that a number of cell wall biogenesis and cell division-related genes are dramatically downregulated in the sgv DZ. Interestingly, both transcriptomic and brassinosteroid (BR) contents detection, as well as quantitative real-time PCR analyses show that these irregularities have resulted from the BR signaling pathway defects. Brassinosteroid defect might also cause the erect leaves and branches as well as the irregular epidermis of the sgv. These results suggest that BR signaling pathway plays critical roles in bamboo internode DZ and leaf development from a mutant perspective and also explain the upstream mechanisms causing the dwarf internode of the sgv bamboo.

Allele-aware chromosome-scale assembly of the allopolyploid genome of hexaploid Ma bamboo (Dendrocalamus latiflorus Munro)

Dendrocalamus latiflorus Munro is a woody clumping bamboo with rapid shoot growth. Both genetic transformation and clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9) gene editing techniques are available for D. latiflorus, enabling reverse genetic approaches. Thus, D. latiflorus has the potential to be a model bamboo species. However, the genome sequence of D. latiflorus has remained unreported due to its polyploidy and large genome size. Here, we sequenced the D. latiflorus genome and assembled it into three allele-aware subgenomes (AABBCC), representing the largest genome of a major bamboo species. We assembled 70 allelic chromosomes (2, 737 Mb) for hexaploid D. latiflorus using both single-molecule sequencing from the Pacific Biosciences (PacBio) Sequel platform and chromosome conformation capture sequencing (Hi-C). Repetitive sequences comprised 52.65% of the D. latiflorus genome. We annotated 135 231 protein-coding genes in the genome based on transcriptomes from eight different tissues. Transcriptome sequencing using RNA-Seq and PacBio single-molecule real-time long-read isoform sequencing revealed highly differential alternative splicing (AS) between non-abortive and abortive shoots, suggesting that AS regulates the abortion rate of bamboo shoots. This high-quality hexaploid genome and comprehensive strand-specific transcriptome datasets for this Poaceae family member will pave the way for bamboo research using D. latiflorus as a model species.

Anatomical Characteristics and Variation Mechanisms on the Thick-Walled and Dwarfed Culm of Shidu Bamboo (Phyllostachys nidularia f. farcta)

Stable culm variants are valuable and important material for the study of culm development in bamboo plants. However, to date, there are few reports on the mechanism of variation of these bamboo variants. Phyllostachys nidularia f. farcta (Shidu bamboo) is a bamboo variant with stable phenotypes such as a dwarf culm with a thickened wall. In this study, we systematically investigated the cytological characteristics and underlying mechanism of morphological variation in culms of this variant using anatomical, mathematical statistical, physiological, and genomic methods. The anatomical observation and statistical results showed that the lateral increase of ground tissue in the inner layer of culm wall and the enlargement of vascular bundles are the anatomical essence of the wall thickening of Shidu bamboo; the limited elongation of fiber cells and the decrease in the number of parenchyma cells longitudinally are probably the main causes of the shortening of its internodes. A number of genes involved in the gibberellin synthesis pathway and in the synthesis of cell wall components are differentially expressed between the variant and its prototype, Ph. nidularia, and may play an important role in determining the phenotype of internode shortening in Shidu bamboo. The decrease in gibberellin content and the content of the major chemical components of the cell wall of Shidu bamboo confirmed the results of the above transcriptome. In addition, the variation in culm morphology in Shidu bamboo had little effect on the volume of the culm wall of individual internodes, suggesting that the decrease in the total number of internodes and the decrease in dry matter content (lignin, cellulose, etc.) may be the main factor for the sharp decline in culm biomass of Shidu bamboo.

The Effect of Guttation on the Growth of Bamboo Shoots

Guttation is the process of exudating droplets from the tips, edges, and adaxial and abaxial surfaces of the undamaged leaves. Guttation is a natural and spontaneous biological phenomenon that occurs in a wide variety of plants. Despite its generally positive effect on plant growth, many aspects of this cryptic process are unknown. In this study, the guttation phenomenon characteristic of bamboo shoots and the anatomical feature of these and culm sheaths were systematically observed. In addition, the water transport pathway and the compounds in guttation droplets of bamboo shoots were analyzed, and the effect of bamboo sheaths’ guttation on the growth of bamboo shoots was assessed. The results revealed that bamboo shoots began to exudate liquid in the evening through to the next morning, during which period the volume of guttation liquid gradually increases and then decreases before sunrise. Many vascular bundles are in bamboo shoots and culm sheaths to facilitate this water transport. The exudate liquid contains organic acids, sugars, and hormones, among other compounds. Our findings suggest that the regular guttation of the sheath blade is crucial to maintain the normal growth of bamboo shoots.

A regulatory network driving shoot lignification in rapidly growing bamboo

Woody bamboo is environmentally friendly, abundant, and an alternative to conventional timber. Degree of lignification and lignin content and deposition affect timber properties. However, the lignification regulatory network in monocots is poorly understood. To elucidate the regulatory mechanism of lignification in moso bamboo (Phyllostachys edulis), we conducted integrated analyses using transcriptome, small RNA, and degradome sequencing followed by experimental verification. The lignification degree and lignin content increased with increased bamboo shoot height, whereas phenylalanine ammonia-lyase and Laccase activities first increased and then decreased with shoot growth. Moreover, we identified 11,504 differentially expressed genes (DEGs) in different portions of the 13th internodes of different height shoots; most DEGs associated with cell wall and lignin biosynthesis were upregulated, whereas some DEGs related to cell growth were downregulated. We identified a total of 1,502 miRNAs, of which 687 were differentially expressed. Additionally, in silico and degradome analyses indicated that 5,756 genes were targeted by 691 miRNAs. We constructed a regulatory network of lignification, including 11 miRNAs, 22 transcription factors, and 36 enzyme genes, in moso bamboo. Furthermore, PeLAC20 overexpression increased lignin content in transgenic Arabidopsis (Arabidopsis thaliana) plants. Finally, we proposed a reliable miRNA-mediated "MYB-PeLAC20" module for lignin monomer polymerization. Our findings provide definite insights into the genetic regulation of bamboo lignification. In addition to providing a platform for understanding related mechanisms in other monocots, these insights could be used to develop strategies to improve bamboo timber properties.

New Genes Interacted with Recent Whole Genome Duplicates in the Fast Stem Growth of Bamboos

As drivers of evolutionary innovations, new genes allow organisms to explore new niches. However, clear examples of this process remain scarce. Bamboos, the unique grass lineage diversifying into the forest, have evolved with a key innovation of fast growth of woody stem, reaching up to 1 m/day. Here, we identify 1,622 bamboo-specific orphan genes that appeared in recent 46 million years, and 19 of them evolved from noncoding ancestral sequences with entire de novo origination process reconstructed. The new genes evolved gradually in exon−intron structure, protein length, expression specificity, and evolutionary constraint. These new genes, whether or not from de novo origination, are dominantly expressed in the rapidly developing shoots, and make transcriptomes of shoots the youngest among various bamboo tissues, rather than reproductive tissue in other plants. Additionally, the particularity of bamboo shoots has also been shaped by recent whole-genome duplicates (WGDs), which evolved divergent expression patterns from ancestral states. New genes and WGDs have been evolutionarily recruited into coexpression networks to underline fast-growing trait of bamboo shoot. Our study highlights the importance of interactions between new genes and genome duplicates in generating morphological innovation.

Integrative lncRNA landscape reveals lncRNA-coding gene networks in the secondary cell wall biosynthesis pathway of moso bamboo (Phyllostachys edulis)

Background

LncRNAs are extensively involved in plant biological processes. However, the lack of a comprehensive lncRNA landscape in moso bamboo has hindered the molecular study of lncRNAs. Moreover, the role of lncRNAs in secondary cell wall (SCW) biosynthesis of moso bamboo is elusive.

Results

For comprehensively identifying lncRNA throughout moso bamboo genome, we collected 231 RNA-Seq datasets, 1 Iso-Seq dataset, and 1 full-length cDNA dataset. We used a machine learning approach to improve the pipeline of lncRNA identification and functional annotation based on previous studies and identified 37,009 lncRNAs in moso bamboo. Then, we established a network of potential lncRNA-coding gene for SCW biosynthesis and identified SCW-related lncRNAs. We also proposed that a mechanism exists in bamboo to direct phenylpropanoid intermediates to lignin or flavonoids biosynthesis through the PAL/4CL/C4H genes. In addition, we identified 4 flavonoids and 1 lignin-preferred genes in the PAL/4CL/C4H gene families, which gained implications in molecular breeding.

Conclusions

We provided a comprehensive landscape of lncRNAs in moso bamboo. Through analyses, we identified SCW-related lncRNAs and improved our understanding of lignin and flavonoids biosynthesis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07953-z.

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.

Characterization of the development dynamics within the linear growth bamboo leaf

The leaf is the main photosynthetic organ in plants, such as bamboo. Leaves from bamboo are used as a food additive. However, according to our investigation, to date there have been no reports concerning the leaf development of bamboo. By measuring over 7500 parenchymal cells, we discovered that the linear leaf growth of Pseudosasa japonica could be divided into three growth sections. The first one is a roughly 1-cm long division zone (DZ), containing about 1580 cells, located at the bottom of the leaf; the second one is an about 3-cm long elongation zone (EZ), with ~1905 cells, located above the DZ; and the last is a mature zone (MZ) in which cell elongation is completed. The cell production rate of the DZ was 25.33-35.81 cells per hour, with an average of 29.73; and the cell division rate was ~0.45 cells per cell every 24 h. PacBio and Illumina transcriptome sequencing found 21 933 unigenes expressed in these zones. Further analysis revealed a dynamic transcriptome, with transcripts for cell division in the DZ changing to transcripts for cell elongation, photosynthetic development, secondary metabolism, stress resistance, and nutrition transport toward the leaf distal. Those transcriptome transformations correlated well with the changes of relative water content, biomass accumulation, and cellulose crystal degree and were supported by quantitative polymerase chain reaction data. These results revealed a developmental gradient of the bamboo linear growth leaf, which offers a foundation to elucidate and engineer leaf development in bamboo, an economically valuable plant.

Refinement of bamboo genome annotations through integrative analyses of transcriptomic and epigenomic data

Bamboo, one of the most crucial nontimber forest resources worldwide, has the capacity for rapid growth. In recent years, the genome of moso bamboo (Phyllostachys edulis) has been decoded, and a large amount of transcriptome data has been published. In this study, we generated the genome-wide profiles of the histone modification H3K4me3 in leaf, stem, and root tissues of bamboo. The trends in the distribution patterns were similar to those in rice. We developed a processing pipeline for predicting novel transcripts to refine the structural annotation of the genome using H3K4me3 ChIP-seq data and 29 RNA-seq datasets. As a result, 12,460 novel transcripts were predicted in the bamboo genome. Compared with the transcripts in the newly released version 2.0 of the bamboo genome, these novel transcripts are tissue-specific and shorter, and most have a single exon. Some representative novel transcripts were validated by semiquantitative RT-PCR and qRT-PCR analyses. Furthermore, we put these novel transcripts back into the ChIP-seq analysis pipeline and discovered that the percentages of H3K4me3 in genic elements were increased. Overall, this work integrated transcriptomic data and epigenomic data to refine the annotation of the genome in order to discover more functional genes and study bamboo growth and development, and the application of this predicted pipeline may help refine the structural annotation of the genome in other species.

Genome-wide identification and functional characterization of the PheE2F/DP gene family in Moso bamboo

BACKGROUND

E2F/DP proteins have been shown to regulate genes implicated in cell cycle control and DNA repair. However, to date, research into the potential role of the Moso bamboo E2F/DP family has been limited.

RESULTS

Here, we identified 23 E2F/DPs in the Moso bamboo genome, including nine E2F genes, six DP genes, eight DEL genes and one gene with a partial E2F domain. An estimation of the divergence time of the paralogous gene pairs suggested that the E2F/DP family expansion primarily occurred through a whole-genome duplication event. A regulatory element and coexpression network analysis indicated that E2F/DP regulated the expression of cell cycle-related genes. A yeast two-hybrid assay and expression analysis based on transcriptome data and in situ hybridization indicated that the PheE2F-PheDP complex played important roles in winter Moso bamboo shoot growth. The qRT-PCR results showed that the PheE2F/DPs exhibited diverse expression patterns in response to drought and salt treatment and diurnal cycles.

CONCLUSION

Our findings provide novel insights into the Moso bamboo E2F/DP family and partial experimental evidence for further functional verification of the PheE2F/DPs.

Multi-analysis of sheath senescence provides new insights into bamboo shoot development at the fast growth stage

Sheath senescence is an important part of bamboo shoot development during the fast growth stage. However, no information has been reported about this distinctive process until now. Using multiple approaches, we found that sheath senescence is a complex process that occurs sequentially with chloroplast corruption, chlorophyll degradation and water loss. Reactive oxygen species (ROS), salicylic acid and abscisic acid also accumulate in the senescing sheath. Transcriptome analysis showed that NAC and WRKY transcription factors, such as NAC2 and WRKY75, as well as their possible downstream target genes, such as those involved in ROS production, proteolysis and nutrition recycling, constitute the gene network of the bamboo sheath senescence process. Furthermore, the initiation of sheath senescence might be triggered by hexokinase genes, such as HXK6, which is localized to the mitochondrion and could promote leaf senescence when overexpressed in Arabidopsis. Sheath senescence occurs after the growth decrease of the internodes, which provides assimilates. The slowing of internode growth possibly results in sugar accumulation, such as glucose, in the sheath, which finally upregulates hexokinase genes and initiates sheath senescence. These findings reveal that sheath senescence is a multilevel regulation process and has a close link to the corresponding internode growth, which provides new insights into the shoot development of bamboo during the fast growth stage.

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.

The acid invertase gene family is involved in internode elongation in Phyllostachys heterocycla cv. pubescens

Acid invertases (INVs) play a pivotal role in both vegetative and reproductive growth of plants. However, their possible functions in fast-growing plants such as bamboo are largely unknown. Here, we report the molecular characterization of acid INVs in Phyllostachys heterocycla cv. pubescens, a fast-growing bamboo species commercially grown worldwide. Nine acid INVs (PhINVs), including seven cell wall INVs (PhCWINV1, PhCWINV2, PhCWINV3, PhCWINV4, PhCWINV5, PhCWINV6 and PhCWINV7) and two vacuolar INVs (PhVINV11 and PhVINV12) were isolated. Bioinformatic analyses demonstrated that they all share high amino acid identity with other INVs from different plant species and contain the motifs typically conserved in acid INV. Enzyme activity assays revealed a significantly higher INV activity in the fast-growing tissues, such as the elongating internodes of stems. Detailed quantitative reverse-transcription PCR analyses showed various expression patterns of PhINVs at different developmental stages of the elongating stems. With the exception of PhCWINV6, all PhINVs were ubiquitously expressed in a developmental-specific manner. Further studies in Arabidopsis exhibited that constitutive expression of PhCWINV1, PhCWINV4 or PhCWINV7 increased the biomass production of transgenic plants, as indicated by augmented plant heights and shoot dry weights than the wild-type plants. All these results suggest that acid INVs play a crucial role in the internode elongation of P. heterocycla cv. pubescens and would provide valuable information for the dissection of their exact biological functions in the fast growth of bamboo.

Physiological Mechanism of Internode Bending Growth After the Excision of Shoot Sheath in Fargesia yunnanensis and Its Implications for Understanding the Rapid Growth of Bamboos

The physiological function of bamboo shoot sheaths is still unclear. In the present study, we investigated the anatomical and physiological influences of bamboo shoot sheaths on internode elongation by longitudinally striping parts of sheaths. The internodes would bend toward the bare sides during night. The results showed that amounts of water leaked at the cut of shoot sheaths during night, which impeded the increase of water, water pressure and assimilate transport rates, and decreased starch and soluble sugar catabolism in the bare side of the internodes. A higher level of water pressure and sugar metabolism increased the vacuole expansion and promoted the cell expansion in the outer sides as compared to the bare sides. The bending growth of internodes was mainly due to the significant differences in cell expansion, which was led by the difference in water pressure and sugar hydrolysis levels between the inner and outer sides. Bamboo internode elongation mainly relied on the increase of water pressure and soluble sugar concentration. Shoot sheaths played an important role in the rapid growth of bamboo shoots as a controller in water and assimilate transportation. This study gave a new insight into understanding the rapid growth mechanism of bamboo plants.