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

Phyllostachys edulis argonaute genes function in the shoot architecture

Argonaute (AGO) proteins are the core components of the RNA-induced silencing complexes (RISC) in the cytoplasm and nucleus, and are necessary for the development of plant shoot meristem, which gives rise to the above-ground plant body. In this study, we identified 23 Phyllostachys edulis AGO genes (PhAGOs) that were distributed unequally on the 14 unmapped scaffolds. Gene collinearity and phylogeny analysis showed that the innovation of PhAGO genes was mainly due to dispersed duplication and whole-genome duplication, which resulted in the enlarged PhAGO family. PhAGO genes were expressed in a temporal-spatial expression pattern, and they encoded proteins differently localized in the cytoplasm and/or nucleus. Overexpression of the PhAGO2 and PhAGO4 genes increased the number of tillers or leaves in Oryza sativa and affected the shoot architecture of Arabidopsis thaliana. These results provided insight into the fact that PhAGO genes play important roles in plant development.

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.

Evolutionary relationship of moso bamboo forms and a multihormone regulatory cascade involving culm shape variation

Moso bamboo (Phyllostachys edulis) known as Mao Zhu (MZ) in Chinese exhibits various forms with distinct morphological characteristics. However, the evolutionary relationship among MZ forms and the mechanisms of culm shape variation are still lacking. Here, the main differences among MZ forms were identified as culm shape variation, which were confirmed by analysing MZ forms (799 bamboo culms) and MZ (458 bamboo culms) populations. To unravel the genetic basis underlying the morphological variations, 20 MZ forms were subjected to whole-genome resequencing. Further analysis yielded 3 230 107 high-quality SNPs and uncovered low genetic diversity and high genotype heterozygosity associated with MZ forms' formation. By integrating the SNP data of 427 MZ individuals representing 15 geographic regions, the origins of eight MZ forms were successfully traced using the phylogenetic tree and the identified common heterozygous loci. Meanwhile, transcriptomic analysis was performed using shoots from MZ and its two forms with culm shape variation. The results, combined with genomic analyses, demonstrated that hormone signalling related genes played crucial roles in culm variation. Co-expression network analysis uncovered genes associated with multiple plant hormone signal transduction, especially auxin and cytokinin were involved in culm shape variation. Furthermore, the regulatory relationships of a specific transcription factor and their target genes associated with auxin and ethylene signalling were validated by yeast one-hybrid, electrophoretic mobility shift assays, and dual-luciferase reporter. Overall, this study provides important insights into the culm shape variation formation in bamboo, which facilitates to breed new varieties with novel culms.

Variations and trade-offs in leaf and culm functional traits among 77 woody bamboo species

BACKGROUND

Woody bamboos are the only diverse large perennial grasses in mesic-wet forests and are widely distributed in the understory and canopy. The functional trait variations and trade-offs in this taxon remain unclear due to woody bamboo syndromes (represented by lignified culm of composed internodes and nodes). Here, we examined the effects of heritable legacy and occurrence site climates on functional trait variations in leaf and culm across 77 woody bamboo species in a common garden. We explored the trade-offs among leaf functional traits, the connection between leaf nitrogen (N), phosphorus (P) concentrations and functional niche traits, and the correlation of functional traits between leaves and culms.

RESULTS

The Bayesian mixed models reveal that the combined effects of heritable legacy (phylogenetic distances and other evolutionary processes) and occurrence site climates accounted for 55.10-90.89% of the total variation among species for each studied trait. The standardized major axis analysis identified trade-offs among leaf functional traits in woody bamboo consistent with the global leaf economics spectrum; however, compared to non-bamboo species, the woody bamboo exhibited lower leaf mass per area but higher N, P concentrations and assimilation, dark respiration rates. The canonical correlation analysis demonstrated a positive correlation (ρ = 0.57, P-value < 0.001) between leaf N, P concentrations and morphophysiology traits. The phylogenetic principal components and trait network analyses indicated that leaf and culm traits were clustered separately, with leaf assimilation and respiration rates associated with culm ground diameter.

CONCLUSION

