Genome-wide analysis of shoot growth-associated alternative splicing in moso bamboo

Comprehensive profiling of epigenetic modifications in fast-growing Moso bamboo shoots

The fast growth of Moso bamboo (Phyllostachys edulis) shoots is caused by the rapid elongation of each internode. However, the key underlying cellular processes and epigenetic mechanisms remain largely unexplored. We used microscopy and multi-omics approaches to investigate two regions (bottom and middle) of the 18th internode from shoots of two different heights (2 and 4 m). We observed that internode cells become longer, and that lignin biosynthesis and glycosyltransferase family 43 (GT43) genes are substantially upregulated with shoot height. Nanopore direct RNA sequencing (DRS) revealed a higher N6-methyladenine (m6A) modification rate in 2-m shoots than in 4-m shoots. In addition, different specific m6A modification sites were enriched at different growth stages. Global DNA methylation profiling indicated that DNA methylation levels are higher in 4-m shoots than in 2-m shoots. We also detected shorter poly(A) tail lengths (PALs) in 4-m shoots compared with 2-m shoots. Genes showing differential PAL were mainly enriched in the functional terms of protein translation and vesicle fusion. An association analysis between PALs and DNA methylation strongly suggested that gene body CG methylation levels are positively associated with PAL. This study provides valuable information to better understand post-transcriptional regulations responsible for fast-growing shoots in Moso bamboo.

Comparative Analysis of Alternative Splicing in Two Contrasting Apple Cultivars Defense against Alternaria alternata Apple Pathotype Infection

Alternaria blotch disease, caused by the Alternaria alternata apple pathotype (A. alternata AP), is one of the most serious fungal diseases in apples. Alternative splicing (AS), one of the pivotal post-transcriptional regulatory mechanisms, plays essential roles in various disease resistance responses. Here, we performed RNA-Seq for two apple cultivars (resistant cultivar 'Jonathan' (J) and susceptible cultivar 'Starking Delicious' (SD)) infected by A. alternata AP to further investigate their AS divergence. In total, 1454, 1780, 1367 and 1698 specifically regulated differential alternative splicing (DAS) events were detected in J36, J72, SD36 and SD72 groups, respectively. Retained intron (RI) was the dominant AS pattern. Conformably, 642, 764, 585 and 742 uniquely regulated differentially spliced genes (DSGs) were found during A. alternata AP infection. Comparative analysis of AS genes in differential splicing and expression levels suggested that only a small proportion of DSGs overlapped with differentially expressed genes (DEGs). Gene ontology (GO) enrichment analysis demonstrated that the DSGs were significantly enriched at multiple levels of gene expression regulation. Briefly, the specific AS was triggered in apple defense against A. alternata AP. Therefore, this study facilitates our understanding on the roles of AS regulation in response to A. alternata AP infection in apples.

MicroRNAs play important roles in regulating the rapid growth of the Phyllostachys edulis culm internode

Moso bamboo (Phyllostachys edulis) is a fast-growing species with uneven growth and lignification from lower to upper segments within one internode. MicroRNAs (miRNAs) play a vital role in post-transcriptional regulation in plants. However, how miRNAs regulate fast growth in bamboo internodes is poorly understood. In this study, one moso bamboo internode was divided during early rapid growth into four segments called F4 (bottom) to F1 (upper) and these were then analysed for transcriptomes, miRNAs and degradomes. The F4 segment had a higher number of actively dividing cells as well as a higher content of auxin (IAA), cytokinin (CK) and gibberellin (GA) compared with the F1 segment. RNA-seq analysis showed DNA replication and cell division-associated genes highly expressed in F4 rather than in F1. In total, 63 miRNAs (DEMs) were identified as differentially expressed between F4 and F1. The degradome and the transcriptome indicated that many downstream transcription factors and hormonal responses genes were modulated by DEMs. Several miR-target interactions were further validated by tobacco co-infiltration. Our findings give new insights into miRNA-mediated regulatory pathways in bamboo, and will contribute to a comprehensive understanding of the molecular mechanisms governing rapid growth.

Whole-genome characterization of chronological age-associated changes in methylome and circular RNAs in moso bamboo (Phyllostachys edulis) from vegetative to floral growth

