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

A Tissue-Specific Landscape of Alternative Polyadenylation, lncRNAs, TFs, and Gene Co-expression Networks in Liriodendron chinense

Liriodendron chinense is an economically and ecologically important deciduous tree species. Although the reference genome has been revealed, alternative polyadenylation (APA), transcription factors (TFs), long non-coding RNAs (lncRNAs), and co-expression networks of tissue-specific genes remain incompletely annotated. In this study, we used the bracts, petals, sepals, stamens, pistils, leaves, and shoot apex of L. chinense as materials for hybrid sequencing. On the one hand, we improved the annotation of the genome. We detected 13,139 novel genes, 7,527 lncRNAs, 1,791 TFs, and 6,721 genes with APA sites. On the other hand, we found that tissue-specific genes play a significant role in maintaining tissue characteristics. In total, 2,040 tissue-specific genes were identified, among which 9.2% of tissue-specific genes were affected by APA, and 1,809 tissue-specific genes were represented in seven specific co-expression modules. We also found that bract-specific hub genes were associated plant defense, leaf-specific hub genes were involved in energy metabolism. Moreover, we also found that a stamen-specific hub TF Lchi25777 may be involved in the determination of stamen identity, and a shoot-apex-specific hub TF Lchi05072 may participate in maintaining meristem characteristic. Our study provides a landscape of APA, lncRNAs, TFs, and tissue-specific gene co-expression networks in L. chinense that will improve genome annotation, strengthen our understanding of transcriptome complexity, and drive further research into the regulatory mechanisms of tissue-specific genes.

A moso bamboo transcription factor, Phehdz1, positively regulates the drought stress response of transgenic rice

78 HD-Zip family genes in Phyllostachys edulis were analyzed. Overexpression of Phehdz1 can improve the drought tolerance of transgenic rice and affect its secondary metabolism. Many studies suggested homeodomain-leucine zipper (HD-Zip) transcription factors are important regulators of plant growth and development, signal transduction, and responses to environmental stresses. In this study, 78 moso bamboo (Phyllostachys edulis) HD-Zip genes were investigated and classified into four subfamilies (HD-Zip I-IV). Additionally, Phehdz1 (HD-Zip I gene) was isolated and confirmed to be highly expressed in the roots. A quantitative real-time PCR analysis indicated Phehdz1 expression was significantly induced by drought, high salinity, and abscisic acid (ABA). A transient expression assay proved that Phehdz1 was localized in the nucleus of tobacco cells. Moreover, it could bind to the core region encoded by the H-box sequence (CAATAATTG) in yeast. In response to mannitol treatments, the Phehdz1-overexpressing transgenic rice had a higher germination rate and longer shoots than the wild-type controls. Moreover, Phehdz1-overexpressing rice plants had a higher survival rate as well as higher relative water and proline contents, but a lower malondialdehyde content, than the WT plants after a 30% polyethylene glycol 6000 treatment. Accordingly, the overexpression of Phehdz1 enhances the drought tolerance of transgenic rice. Many of the differentially expressed genes identified by a transcriptome analysis are involved in MAPK signal transduction and the biosynthesis of secondary metabolites. Thus, the overexpression of Phehdz1 enhances the drought stress tolerance of transgenic rice, while also potentially modulating the expression of metabolism-related genes.

Genome-Wide Identification of the Expansin Gene Family and Its Potential Association with Drought Stress in Moso Bamboo

Expansins, a group of cell wall-loosening proteins, are involved in cell-wall loosening and cell enlargement in a pH-dependent manner. According to previous study, they were involved in plant growth and abiotic stress responses. However, information on the biological function of the expansin gene in moso bamboo is still limited. In this study, we identified a total of 82 expansin genes in moso bamboo, clustered into four subfamilies (α-expansin (EXPA), β-expansin (EXPB), expansin-like A (EXLA) and expansin-like B (EXPB)). Subsequently, the molecular structure, chromosomal location and phylogenetic relationship of the expansin genes of Phyllostachys edulis (PeEXs) were further characterized. A total of 14 pairs of tandem duplication genes and 31 pairs of segmented duplication genes were also identified, which may promote the expansion of the expansin gene family. Promoter analysis found many cis-acting elements related to growth and development and stress response, especially abscisic acid response element (ABRE). Expression pattern revealed that most PeEXs have tissue expression specificity. Meanwhile, the expression of some selected PeEXs was significantly upregulated mostly under abscisic acid (ABA) and polyethylene glycol (PEG) treatment, which implied that these genes actively respond to expression under abiotic stress. This study provided new insights into the structure, evolution and function prediction of the expansin gene family in moso bamboo.

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

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