Our study confirms the applicability of the leaf economics spectrum and the biogeochemical niche in woody bamboo taxa, improves the understanding of woody bamboo leaf and culm functional trait variations and trade-offs, and broadens the taxonomic units considered in plant functional trait studies, which contributes to our comprehensive understanding of terrestrial forest ecosystems.

Genome assemblies of 11 bamboo species highlight diversification induced by dynamic subgenome dominance

Polyploidy (genome duplication) is a pivotal force in evolution. However, the interactions between parental genomes in a polyploid nucleus, frequently involving subgenome dominance, are poorly understood. Here we showcase analyses of a bamboo system (Poaceae: Bambusoideae) comprising a series of lineages from diploid (herbaceous) to tetraploid and hexaploid (woody), with 11 chromosome-level de novo genome assemblies and 476 transcriptome samples. We find that woody bamboo subgenomes exhibit stunning karyotype stability, with parallel subgenome dominance in the two tetraploid clades and a gradual shift of dominance in the hexaploid clade. Allopolyploidization and subgenome dominance have shaped the evolution of tree-like lignified culms, rapid growth and synchronous flowering characteristic of woody bamboos as large grasses. Our work provides insights into genome dominance in a remarkable polyploid system, including its dependence on genomic context and its ability to switch which subgenomes are dominant over evolutionary time.

Xylan-directed cell wall assembly in grasses

Xylan is the most abundant hemicellulosic polysaccharide in the cell walls of grasses and is pivotal for the assembly of distinct cell wall structures that govern various cellular functions. Xylan also plays a crucial role in regulating biomass recalcitrance, ultimately affecting the utilization potential of lignocellulosic materials. Over the past decades, our understanding of the xylan biosynthetic machinery and cell wall organization has substantially improved due to the innovative application of multiple state-of-the-art techniques. Notably, novel xylan-based nanostructures have been revealed in the cell walls of xylem vessels, promoting a more extensive exploration of the role of xylan in the formation of cell wall structures. This Update summarizes recent achievements in understanding xylan biosynthesis, modification, modeling, and compartmentalization in grasses, providing a brief overview of cell wall assembly regarding xylan. We also discuss the potential for tailoring xylan to facilitate the breeding of elite energy and feed crops.

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 identification and expression analysis of the HAK/KUP/KT gene family in Moso bamboo

The K+ uptake permease/high-affinity K+/K+ transporter (KUP/HAK/KT) family is the most prominent group of potassium (K+) transporters, playing a key role in K+ uptake, transport, plant growth and development, and stress tolerance. However, the presence and functions of the KUP/HAK/KT family in Moso bamboo (Phyllostachys edulis (Carriere) J. Houzeau), the fastest-growing plant, have not been studied. In this study, we identified 41 KUP/HAK/KT genes (PeHAKs) distributed across 18 chromosomal scaffolds of the Moso bamboo genome. PeHAK is a typical membrane protein with a conserved structural domain and motifs. Phylogenetic tree analysis classified PeHAKs into four distinct clusters, while collinearity analysis revealed gene duplications resulting from purifying selection, including both tandem and segmental duplications. Enrichment analysis of promoter cis-acting elements suggested their plausible role in abiotic stress response and hormone induction. Transcriptomic data and STEM analyses indicated that PeHAKs were involved in tissue and organ development, rapid growth, and responded to different abiotic stress conditions. Subcellular localization analysis demonstrated that PeHAKs are predominantly expressed at the cell membrane. In-situ PCR experiments confirmed that PeHAK was mainly expressed in the lateral root primordia. Furthermore, the involvement of PeHAKs in potassium ion transport was confirmed by studying the potassium ion transport properties of a yeast mutant. Additionally, through homology modeling, we revealed the structural properties of HAK as a transmembrane protein associated with potassium ion transport. This research provides a solid basis for understanding the classification, characterization, and functional analysis of the PeHAK family in Moso bamboo.

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.

Low phosphorus impact on Moso bamboo (Phyllostachys edulis) root morphological polymorphism and expression pattern of the related genes

