Chromosome-level reference genome and alternative splicing atlas of moso bamboo (Phyllostachys edulis)

Uncovering PheCLE1 and PheCLE10 Promoting Root Development Based on Genome-Wide Analysis

Moso bamboo (Phyllostachys edulis), renowned for its rapid growth, is attributed to the dynamic changes in its apical meristem. The CLAVATA3/EMBRYO SURROUNDING REGION-RELATED (CLE) family genes are known to play crucial roles in regulating meristem and organ formation in model plants, but their functions in Moso bamboo remain unclear. Here, we conducted a genome-wide identification of the CLE gene family of Moso bamboo and investigated their gene structure, chromosomal localization, evolutionary relationships, and expression patterns. A total of 11 PheCLE genes were identified, all of which contained a conserved CLE peptide core functional motif (Motif 1) at their C-termini. Based on Arabidopsis classification criteria, these genes were predominantly distributed in Groups A-C. Collinearity analysis unveiled significant synteny among CLE genes in Moso bamboo, rice, and maize, implying potential functional conservation during monocot evolution. Transcriptomic analysis showed significant expression of these genes in the apical tissues of Moso bamboo, including root tips, shoot tips, rhizome buds, and flower buds. Particularly, single-cell transcriptomic data and in situ hybridization further corroborated the heightened expression of PheCLE1 and PheCLE10 in the apical tissue of basal roots. Additionally, the overexpression of PheCLE1 and PheCLE10 in rice markedly promoted root growth. PheCLE1 and PheCLE10 were both located on the cell membrane. Furthermore, the upstream transcription factors NAC9 and NAC6 exhibited binding affinity toward the promoters of PheCLE1 and PheCLE10, thereby facilitating their transcriptional activation. In summary, this study not only systematically identified the CLE gene family in Moso bamboo for the first time but also emphasized their central roles in apical tissue development. This provides a valuable theoretical foundation for the further exploration of functional peptides and their signaling regulatory networks in bamboo species.

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.

Phylogenomic analyses reveal reticulate evolution between Neomicrocalamus and Temochloa (Poaceae: Bambusoideae)

Neomicrocalamus and Temochloa are closely related to bamboo genera. However, when considered with newly discovered and morphologically similar material from China and Vietnam, the phylogenetic relationship among these three groups was ambiguous in the analyses based on DNA regions. Here, as a means of investigating the relationships among the three bamboo groups and exploring potential sources of genomic conflicts, we present a phylogenomic examination based on the whole plastome, single-nucleotide polymorphism (SNP), and single-copy nuclear (SCN) gene datasets. Three different phylogenetic hypotheses were found. The inconsistency is attributed to the combination of incomplete lineage sorting and introgression. The origin of newly discovered bamboos is from introgressive hybridization between Temochloa liliana (which contributed 80.7% of the genome) and Neomicrocalamus prainii (19.3%), indicating that the newly discovered bamboos are closer to T. liliana in genetics. The more similar morphology and closer distribution elevation also imply a closer relationship between Temochloa and newly discovered bamboos.

Chromosome-scale genomes of commercial timber trees (Ochroma pyramidale, Mesua ferrea, and Tectona grandis)

Wood is the most important natural and endlessly renewable source of energy. Despite the ecological and economic importance of wood, many aspects of its formation have not yet been investigated. We performed chromosome-scale genome assemblies of three timber trees (Ochroma pyramidale, Mesua ferrea, and Tectona grandis) which exhibit different wood properties such as wood density, hardness, growth rate, and fiber cell wall thickness. The combination of 10X, stLFR, Hi-Fi sequencing and HiC data led us to assemble high-quality genomes evident by scaffold N50 length of 55.97 Mb (O. pyramidale), 22.37 Mb (M. ferrea) and 14.55 Mb (T. grandis) with >97% BUSCO completeness of the assemblies. A total of 35774, 24027, and 44813 protein-coding genes were identified in M. ferrea, T. grandis and O. pyramidale, respectively. The data generated in this study is anticipated to serve as a valuable genetic resource and will promote comparative genomic analyses, and it is of practical importance in gaining a further understanding of the wood properties in non-model woody species.

Genome-Wide Identification of SPX Family Genes and Functional Characterization of PeSPX6 and PeSPX-MFS2 in Response to Low Phosphorus in Phyllostachys edulis

Moso bamboo (Phyllostachys edulis) is the most widely distributed bamboo species in the subtropical regions of China. Due to the fast-growing characteristics of P. edulis, its growth requires high nutrients, including phosphorus. Previous studies have shown that SPX proteins play key roles in phosphorus signaling and homeostasis. However, the systematic identification, molecular characterization, and functional characterization of the SPX gene family have rarely been reported in P. edulis. In this study, 23 SPXs were identified and phylogenetic analysis showed that they were classified into three groups and distributed on 13 chromosomes. The analysis of conserved domains indicated that there was a high similarity between PeSPXs among SPX proteins in other species. RNA sequencing and qRT-PCR analysis indicated that PeSPX6 and PeSPX-MFS2, which were highly expressed in roots, were clearly upregulated under low phosphorus. Co-expression network analysis and a dual luciferase experiment in tobacco showed that PeWRKY6 positively regulated the PeSPX6 expression, while PeCIGR1-2, PeMYB20, PeWRKY6, and PeWRKY53 positively regulated the PeSPX-MFS2 expression. Overall, these results provide a basis for the identification of SPX genes in P. edulis and further exploration of their functions in mediating low phosphorus responses.

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.

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.

Overexpression of <italic>PvSVP1</italic>, an <italic>SVP</italic>-like gene of bamboo, causes early flowering and abnormal floral organs in <italic>Arabidopsis</italic> and rice

<p indent="0mm">Bamboo is a nontimber woody plant featuring a long vegetative stage and uncertain flowering time. Therefore, the genes belonging to flowering repressors might be essential in regulating the transition from the vegetative to reproductive stage in bamboo. The <italic>Short Vegetative Phase</italic> ( <italic>SVP</italic>) gene plays a pivotal role in floral transition and development. However, little is known about the bamboo <italic>SVP</italic> homologues. In this study, <italic>Phyllostachys violascens</italic> <italic>PvSVP1</italic> is isolated by analysis of the <italic>P</italic>. <italic>edulis</italic> transcriptome database. Phylogenetic analysis shows that <italic>PvSVP1</italic> is closely related to <italic>OsMADS55</italic> (rice <italic>SVP</italic> homolog). <italic>PvSVP1</italic> is ubiquitously expressed in various tissues, predominantly in vegetative tissues. To investigate the function of <italic>PvSVP1</italic>, <italic>PvSVP1</italic> is overexpressed in <italic>Arabidopsis</italic> and rice under the influence of the 35S promoter. Overexpression of <italic>PvSVP1</italic> in <italic>Arabidopsis</italic> causes early flowering and produces abnormal petals and sepals. Quantitative real-time PCR reveals that overexpression in <italic>Arabidopsis</italic> produces an early flowering phenotype by downregulating <italic>FLC</italic> and upregulating <italic>FT</italic> and produces abnormal floral organs by upregulating <italic>AP1</italic>, <italic>AP3</italic> and <italic>PI</italic> expressions. Simultaneously, overexpression of <italic>PvSVP1</italic> in rice alters the expressions of flowering-related genes such as <italic>Hd3a</italic>, <italic>RFT1</italic>, <italic>OsMADS56</italic> and <italic>Ghd7</italic> and promotes flowering under field conditions. In addition, PvSVP1 may be a nuclear protein which interacts with PvVRN1 and PvMADS56 on the yeast two-hybrid and BiFC systems. Our study suggests that <italic>PvSVP1</italic> may play a vital role in flowering time and development by interacting with PvVRN1 and PvMADS56 in the nucleus. Furthermore, this study paves the way toward understanding the complex flowering process of bamboo. </p>.