In mammals, DNA methylation is associated with aging. However, age-related DNA methylation changes during phase transitions largely remain unstudied in plants. Moso bamboo (Phyllostachys edulis) requires a very long time to transition from the vegetative to the floral phase. To comprehensively investigate the association of DNA methylation with aging, we present here single-base-resolution DNA methylation profiles using both high-throughput bisulfite sequencing and single-molecule nanopore-based DNA sequencing, covering the long period of vegetative growth and transition to flowering in moso bamboo. We discovered that CHH methylation gradually accumulates from vegetative to reproductive growth in a time-dependent fashion. Differentially methylated regions, correlating with chronological aging, occurred preferentially at both transcription start sites and transcription termination sites. Genes with CG methylation changes showed an enrichment of Gene Ontology (GO) categories in 'vegetative to reproductive phase transition of meristem'. Combining methylation data with mRNA sequencing revealed that DNA methylation in promoters, introns and exons may have different roles in regulating gene expression. Finally, circular RNA (circRNA) sequencing revealed that the flanking introns of circRNAs are hypermethylated and enriched in long terminal repeat (LTR) retrotransposons. Together, the observations in this study provide insights into the dynamic DNA methylation and circRNA landscapes, correlating with chronological age, which paves the way to study further the impact of epigenetic factors on flowering in moso bamboo.

Profiling of circular RNA N6 -methyladenosine in moso bamboo (Phyllostachys edulis) using nanopore-based direct RNA sequencing

N⁶ -methyladenosine (m⁶ A) is a prevalent modification in messenger RNAs and circular RNAs that play important roles in regulating various aspects of RNA metabolism. However, the occurrence of the m⁶ A modification in plant circular RNAs has not been reported. A widely used method to identify m⁶ A modifications relies on m⁶ A-specific antibodies followed by next-generation sequencing of precipitated RNAs (MeRIP-Seq). However, one limitation of MeRIP-Seq is that it does not provide the precise location of m⁶ A at single-nucleotide resolution. Although more recent sequencing techniques such as Nanopore-based direct RNA sequencing (DRS) can overcome such limitations, the technology does not allow sequencing of circular RNAs, as these molecules lack a poly(A) tail. Here, we developed a novel method to detect the precise location of m⁶ A modifications in circular RNAs using Nanopore DRS. We first enriched our samples for circular RNAs, which we then fragmented and sequenced on the Nanopore platform with a customized protocol. Using this method, we identified 470 unique circular RNAs from DRS reads based on the back-spliced junction region. Among exonic circular RNAs, about 10% contained m⁶ A sites, which mainly occurred around acceptor and donor splice sites. This study demonstrates the utility of our antibody-independent method in identifying total and methylated circular RNAs using Nanopore DRS. This method has the additional advantage of providing the exact location of m⁶ A sites at single-base resolution in circular RNAs or linear transcripts from non-coding RNA without poly(A) tails.

Multi-omics analysis of cellular pathways involved in different rapid growth stages of moso bamboo

Moso bamboo (Phyllostachys edulis (Carriere) J. Houzeau) is a rapidly growing grass of industrial and ecological importance. However, the molecular mechanisms of its remarkable growth are not well understood. In this study, we investigated the early-stage growth of moso bamboo shoots and defined three different growth stages based on histological and biochemical analyses, namely, starting of cell division (SD), rapid division (RD) and rapid elongation (RE). Further analyses on potentially relevant cellular pathways in these growth stages using multi-omics approaches such as transcriptomics and proteomics revealed the involvement of multiple cellular pathways, including DNA replication, repair and ribosome biogenesis. A total of 8045 differentially expressed genes (DEGs) and 1053 differentially expressed proteins (DEPs) were identified in our analyses. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses of detected DEGs identified several key biological pathways such as phytohormone metabolism, signal transduction, cell wall development and carbohydrate metabolism. The comparative analysis of proteins displayed that a total of 213 DEPs corresponded with DEGs and 3 significant expression profiles that could be promoting the fast growth of bamboo internodes. Moreover, protein-protein interaction network prediction analysis is suggestive of the involvement of five major proteins of signal transduction, DNA synthesis and RNA transcription, and may act as key elements responsible for the rapid shoot growth. Our work exploits multi-omics and bioinformatic approaches to unfurl the complexity of molecular networks involved in the rapid growth of moso bamboo and opens up questions related to the interactions between the functions played by individual molecular pathway.

The interplay between microRNA and alternative splicing of linear and circular RNAs in eleven plant species

MOTIVATION

MicroRNA (miRNA) and alternative splicing (AS)-mediated post-transcriptional regulation has been extensively studied in most eukaryotes. However, the interplay between AS and miRNAs has not been explored in plants. To our knowledge, the overall profile of miRNA target sites in circular RNAs (circRNA) generated by alternative back splicing has never been reported previously. To address the challenge, we identified miRNA target sites located in alternatively spliced regions of the linear and circular splice isoforms using the up-to-date single-molecule real-time (SMRT) isoform sequencing (Iso-Seq) and Illumina sequencing data in eleven plant species.