Moso bamboo typically grows in phosphorus (P)-deficient soil that limits its growth and development. In this study, 10 Moso bamboo genotypes (Ph-1 to Ph-10) were evaluated for their responses to P deficiency during the seedling stage by growing them in both P-sufficient and P-deficient conditions. Adaptive responses to low P (LP) conditions were observed in the majority of genotypes. Under P deficiency conditions, the total biomass decreased in several genotypes, but at the same time, the root-to-shoot ratio increased. Principal component analysis identified two main comprehensive traits (PC1 and PC2) related to the root volume and surface area and P concentration and accumulation. Based on the analysis, two genotypes (Ph-6 and Ph-10) were identified with significantly different levels of tolerance to P deficiency. The results revealed that the genotype Ph-10 responded to P deficiency by significantly increasing the root surface area and volume, while simultaneously reducing the number of root cortex cells when compared with the genotype Ph-6, which showed the lowest tolerance (intolerant). The genotype Ph-10 exhibited a robust response to external LP conditions, marked by elevated expression levels of PHOSPHATE TRANSPORTERs and SYG1/PHO81/XPR1s. In situ Polymerase Chain Reaction (PCR) analysis also revealed distinct tissue-specific expression patterns of the genes in the roots, particularly highlighting the differences between Ph-6 and Ph-10. The results provide a foundation for elucidating the mechanism of LP tolerance, thus potentially contributing to developing high P-use efficiency in Moso bamboo species.

Auxin response factors fine-tune lignin biosynthesis in response to mechanical bending in bamboo

Lignin contributes to plant mechanical properties during bending loads. Meanwhile, phytohormone auxin controls various plant biological processes. However, the mechanism of auxin's role in bending-induced lignin biosynthesis was unclear, especially in bamboo, celebrated for its excellent deformation stability. Here, we reported that auxin response factors (ARF) 3 and ARF6 from Moso bamboo (Phyllostachys edulis (Carrière) J. Houz) directly regulate lignin biosynthesis pathway genes, and affect lignin biosynthesis in bamboo. Auxin and lignin exhibited asymmetric distribution patterns, and auxin promoted lignin biosynthesis in response to bending loads in bamboo. Employing transcriptomic and weighted gene co-expression network analysis approach, we discovered that expression patterns of ARF3 and ARF6 strongly correlated with lignin biosynthesis genes. ARF3 and ARF6 directly bind to the promoter regions of 4-coumarate: coenzyme A ligase (4CL3, 4CL7, and 4CL9) or caffeoyl-CoA O-methyltransferase (CCoAOMT2) genes, pivotal to lignin biosynthesis, and activate their expressions. Notably, the efficacy of this binding hinges on auxin levels. Alternation of the expressions of ARF3 and ARF6 substantially altered lignin accumulation in transgenic bamboo. Collectively, our study shed light on bamboo lignification genetics. Auxin signaling could directly modulate lignin biosynthesis genes to impact plant lignin content.

Dynamics of rhizosphere microbial structure and function associated with the biennial bearing of moso bamboo

Moso bamboo (Phyllostachys edulis) is a valuable nontimber forestry product with a biennial cycle, producing abundant bamboo shoots within one year (on-year) and few shoots within the following year (off-year). Moso bamboo plants undergo clonal reproduction, resulting in similar genetic backgrounds. However, the number of moso bamboo shoots produced each year varies. Despite this variation, the impact of soil nutrients and the root microbiome on the biennial bearing of moso bamboo is poorly understood. We collected 139 soil samples and determined 14 major physicochemical properties of the rhizosphere, rhizoplane, and bulk soil in different seasons (i.e., the growing and deciduous seasons) and different years (i.e., on- and off-years). Based on 16S rRNA and metagenomic sequencing, major variations were found in the rhizospheric microbial composition during different seasons and years in the moso bamboo forest. Environmental driver analysis revealed that essential nutrients (i.e., SOC, TOC, TN, P, and NH4+) were the main drivers of the soil microbial community composition and were correlated with the on- and off-year cycles. Moreover, 19 MAGs were identified as important biomarkers that could distinguish on- and off-years. We found that both season and year influenced both the microbial community structure and functional pathways through the biosynthesis of nutrients that potentially interact with the moso bamboo growth rhythm, especially the on-year root-associated microbiome, which had a greater abundance of specific nutrients such as gibberellins and vitamin B6. This work provides a dynamic perspective of the differential responses of various on- and off-year microbial communities and enhances our understanding of bamboo soil microbiome biodiversity and stability.

