Transcriptome and miRNAome analysis reveals components regulating tissue differentiation of bamboo shoots

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.

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.

The Construction of lncRNA/circRNA-miRNA-mRNA Networks Reveals Functional Genes Related to Growth Traits in Schima superba

Schima superba is a precious timber and fire-resistant tree species widely distributed in southern China. Currently, there is little knowledge related to its growth traits, especially with respect to molecular breeding. The lack of relevant information has delayed the development of modern breeding. The purpose is to identify probable functional genes involved in S. superba growth through whole transcriptome sequencing. In this study, a total of 32,711 mRNAs, 525 miRNAs, 54,312 lncRNAs, and 1522 circRNAs were identified from 10 S. superba individuals containing different volumes of wood. Four possible regulators, comprising three lncRNAs, one circRNA, and eleven key miRNAs, were identified from the regulatory networks of lncRNA-miRNA-mRNA and circRNA-miRNA-mRNA to supply information on ncRNAs. Several candidate genes involved in phenylpropane and cellulose biosynthesis pathways, including Ss4CL2, SsCSL1, and SsCSL2, and transcription factors, including SsDELLA2 (SsSLR), SsDELLA3 (SsSLN), SsDELLA5 (SsGAI-like2), and SsNAM1, were identified to reveal the molecular regulatory mechanisms regulating the growth traits of S. superba. The results not merely provide candidate functional genes related to S. superba growth trait and will be useful to carry out molecular breeding, but the strategy and method also provide scientists with an effective approach to revealing mechanisms behind important economic traits in other species.

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.

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.

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.

Identification and characterization of CircRNA-associated CeRNA networks in moso bamboo under nitrogen stress

BACKGROUND

Nitrogen is a macronutrient element for plant growth and development. Circular RNAs (circRNAs) serve as pivotal regulators for the coordination between nutrient supply and plant demand. Moso bamboo (Phyllostachys edulis) is an excellent plant with fast growth, and the mechanism of the circRNA-target module in response to nitrogen remains unclear.

RESULTS

Deep small RNA sequencing results of moso bamboo seedlings under different concentrations of KNO₃ (N0 = 0 mM, N6 = 6 mM, N18 = 18 mM) were used to identify circRNAs. A total of 549 circRNAs were obtained, of which 309 were generated from corresponding parental coding genes including 66 new ones. A total of 536 circRNA-parent genes were unevenly distributed in 24 scaffolds and were associated with root growth and development. Furthermore, 52 differentially expressed circRNAs (DECs) were obtained, including 24, 33 and 15 DECs from three comparisons of N0 vs. N6, N0 vs. N18 and N6 vs. N18, respectively. Based on integrative analyses of the identified DECs, differentially expressed mRNAs (DEGs), and miRNAs (DEMs), a competitive endogenous RNA (ceRNA) network was constructed, including five DECs, eight DEMs and 32 DEGs. A regulatory module of PeSca_6:12,316,320|12,372,905-novel_miR156-PH02Gene35622 was further verified by qPCR and dual-luciferase reporter assays.

CONCLUSION

The results indicated that circRNAs could participate in multiple biological processes as miRNA sponges, including organ nitrogen compound biosynthesis and metabolic process regulation in moso bamboo. Our results provide valuable information for further study of circRNAs in moso bamboo under fluctuating nitrogen conditions.

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.

Development of an efficient expression system with large cargo capacity for interrogation of gene function in bamboo based on bamboo mosaic virus

Bamboo is one of the fastest growing plants among monocotyledonous species and is grown extensively in subtropical regions. Although bamboo has high economic value and produces much biomass quickly, gene functional research is hindered by the low efficiency of genetic transformation in this species. We therefore explored the potential of a bamboo mosaic virus (BaMV)-mediated expression system to investigate genotype-phenotype associations. We determined that the sites between the triple gene block proteins (TGBps) and the coat protein (CP) of BaMV are the most efficient insertion sites for the expression of exogenous genes in both monopodial and sympodial bamboo species. Moreover, we validated this system by individually overexpressing the two endogenous genes ACE1 and DEC1, which resulted in the promotion and suppression of internode elongation, respectively. In particular, this system was able to drive the expression of three 2A-linked betalain biosynthesis genes (more than 4 kb in length) to produce betalain, indicating that it has high cargo capacity and may provide the prerequisite basis for the development of a DNA-free bamboo genome editing platform in the future. Since BaMV can infect multiple bamboo species, we anticipate that the system described in this study will greatly contribute to gene function research and further promote the molecular breeding of bamboo.

Revealing the mechanisms of the bioactive ingredients accumulation in Polygonatum cyrtonema by multiomics analyses

