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

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

Analyzing lignin biosynthesis pathways in rattan using improved co-expression networks of NACs and MYBs.

BACKGROUND: The rattan is a valuable plant resource with multiple applications in tropical forests. Calamus simplicifolius and Daemonorops jenkinsiana are the two most representative rattan species, supplying over 95% of the raw materials for the rattan industry. Hence, the wood properties of both rattans have always attracted researchers' attention. RESULTS: We re-annotated the genomes, obtained 81 RNA-Seq datasets, and developed an improved pipeline to increase the reliability of co-expression networks of both rattans. Based on the data and pipeline, co-expression relationships were detected in 11 NACs, 49 MYBs, and 86 lignin biosynthesis genes in C. simplicifolius and four NACs, 59 MYBs, and 76 lignin biosynthesis genes in D. jenkinsiana, respectively. Among these co-expression pairs, several genes had a close relationship to the development of wood properties. Additionally, we detected the enzyme gene on the lignin biosynthesis pathway was regulated by either NAC or MYB, while LACCASES was regulated by both NAC and MYB. For D. jenkinsiana, the lignin biosynthesis regulatory network was characterized by positive regulation, and MYB possible negatively regulate non-expressed lignin biosynthesis genes in stem tissues. For C. simplicifolius, NAC may positively regulate highly expressed genes and negatively regulate non-expressed lignin biosynthesis genes in stem tissues. Furthermore, we established core regulatory networks of NAC and MYB for both rattans. CONCLUSIONS: This work improved the accuracy of rattan gene annotation by integrating an efficient co-expression network analysis pipeline, enhancing gene coverage and accuracy of the constructed network, and facilitating an understanding of co-expression relationships among NAC, MYB, and lignin biosynthesis genes in rattan and other plants.

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.

Targeted enrichment of novel chloroplast-based probes reveals a large-scale phylogeny of 412 bamboos.