Transcriptome analysis of lateral buds from Phyllostachys edulis rhizome during germination and early shoot stages

BACKGROUND

The vegetative growth is an important stage for plants when they conduct photosynthesis, accumulate and collect all resources needed and prepare for reproduction stage. Bamboo is one of the fastest growing plant species. The rapid growth of Phyllostachys edulis results from the expansion of intercalary meristem at the basal part of nodes, which are differentiated from the apical meristem of rhizome lateral buds. However, little is known about the major signaling pathways and players involved during this rapid development stage of bamboo. To study this question, we adopted the high-throughput sequencing technology and compared the transcriptomes of Moso bamboo rhizome buds in germination stage and late development stage.

RESULTS

We found that the development of Moso bamboo rhizome lateral buds was coordinated by multiple pathways, including meristem development, sugar metabolism and phytohormone signaling. Phytohormones have fundamental impacts on the plant development. We found the evidence of several major hormones participating in the development of Moso bamboo rhizome lateral bud. Furthermore, we showed direct evidence that Gibberellic Acids (GA) signaling participated in the Moso bamboo stem elongation.

CONCLUSION

Significant changes occur in various signaling pathways during the development of rhizome lateral buds. It is crucial to understand how these changes are translated to Phyllostachys edulis fast growth. These results expand our knowledge on the Moso bamboo internodes fast growth and provide research basis for further study.

Genome-Wide Identification and Characterization of Hexokinase Genes in Moso Bamboo (Phyllostachys edulis)

Plant hexokinases (HXKs) are a class of multifunctional proteins that not only act as the enzymes required for hexose phosphorylation but also serve as sugar sensors that repress the expression of some photosynthetic genes when internal glucose level increases and regulators of cell metabolism and some sugar-related signaling pathways independent on their catalytic actives. The HXKs have been studied in many plants; however, limited information is available on HXKs of moso bamboo (Phyllostachys edulis). In this study, we identified and characterized 12 hexokinase genes in moso bamboo. Phylogenetic analysis revealed that the moso bamboo hexokinases (PeHXKs) were classifiable into five subfamilies which represented the three types of hexokinases in plants. Gene structure and conserved motif analysis showed that the PeHXK genes contained diverse numbers of introns and exons and that the encoded proteins showed similar motif organization within each subfamily. Multiple sequence alignment revealed that the PeHXK proteins contained conserved domains, such as phosphate 1 (P1), phosphate 2 (P2), adenosine, and a sugar-binding domain. Evolutionary divergence analysis indicated that the PeHXK, OsHXK, and BdHXK families underwent negative selection and experienced a large-scale duplication event approximately 19-319 million years ago. Expression analysis of the PeHXK genes in the leaf, stem, root, and rhizome of moso bamboo seedlings indicated that the PeHXKs perform pivotal functions in the development of moso bamboo. A protein subcellular localization assay showed that PeHXK5a, PeHXK8, and PeHXK3b were predominantly localized in mitochondria, and PeHXK8 protein was also detected in the nucleus. The HXK activity of the PeHXK5a, PeHXK8, and PeHXK3b was verified by a functional complementation assay using the HXK-deficient triple-mutant yeast strain YSH7.4-3C (hxk1, hxk2, and glk1), and the results showed that the three PeHXKs had the plant HXK-specific enzyme traits. The present findings would provide a foundation for further functional analysis of the PeHXK gene family.

Identification and Evolution of the WUSCHEL-Related Homeobox Protein Family in Bambusoideae