Exploring N6-methyladenosine (m6A) modification in tree species: opportunities and challenges

N 6-methyladenosine (m6A) in eukaryotes is the most common and widespread internal modification in mRNA. The modification regulates mRNA stability, translation efficiency, and splicing, thereby fine-tuning gene regulation. In plants, m6A is dynamic and critical for various growth stages, embryonic development, morphogenesis, flowering, stress response, crop yield, and biomass. Although recent high-throughput sequencing approaches have enabled the rapid identification of m6A modification sites, the site-specific mechanism of this modification remains unclear in trees. In this review, we discuss the functional significance of m6A in trees under different stress conditions and discuss recent advancements in the quantification of m6A. Quantitative and functional insights into the dynamic aspect of m6A modification could assist researchers in engineering tree crops for better productivity and resistance to various stress conditions.

The scaling of elemental stoichiometry and growth rate over the course of bamboo ontogeny

Stoichiometric rules may explain the allometric scaling among biological traits and body size, a fundamental law of nature. However, testing the scaling of elemental stoichiometry and growth to size over the course of plant ontogeny is challenging. Here, we used a fast-growing bamboo species to examine how the concentrations and contents of carbon (C), nitrogen (N) and phosphorus (P), relative growth rate (G), and nutrient productivity scale with whole-plant mass (M) at the culm elongation and maturation stages. The whole-plant C content vs M and N content vs P content scaled isometrically, and the N or P content vs M scaled as a general 3/4 power function across both growth stages. The scaling exponents of G vs M and N (and P) productivity in newly grown mass vs M relationships across the whole growth stages decreased as a -1 power function. These findings reveal the previously undocumented generality of stoichiometric allometries over the course of plant ontogeny and provide new insights for understanding the origin of ubiquitous quarter-power scaling laws in the biosphere.

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.

Emerging Insights into the Roles of the Rhizome-Culm System in Bamboo Shoot Development through Analysis of Non-Structural Carbohydrate Changes

Non-structural carbohydrates (NSCs) required for bamboo shoot development, the critical stage that determines the yield of a bamboo stand, originate from the parent bamboo with the complex underground system. However, the metabolic mechanism of NSCs in the rhizome-culm system during bamboo shoot development remains unclear. In this study, we focused on the changes of NSCs in the rhizome-culm system and used anatomical, physiological, and biochemical methods to investigate the metabolism of NSCs in bamboo shoots of Phyllostachys edulis and the role of NSCs supply in the parent bamboo at different ages. The results showed that NSCs were accumulated and consumed from the bottom to the top in a bamboo shoot, which was consistent with the developmental pattern. The starch granules were stored in advance. The bamboo sheath stored starch from the dormant stage of shoot buds. The functions of culms and rhizomes showed age-dependent differences. Adult culms showed the highest capacity to provide NSCs, with more stored NSCs and higher β-amylase activity. Conversely, young culms seemed to prefer their growth, while old culms tended to store starch. Accordingly, adult rhizomes preferred sugar transport due to the lowest starch storage, lower ADP-glucose pyrophosphorylase (AGPase) activity, and higher β-amylase activity, while young and old rhizomes tended to prefer starch storage. These results provide a basis for further understanding of nutrient metabolism in bamboo stands.

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.

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.

Effects of different phosphorus levels on tiller bud development in hydroponic Phyllostachys edulis seedlings