Comprehensive identification and expression analysis of CAMTA gene family in Phyllostachys edulis under abiotic stress

Background

Calmodulin-binding transcription factor (CAMTA) is a major transcription factor regulated by calmodulin (CaM) that plays an essential role in plant growth, development and response to biotic and abiotic stresses. The CAMTA gene family has been identified in Arabidopsis thaliana, rice (Oryza sativa) and other model plants, and its gene function in moso bamboo (Phyllostachys edulis) has not been identified.

Results

In this study, a total of 11 CAMTA genes were identified in P. edulis genome. Conserved domain and multiplex sequence alignment analysis showed that the structure between these genes was highly similar, with all members having CG-1 domains and some members having TIG and IQ domains. Phylogenetic relationship analysis showed that the CAMTA genes were divided into five subfamilies, and gene fragment replication promoted the evolution of this gene family. Promoter analysis revealed a large number of drought stress-related cis-acting elements in PeCAMTAs, and similarly high expression of the CAMTA gene family was found in drought stress response experiments, indicating the involvement of this gene family in drought stress. Gene expression pattern according to transcriptome data revealed participation of the PeCAMTA genes in tissue development.

Conclusions

Our results present new findings for the P. edulis CAMTA gene family and provide partial experimental evidence for further validation of the function of PeCAMTAs.

The identity of Sasaoblongula C.H.Hu (Poaceae, Bambusoideae, Arundinarieae): evidence from morphology and molecular data

Sasaoblongula was described in 1987 based on a cultivated plant at the bamboo garden of Sun Yat-sen University. This species has two or three branches at the upper nodes, which differ from the rest of Sasa species that have a single branch per node. During the field trip to Baishi Town, Yunfu City, Guangdong Province in July 2021, one bamboo species with oblong foliage leaves was collected and matches the isotype. Then, our question was to test the identity of S.oblongula concerning other Sasa species based on morphology and molecular data. To do that, we sequenced the whole chloroplast genome of S.oblongula and did a phylogenetic analysis. Our morphological results indicate that the new collection is S.oblongula. The phylogenetic tree showed that S.oblongula is close to Pseudosasa, instead of Sasa species. Therefore, we transferred it to the genus Pseudosasa, and a revised description of P.oblongula is provided here.

The pattern of DNA methylation alteration, and its association with the expression changes of non-coding RNAs and mRNAs in Moso bamboo under abiotic stress

Epigenetic changes play an important role in plant growth and development and in stress response. However, DNA methylation pattern and its relationship with the expression changes of non-coding RNAs and mRNAs of Moso bamboo in response to abiotic stress is still largely unknown. In this work, we used whole-genome bisulfite sequencing in combination with whole-transcriptome sequencing to analyze the DNA methylation and transcription patterns of mRNAs and non-coding RNAs in Moso bamboo under abiotic stresses such as cold, heat, ultraviolet (UV) and salinity. We found that CHH methylation in the promoter region was positively correlated with gene expression, while CHG and CHH methylations in the gene body regions were negatively associated with gene expression. Moreover, CG and CHG methylations in the promoter regions were negatively correlated with the transcript abundance of long non-coding RNAs (lncRNAs), microRNAs (miRNAs) and circular RNAs (circRNAs). Similarly, the methylation levels of three contexts in the genic regions were negatively correlated with the transcript abundance of lncRNAs and miRNAs but positively correlated with that of circRNAs. In addition, we suggested that the reduction of 21-nt and 24-nt small interfering RNA (siRNA) expression tended to increase methylation levels in the genic regions. We found that stress-responsive genes such as CRPK1, HSFB2A and CIPK were differentially methylated and expressed. Our results also proposed that DNA methylation may regulate the expression of the transcription factors (TFs) and plant hormone signalling genes such as IAA9, MYC2 and ERF110 in response to abiotic stress. This study firstly reports the abiotic stress-responsive DNA methylation pattern and its involvement of expression of coding RNAs and non-coding RNAs in Moso bamboo. The results expand the knowledge of epigenetic mechanisms in Moso bamboo under abiotic stress and support in-depth deciphering of the function of specific non-coding RNAs in future studies.

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.

Full-length fruit transcriptomes of southern highbush (Vaccinium sp.) and rabbiteye (V. virgatum Ait.) blueberry

Background

Blueberries (Vaccinium sp.) are native to North America and breeding efforts to improve blueberry fruit quality are focused on improving traits such as increased firmness, enhanced flavor and greater shelf-life. Such efforts require additional genomic resources, especially in southern highbush and rabbiteye blueberries.

Results

We generated the first full-length fruit transcriptome for the southern highbush and rabbiteye blueberry using the cultivars, Suziblue and Powderblue, respectively. The transcriptome was generated using the Pacific Biosciences single-molecule long-read isoform sequencing platform with cDNA pooled from seven stages during fruit development and postharvest storage. Raw reads were processed through the Isoseq pipeline and full-length transcripts were mapped to the ‘Draper’ genome with unmapped reads collapsed using Cogent. Finally, we identified 16,299 and 15,882 non-redundant transcripts in ‘Suziblue’ and ‘Powderblue’ respectively by combining the reads mapped to Northern Highbush blueberry ‘Draper’ genome and Cogent analysis. In both cultivars, > 80% of sequences were longer than 1,000 nt, with the median transcript length around 1,700 nt. Functionally annotated transcripts using Blast2GO were > 92% in both ‘Suziblue’ and ‘Powderblue’ with overall equal distribution of gene ontology (GO) terms in the two cultivars. Analyses of alternative splicing events indicated that around 40% non-redundant sequences exhibited more than one isoform. Additionally, long non-coding RNAs were predicted to represent 5.6% and 7% of the transcriptomes in ‘Suziblue’ and ‘Powderblue’, respectively. Fruit ripening is regulated by several hormone-related genes and transcription factors. Among transcripts associated with phytohormone metabolism/signaling, the highest number of transcripts were related to abscisic acid (ABA) and auxin metabolism followed by those for brassinosteroid, jasmonic acid and ethylene metabolism. Among transcription factor-associated transcripts, those belonging to ripening-related APETALA2/ethylene-responsive element-binding factor (AP2/ERF), NAC (NAM, ATAF1/2 and CUC2), leucine zipper (HB-zip), basic helix-loop-helix (bHLH), MYB (v-MYB, discovered in avian myeloblastosis virus genome) and MADS-Box gene families, were abundant.

Further we measured three fruit ripening quality traits and indicators [ABA, and anthocyanin concentration, and texture] during fruit development and ripening. ABA concentration increased during the initial stages of fruit ripening and then declined at the Ripe stage, whereas anthocyanin content increased during the final stages of fruit ripening in both cultivars. Fruit firmness declined during ripening in ‘Powderblue’. Genes associated with the above parameters were identified using the full-length transcriptome. Transcript abundance patterns of these genes were consistent with changes in the fruit ripening and quality-related characteristics.