Bamboos (Bambusoideae) are fast-growing species due to their rapid growth rate and ability to reproduce annually via cloned buds produced on the rhizome. WUSCHEL-related homeobox (WOX) genes have been reported to regulate shoot apical meristem organization, lateral organ formation, cambium and vascular proliferation, and so on, but have rarely been studied in bamboos. In this study, the WOXs of both herbaceous bamboo species (12 OlaWOXs and nine RguWOXs) and woody bamboo species (18 GanWOXs, 27 PheWOXs, and 26 BamWOXs) were identified and categorized into three clades based on their phylogenetic relationship-ancient, intermediate, or WUS clade. Polyploidy is the major driver of the expansion of the bamboo WOX family. Eight conserved domains, besides the homeodomain, were identified by comparatively analyzing the WOXs of dicot and monocot species. Intensive purifying selection pressure in the coding region of specific domains explained the functional similarity of WOXs between different species. For Bambusoideae WOXs, polyploidy is the major driver of the expansion of the WOX family. Stronger purifying selection was found in orthologous WOXs of Bambusoideae, especially for WOX4s and WOX5s, which are conserved not only at the translational levels, but also at the genome level. Several conserved cis-acting elements were discovered at similar position in the promoters of the orthologous WOXs. For example, AP2/ERF protein-binding elements and B3 protein-binding elements were found in the promoters of the bamboo WOX4, while MYB protein-binding elements and Dof protein-binding elements were found in the promoters of bamboo WOX5, and MADS protein-binding sites was found in the promoters of bamboo WUS, WOX3, and WOX9. These conserved positions may play an important role in regulating the expression of bamboo WOXs. Our work provides insight into the origin and evolution of bamboo WOXs, and will facilitate functional investigations of the clonal propagation of bamboos.

Overexpression of PheNAC3 from moso bamboo promotes leaf senescence and enhances abiotic stress tolerance in Arabidopsis

The NAC family is one of the largest transcription factor families unique to plants, which regulates the growth and development, biotic and abiotic stress responses, and maturation and senescence in plants. In this study, PheNAC3, a NAC gene, was isolated and characterized from moso bamboo (Phyllostachys edulis). PheNAC3 belong to the NAC1 subgroup and has a conserved NAC domain on the N-terminus, which with 88.74% similarity to ONAC011 protein. PheNAC3 localized in the nucleus and exhibited transactivation activity. PheNAC3 was upregulated during the process of senescence of leaves and detected shoots. PheNAC3 was also induced by ABA, MeJA, NaCl and darkness, but it had no remarkable response to PEG and SA treatments. Overexpression of PheNAC3 could cause precocious senescence in Arabidopsis. Transgenic Arabidopsis displayed faster seed germination, better seedling growth, and a higher survival rate than the wild-type under salt or drought stress conditions. Moreover, AtSAG12 associated with senescence and AtRD29A and AtRD29b related to ABA were upregulated by PheNAC3 overexpression, but AtCAB was inhibited. These findings show that PheNAC3 may participate in leaf senescence and play critical roles in the salt and drought stress response.

Transcriptome profiling reveals the crucial biological pathways involved in cold response in Moso bamboo (Phyllostachys edulis)

Most bamboo species including Moso bamboo (Phyllostachys edulis) are tropical or subtropical plants that greatly contribute to human well-being. Low temperature is one of the main environmental factors restricting bamboo growth and geographic distribution. Our knowledge of the molecular changes during bamboo adaption to cold stress remains limited. Here, we provided a general overview of the cold-responsive transcriptional profiles in Moso bamboo by systematically analyzing its transcriptomic response under cold stress. Our results showed that low temperature induced strong morphological and biochemical alternations in Moso bamboo. To examine the global gene expression changes in response to cold, 12 libraries (non-treated, cold-treated 0.5, 1 and 24 h at -2 °C) were sequenced using an Illumina sequencing platform. Only a few differentially expressed genes (DEGs) were identified at early stage, while a large number of DEGs were identified at late stage in this study, suggesting that the majority of cold response genes in bamboo are late-responsive genes. A total of 222 transcription factors from 24 different families were differentially expressed during 24-h cold treatment, and the expressions of several well-known C-repeat/dehydration responsive element-binding factor negative regulators were significantly upregulated in response to cold, indicating the existence of special cold response networks. Our data also revealed that the expression of genes related to cell wall and the biosynthesis of fatty acids were altered in response to cold stress, indicating their potential roles in the acquisition of bamboo cold tolerance. In summary, our studies showed that both plant kingdom-conserved and species-specific cold response pathways exist in Moso bamboo, which lays the foundation for studying the regulatory mechanisms underlying bamboo cold stress response and provides useful gene resources for the construction of cold-tolerant bamboo through genetic engineering in the future.

Alternative Splicing Enhances the Transcriptome Complexity of Liriodendron chinense