RESULTS

In total, we identified 399 401 and 114 574 AS events from linear and circular RNAs, respectively. Among them, there were 64 781 and 41 146 miRNA target sites located in linear and circular AS region, respectively. In addition, we found 38 913 circRNAs to be overlapping with 45 648 AS events of its own parent isoforms, suggesting circRNA regulation of AS of linear RNAs by forming R-loop with the genomic locus. Here, we present a comprehensive database of miRNA targets in alternatively spliced linear and circRNAs (ASmiR) and a web server for deposition and identification of miRNA target sites located in the alternatively spliced region of linear and circular RNAs. This database is accompanied by an easy-to-use web query interface for meaningful downstream analysis. Plant research community can submit user-defined datasets to the web service to search AS regions harboring small RNA target sites. In conclusion, this study provides an unprecedented resource to understand regulatory relationships between miRNAs and AS in both gymnosperms and angiosperms.

AVAILABILITY AND IMPLEMENTATION

The readily accessible database and web-based tools are available at http://forestry.fafu.edu.cn/bioinfor/db/ASmiR.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

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.

Global identification of alternative splicing in Shiraia bambusicola and analysis of its regulation in hypocrellin biosynthesis

Hypocrellins, as natural pigments from Shiraia bambusicola, have extensive applications in the agricultural, cosmetic, food, and feed industries, and play a vital role in photodynamic therapy for anticancer and antiviral treatments. However, environmental stresses are always the bottlenecks for increasing hypocrellin yield during the process of fermentation. Pre-mRNA alternative splicing (AS) is an essential mechanism in the defense of abiotic stresses in the animal and plant kingdom, but is seldom involved in fungi. In this study, AS from genome-wide sequencing and RNA-seq data for S. bambusicola was analyzed for the first time. Interestingly, the proportion of AS in S. bambusicola was 38.44% (most of them participated in metabolic processes, covering pigmentation and response to stimulus), a much higher ratio than seen in that of other fungal species (1.3-18%). Here, we identified the relationship of AS and secondary metabolic (SM) biosynthesis under a series of abiotic stresses. Suitable fungicides suppressed hypocrellin production significantly, and AS occurred in key functional genes (sbFLO, sbMFS, sbPKS) of hypocrellin biosynthesis. In contrast, H₂O₂ improved the yield of hypocrellins, but AS were not found in the corresponding gene cluster. A further study showed that overexpressing an isoform of sbPKS (sbPKSa) in Shiraia bambusicola could dramatically down-regulate the expression of the original gene sbPKS and nearly inhibit the production of hypocrellins. Altogether, our study strongly supported the hypothesis that AS had a vital role in the regulation of hypocrellin biosynthesis under stresses, and initially explored whether SM functional genes were relevant for fungi.

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.

Expression Analysis and Regulation Network Identification of the CONSTANS-Like Gene Family in Moso Bamboo (Phyllostachys edulis) Under Photoperiod Treatments

CONSTANS (CO)/CONSTANS-like (COL) genes that have been studied in annual model plants such as Arabidopsis thaliana and Oryza sativa play key roles in the photoperiodic flowering pathway. Moso bamboo is a perennial plant characterized by a long vegetative stage and flowers synchronously followed by widespread death. However, the characteristics of COL in moso bamboo remain unclear. In view of this, we performed a genome-wide identification and expression analysis of the COL gene family in moso bamboo. Fourteen nonredundant PheCOL genes were identified, and we analyzed gene structures, phylogeny, and subcellular location predictions. Phylogenetic analyses indicated that 14 PheCOLs could be clustered into three groups, and each clade was well supported by conserved intron/exon structures and motifs. A number of light-related and tissue-specific cis-elements were randomly distributed within the promoter sequences of the PheCOLs. The expression profiling of PheCOL genes in various tissues and developmental stages revealed that most of PheCOL genes were most highly expressed in the leaves and took part in moso bamboo flower development and rapid shoot growth. In addition, the transcription of PheCOLs exhibited a clear diurnal oscillation in both long-day and short-day conditions. Most of the PheCOL genes were inhibited under light treatment and upregulated in dark conditions. PheCOLs can interact with each other. Subcellular localization result showed that PheCOL14 encoded a cell membrane protein, and it bound to the promoter of PheCOL3. Taken together, the results of this study will be useful not only as they contribute to comprehensive information for further analyses of the molecular functions of the PheCOL gene family, but also will provide a theoretical foundation for the further construction of moso bamboo photoperiod regulation networks.