Polygonatum cyrtonema is a medicinal and edible herb rich in polysaccharides, steroidal saponins, and flavonoids that has been widely used as a food, vegetable, and medicine over the years. Although previous studies have preliminarily explored the metabolic and transcriptional regulatory mechanisms of the main secondary metabolites in P. cyrtonema, the complex mechanism of microRNA (miRNA)-mediated posttranscriptional regulation remains unclear. Metabolome analysis showed that iso-ophiopogonanone B, (25S)-pratioside D1, disporopsin, and isodiosgenin-Glc-Glc, which are associated with intermediates in the flavonoids and saponins pathways, were significantly upregulated in the stem and leaf compared with the rhizome, and most saccharides, including arabinose, cellobiose, maltotetraose, and panose, showed the opposite trend, suggesting that they may contribute to the formation and accumulation of the main active ingredients in P. cyrtonema. We found that 4-hydroxymandelonitrile have a relatively good inhibitory effect on α-glucosidase, indicating that it may play a role in hypoglycemic functions. Transcriptome and weighted gene coexpression network analysis (WGCNA) were combined to reveal several candidate genes involved in the accumulation of polysaccharides, saponins, and flavonoids, including PcSQLE, PcCYP71A1, PcSUS, PcFK, and PcMYB102. Integrated analyses of miRNAs and messengerRNAs (mRNAs) showed that novel_miR14, novel_miR49, novel_miR75, and aof_miR164 were negatively correlated with alpha-linolenic acid metabolism and the mitogen activated protein kinase (MAPK) signaling pathway, including PcAOS, PcSPLA2, PcFRK1, and PcDELLA, indicating that these miRNAs may coordinately regulate the biosynthesis of other secondary metabolites in P. cyrtonema. These findings will facilitate in-depth research on the functions of these miRNAs and mRNAs related to the main active substances for pathological and biological regulation, which will be beneficial to provide theoretical guidance for the molecular breeding of P. cyrtonema.

Integrative analyses of morphology, physiology, and transcriptional expression profiling reveal miRNAs involved in culm color in bamboo

Culm color variation is an interesting phenomenon that contributes to the breeding of new varieties of ornamental plants during domestication. De-domesticated variation is considered ideal for identifying and interpreting the molecular mechanisms of plant mutations. However, the variation in culm color of bamboo remains unknown. In the present study, yellow and green culms generated from the same rhizome of Phyllostachys vivax cv. Aureocaulis (P. vivax) were used to elucidate the molecular mechanism of culm color formation. Phenotypic and physiological data showed that environmental suitability was higher in green culms than in yellow culms. High-throughput sequencing analysis showed 295 differentially expressed genes (DEGs) and 22 differentially expressed miRNAs (DEMs) in two different colored bamboo culms. There were 103 DEM-DEG interaction pairs, of which a representative "miRNA-mRNA" regulatory module involved in photosynthesis and pigment metabolism was formed by 14 DEM-DEG pairs. The interaction of the three key pairs was validated by qPCR and dual-luciferase assays. This study provides new insights into the molecular mechanism of miRNAs involved in P. vivax culm color formation, which provides evidence for plant de-domestication and is helpful for revealing the evolutionary mechanism of bamboo.

Transcriptome Analysis to Identify Genes Related to Flowering Reversion in Tomato

Flowering reversion is a common phenomenon in plant development in which differentiated floral organs switch from reproductive growth to vegetative growth and ultimately form abnormal floral organs or vegetative organs. This greatly reduces tomato yield and quality. Research on this phenomenon has recently increased, but there is a lack of research at the molecular and gene expression levels. Here, transcriptomic analyses of the inflorescence meristem were performed in two kinds of materials at different developmental stages, and a total of 3223 differentially expressed genes (DEGs) were screened according to the different developmental stages and trajectories of the two materials. The analysis of database annotations showed that these DEGs were closely related to starch and sucrose metabolism, DNA replication and modification, plant hormone synthesis and signal transduction. It was further speculated that tomato flowering reversion may be related to various biological processes, such as cell signal transduction, energy metabolism and protein post-transcriptional regulation. Combined with the results of previous studies, our work showed that the gene expression levels of CLE9, FA, PUCHI, UF, CLV3, LOB30, SFT, S-WOX9 and SVP were significantly different in the two materials. Endogenous hormone analysis and exogenous hormone treatment revealed a variety of plant hormones involved in flowering reversion in tomato. Thus, tomato flowering reversion was studied comprehensively by transcriptome analysis for the first time, providing new insights for the study of flower development regulation in tomato and other plants.

The miR166–mRNA network regulates vascular tissue differentiation in Moso bamboo

miR166s play an important role in plant tissue differentiation. However, the functions of miR166s in the differentiation of vascular tissue in bamboo have not yet been elucidated. Here, we showed that five miR166s are overexpressed (tags per million reads &gt; 2,000) in underground shoot samples of wild-type (WT) Moso bamboo (Phyllostachys edulis) and a thick-walled variant (P. edulis “Pachyloen”) throughout the developmental process. Potential targets of these miR166s include some genes encoding homeodomain-leucine zipper (HD-ZIP) transcription factors and protein kinases. Cleavage sites for miR166s were identified in seven PeHD-ZIP homologs and a protein kinase gene via degradome sequencing (p &lt; 0.05). Dual-luciferase and transient expression assays confirmed the binding of miR166s to PeHOXs. Fluorescence in situ hybridization revealed that miR166s were localized to the xylem of the leaf, root, and internode of 2-month-old pot seedlings of WT Moso bamboo. Overall, these findings reveal that miR166s are regulators of vascular tissue differentiation in bamboo. The miR166s identified in our study provide novel targets for bamboo breeding.