Alternative splicing (AS) plays pivotal roles in regulating plant growth and development, flowering, biological rhythms, signal transduction, and stress responses. However, no studies on AS have been performed in Liriodendron chinense, a deciduous tree species that has high economic and ecological value. In this study, we used multiple tools and algorithms to analyze transcriptome data derived from seven tissues via hybrid sequencing. Although only 17.56% (8,503/48,408) of genes in L. chinense were alternatively spliced, these AS genes occurred in 37,844 AS events. Among these events, intron retention was the most frequent AS event, producing 1,656 PTC-containing and 3,310 non-PTC-containing transcripts. Moreover, 183 long noncoding RNAs (lncRNAs) also underwent AS events. Furthermore, weighted gene coexpression network analysis (WGCNA) revealed that there were great differences in the activities of transcription and post-transcriptional regulation between pistils and leaves, and AS had an impact on many physiological and biochemical processes in L. chinense, such as photosynthesis, sphingolipid metabolism, fatty acid biosynthesis and metabolism. Moreover, our analysis showed that the features of genes may affect AS, as AS genes and non-AS genes had differences in the exon/intron length, transcript length, and number of exons/introns. In addition, the structure of AS genes may impact the frequencies and types of AS because AS genes with more exons or introns tended to exhibit more AS events, and shorter introns tended to be retained, whereas shorter exons tended to be skipped. Furthermore, eight AS genes were verified, and the results were consistent with our analysis. Overall, this study reveals that AS and gene interaction are mutual-on one hand, AS can affect gene expression and translation, while on the other hand, the structural characteristics of the gene can also affect AS. This work is the first to comprehensively report on AS in L. chinense, and it can provide a reference for further research on AS in L. chinense.

Development and Deployment of High-Throughput Retrotransposon-Based Markers Reveal Genetic Diversity and Population Structure of Asian Bamboo

Bamboo, a non-timber grass species, known for exceptionally fast growth is a commercially viable crop. Long terminal repeat (LTR) retrotransposons, the main class I mobile genetic elements in plant genomes, are highly abundant (46%) in bamboo, contributing to genome diversity. They play significant roles in the regulation of gene expression, chromosome size and structure as well as in genome integrity. Due to their random insertion behavior, interspaces of retrotransposons can vary significantly among bamboo genotypes. Capitalizing this feature, inter-retrotransposon amplified polymorphism (IRAP) is a high-throughput marker system to study the genetic diversity of plant species. To date, there are no transposon based markers reported from the bamboo genome and particularly using IRAP markers on genetic diversity. Phyllostachys genus of Asian bamboo is the largest of the Bambusoideae subfamily, with great economic importance. We report structure-based analysis of bamboo genome for the LTR-retrotransposon superfamilies, Ty3-gypsy and Ty1-copia, which revealed a total of 98,850 retrotransposons with intact LTR sequences at both the ends. Grouped into 64,281 clusters/scaffold using CD-HIT-EST software, only 13 clusters of retroelements were found with more than 30 LTR sequences and with at least one copy having all intact protein domains such as gag and polyprotein. A total of 16 IRAP primers were synthesized, based on the high copy numbers of conserved LTR sequences. A study using these IRAP markers on genetic diversity and population structure of 58 Asian bamboo accessions belonging to the genus Phyllostachys revealed 3340 amplicons with an average of 98% polymorphism. The bamboo accessions were collected from nine different provinces of China, as well as from Italy and America. A three phased approach using hierarchical clustering, principal components and a model based population structure divided the bamboo accessions into four sub-populations, PhSP1, PhSP2, PhSP3 and PhSP4. All the three analyses produced significant sub-population wise consensus. Further, all the sub-populations revealed admixture of alleles. The analysis of molecular variance (AMOVA) among the sub-populations revealed high intra-population genetic variation (75%) than inter-population. The results suggest that Phyllostachys bamboos are not well evolutionarily diversified, although geographic speciation could have occurred at a limited level. This study highlights the usability of IRAP markers in determining the inter-species variability of Asian bamboos.

Affinities of Terminal Inverted Repeats to DNA Binding Domain of Transposase Affect the Transposition Activity of Bamboo Ppmar2 Mariner-Like Element