An appropriate amount of phosphate fertilizer can improve the germination rate of bamboo buds and increase the bamboo shoot output. However, the underlying biological mechanisms of phosphate fertilizer in bamboo shoot development have not been systematically reported. Herein, the effects of low (LP, 1 μM), normal (NP, 50 μM) and high (HP, 1000 μM) phosphorus (P) on the growth and development of moso bamboo (Phyllostachys edulis) tiller buds were first investigated. Phenotypically, the seedling biomass, average number of tiller buds and bud height growth rate under the LP and HP treatments were significantly lower than those under the NP treatment. Next, the microstructure difference of tiller buds in the late development stage (S4) at three P levels was analyzed. The number of internode cells and vascular bundles were significantly lower in the LP treatments than in the NP treatments. The relative expression levels of eight P transport genes, eight hormone-related genes and four bud development genes at the tiller bud developmental stage (S2-S4) and the tiller bud re-tillering stage were analyzed with real-time polymerase chain reaction. The results showed that the expression trends for most P transport genes, hormone-related genes and bud development genes from S2 to S4 were diversified at different P levels, and the expression levels were also different at different P levels. In the tiller bud re-tillering stage, the expression levels of seven P transport genes and six hormone-related genes showed a downward trend with increasing P level. REV expression level decreased under LP and HP conditions. TB1 expression level increased under HP condition. Therefore, we conclude that P deficiency inhibits tiller bud development and re-tillering, and that P depends on the expression of REV and TB1 genes and auxin, cytokinin and strigolactones synthesis and transporter genes to mediate tiller bud development and re-tillering.

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.

Sources of carbon supporting the fast growth of developing immature moso bamboo (Phyllostachys edulis) culms: inference from carbon isotopes and anatomy

Phyllostachys edulis is a spectacularly fast-growing species that completes its height growth within 2 months after the shoot emerges without producing leaves (fast-growing period, FGP). This phase was considered heterotrophic, with the carbon necessary for the growth being transferred from the mature culms via the rhizomes, although previous studies observed key enzymes and anatomical features related to C₄-carbon fixation in developing culms. We tested whether C₄-photosynthesis or dark-CO₂ fixation through anaplerotic reactions significantly contributes to the FGP, resulting in differences in the natural abundance of δ¹³C in bulk organic matter and organic compounds. Further, pulse-¹³CO₂-labelling was performed on developing culms, either from the surface or from the internal hollow, to ascertain whether significant CO₂ fixation occurs in developing culms. δ¹³C of young shoots and developing culms were higher (-26.3 to -26.9 ‰) compared to all organs of mature bamboos (-28.4 to -30.1 ‰). Developing culms contained chlorophylls, most observed in the skin tissues. After pulse-¹³CO₂-labelling, the polar fraction extracted from the skin tissues was slightly enriched in ¹³C, and only a weak ¹³C enrichment was observed in inner tissues. Main carbon source sustaining the FGP was not assimilated by the developing culm, while a limited anaplerotic fixation of respired CO₂ cannot be excluded and is more likely than C₄-photosynthetic carbon fixation.

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.

Genome-Wide Identification and Expression Analysis of Dendrocalamus farinosus CCoAOMT Gene Family and the Role of DfCCoAOMT14 Involved in Lignin Synthesis

As the main component of plant cell walls, lignin can not only provide mechanical strength and physical defense for plants, but can also be an important indicator affecting the properties and quality of wood and bamboo. Dendrocalamus farinosus is an important economic bamboo species for both shoots and timber in southwest China, with the advantages of fast growth, high yield and slender fiber. Caffeoyl-coenzyme A-O-methyltransferase (CCoAOMT) is a key rate-limiting enzyme in the lignin biosynthesis pathway, but little is known about it in D. farinosus. Here, a total of 17 DfCCoAOMT genes were identified based on the D. farinosus whole genome. DfCCoAOMT1/14/15/16 were homologs of AtCCoAOMT1. DfCCoAOMT6/9/14/15/16 were highly expressed in stems of D. farinosus; this is consistent with the trend of lignin accumulation during bamboo shoot elongation, especially DfCCoAOMT14. The analysis of promoter cis-acting elements suggested that DfCCoAOMTs might be important for photosynthesis, ABA/MeJA responses, drought stress and lignin synthesis. We then confirmed that the expression levels of DfCCoAOMT2/5/6/8/9/14/15 were regulated by ABA/MeJA signaling. In addition, overexpression of DfCCoAOMT14 in transgenic plants significantly increased the lignin content, xylem thickness and drought resistance of plants. Our findings revealed that DfCCoAOMT14 can be a candidate gene that is involved in the drought response and lignin synthesis pathway in plants, which could contribute to the genetic improvement of many important traits in D. farinosus and other species.