Conclusions

A full-length, well-annotated fruit transcriptome was generated for two blueberry species commonly cultivated in the southeastern United States. The robustness of the transcriptome was verified by the identification and expression analyses of multiple fruit ripening and quality–regulating genes. The full-length transcriptome is a valuable addition to the blueberry genomic resources and will aid in further improving the annotation. It will also provide a useful resource for the investigation of molecular aspects of ripening and postharvest processes.

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.

GGDB: A Grameneae genome alignment database of homologous genes hierarchically related to evolutionary events

The genomes of Gramineae plants have been preferentially sequenced owing to their economic value. These genomes are often quite complex, for example harboring many duplicated genes, and are the main source of genetic innovation and often the result of recurrent polyploidization. Deciphering these complex genome structures and linking duplicated genes to specific polyploidization events are important for understanding the biology and evolution of plants. However, efforts have been hampered by the complexity of analyzing these genomes. Here, we analyzed 29 well-assembled and up-to-date Gramineae genome sequences by hierarchically relating duplicated genes in collinear regions to specific polyploidization or speciation events. We separated duplicated genes produced by each event, established lists of paralogous and orthologous genes, and ultimately constructed an online database, GGDB (http://www.grassgenome.com/). Homologous gene lists from each plant and between plants can be displayed, searched, and downloaded from the database. Interactive comparison tools are deployed to demonstrate homology among user-selected plants and to draw genome-scale or local alignment figures and gene-based phylogenetic trees corrected by exploiting gene collinearity. Using these tools and figures, users can easily detect structural changes in genomes and explore the effects of paleo-polyploidy on crop genome structure and function. The GGDB will provide a useful platform for improving our understanding of genome changes and functional innovation in Gramineae plants.

Comprehensive Analyses of Simple Sequence Repeat (SSR) in Bamboo Genomes and Development of SSR Markers with Peroxidase Genes

Simple sequence repeats (SSRs) are one of the most important molecular markers, which are widespread in plants. Bamboos are important forest resources worldwide. Here, the comprehensive identification and comparative analysis of SSRs were performed in three woody and two herbaceous bamboo species. Altogether 567,175 perfect SSRs and 71,141 compound SSRs were identified from 5737.8 Mb genome sequences of five bamboo species. Di-nucleotide SSRs were the most predominant type, with an average of ~50,152.2 per species. Most SSRs were located in intergenic regions, while those located in genic regions were relatively less. Moreover, the results of annotation distribution indicated that terms with P450, peroxidase and ATP-binding cassette transporter related to lignin biosynthesis might play important roles in woody and herbaceous bamboos under the mediation of SSRs. Furthermore, the peroxidase gene family consisted of a large number of genes containing SSRs was selected for the evolutionary relationship analysis and SSR markers development. Fifteen SSR markers derived from peroxidase family genes of Phyllostachys edulis were identified as polymorphic in 34 accessions belonging to seven genera in Bambusoideae. These results provided a comprehensive insight of SSR markers into bamboo genomes, which would facilitate bamboo research related to comparative genomics, evolution and marker-assisted selection.

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.

Study of genetic variation and its association with tensile strength among bamboo species through whole genome resequencing

Moso bamboo (Phyllostachys edulis) is a versatile plant species that is widely used as a construction material by many low-income countries due to the lack of major construction materials such as steel and reinforced concrete. It is also widely used in China. Bamboo is an economically sustainable material that behaves exceptionally in natural disasters such as earthquakes and it can offer viable solutions for contemporary engineering challenges. Despite bamboo's potential in the engineering sector, biological features such as its long generation time, its large genome size, and its polyploidy are constraining factors for genetic and genomic studies that potentially can assist the breeding efforts. This study re-sequenced 8 Phyllostachys species and 18 natural accessions of Ph. edulis, generating a large set of functionally annotated molecular markers (Single Nucleotide Polymorphisms (SNPs) and InDels) providing key genomic resource information. Moreover, all this genomic information was used to carry out a preliminary genome-wide association analysis and several candidate genes were identified to be correlated with a mechanical property that is of high interest to structural engineers: its tensile strength normal to its fibers (i.e., splitting).

The Original Form of C4-Photosynthetic Phosphoenolpyruvate Carboxylase Is Retained in Pooids but Lost in Rice

Poaceae is the most prominent monocot family that contains the primary cereal crops wheat, rice, and maize. These cereal species exhibit physiological diversity, such as different photosynthetic systems and environmental stress tolerance. Phosphoenolpyruvate carboxylase (PEPC) in Poaceae is encoded by a small multigene family and plays a central role in C4-photosynthesis and dicarboxylic acid metabolism. Here, to better understand the molecular basis of the cereal species diversity, we analyzed the PEPC gene family in wheat together with other grass species. We could designate seven plant-type and one bacterial-type grass PEPC groups, ppc1a, ppc1b, ppc2a, ppc2b, ppc3, ppc4, ppcC4, and ppc-b, respectively, among which ppc1b is an uncharacterized type of PEPC. Evolutionary inference revealed that these PEPCs were derived from five types of ancient PEPCs (ppc1, ppc2, ppc3, ppc4, and ppc-b) in three chromosomal blocks of the ancestral Poaceae genome. C4-photosynthetic PEPC (ppcC4 ) had evolved from ppc1b, which seemed to be arisen by a chromosomal duplication event. We observed that ppc1b was lost in many Oryza species but preserved in Pooideae after natural selection. In silico analysis of cereal RNA-Seq data highlighted the preferential expression of ppc1b in upper ground organs, selective up-regulation of ppc1b under osmotic stress conditions, and nitrogen response of ppc1b. Characterization of wheat ppc1b showed high levels of gene expression in young leaves, transcriptional responses under nitrogen and abiotic stress, and the presence of a Dof1 binding site, similar to ppcC4 in maize. Our results indicate the evolving status of Poaceae PEPCs and suggest the functional association of ppc1-derivatives with adaptation to environmental changes.

PheGRF4e initiated auxin signaling during moso bamboo shoot development

BACKGROUND

As a ubiquitous acid-regulating protein family in eukaryotes, general regulatory factors (GRFs) are active in various life activities of plants. However, detailed investigations of the GRFs gene family in moso bamboo are scarce.

METHODS AND RESULTS

Genome-wide characteristics of the GRF gene family in moso bamboo were analyzed using the moso bamboo genome. GRF phylogeny, gene structure, conserved domains, cis-element promoters, and gene expression were systematically analyzed. A total of 20 GRF gene family members were identified in the moso bamboo genome. These genes were divided into ε and non-ε groups. qRT-PCR (real-time quantitative reverse transcription polymerase chain reaction) showed that PheGRF genes responded to auxin and gibberellin treatment. To further study PheGRF gene functions, a yeast two-hybrid experiment was performed and verified by a bimolecular fluorescence complementation experiment. The results showed that PheGRF4e could interact with PheIAA30 (auxin/indole-3-acetic acid, an Aux/IAA family gene), and both were found to act mainly on the root tip meristem and vascular bundle cells of developing shoots by in situ hybridization assay.

CONCLUSIONS

This study revealed that PheGRF genes were involved in hormone response during moso bamboo shoot development, and the possible regulatory functions of PheGRF genes were enriched by the fact that PheGRF4e initiated auxin signaling by binding to PheIAA30.

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

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

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

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

Identification of TCP family in moso bamboo (Phyllostachys edulis) and salt tolerance analysis of PheTCP9 in transgenic Arabidopsis

Bioinformatic analysis of moso bamboo TEOSINTE BRANCHED 1, CYCLOIDEA, and PROLIFERATING CELL FACTORS (TCP) transcription factors reveals their conservation and variation as well as the probable biological functions in abiotic stress response. Overexpressing PheTCP9 in Arabidopsis thaliana illustrates it may exhibit a new vision in different aspects of response to salt stress. Plant specific TCPs play important roles in plant growth, development and stress response, but studies of TCP in moso bamboo are limited. Therefore, in this study, a total of 40 TCP genes (PheTCP1 ~ 40) were identified and characterized from moso bamboo genome and divided into three different subfamilies, namely, 7 in TEOSINTE BRANCHED 1 / CYCLOIDEA (TB1/CYC), 14 in CINCINNATA (CIN) and 19 in PROLIFERATING CELL FACTOR (PCF). Subsequently, we analyzed the gene structures and conserved domain of these genes and found that the members from the same subfamilies exhibited similar exon/intron distribution patterns. Selection pressure and gene duplication analysis results indicated that PheTCP genes underwent strong purification selection during evolution. There were many cis-elements related to phytohermone and stress responsive existing in the upstream promoter regions of PheTCP genes, such as ABRE, CGTCA-motif and ARE. Subcellular localization experiments showed that PheTCP9 was a nuclear localized protein. As shown by β-glucuronidase (GUS) activity, the promoter of PheTCP9 was significantly indicated by salt stress. PheTCP9 was significantly induced in the roots, stems and leaves of moso bamboo. It was also significantly induced by NaCl solution. Overexpressing PheTCP9 increased the salt tolerance of transgenic Arabidopsis. Meanwhile, H₂O₂ and malondialdehyde (MDA) contents were significantly lower in PheTCP9 over expression (OE) transgenic Arabidopsis than WT. Catalase (CAT) activity, K⁺/Na⁺ ratio as well as CAT2 expression level was also much improved in transgenic Arabidopsis than WT under salt conditions. In addition, PheTCP9 OE transgenic Arabidopsis held higher survival rates of seedlings than WT under NaCl conditions. These results showed the positive regulation functions of PheTCP9 in plants under salt conditions.

An Integrated Regulatory Network of mRNAs, microRNAs, and lncRNAs Involved in Nitrogen Metabolism of Moso Bamboo

Nitrogen is a key macronutrient essential for plant growth and development, and its availability has a strong influence on biological processes. Nitrogen fertilizer has been widely applied in bamboo forests in recent decades; however, the mechanism of nitrogen metabolism in bamboo is not fully elucidated. Here, we characterized the morphological, physiological, and transcriptome changes of moso bamboo in response to different schemes for nitrogen addition to illuminate the regulation mechanism of nitrogen metabolism. The appropriate addition of nitrogen improved the chlorophyll content and Pn (net photosynthetic rate) of leaves, the nitrogen and ammonium contents of the seedling roots, the biomass of the whole seedling, the number of lateral roots, and the activity of enzymes involved in nitrogen metabolism in the roots. Based on the whole transcriptome data of the roots, a total of 8,632 differentially expressed mRNAs (DEGs) were identified under different nitrogen additions, such as 52 nitrate transporter genes, 6 nitrate reductase genes, 2 nitrite reductase genes, 2 glutamine synthase genes, 2 glutamate synthase genes (GOGAT), 3 glutamate dehydrogenase genes, and 431 TFs belonging to 23 families. Meanwhile, 123 differentially expressed miRNAs (DEMs) and 396 differentially expressed lncRNAs (DELs) were characterized as nitrogen responsive, respectively. Furthermore, 94 DEM-DEG pairs and 23 DEL-DEG pairs involved in nitrogen metabolism were identified. Finally, a predicted regulatory network of nitrogen metabolism was initially constructed, which included 17 nitrogen metabolic pathway genes, 15 TFs, 4 miRNAs, and 10 lncRNAs by conjoint analysis of DEGs, DEMs, and DELs and their regulatory relationships, which was supported by RNA-seq data and qPCR results. The lncRNA-miRNA-mRNA network provides new insights into the regulation mechanism of nitrogen metabolism in bamboo, which facilitates further genetic improvement for bamboo to adapt to the fluctuating nitrogen environment.

Genome-wide identification and expression characterization of the DoG gene family of moso bamboo (Phyllostachys edulis)

BACKGROUND

The DoG (Delay of Germination1) family plays a key regulatory role in seed dormancy and germination. However, to date, there is no complete genomic overview of the DoG gene family of any economically valuable crop, including moso bamboo (Phyllostachys edulis), and no studies have been conducted to characterize its expression profile. To identify the DoG gene members of moso bamboo (PeDoG) and to investigate their family structural features and tissue expression profile characteristics, a study was conducted. Based on the whole genome and differential transcriptome data, in this investigation, we have scrutinized the physicochemical properties, gene structure, cis-acting elements, phylogenetic relationships, conserved structural (CS) domains, CS motifs and expression patterns of the PeDoG1 family of moso bamboo.

RESULTS

The DoG family genes of moso bamboo were found distributed across 16 chromosomal scaffolds with 24 members. All members were found to carry DoG1 structural domains, while 23 members additionally possessed basic leucine zipper (bZIP) structural domains. We could divide the PeDoG genes into three subfamilies based on phylogenetic relationships. Covariance analysis revealed that tandem duplication was the main driver of amplification of the PeDoG genes. The upstream promoter of these genes containing several cis-acting elements indicates a plausible role in abiotic stress and hormone induction. Gene expression pattern according to transcriptome data revealed participation of the PeDoG genes in tissue and organ development. Analysis using Short Time-series Expression Miner (STEM) tool revealed that the PeDoG gene family is also associated with rapid early shoot growth. Gene ontology (GO) and KEGG analyses showed a dual role of the PeDoG genes. We found that PeDoGs has a possible role as bZIP transcription factors by regulating Polar like1 (PL1) gene expression, and thereby playing a disease response role in moso bamboo. Quantitative gene expression of the PeDoG genes revealed that they were abundantly expressed in roots and leaves, and could be induced in response to gibberellin (GA).

CONCLUSION

In this study, we found that the PeDoG genes are involved in a wide range of activities such as growth and development, stress response and transcription. This forms the first report of PeDoG genes and their potential roles in moso bamboo.

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

Primary thickening determines bamboo yield and wood property. However, little is known about the regulatory networks involved in this process. This study identified a total of 58,652 genes and 150 miRNAs via transcriptome and small RNA sequencing using the underground thickening shoot samples of wild-type (WT) Moso bamboo (Phyllostachys edulis) and a thick wall (TW) variant (P. edulis "Pachyloen") at five developmental stages (WTS1/TWS1-WTS5/TWS5). A total of 14,029 (65.17%) differentially expressed genes and 68 (45.33%) differentially expressed miRNAs were identified from the WT, TW, and WTTW groups. The first two groups were composed of four pairwise combinations, each between two successive stages (WTS2/TWS2_versus_WTS1/TWS1, WTS3/TWS3_versus_WTS2/TWS2, WTS4/TWS4_versus_WTS3/TWS3, and WTS5/TWS5_versus_WTS4/TWS4), and the WTTW group was composed of five combinations, each between two relative stages (TWS1-5_versus_WTS1-5). Additionally, among the phytohormones, zeatin showed more remarkable changes in concentrations than indole-3-acetic acid, gibberellic acid, and abscisic acid throughout the five stages in the WT and the TW groups. Moreover, 125 cleavage sites were identified for 387 miRNA-mRNA pairs via degradome sequencing (P < 0.05). The dual-luciferase reporter assay confirmed that 13 miRNAs bound to 12 targets. Fluorescence in situ hybridization localized miR166 and miR160 in the shoot apical meristem and the procambium of Moso bamboo shoots at the S1 stage. Thus, primary thickening is a complex process regulated by miRNA-gene-phytohormone networks, and the miRNAome and transcriptome dynamics regulate phenotypic plasticity. These findings provide insights into the molecular mechanisms underlying wood formation and properties and propose targets for bamboo breeding.

Identification and Characterization of circRNAs under Drought Stress in Moso Bamboo (Phyllostachys edulis)

Circular RNAs (circRNAs) are a class of endogenous noncoding RNAs formed by 3′-5′ ligation during splicing. They play an important role in the regulation of transcription and miRNA in eukaryotes. Drought is one of the detrimental abiotic stresses that limit plant growth and productivity. How circRNAs influence the response to drought stress in moso bamboo (Phyllostachys edulis) remains elusive. In this study, we investigate the expression pattern of circRNAs in moso bamboo at 6 h, 12 h, 24 h and 48 h after drought treatment by deep sequencing and bioinformatics analysis and identify 4931 circRNAs, 52 of which are differentially expressed (DEcircRNAs) in drought-treated and untreated moso bamboo. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses of the host genes that generate the DEcircRNAs indcate that these DEcircRNAs are predicted to be involved in biochemical processes in response to drought, such as ubiquitin-mediated proteolysis, calcium-dependent protein kinase phosphorylation, amino acid biosynthesis and plant hormone signal transduction including abscisic acid. In addition, some circRNAs are shown to act as sponges for 291 miRNAs. Taken together, our results characterize the transcriptome profiles of circRNAs in drought responses and provide new insights into resistance breeding of moso bamboo.

An Efficient Genetic Transformation and CRISPR/Cas9-Based Genome Editing System for Moso Bamboo (Phyllostachys edulis)

Moso bamboo (Phyllostachys edulis) is the most important monopodial bamboo species worldwide. Without a genetic transformation system, it is difficult to verify the functions of genes controlling important traits and conduct molecular breeding in moso bamboo. Here, we established a plant regeneration system from immature embryos. Calli were induced on MS medium added 4-6 mg⋅L-1 2,4-dichlorophenoxyacetic acid (2,4-D) with high efficiency (>60%). A plant growth regulator combination of 0.5 mg⋅L-1 1-naphthylacetic acid (NAA), 2.0 mg⋅L-1 6-benzylaminopurine (BAP), and 3.0 mg⋅L-1 zeatin (ZT) was suitable for shoot differentiation, and the shoot induction frequency was increased to 43% after 0.5 mg⋅L-1 abscisic acid (ABA) pretreatment. An effective antibiotic screening concentration was determined by hygromycin sensitivity test. We further optimized the Agrobacterium concentration and added vacuum infiltration for infection, which improves the transient expression efficiency. A genetic transformation system was established for the first time in moso bamboo, with the transformation efficiency of approximately 5%. To optimize genome editing, two endogenous U3 small nuclear RNA (snRNA) promoters were isolated and used to drive small guide RNA (sgRNA) expression. The results showed that the PeU3.1 promoter exhibited higher efficiency, and it was used for subsequent genome editing. Finally, homozygous pds1pds2 mutants were obtained by an efficient CRISPR/Cas9 genome-editing system. These technical systems will be conducive to gene functional validation and accelerate the molecular breeding process of moso bamboo.

Genome-Wide Identification of JRL Genes in Moso Bamboo and Their Expression Profiles in Response to Multiple Hormones and Abiotic Stresses

Jacalin-related lectins (JRLs) are a new subfamily of plant lectins that has recently been recognized and plays an important role in plant growth, development, and abiotic stress response. Although moso bamboo (Phyllostachys edulis) is an economically and industrially important bamboo worldwide, there has been no systematic identification of JRLs in this species. Here, we identified 25 JRL genes in moso bamboo, and these genes are unequally distributed among 10 genome scaffolds. Phylogenetic analysis showed that the moso bamboo JRLs were clustered into four JRL subgroups: I, II, V, and VII. Numerous stress-responsive and hormone-regulated cis-elements were detected in the upstream promoter regions of the JRLs. Genome collinearity analyses showed that the JRL genes of moso bamboo are more closely related to those of Brachypodium distachyon than to those of Oryza sativa and Zea mays. Sixty-four percent of the PeJRL genes are present as segmental and tandem duplicates. qRT-PCR expression analysis showed that JRL genes in the same subgroup were significantly downregulated in response to salicylic acid (SA), abscisic acid (ABA), and methyl jasmonate (MeJA) treatments and significantly upregulated under low temperature, drought, and salt stress; they also exhibited tissue-specific expression patterns. Subcellular localization experiments revealed that PeJRL04 and PeJRL13 were localized to the cell membrane, nucleus, and cytoplasm. Three dimensional structure prediction and yeast two-hybrid assays were used to verify that PeJRL13 exists as a self-interacting homodimer in vivo. These findings provide an important reference for understanding the functions of specific moso bamboo JRL genes and for the effective selection of stress-related genes.

Genome-wide identification and evolution of WNK kinases in Bambusoideae and transcriptional profiling during abiotic stress in Phyllostachys edulis

With-no-lysine (WNK) kinases play vital roles in abiotic stress response, circadian rhythms, and regulation of flowering time in rice, Arabidopsis, and Glycine max. However, there are no previous reports of WNKs in the Bambusoideae, although genome sequences are available for diploid, tetraploid, and hexaploid bamboo species. In the present study, we identified 41 WNK genes in five bamboo species and analysed gene evolution, phylogenetic relationship, physical and chemical properties, cis-elements, and conserved motifs. We predicted the structure of PeWNK proteins of moso bamboo and determined the exposed, buried, structural and functional amino acids. Real-time qPCR analysis revealed that PeWNK5, PeWNK7, PeWNK8, and PeWNK11 genes are involved in circadian rhythms. Analysis of gene expression of different organs at different developmental stages revealed that PeWNK genes are tissue-specific. Analysis of various abiotic stress transcriptome data (drought, salt, SA, and ABA) revealed significant gene expression levels in all PeWNKs except PeWNK11. In particular, PeWNK8 and PeWNK9 were significantly down- and up-regulated, respectively, after abiotic stress treatment. A co-expression network of PeWNK genes also showed that PeWNK2, PeWNK4, PeWNK7, and PeWNK8 were co-expressed with transcriptional regulators related to abiotic stress. In conclusion, our study identified the PeWNKs of moso bamboo involved in circadian rhythms and abiotic stress response. In addition, this study serves as a guide for future functional genomic studies of the WNK genes of the Bambusoideae.

Entailing the Next-Generation Sequencing and Metabolome for Sustainable Agriculture by Improving Plant Tolerance

Crop production is a serious challenge to provide food for the 10 billion individuals forecasted to live across the globe in 2050. The scientists' emphasize establishing an equilibrium among diversity and quality of crops by enhancing yield to fulfill the increasing demand for food supply sustainably. The exploitation of genetic resources using genomics and metabolomics strategies can help generate resilient plants against stressors in the future. The innovation of the next-generation sequencing (NGS) strategies laid the foundation to unveil various plants' genetic potential and help us to understand the domestication process to unmask the genetic potential among wild-type plants to utilize for crop improvement. Nowadays, NGS is generating massive genomic resources using wild-type and domesticated plants grown under normal and harsh environments to explore the stress regulatory factors and determine the key metabolites. Improved food nutritional value is also the key to eradicating malnutrition problems around the globe, which could be attained by employing the knowledge gained through NGS and metabolomics to achieve suitability in crop yield. Advanced technologies can further enhance our understanding in defining the strategy to obtain a specific phenotype of a crop. Integration among bioinformatic tools and molecular techniques, such as marker-assisted, QTLs mapping, creation of reference genome, de novo genome assembly, pan- and/or super-pan-genomes, etc., will boost breeding programs. The current article provides sequential progress in NGS technologies, a broad application of NGS, enhancement of genetic manipulation resources, and understanding the crop response to stress by producing plant metabolites. The NGS and metabolomics utilization in generating stress-tolerant plants/crops without deteriorating a natural ecosystem is considered a sustainable way to improve agriculture production. This highlighted knowledge also provides useful research that explores the suitable resources for agriculture sustainability.

Genome-Wide Identification and Expression Analysis of the Plant U-Box Protein Gene Family in Phyllostachys edulis

The U-box gene encodes a ubiquitin ligase that contains a U-box domain. The plant U-box (PUB) protein plays an important role in the plant stress response; however, very few studies have investigated the role of these proteins in Moso bamboo (Phyllostachys edulis). Thus, more research on PUB proteins is necessary to understand the mechanisms of stress tolerance in P. edulis. In this study, we identified 121 members of the PUB family in P. edulis (PePUB), using bioinformatics based on the P. edulis V2 genome build. The U-box genes of P. edulis showed an uneven distribution among the chromosomes. Phylogenetic analysis of the U-box genes between P. edulis and Arabidopsis thaliana suggested that these genes can be classified into eight subgroups (Groups I–VIII) based on their structural and phylogenetic features. All U-box genes and the structure of their encoded proteins were identified in P. edulis. We further investigated the expression pattern of PePUB genes in different tissues, including the leaves, panicles, rhizomes, roots, and shoots. The qRT-PCR results showed that expression of three genes, PePUB15, PePUB92, and PePUB120, was upregulated at low temperatures compared to that at 25°C. The expression levels of two PePUBs, PePUB60 and PePUB120, were upregulated under drought stress. These results suggest that the PePUB genes play an important role in resistance to low temperatures and drought in P. edulis. This research provides new insight into the function, diversity, and characterization of PUB genes in P. edulis and provides a basis for understanding their biological roles and molecular mechanisms.

Long terminal repeats (LTR) and transcription factors regulate PHRE1 and PHRE2 activity in Moso bamboo under heat stress

BACKGROUND

LTR retrotransposons play a significant role in plant growth, genome evolution, and environmental stress response, but their regulatory response to heat stress remains unclear. We have investigated the activities of two LTR retrotransposons, PHRE1 and PHRE2, of moso bamboo (Phyllostachys edulis) in response to heat stress.

RESULTS

The differential overexpression of PHRE1 and PHRE2 with or without CaMV35s promoter showed enhanced expression under heat stress in transgenic plants. The transcriptional activity studies showed an increase in transposition activity and copy number among moso bamboo wild type and Arabidopsis transgenic plants under heat stress. Comparison of promoter activity in transgenic plants indicated that 5'LTR promoter activity was higher than CaMV35s promoter. Additionally, yeast one-hybrid (Y1H) system and in planta biomolecular fluorescence complementation (BiFC) assay revealed interactions of heat-dependent transcription factors (TFs) with 5'LTR sequence and direct interactions of TFs with pol and gag.

CONCLUSIONS

Our results conclude that the 5'LTR acts as a promoter and could regulate the LTR retrotransposons in moso bamboo under heat stress.

Identification and expression analysis of the glycosyltransferase GT43 family members in bamboo reveal their potential function in xylan biosynthesis during rapid growth

Background

Xylan is one of the most abundant hemicelluloses and can crosslink cellulose and lignin to increase the stability of cell walls. A number of genes encoding glycosyltransferases play vital roles in xylan biosynthesis in plants, such as those of the GT43 family. However, little is known about glycosyltransferases in bamboo, especially woody bamboo which is a good substitute for timber.

Results

A total of 17 GT43 genes (PeGT43–1 ~ PeGT43–17) were identified in the genome of moso bamboo (Phyllostachys edulis), which belong to three subfamilies with specific motifs. The phylogenetic and collinearity analyses showed that PeGT43s may have undergone gene duplication, as a result of collinearity found in 12 pairs of PeGT43s, and between 17 PeGT43s and 10 OsGT43s. A set of cis-acting elements such as hormones, abiotic stress response and MYB binding elements were found in the promoter of PeGT43s. PeGT43s were expressed differently in 26 tissues, among which the highest expression level was found in the shoots, especially in the rapid elongation zone and nodes. The genes coexpressed with PeGT43s were annotated as associated with polysaccharide metabolism and cell wall biosynthesis. qRT–PCR results showed that the coexpressed genes had similar expression patterns with a significant increase in 4.0 m shoots and a peak in 6.0 m shoots during fast growth. In addition, the xylan content and structural polysaccharide staining intensity in bamboo shoots showed a strong positive correlation with the expression of PeGT43s. Yeast one-hybrid assays demonstrated that PeMYB35 could recognize the 5′ UTR/promoter of PeGT43–5 by binding to the SMRE cis-elements.

Conclusions

PeGT43s were found to be adapted to the requirement of xylan biosynthesis during rapid cell elongation and cell wall accumulation, as evidenced by the expression profile of PeGT43s and the rate of xylan accumulation in bamboo shoots. Yeast one-hybrid analysis suggested that PeMYB35 might be involved in xylan biosynthesis by regulating the expression of PeGT43–5 by binding to its 5′ UTR/promoter. Our study provides a comprehensive understanding of PeGT43s in moso bamboo and lays a foundation for further functional analysis of PeGT43s for xylan biosynthesis during rapid growth.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08192-y.

Homo- and Hetero-Dimers of CAD Enzymes Regulate Lignification and Abiotic Stress Response in Moso Bamboo

Lignin biosynthesis enzymes form complexes for metabolic channelling during lignification and these enzymes also play an essential role in biotic and abiotic stress response. Cinnamyl alcohol dehydrogenase (CAD) is a vital enzyme that catalyses the reduction of aldehydes to alcohols, which is the final step in the lignin biosynthesis pathway. In the present study, we identified 49 CAD enzymes in five Bambusoideae species and analysed their phylogenetic relationships and conserved domains. Expression analysis of Moso bamboo PheCAD genes in several developmental tissues and stages revealed that among the PheCAD genes, PheCAD2 has the highest expression level and is expressed in many tissues and PheCAD1, PheCAD6, PheCAD8 and PheCAD12 were also expressed in most of the tissues studied. Co-expression analysis identified that the PheCAD2 positively correlates with most lignin biosynthesis enzymes, indicating that PheCAD2 might be the key enzyme involved in lignin biosynthesis. Further, more than 35% of the co-expressed genes with PheCADs were involved in biotic or abiotic stress responses. Abiotic stress transcriptomic data (SA, ABA, drought, and salt) analysis identified that PheCAD2, PheCAD3 and PheCAD5 genes were highly upregulated, confirming their involvement in abiotic stress response. Through yeast two-hybrid analysis, we found that PheCAD1, PheCAD2 and PheCAD8 form homo-dimers. Interestingly, BiFC and pull-down experiments identified that these enzymes form both homo- and hetero- dimers. These data suggest that PheCAD genes are involved in abiotic stress response and PheCAD2 might be a key lignin biosynthesis pathway enzyme. Moreover, this is the first report to show that three PheCAD enzymes form complexes and that the formation of PheCAD homo- and hetero- dimers might be tissue specific.

Identification and Expression Analyses of Invertase Genes in Moso Bamboo Reveal Their Potential Drought Stress Functions

Invertases (INVs) can irreversibly hydrolyze sucrose into fructose and glucose, which play principal roles in carbon metabolism and responses to various stresses in plants. However, little is known about the INV family in bamboos, especially their potential function in drought stress. In this study, 29 PeINVs were identified in moso bamboo (Phyllostachys edulis). They were clustered into alkaline/neutral invertase (NINV) and acid invertase (AINV) groups based on the gene structures, conserved motifs, and phylogenetic analysis results. The collinearity analysis showed nine segmental duplication pairs within PeINVs, and 25 pairs were detected between PeINVs and OsINVs. PeINVs may have undergone strong purification selection during evolution, and a variety of stress and phytohormone-related regulatory elements were found in the promoters of PeINVs. The tissue-specific expression analysis showed that PeINVs were differentially expressed in various moso bamboo tissues, which suggested that they showed functional diversity. Both the RNA-seq and quantitative real-time PCR results indicated that four PeINVs were significantly upregulated under drought stress. Co-expression network and Pearson’s correlation coefficient analyses showed that these PeINVs co-expressed positively with sugar and water transport genes (SWTGs), and the changes were consistent with sugar content. Overall, we speculate that the identified PeINVs are spatiotemporally expressed, which enables them to participate in moso bamboo growth and development. Furthermore, PeINVs, together with SWTGs, also seem to play vital roles in the response to drought stress. These results provide a comprehensive information resource for PeINVs, which will facilitate further study of the molecular mechanism underlying PeINVs involvement in the response to drought stress in moso bamboo.

Analysis of 427 genomes reveals moso bamboo population structure and genetic basis of property traits

Moso bamboo (Phyllostachys edulis) is an economically and ecologically important nontimber forestry species. Further development of this species as a sustainable bamboo resource has been hindered by a lack of population genome information. Here, we report a moso bamboo genomic variation atlas of 5.45 million single-nucleotide polymorphisms (SNPs) from whole-genome resequencing of 427 individuals covering 15 representative geographic areas. We uncover low genetic diversity, high genotype heterozygosity, and genes under balancing selection underlying moso bamboo population adaptation. We infer its demographic history with one bottleneck and its recently small population without a rebound. We define five phylogenetic groups and infer that one group probably originated by a single-origin event from East China. Finally, we conduct genome-wide association analysis of nine important property-related traits to identify candidate genes, many of which are involved in cell wall, carbohydrate metabolism, and environmental adaptation. These results provide a foundation and resources for understanding moso bamboo evolution and the genetic mechanisms of agriculturally important traits.

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.

Genome-wide identification and analysis of the heat shock transcription factor family in moso bamboo (Phyllostachys edulis)

Heat shock transcription factors (HSFs) are central elements in the regulatory network that controls plant heat stress response. They are involved in multiple transcriptional regulatory pathways and play important roles in heat stress signaling and responses to a variety of other stresses. We identified 41 members of the HSF gene family in moso bamboo, which were distributed non-uniformly across its 19 chromosomes. Phylogenetic analysis showed that the moso bamboo HSF genes could be divided into three major subfamilies; HSFs from the same subfamily shared relatively conserved gene structures and sequences and encoded similar amino acids. All HSF genes contained HSF signature domains. Subcellular localization prediction indicated that about 80% of the HSF proteins were located in the nucleus, consistent with the results of GO enrichment analysis. A large number of stress response-associated cis-regulatory elements were identified in the HSF upstream promoter sequences. Synteny analysis indicated that the HSFs in the moso bamboo genome had greater collinearity with those of rice and maize than with those of Arabidopsis and pepper. Numerous segmental duplicates were found in the moso bamboo HSF gene family. Transcriptome data indicated that the expression of a number of PeHsfs differed in response to exogenous gibberellin (GA) and naphthalene acetic acid (NAA). A number of HSF genes were highly expressed in the panicles and in young shoots, suggesting that they may have functions in reproductive growth and the early development of rapidly-growing shoots. This study provides fundamental information on members of the bamboo HSF gene family and lays a foundation for further study of their biological functions in the regulation of plant responses to adversity.

Genome-wide identification and analysis of the YABBY gene family in Moso Bamboo (Phyllostachys edulis (Carrière) J. Houz)

BACKGROUND

The YABBY gene family is a family of small zinc finger transcription factors associated with plant morphogenesis, growth, and development. In particular, it is closely related to the development of polarity in the lateral organs of plants. Despite being studied extensively in many plant species, there is little information on genome-wide characterization of this gene family in Moso bamboo.

METHODS

In the present study, we identified 16 PeYABBY genes, which were unequally distributed on 11 chromosomes, through genome-wide analysis of high-quality genome sequences of M oso bamboo by bioinformatics tools and biotechnological tools. Gene expression under hormone stress conditions was verified by quantitative real-time PCR (qRT-PCR) experiments.

RESULTS

Based on peptide sequences and similarity of exon-intron structures, we classified the PeYABBY genes into four subfamilies. Analysis of putative cis-acting elements in promoters of these genes revealed that PeYABBYs contained a large number of hormone-responsive and stress-responsive elements. Expression analysis showed that they were expressed at a high level in Moso bamboo panicles, rhizomes, and leaves. Expression patterns of putative PeYABBY genes in different organs and hormone-treated were analyzed using RNA-seq data, results showed that some PeYABBY genes were responsive to gibberellin (GA) and abscisic acid (ABA), indicating that they may play an important role in plant hormone responses. Gene Ontology (GO) analyses of YABBY proteins indicated that they may be involved in many developmental processes, particularly high level of enrichment seen in plant leaf development. In summary, our results provide a comprehensive genome-wide study of the YABBY gene family in bamboos, which could be useful for further detailed studies of the function and evolution of the YABBY genes, and to provide a fundamental basis for the study of YABBY in Gramineae for resistance to stress and hormonal stress.

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

BACKGROUND

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

RESULT

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

CONCLUSIONS

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

A bamboo leaf-specific aquaporin gene PePIP2;7 is involved in abiotic stress response

PePIP2;7, a leaf-specific aquaporin gene in bamboo, is upregulated under abiotic stresses. Overexpressing PePIP2;7 confers abiotic stresses tolerance in transgenic Arabidopsis plant and yeast. Aquaporins (AQPs) participate in the regulation of water balance in plants. However, the function of AQPs in bamboo remains unclear. Here, PePIP2;7 was identified as a leaf-specific aquaporin gene in moso bamboo based on the expression analysis of transcriptome data and PCR. In situ hybridization further indicated that PePIP2;7 was mainly expressed in mesophyll cells of mature leaves, while in immature leaves it was dominant in blade edge cells followed by mesophyll cells. Interestingly, PePIP2;7 was strongly expressed in the mesophyll cells near bulliform cells of immature leaves, suggesting that PePIP2;7 might function in water transport and contribute to leaf unfolding. The transient expression assay showed that PePIP2;7 was a plasma membrane intrinsic protein. Furthermore, PePIP2;7 was upregulated under abiotic stresses such as high light, drought, and NaCl. Compared with Col-0, transgenic Arabidopsis plants overexpressing PePIP2;7 had better seed germination rate, longer taproot length, higher SOD activity, and lower MDA content under abiotic stresses. Besides, yeasts expressing PePIP2;7 also had higher tolerance to stress compared to the control. Taken together, our results show that PePIP2;7 is leaf-specific and involved in stress response, which provides new insights into aquaporin function in bamboo.

Genome-wide identification and expression analysis of LBD transcription factor genes in Moso bamboo (Phyllostachys edulis)

BACKGROUND

Moso bamboo, the fastest growing plant on earth, is an important source for income in large areas of Asia, mainly cultivated in China. Lateral organ boundaries domain (LBD) proteins, a family of transcription factors unique to plants, are involved in multiple transcriptional regulatory pathways and play important roles in lateral organ development, pathogen response, secondary growth, and hormone response. The LBD gene family has not previously been characterized in moso bamboo (Phyllostachys edulis).

RESULTS

In this study, we identified 55 members of the LBD gene family from moso bamboo and found that they were distributed non-uniformly across its 18 chromosomes. Phylogenetic analysis showed that the moso bamboo LBD genes could be divided into two classes. LBDs from the same class share relatively conserved gene structures and sequences encoding similar amino acids. A large number of hormone response-associated cis-regulatory elements were identified in the LBD upstream promoter sequences. Synteny analysis indicated that LBDs in the moso bamboo genome showed greater collinearity with those of O. sativa (rice) and Zea mays (maize) than with those of Arabidopsis and Capsicum annuum (pepper). Numerous segmental duplicates were found in the moso bamboo LBD gene family. Gene expression profiles in four tissues showed that the LBD genes had different spatial expression patterns. qRT-PCR assays with the Short Time-series Expression Miner (STEM) temporal expression analysis demonstrated that six genes (PeLBD20, PeLBD29, PeLBD46, PeLBD10, PeLBD38, and PeLBD06) were consistently up-regulated during the rapid growth and development of bamboo shoots. In addition, 248 candidate target genes that function in a variety of pathways were identified based on consensus LBD binding motifs.

CONCLUSIONS

In the current study, we identified 55 members of the moso bamboo transcription factor LBD and characterized for the first time. Based on the short-time sequence expression software and RNA-seq data, the PeLBD gene expression was analyzed. We also investigated the functional annotation of all PeLBDs, including PPI network, GO, and KEGG enrichment based on String database. These results provide a theoretical basis and candidate genes for studying the molecular breeding mechanism of rapid growth of moso bamboo.

The SOC1-like gene BoMADS50 is associated with the flowering of Bambusa oldhamii

Bamboo is known for its edible shoots and beautiful texture and has considerable economic and ornamental value. Unique among traditional flowering plants, many bamboo plants undergo extensive synchronized flowering followed by large-scale death, seriously affecting the productivity and application of bamboo forests. To date, the molecular mechanism of bamboo flowering characteristics has remained unknown. In this study, a SUPPRESSOR OF OVEREXPRESSION OF CONSTANS1 (SOC1)-like gene, BoMADS50, was identified from Bambusa oldhamii. BoMADS50 was highly expressed in mature leaves and the floral primordium formation period during B. oldhamii flowering and overexpression of BoMADS50 caused early flowering in transgenic rice. Moreover, BoMADS50 could interact with APETALA1/FRUITFULL (AP1/FUL)-like proteins (BoMADS14-1/2, BoMADS15-1/2) in vivo, and the expression of BoMADS50 was significantly promoted by BoMADS14-1, further indicating a synergistic effect between BoMADS50 and BoAP1/FUL-like proteins in regulating B. oldhamii flowering. We also identified four additional transcripts of BoMADS50 (BoMADS50-1/2/3/4) with different nucleotide variations. Although the protein-CDS were polymorphic, they had flowering activation functions similar to those of BoMADS50. Yeast one-hybrid and transient expression assays subsequently showed that both BoMADS50 and BoMADS50-1 bind to the promoter fragment of itself and the SHORT VEGETATIVE PHASE (SVP)-like gene BoSVP, but only BoMADS50-1 can positively induce their transcription. Therefore, nucleotide variations likely endow BoMADS50-1 with strong regulatory activity. Thus, BoMADS50 and BoMADS50-1/2/3/4 are probably important positive flowering regulators in B. oldhamii. Moreover, the functional conservatism and specificity of BoMADS50 and BoMADS50-1 might be related to the synchronized and sporadic flowering characteristics of B. oldhamii.

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.

Genome-wide identification and expression analysis of dirigent-jacalin genes from plant chimeric lectins in Moso bamboo (Phyllostachys edulis)

Dirigent-jacalin (D-J) genes belong to the plant chimeric lectin family, and play vital roles in plant growth and resistance to abiotic and biotic stresses. To explore the functions of the D-J family in the growth and development of Moso bamboo (Phyllostachys edulis), their physicochemical properties, phylogenetic relationships, gene and protein structures, and expression patterns were analyzed in detail. Four putative PeD-J genes were identified in the Moso bamboo genome, and microsynteny and phylogenetic analyses indicated that they represent a new branch in the evolution of plant lectins. PeD-J proteins were found to be composed of a dirigent domain and a jacalin-related lectin domain, each of which contained two different motifs. Multiple sequence alignment and homologous modeling analysis indicated that the three-dimensional structure of the PeD-J proteins was significantly different compared to other plant lectins, primarily due to the tandem dirigent and jacalin domains. We surveyed the upstream putative promoter regions of the PeD-Js and found that they mainly contained cis-acting elements related to hormone and abiotic stress response. An analysis of the expression patterns of root, leaf, rhizome and panicle revealed that four PeD-J genes were highly expressed in the panicle, indicating that they may be required during the formation and development of several different tissue types in Moso bamboo. Moreover, PeD-J genes were shown to be involved in the rapid growth and development of bamboo shoots. Quantitative Real-time PCR (qRT PCR) assays further verified that D-J family genes were responsive to hormones and stresses. The results of this study will help to elucidate the biological functions of PeD-Js during bamboo growth, development and stress response.

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

BACKGROUND

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

RESULT

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

CONCLUSION

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