Mariner-like elements (MLE) are a super-family of DNA transposons widespread in animal and plant genomes. Based on their transposition characteristics, such as random insertions and high-frequency heterogeneous transpositions, several MLEs have been developed to be used as tools in gene tagging and gene therapy. Two active MLEs, Ppmar1 and Ppmar2, have previously been identified in moso bamboo (Phyllostachys edulis). Both of these have a preferential insertion affinity to AT-rich region and their insertion sites are close to random in the host genome. In Ppmar2 element, we studied the affinities of terminal inverted repeats (TIRs) to DNA binding domain (DBD) and their influence on the transposition activity. We could identify two putative boxes in the TIRs which play a significant role in defining the TIR's affinities to the DBD. Seven mutated TIRs were constructed, differing in affinities based on similarities with those of other plant MLEs. Gel mobility shift assays showed that the TIR mutants with mutation sites G669A-C671A had significantly higher affinities than the mutants with mutation sites C657T-A660T. The high-affinity TIRs indicated that their transposition frequency was 1.5-2.0 times higher than that of the wild type TIRs in yeast transposition assays. The MLE mutants with low-affinity TIRs had relatively lower transposition frequency from that of wild types. We conclude that TIR affinity to DBD significantly affects the transposition activity of Ppmar2. The mutant MLEs highly active TIRs constructed in this study can be used as a tool for bamboo genetic studies.

Large Scale Profiling of Protein Isoforms Using Label-Free Quantitative Proteomics Revealed the Regulation of Nonsense-Mediated Decay in Moso Bamboo (Phyllostachys edulis)

Moso bamboo is an important forest species with a variety of ecological, economic, and cultural values. However, the gene annotation information of moso bamboo is only based on the transcriptome sequencing, lacking the evidence of proteome. The lignification and fiber in moso bamboo leads to a difficulty in the extraction of protein using conventional methods, which seriously hinders research on the proteomics of moso bamboo. The purpose of this study is to establish efficient methods for extracting the total proteins from moso bamboo for following mass spectrometry-based quantitative proteome identification. Here, we have successfully established a set of efficient methods for extracting total proteins of moso bamboo followed by mass spectrometry-based label-free quantitative proteome identification, which further improved the protein annotation of moso bamboo genes. In this study, 10,376 predicted coding genes were confirmed by quantitative proteomics, accounting for 35.8% of all annotated protein-coding genes. Proteome analysis also revealed the protein-coding potential of 1015 predicted long noncoding RNA (lncRNA), accounting for 51.03% of annotated lncRNAs. Thus, mass spectrometry-based proteomics provides a reliable method for gene annotation. Especially, quantitative proteomics revealed the translation patterns of proteins in moso bamboo. In addition, the 3284 transcript isoforms from 2663 genes identified by Pacific BioSciences (PacBio) single-molecule real-time long-read isoform sequencing (Iso-Seq) was confirmed on the protein level by mass spectrometry. Furthermore, domain analysis of mass spectrometry-identified proteins encoded in the same genomic locus revealed variations in domain composition pointing towards a functional diversification of protein isoform. Finally, we found that part transcripts targeted by nonsense-mediated mRNA decay (NMD) could also be translated into proteins. In summary, proteomic analysis in this study improves the proteomics-assisted genome annotation of moso bamboo and is valuable to the large-scale research of functional genomics in moso bamboo. In summary, this study provided a theoretical basis and technical support for directional gene function analysis at the proteomics level in moso bamboo.

Analysis of Transcriptome and Epitranscriptome in Plants Using PacBio Iso-Seq and Nanopore-Based Direct RNA Sequencing

Nanopore sequencing from Oxford Nanopore Technologies (ONT) and Pacific BioSciences (PacBio) single-molecule real-time (SMRT) long-read isoform sequencing (Iso-Seq) are revolutionizing the way transcriptomes are analyzed. These methods offer many advantages over most widely used high-throughput short-read RNA sequencing (RNA-Seq) approaches and allow a comprehensive analysis of transcriptomes in identifying full-length splice isoforms and several other post-transcriptional events. In addition, direct RNA-Seq provides valuable information about RNA modifications, which are lost during the PCR amplification step in other methods. Here, we present a comprehensive summary of important applications of these technologies in plants, including identification of complex alternative splicing (AS), full-length splice variants, fusion transcripts, and alternative polyadenylation (APA) events. Furthermore, we discuss the impact of the newly developed nanopore direct RNA-Seq in advancing epitranscriptome research in plants. Additionally, we summarize computational tools for identifying and quantifying full-length isoforms and other co/post-transcriptional events and discussed some of the limitations with these methods. Sequencing of transcriptomes using these new single-molecule long-read methods will unravel many aspects of transcriptome complexity in unprecedented ways as compared to previous short-read sequencing approaches. Analysis of plant transcriptomes with these new powerful methods that require minimum sample processing is likely to become the norm and is expected to uncover novel co/post-transcriptional gene regulatory mechanisms that control biological outcomes during plant development and in response to various stresses.