Morphological and Anatomical Analysis of the Internodes of a New Dwarf Variant of Moso Bamboo, Phyllostachys edulis f. exaurita

The lack of mutants due to the long periods between flowering of bamboo plants is one of the limiting factors inhibiting research progress in the culm development of bamboo plants. In this study, a stable new dwarf variant of Phyllostachys edulis (Moso bamboo), Phyllostachys edulis f. exaurita T. G. Chen, was discovered and was characterized morphologically, anatomically, and physiologically. The height, diameter at breast height, number of internodes, length and wall thickness of internodes, length, width and number of parenchyma cells of internodes, and morphology of the wide-type (WT) and dwarf variant vascular bundles were compared. The height of the variant was only 49% that of the WT Moso bamboo. It was concluded that the decrease in internode number and length was the cause of dwarfism in P. edulis f. exaurita. The decreased internode length was the result of a decrease in cell number and cell length in the internode. In addition, the laws of change of internode length, internode thickness, cell length, and cell number differed between the WT Moso bamboo and the variant. Furthermore, lower IAA and zeatin concentrations were detected in the buds of the variant. These results suggest that P. edulis f. exaurita is a variant with inhibited primary thickening growth, which is valuable for interpretating the molecular mechanisms underlying the primary thickening growth of bamboo that are still largely unknown.

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 Characterization of PRE Genes in Moso Bamboo (Phyllostachys edulis)

Basic helix-loop-helix (bHLH)/HLH transcription factors are involved in various aspects of the growth and development of plants. Here, we identified four HLH genes, PePRE1-4, in moso bamboo plants that are homologous to Arabidopsis PRE genes. In bamboo seedlings, PePRE1/3 were found to be highly expressed in the internode and lamina joint by using quantitative RT-PCR analysis. In the elongating internode of bamboo shoots, PePRE genes are expressed at higher levels in the basal segment than in the mature top segment. Overexpression of PePREs (PePREs-OX) in Arabidopsis showed longer petioles and hypocotyls, as well as earlier flowering. PePRE1 overexpression restored the phenotype due to the deficiency of AtPRE genes caused by artificial micro-RNA. PePRE1-OX plants showed hypersensitivity to propiconazole treatment compared with the wild type. In addition, PePRE1/3 but not PePRE2/4 proteins accumulated as punctate structures in the cytosol, which was disrupted by the vesicle recycling inhibitor brefeldin A (BFA). PePRE genes have a positive function in the internode elongation of moso bamboo shoots, and overexpression of PePREs genes promotes flowering and growth in Arabidopsis. Our findings provided new insights about the fast-growing mechanism of bamboo shoots and the application of PRE genes from bamboo.

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.

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.

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.

Dynamics of stand productivity in Moso bamboo forest after strip cutting

Strip cutting can effectively reduce the cutting cost of bamboo forests and promote the transformation and upgradation of bamboo forests through mechanization and modernization. Despite the rapid accumulation of Moso bamboo biomass, the dynamics of five years changes in stand characteristics and productivity after cutting remain unclear. This is critical for formulating efficient bamboo forest management measures. In this paper, plots with an 8 m width strip cut (SC) and respective reserved belts (RB) were selected as the research object, and the traditional management forest (CK) as control. The dynamic characteristics of stand, biomass distribution pattern, and productivity change in the different treatment plots were studied for 5 years after cutting. The results showed that cutting increased the number of shoots and new bamboo, and decreased the diameter at breast height, height to crown base, and height of new bamboo (p<0.05). Cutting reduces the productivity of both SC and RB, and allocates more biomass to the bamboo leaves to capture light in SC (p<0.05). Over time, the characteristics of new bamboo in SC reached the level of CK, and the density of standing bamboo, and productivity, were higher than those in CK. However, the number and productivity of new bamboo decreased significantly in the RB (p<0.05), which reflected the density restriction effect of bamboo forest. Further analysis showed that the increase in productivity in SC and CK was mainly from Moso bamboo at II and III "du", which positively correlated with the soil contents of total nitrogen, total phosphorus, and available phosphorus. It was suggested that after three On-year restorations, the SC could reach the level of CK, however it is necessary to density manage RB from the second On-year after cutting.

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.