BACKGROUND: The subfamily Bambusoideae belongs to the grass family Poaceae and has significant roles in culture, economy, and ecology. However, the phylogenetic relationships based on large-scale chloroplast genomes (CpGenomes) were elusive. Moreover, most of the chloroplast DNA sequencing methods cannot meet the requirements of large-scale CpGenome sequencing, which greatly limits and impedes the in-depth research of plant genetics and evolution. RESULTS: To develop a set of bamboo probes, we used 99 high-quality CpGenomes with 6 bamboo CpGenomes as representative species for the probe design, and assembled 15 M unique sequences as the final pan-chloroplast genome. A total of 180,519 probes for chloroplast DNA fragments were designed and synthesized by a novel hybridization-based targeted enrichment approach. Another 468 CpGenomes were selected as test data to verify the quality of the newly synthesized probes and the efficiency of the probes for chloroplast capture. We then successfully applied the probes to synthesize, enrich, and assemble 358 non-redundant CpGenomes of woody bamboo in China. Evaluation analysis showed the probes may be applicable to chloroplasts in Magnoliales, Pinales, Poales et al. Moreover, we reconstructed a phylogenetic tree of 412 bamboos (358 in-house and 54 published), supporting a non-monophyletic lineage of the genus Phyllostachys. Additionally, we shared our data by uploading a dataset of bamboo CpGenome into CNGB ( https://db.cngb.org/search/project/CNP0000502/ ) to enrich resources and promote the development of bamboo phylogenetics. CONCLUSIONS: The development of the CpGenome enrichment pipeline and its performance on bamboos recommended an inexpensive, high-throughput, time-saving and efficient CpGenome sequencing strategy, which can be applied to facilitate the phylogenetics analysis of most green plants.

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

KEY MESSAGE: 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 brassinosteroid action-related genes during the shoot growth of moso bamboo.

Brassinosteroids (BRs) are a group of plant steroid hormones that play crucial roles in a range of plant growth and development processes. BR action includes active BR formation by a complex biosynthesis process and driving BR biological function through signal transduction. Although the characterization of several BR action-related genes has been conducted in a few model plants, systematic information about these genes in bamboo is still lacking. We identified 64 genes related to BR action from the genome of moso bamboo (Phyllostachys edulis), including twenty that participated in BR biosynthesis and forty-four involved in BR signal transduction. The characteristics of all these candidate genes were identified by bioinformatics methods, including the gene structures, basic physical and chemical properties of proteins, conserved domains and evolutionary relationships. Based on the transcriptome data, the candidate genes demonstrated different expression patterns, which were further validated by qRT-PCR using templates from bamboo shoots with different heights. Thirty-four positive and three negative co-expression modules were identified by 44 candidate genes in the newly emerging bamboo shoot. The gene expression patterns and co-expression modules of BR action-related genes in bamboo shoots indicated that they might function to promote bamboo growth through BR biosynthesis and signal transduction processes. This study provides the first step towards the cloning and functional dissection of the role of BR action-related genes in moso bamboo, which also presents an excellent opportunity for genetic engineering using the candidate genes to improve bamboo quantity and quality.

Expression and functional analysis of two PsbS genes in bamboo (Phyllostachys edulis).

Higher plants have an array of photoprotection mechanisms alleviating the harmful effects of light. Non-photochemical quenching (NPQ) is one of the photoprotective mechanisms, which dissipates the excess of light energy absorbed in the light-harvesting complexes (LHCs) into thermal energy. The photosystem II subunit S (PsbS), a member of the LHC family thought to be present exclusively in higher plants, is supposed to activate NPQ through interactions with antenna proteins. However, the roles of PsbS in bamboo remain unclear. Here, two genes of bamboo (Phyllostachys edulis), PePsbS1 and PePsbS2, are investigated and functionally analyzed. PePsbS1 and PePsbS2 have a similar gene structure with three introns separated by two exons, which encode 269 and 268 amino acid residues, respectively. Tissue-specific analysis showed that PePsbS1 and PePsbS2 are highly expressed in leaf blade. Besides, they are both upregulated in the leaf blade when plantlets are submitted to an increased and prolonged light intensity, suggesting that they are light-induced. Western blot analysis indicated that the accumulation level of total PePsbSs is consistent with what obtained by quantitative real-time polymerase chain reaction for PePsbS1 and PePsbS2. Transgenic Arabidopsis plants overexpressing PePsbS1 and PePsbS2 both displayed an enhanced photoprotection. Moreover, the expression of PePsbS1 and PePsbS2 could both rescue the NPQ of Arabidopsis npq4 mutant, indicating that the PsbSs are functionally conserved between monocots and dicots. These results indicated that both PePsbS1 and PePsbS2 could circumvent photoinhibition and enhance photoprotection, which are key factors for bamboo's adaptation to different light environment.

The bamboo aquaporin gene PeTIP4;1-1 confers drought and salinity tolerance in transgenic Arabidopsis.

KEY MESSAGE: PeTIP4;1-1, an aquaporin gene involved in bamboo shoot growth, is regulated by abiotic stresses. Overexpression of PeTIP4;1-1 confers drought and salinity tolerance in transgenic Arabidopsis. Aquaporins play a central role in numerous physiological processes throughout plant growth and development. PeTIP4;1-1, an aquaporin gene isolated from moso bamboo (Phyllostachys edulis), comprises an open reading frame (ORF) of 756 bp encoding a peptide of 251 amino acids. The genomic sequence corresponding to the ORF of PeTIP4;1-1 was 1777 bp and contained three exons separated by two introns. PeTIP4;1-1 was constitutively expressed at the highest level in culms, and the expression level was elevated with increasing height of the bamboo shoot. PeTIP4;1-1 was significantly up-regulated in response to drought and salinity stresses in bamboo roots and leaves. To investigate the role of PeTIP4;1-1 in response to drought and salinity stresses, transgenic Arabidopsis plants overexpressing PeTIP4;1-1 under the control of CaMV 35S promoter were generated and subjected to morphological and physiological assays. Compared with Col-0, the transgenic plants showed enhanced tolerance to drought and salinity stresses and produced longer taproots, which had more green leaves, higher F v/F m and NPQ values, higher activities of SOD, POD and CAT, lower MDA concentration and higher water content. Transcript levels of three stress-related genes (AtP5CS, AtNHX1 and AtLEA) were enhanced. These results indicated that PeTIP4;1-1 might play an important function in response to drought and salinity stresses, and is a candidate gene for breeding of stress tolerance in other crops through genetic engineering.

Transcriptome and comparative gene expression analysis of Phyllostachys edulis in response to high light.

BACKGROUND: Photosynthesis plays a vital role as an energy source for plant metabolism, and its efficiency may be drastically reduced owing to abiotic stresses. Moso bamboo (Phyllostachys edulis), is a renewable and versatile resource with significant ecological and economic value, which encounters high light stress with large amplitude in natural environment. However, the gene expression profiles in response to high light were elusive in bamboo. RESULTS: We firstly performed physiological experiments on moso bamboo leaves treated with high light (1200 µmol · m(-2) · s(-1)). Based on the physiological results, three samples of leaves treated with high light for 0 h (CK), 0.5 h (0.5H), and 8 h (8H) were selected to perform further high-throughput RNA sequencing (RNA-Seq), respectively. Then, the transcriptomic result demonstrated that the most genes were expressed at a statistically significant value (FPKM ≥ 1) and the RNA-Seq data were validated via quantitative real time PCR. Moreover, some significant gene expression changes were detected. For instance, 154 differentially expressed genes were detected in 0.5H vs. CK, those in 8H vs. CK were 710, and 429 differentially expressed genes were also identified in 0.5H vs.8 H. Besides, 47 gene annotations closely related to photosynthesis were refined, including 35 genes annotated as light-harvesting chlorophyll a/b-binding (LHC) proteins, 9 LHC-like proteins and 3 PsbSs. Furthermore, the pathway of reactive oxygen species (ROS) in photosynthesis was further analyzed. A total of 171 genes associated with ROS-scavenging were identified. Some up-regulated transcript factors, such as NAC, WRKY, AR2/ERF, and bHLH, mainly concentrated in short-term response, while C2H2, HSF, bZIP, and MYB were largely involved in short and middle terms response to high light. CONCLUSION: Based on the gene expression analysis of moso bamboo in response to high light, we thus identified the global gene expression patterns, refined the annotations of LHC protein, LHC-like protein and PsbS, detected the pathway of ROS as well as identified ROS-scavenging genes and transcription factors in the regulation of photosynthetic and related metabolisms. These findings maybe provide a starting point to interpret the molecular mechanism of photosynthesis in moso bamboo under high light stress.

Genome-wide identification and characterization of aquaporin gene family in moso bamboo (Phyllostachys edulis).

Aquaporins (AQPs) are known to play a major role in maintaining water and hydraulic conductivity balance in the plant system. Numerous studies have showed AQPs execute multi-function throughout plant growth and development, including water transport, nitrogen, carbon, and micronutrient acquisition etc. However, little information on AQPs is known in bamboo. In this study, we present the first genome-wide identification and characterization of AQP genes in moso bamboo (Phyllostachys edulis) using bioinformatics. In total, 26 AQP genes were identified by homologous analysis, which were divided into four groups (PIPs, TIPs, NIPs, and SIPs) based on the phylogenetic analysis. All the genes were located on 26 different scaffolds respectively on basis of the gene mapped to bamboo genome. Evolutionary analysis indicated that Ph. edulis was more close to Oryza sativa than Zea mays in the genetic relationship. Besides, qRT-PCR was used to analyze gene expression profiles, which revealed that AQP genes were expressed constitutively in all the detected tissues, and were all responsive to the environmental cues such as drought, water, and NaCl stresses. This data suggested that AQPs may play fundamental roles in maintaining normal growth and development of bamboo, which would contribute to better understanding for the complex regulation mechanism involved in the fast-growing process of bamboo. Furthermore, the result could provide valuable information for further research on bamboo functional genomics.

Characterization and primary functional analysis of a bamboo NAC gene targeted by miR164b.

KEY MESSAGE: PeSNAC1 , a stress-related NAC1 from Phyllostachys edulis , was characterized. Ectopic expression in Arabidopsis indicated that PeSNAC1 together with ped -miR164b participated in the regulation of organ boundaries and stress tolerance. NAC (NAM, ATAF1/2 and CUC2) participates in many different processes regulating plant growth, development, and stress response. A total of 125 NAC genes have been predicted in moso bamboo (Phyllostachys edulis), but their roles are poorly understood. PeSNAC1 targeted by ped-miR164b was focused for further study. The cleavage of PeSNAC1 mRNA guided by ped-miR164b was validated using RLM-5' RACE. Tissue-specific expression analysis demonstrated that ped-miR164b had a declining trend from root, sheath, leaf, to that of stem, which was opposite to that of PeSNAC1. Transgenic Arabidopsis plants overexpressing either PeSNAC1 (OX-PeSNAC1) or, ped-miR164b (OX-ped-miR164b) driven by the CaMV35S promoter were generated. OX-ped-miR164b plants showed similar phenotype of cuc2 mutants whose growth was seriously suppressed. Compared with Col-0, sense OX-PeSNAC1 plants grew rapidly and flowered earlier, whereas antisense plants grew slowly and exhibited delayed flowering. Sense OX-PeSNAC1 plants had the greatest number of lateral roots, while antisense OX-PeSNAC1 and OX-ped-miR164b plants had fewer lateral roots than Col-0. Under NaCl and PEG6000 stresses, survival rates were higher and F v/F m values declined more slowly in sense OX-PeSNAC1 plants than in Col-0, with lower survival rates and a more rapid decrease in F v/F m values conversely observed in antisense OX-PeSNAC1 and OX-ped-miR164b plants. These findings indicated that ped-miR164b-targeted PeSNAC1 may play key roles in plant development and tolerance to salinity and drought stresses.

Molecular characterization and primary functional analysis of PeMPEC, a magnesium-protoporphyrin IX monomethyl ester cyclase gene of bamboo (Phyllostachys edulis).

KEY MESSAGE: A homologous gene of MPEC from Phyllostachys edulis was isolated and characterized. Its overexpression in Arabidopsis thaliana increased chlorophyll concentration and photosynthesis efficiency, indicating it is involved in chlorophyll biosynthesis. Magnesium-protoporphyrin IX monomethyl ester cyclase (MPEC) is an essential enzyme in the biosynthesis of chlorophyll, which plays an important role in photosynthesis. However, limited information is available on the roles of MPEC gene in bamboo. A homologous gene, PeMPEC was identified from Phyllostachys edulis, which comprised 1474 bp and contained an open reading frame encoding 415 amino acids. PeMPEC was transcribed abundantly in leaf blade where photosynthesis predominantly occurs, which agreed with the protein accumulation pattern confirmed by Western blotting. The PeMPEC transcription was promoted by continuous darkness for 24 h, and was suppressed by increasing light intensity (100-1500 µmo1 m(-2) s(-1)) and high temperature (42 °C). However, transcription was induced within 0.5 h and thereafter declined with prolonged treatment (up to 12 h) under low temperature (4 °C). Although PeMPEC expressed weakly in etiolated leaves, transcript levels increased gradually with subsequent light treatment (200 µmol m(-2) s(-1)). Overexpression of PeMPEC in Arabidopsis thaliana resulted in increased chlorophyll concentration and photosynthesis efficiency in sense transgenic plants compared with a reduction in antisense transgenic plants. These changes were consistent with the transcript levels of PeMPEC. These results indicated that PeMPEC might be involved in chlorophyll biosynthesis and play important roles in maintaining the stability of photosystems, and provide a basis for the study of chlorophyll biosynthesis and dissection of photosynthesis in bamboo.

Haplotype-based pangenomes reveal genetic variations and climate adaptations in moso bamboo populations.

Moso bamboo (Phyllostachys edulis), an ecologically and economically important forest species in East Asia, plays vital roles in carbon sequestration and climate change mitigation. However, intensifying climate change threatens moso bamboo survival. Here we generate high-quality haplotype-based pangenome assemblies for 16 representative moso bamboo accessions and integrated these assemblies with 427 previously resequenced accessions. Characterization of the haplotype-based pangenome reveals extensive genetic variation, predominantly between haplotypes rather than within accessions. Many genes with allele-specific expression patterns are implicated in climate responses. Integrating spatiotemporal climate data reveals more than 1050 variations associated with pivotal climate factors, including temperature and precipitation. Climate-associated variations enable the prediction of increased genetic risk across the northern and western regions of China under future emissions scenarios, underscoring the threats posed by rising temperatures. Our integrated haplotype-based pangenome elucidates moso bamboo's local climate adaptation mechanisms and provides critical genomic resources for addressing intensifying climate pressures on this essential bamboo. More broadly, this study demonstrates the power of long-read sequencing in dissecting adaptive traits in climate-sensitive species, advancing evolutionary knowledge to support conservation.

The Global Atlas of Bamboo and Rattan (GABR) Phase II: new resources for sustainable development.

Bamboo, the fast-growing grass plant, and rattan, the spiky climbing palm, are both essential forest resources that have been closely linked with human lives, livelihoods and material culture since ancient times. To promote genetic and genomic research in bamboo and rattan, a comprehensive and coordinated international project, the Genome Atlas of Bamboo and Rattan (GABR), was launched in 2017. GABR achieved great success during Phase I (2017-2022). We will focus on investigating and protecting bamboo and rattan germplasm resources in Phase II ( 2022-2027). Here, we briefly review the achievements of Phase I and introduce the goals of Phase II.

The flying spider-monkey tree fern genome provides insights into fern evolution and arborescence.

To date, little is known about the evolution of fern genomes, with only two small genomes published from the heterosporous Salviniales. Here we assembled the genome of Alsophila spinulosa, known as the flying spider-monkey tree fern, onto 69 pseudochromosomes. The remarkable preservation of synteny, despite resulting from an ancient whole-genome duplication over 100 million years ago, is unprecedented in plants and probably speaks to the uniqueness of tree ferns. Our detailed investigations into stem anatomy and lignin biosynthesis shed new light on the evolution of stem formation in tree ferns. We identified a phenolic compound, alsophilin, that is abundant in xylem, and we provided the molecular basis for its biosynthesis. Finally, analysis of demographic history revealed two genetic bottlenecks, resulting in rapid demographic declines of A. spinulosa. The A. spinulosa genome fills a crucial gap in the plant genomic landscape and helps elucidate many unique aspects of tree fern biology.

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.

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.

Coexpression Analysis Reveals Dynamic Modules Regulating the Growth and Development of Cirri in the Rattans (Calamus simplicifolius and Daemonorops jenkinsiana).

Rattan is regarded as one of the major non-timber forest products, second only to wood and bamboo, worldwide. Although the published genomes of Calamus simplicifolius and Daemonorops jenkinsiana have facilitated genome-wide gene functional analyses, coexpression networks (CENs) provide more comprehensive and complete annotations of gene function at the transcriptome level. Thus, we analyzed the CENs of the two rattans, C. simplicifolius and D. jenkinsiana, by integrating the genome sequences and analyzing in-house transcriptome data from different development stages of their cirri using a well-developed strategy. A total of 3,504 and 3,027 functional modules were identified in C. simplicifolius and D. jenkinsiana, respectively, based on a combination of CENs, gene family classification, and function enrichment tools. These modules covered the major developmental processes, including photosynthesis, lignin biosynthesis, flavonoid biosynthesis, and phenylpropanoid biosynthesis. Reference annotations were refined using CENs and functional modules. Moreover, we obtained novel insights into the regulation of cirrus growth and development in rattans. Furthermore, Rattan-NET (http://rattan.bamboogdb.org/), an online database with analysis tools for gene set enrichment analysis, module enrichment, network comparison analysis, and cis-element analysis, was constructed for the easy analysis of gene function and regulation modules involved in the growth and development of cirri in rattans.

Genome-wide identification and expression analysis of the MYB transcription factor in moso bamboo (Phyllostachys edulis).

The MYB family, one of the largest transcription factor (TF) families in the plant kingdom, plays vital roles in cell formation, morphogenesis and signal transduction, as well as responses to biotic and abiotic stresses. However, the underlying function of bamboo MYB TFs remains unclear. To gain insight into the status of these proteins, a total of 85 PeMYBs, which were further divided into 11 subgroups, were identified in moso bamboo (Phyllostachys edulis) by using a genome-wide search strategy. Gene structure analysis showed that PeMYBs were significantly different, with exon numbers varying from 4 to 13. Phylogenetic analysis indicated that PeMYBs clustered into 27 clades, of which the function of 18 clades has been predicted. In addition, almost all of the PeMYBs were differently expressed in leaves, panicles, rhizomes and shoots based on RNA-seq data. Furthermore, qRT-PCR analysis showed that 12 PeMYBs related to the biosynthesis and deposition of the secondary cell wall (SCW) were constitutively expressed, and their transcript abundance levels have changed significantly with increasing height of the bamboo shoots, for which the degree of lignification continuously increased. This result indicated that these PeMYBs might play fundamental roles in SCW thickening and bamboo shoot lignification. The present comprehensive and systematic study on the members of the MYB family provided a reference and solid foundation for further functional analysis of MYB TFs in moso bamboo.