GSK3/shaggy-like kinase 1 ubiquitously regulates cell growth from Arabidopsis to Moso bamboo (Phyllostachys edulis)

Moso bamboo (Phyllostachys edulis) is one of the fastest growing species with a maximum growth rate of 1 m/day. However, the regulator genes for this explosive growth phenomenon have not been functionally studied. Here, we found that Moso bamboo GSK3/shaggy-like kinase 1 (PeGSK1) acts as a negative regulator of cell growth. Over-expression of PeGSK1 in Arabidopsis showed significant growth arrest phenotypes, including dwarfism, small leaves, reduced cell length, and disturbed cell elongation of petiole. Furthermore, Overexpression of PeGSK1 fully inhibited the longer hypocotyl phenotype of Arabidopsis atgsk1 mutants. In addition, PeGSK1-overexpressing lines were resistant to exogenous BR treatment and PeGSK1 interacted with the brassinosteroid signal transduction key regulator BZR1. The BZR1-dependent cell growth genes were down-regulated in PeGSK1-overexpressing lines. These results indicated that PeGSK1 is functionally similar to AtGSK1 and inhibited cell growth via the brassinosteroid signaling pathway. Importantly, PeGSK1 also interacted with PeBZR1, and the expression pattern of PeGSK1 was negatively correlated with the internode elongation of bamboo, indicating that PeGSK1 is involved in the cell growth of bamboo. In summary, our results provide insight into the role of brassinosteroids in the rapid-growth of bamboo culms and identifying target genes for the genetic manipulation of plant height.

Genome-Wide Profiling of Circular RNAs in the Rapidly Growing Shoots of Moso Bamboo (Phyllostachys edulis)

Circular RNAs, including circular exonic RNAs (circRNA), circular intronic RNAs (ciRNA) and exon-intron circRNAs (EIciRNAs), are a new type of noncoding RNAs. Growing shoots of moso bamboo (Phyllostachys edulis) represent an excellent model of fast growth and their circular RNAs have not been studied yet. To understand the potential regulation of circular RNAs, we systematically characterized circular RNAs from eight different developmental stages of rapidly growing shoots. Here, we identified 895 circular RNAs including a subset of mutually inclusive circRNA. These circular RNAs were generated from 759 corresponding parental coding genes involved in cellulose, hemicellulose and lignin biosynthetic process. Gene co-expression analysis revealed that hub genes, such as DEFECTIVE IN RNA-DIRECTED DNA METHYLATION 1 (DRD1), MAINTENANCE OF METHYLATION (MOM), dicer-like 3 (DCL3) and ARGONAUTE 1 (AGO1), were significantly enriched giving rise to circular RNAs. The expression level of these circular RNAs presented correlation with its linear counterpart according to transcriptome sequencing. Further protoplast transformation experiments indicated that overexpressing circ-bHLH93 generating from transcription factor decreased its linear transcript. Finally, the expression profiles suggested that circular RNAs may have interplay with miRNAs to regulate their cognate linear mRNAs, which was further supported by overexpressing miRNA156 decreasing the transcript of circ-TRF-1 and linear transcripts of TRF-1. Taken together, the overall profile of circular RNAs provided new insight into an unexplored category of long noncoding RNA regulation in moso 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.

The R2R3MYB Gene Family in Phyllostachys edulis: Genome-Wide Analysis and Identification of Stress or Development-Related R2R3MYBs

The MYB transcription factor (TF) is one of the largest gene families in plants and involved to multiple biological processes. However, little is known about the MYB family and its functional role in the genome of moso bamboo. In the present study, a total of 114 R2R3MYB genes were first identified from moso bamboo genome and full-length non-chimeric (FLNC) reads. Phylogenetic analysis coupled with gene structure analysis and motif determination resulted in the division of these PheR2R3MYBs into 17 subgroups. The position of eight proteins along an external branch in the phylogenetic tree suggested their relatively ancient origin. The genes in this group were all substituted by (Met, M)/(Arg, R) at conservative W residues in both R2 and R3 repeats, and half were found to possess no transcriptional activation activity. The analysis of evolutionary patterns and divergence suggests that the expansion of PheMYBs was mainly attributable to whole genome duplication (WGD) under different selection pressures. Expressional analysis based on microarray and qRT-PCR data performed diverse expression patterns of R2R3MYBs in response to both various abiotic stimuli and flower development. Furthermore, the co-expression analysis of R2R3MYBs suggested an intricate interplay of growth- and stress-related responses. Finally, we found a hub gene, PheMYB4, was involved in a complex proteins interaction network. Further functional analysis indicated that ectopic overexpression of its homologous gene, PheMYB4-1, could increase tolerance to cold treatment and sensitivity to drought and salt treatment of transgenic Arabidopsis seedlings. These findings provide comprehensive insights into the MYB family members in moso bamboo and offer candidate MYB genes for further studies on their roles in stress resistance.