Co-expression Gene Network Analysis and Functional Module Identification in Bamboo Growth and Development

Bamboo is one of the fastest-growing non-timber forest plants. Moso bamboo (Phyllostachys edulis) is the most economically valuable bamboo in Asia, especially in China. With the release of the whole-genome sequence of moso bamboo, there are increasing demands for refined annotation of bamboo genes. Recently, large amounts of bamboo transcriptome data have become available, including data on the multiple growth stages of tissues. It is now feasible for us to construct co-expression networks to improve bamboo gene annotation and reveal the relationships between gene expression and growth traits. We integrated the genome sequence of moso bamboo and 78 transcriptome data sets to build genome-wide global and conditional co-expression networks. We overlaid the gene expression results onto the network with multiple dimensions (different development stages). Through combining the co-expression network, module classification and function enrichment tools, we identified 1,896 functional modules related to bamboo development, which covered functions such as photosynthesis, hormone biosynthesis, signal transduction, and secondary cell wall biosynthesis. Furthermore, an online database (http://bioinformatics.cau.edu.cn/bamboo) was built for searching the moso bamboo co-expression network and module enrichment analysis. Our database also includes cis-element analysis, gene set enrichment analysis, and other tools. In summary, we integrated public and in-house bamboo transcriptome data sets and carried out co-expression network analysis and functional module identification. Through data mining, we have yielded some novel insights into the regulation of growth and development. Our established online database might be convenient for the bamboo research community to identify functional genes or modules with important traits.

Whole-Genome and Expression Analyses of Bamboo Aquaporin Genes Reveal Their Functions Involved in Maintaining Diurnal Water Balance in Bamboo Shoots

Water supply is essential for maintaining normal physiological function during the rapid growth of bamboo. Aquaporins (AQPs) play crucial roles in water transport for plant growth and development. Although 26 PeAQPs in bamboo have been reported, the aquaporin-led mechanism of maintaining diurnal water balance in bamboo shoots remains unclear. In this study, a total of 63 PeAQPs were identified, based on the updated genome of moso bamboo (Phyllostachys edulis), including 22 PePIPs, 20 PeTIPs, 17 PeNIPs, and 4 PeSIPs. All of the PeAQPs were differently expressed in 26 different tissues of moso bamboo, based on RNA sequencing (RNA-seq) data. The root pressure in shoots showed circadian rhythm changes, with positive values at night and negative values in the daytime. The quantitative real-time PCR (qRT-PCR) result showed that 25 PeAQPs were detected in the base part of the shoots, and most of them demonstrated diurnal rhythm changes. The expression levels of some PeAQPs were significantly correlated with the root pressure. Of the 86 sugar transport genes, 33 had positive co-expression relationships with 27 PeAQPs. Two root pressure-correlated PeAQPs, PeTIP4;1 and PeTIP4;2, were confirmed to be highly expressed in the parenchyma and epidermal cells of bamboo culm, and in the epidermis, pith, and primary xylem of bamboo roots by in situ hybridization. The authors’ findings provide new insights and a possible aquaporin-led mechanism for bamboo fast growth.

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

Background

Bamboo is one of the most important nontimber forestry products worldwide. However, a chromosome-level reference genome is lacking, and an evolutionary view of alternative splicing (AS) in bamboo remains unclear despite emerging omics data and improved technologies.

Results

Here, we provide a chromosome-level de novo genome assembly of moso bamboo (Phyllostachys edulis) using additional abundance sequencing data and a Hi-C scaffolding strategy. The significantly improved genome is a scaffold N50 of 79.90 Mb, approximately 243 times longer than the previous version. A total of 51,074 high-quality protein-coding loci with intact structures were identified using single-molecule real-time sequencing and manual verification. Moreover, we provide a comprehensive AS profile based on the identification of 266,711 unique AS events in 25,225 AS genes by large-scale transcriptomic sequencing of 26 representative bamboo tissues using both the Illumina and Pacific Biosciences sequencing platforms. Through comparisons with orthologous genes in related plant species, we observed that the AS genes are concentrated among more conserved genes that tend to accumulate higher transcript levels and share less tissue specificity. Furthermore, gene family expansion, abundant AS, and positive selection were identified in crucial genes involved in the lignin biosynthetic pathway of moso bamboo.

Conclusions

These fundamental studies provide useful information for future in-depth analyses of comparative genome and AS features. Additionally, our results highlight a global perspective of AS during evolution and diversification in bamboo.