Transcriptome characterization of moso bamboo (Phyllostachys edulis) seedlings in response to exogenous gibberellin applications

BACKGROUND

Moso bamboo (Phyllostachys edulis) is a well-known bamboo species of high economic value in the textile industry due to its rapid growth. Phytohormones, which are master regulators of growth and development, serve as important endogenous signals. However, the mechanisms through which phytohormones regulate growth in moso bamboo remain unknown to date.

RESULTS

Here, we reported that exogenous gibberellins (GA) applications resulted in a significantly increased internode length and lignin condensation. Transcriptome sequencing revealed that photosynthesis-related genes were enriched in the GA-repressed gene class, which was consistent with the decrease in leaf chlorophyll concentrations and the lower rate of photosynthesis following GA treatment. Exogenous GA applications on seedlings are relatively easy to perform, thus we used 4-week-old whole seedlings of bamboo for GA- treatment followed by high throughput sequencing. In this study, we identified 932 cis-nature antisense transcripts (cis-NATs), and 22,196 alternative splicing (AS) events in total. Among them, 42 cis-nature antisense transcripts (cis-NATs) and 442 AS events were differentially expressed upon exposure to exogenous GA3, suggesting that post-transcriptional regulation might be also involved in the GA3 response. Targets of differential expression of cis-NATs included genes involved in hormone receptor, photosynthesis and cell wall biogenesis. For example, LAC4 and its corresponding cis-NATs were GA3-induced, and may be involved in the accumulation of lignin, thus affecting cell wall composition.

CONCLUSIONS

This study provides novel insights illustrating how GA alters post-transcriptional regulation and will shed light on the underlying mechanism of growth modulated by GA in moso bamboo.

The association of hormone signalling genes, transcription and changes in shoot anatomy during moso bamboo growth

Moso bamboo is a large, woody bamboo with the highest ecological, economic and cultural value of all the bamboo types and accounts for up to 70% of the total area of bamboo grown. However, the spatiotemporal variation role of moso bamboo shoot during growth period is still unclear. We found that the bamboo shoot growth can be divided into three distinct periods, including winter growth, early growth and late growth based on gene expression and anatomy. In the early growth period, lateral buds germinated from the top of the bamboo joint in the shoot tip. Intercalary meristems grew vigorously during the winter growth period and early growth period, but in the late growth period, mitosis in the intercalary meristems decreased. The expression of cell cycle-associated genes and the quantity of differentially expressed genes were higher in early growth than those in late growth, appearing to be influenced by hormonal concentrations. Gene expression analysis indicates that hormone signalling genes play key roles in shoot growth, while auxin signalling genes play a central role. In situ hybridization analyses illustrate how auxin signalling genes regulate apical dominance, meristem maintenance and lateral bud development. Our study provides a vivid picture of the dynamic changes in anatomy and gene expression during shoot growth in moso bamboo, and how hormone signalling-associated genes participate in moso bamboo shoot growth.

Comprehensive profiling of rhizome-associated alternative splicing and alternative polyadenylation in moso bamboo (Phyllostachys edulis)

Moso bamboo (Phyllostachys edulis) represents one of the fastest-spreading plants in the world, due in part to its well-developed rhizome system. However, the post-transcriptional mechanism for the development of the rhizome system in bamboo has not been comprehensively studied. We therefore used a combination of single-molecule long-read sequencing technology and polyadenylation site sequencing (PAS-seq) to re-annotate the bamboo genome, and identify genome-wide alternative splicing (AS) and alternative polyadenylation (APA) in the rhizome system. In total, 145 522 mapped full-length non-chimeric (FLNC) reads were analyzed, resulting in the correction of 2241 mis-annotated genes and the identification of 8091 previously unannotated loci. Notably, more than 42 280 distinct splicing isoforms were derived from 128 667 intron-containing full-length FLNC reads, including a large number of AS events associated with rhizome systems. In addition, we characterized 25 069 polyadenylation sites from 11 450 genes, 6311 of which have APA sites. Further analysis of intronic polyadenylation revealed that LTR/Gypsy and LTR/Copia were two major transposable elements within the intronic polyadenylation region. Furthermore, this study provided a quantitative atlas of poly(A) usage. Several hundred differential poly(A) sites in the rhizome-root system were identified. Taken together, these results suggest that post-transcriptional regulation may potentially have a vital role in the underground rhizome-root system.