Huayu Sun

Developing genome-wide microsatellite markers of bamboo and their applications on molecular marker assisted taxonomy for accessions in the genus Phyllostachys

Morphology-based taxonomy via exiguously reproductive organ has severely limitation on bamboo taxonomy, mainly owing to infrequent and unpredictable flowering events of bamboo. Here, we present the first genome-wide analysis and application of microsatellites based on the genome of moso bamboo (Phyllostachys edulis) to assist bamboo taxonomy. Of identified 127,593 microsatellite repeat-motifs, the primers of 1,451 microsatellites were designed and 1,098 markers were physically mapped on the genome of moso bamboo. A total of 917 markers were successfully validated in 9 accessions with ~39.8% polymorphic potential. Retrieved from validated microsatellite markers, 23 markers were selected for polymorphic analysis among 78 accessions and 64 alleles were detected with an average of 2.78 alleles per primers. The cluster result indicated the majority of the accessions were consistent with their current taxonomic classification, confirming the suitability and effectiveness of the developed microsatellite markers. The variations of microsatellite marker in different species were confirmed by sequencing and in silico comparative genome mapping were investigated. Lastly, a bamboo microsatellites database (http://www.bamboogdb.org/ssr) was implemented to browse and search large information of bamboo microsatellites. Consequently, our results of microsatellite marker development are valuable for assisting bamboo taxonomy and investigating genomic studies in bamboo and related grass species.

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.

Comprehensive analysis of multi-tissue transcriptome data and the genome-wide investigation of GRAS family in Phyllostachys edulis

GRAS family is one of plant specific transcription factors and plays diverse roles in the regulation of plant growth and development as well as in the plant disease resistance and abiotic stress responses. However, the investigation of GRAS family and multi-tissue gene expression profiles still remains unavailable in bamboo (Phyllostachys edulis). Here, we applied RNA-Seq analysis to monitor global transcriptional changes and investigate expression patterns in the five tissues of Ph. edulis, and analyzed a large-scale transcriptional events and patterns. Moreover, the tissue-specific genes and DEGs in different tissues were detected. For example, DEGs in panicle and leaf tissues were abundant in photosynthesis, glutathione, porphyrin and chlorophyll metabolism, whereas those in shoot and rhizome were majority in glycerophospholipid metabolism. In the portion of Ph. edulis GRAS (PeGRAS) analyses, we performed the analysis of phylogenetic, gene structure, conserved motifs, and analyzed the expression profiles of PeGRASs in response to high light and made a co-expression analysis. Additionally, the expression profiles of PeGRASs were validated using quantitative real-time PCR. Thus, PeGRASs based on dynamics profiles of gene expression is helpful in uncovering the specific biological functions which might be of critical values for bioengineering to improve bamboo breeding in future.

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

Key messagePeSNAC1, a stress-related NAC1 fromPhyllostachys edulis, was characterized. Ectopic expression in Arabidopsis indicated thatPeSNAC1together withped-miR164b participated in the regulation of organ boundaries and stress tolerance.AbstractNAC (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 Fv/Fm values declined more slowly in sense OX-PeSNAC1 plants than in Col-0, with lower survival rates and a more rapid decrease in Fv/Fm 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.

Discovery and comparative profiling of microRNAs in representative monopodial bamboo (Phyllostachys edulis) and sympodial bamboo (Dendrocalamus latiflorus).

Background According to the growth pattern of bamboo, sympodial bamboo and monopodial bamboo are considered as two mainly kinds of bamboo. They have different phenotypes and different characteristics in developmental stage. Much attention had been paid on the study of bamboo cultivation, processing, physiology, biochemistry and molecular biology, which had made great progresses in the last decade, especially for the highlighted achievement of the bamboo genomics. However, there is no information available on concerning comparative profiling of miRNAs between sympodial bamboo and monopodial bamboo, which might play important roles in the regulation of bamboo development. Methodology/Principal Findings We identified the profiles of small RNAs using leaf tissues from one sympodial bamboo i.e. moso bamboo (Phyllostachys edulis) and another monopodial bamboo i.e. ma bamboo (Dendrocalamus latiflorus). The result showed that there were 19,295,759 and 11,513,888 raw sequence reads, in which 92 and 69 conserved miRNAs, as well as 95 and 62 novel miRNAs were identified in moso bamboo and ma bamboo, respectively. The ratio of high conserved miRNA families in ma bamboo is more than that in moso bamboo. In addition, a total of 49 and 106 potential targets were predicted in moso bamboo and ma bamboo, respectively, in which several targets for novel miRNAs are transcription factors that play important roles in plant development. More importantly, annotation of differentially expressed target genes was performed based on the analysis of pathway and gene ontology terms enrichment. Conclusions/Significance This study provides the first large-scale sight of discovery and comparative characterization of miRNAomes between two representative bamboos belonged to sympodial bamboo and monopodial bamboo, respectively. Although it will be necessary to validate the function of miRNAs through more experimental research in further, these results lay a foundation for unraveling the miRNA-mediated molecular processes in different kinds of bamboo.

Transcriptome-based investigation of cirrus development and identifying microsatellite markers in rattan (Daemonorops jenkinsiana)

Rattan is an important group of regenerating non-wood climbing palm in tropical forests. The cirrus is an essential climbing organ and provides morphological evidence for evolutionary and taxonomic studies. However, limited data are available on the molecular mechanisms underlying the development of the cirrus. Thus, we performed in-depth transcriptomic sequencing analyses to characterize the cirrus development at different developmental stages of Daemonorops jenkinsiana. The result showed 404,875 transcripts were assembled, including 61,569 high-quality unigenes were identified, of which approximately 76.16% were annotated and classified by seven authorized databases. Moreover, a comprehensive analysis of the gene expression profiles identified differentially expressed genes (DEGs) concentrated in developmental pathways, cell wall metabolism, and hook formation between the different stages of the cirri. Among them, 37 DEGs were validated by qRT-PCR. Furthermore, 14,693 transcriptome-based microsatellites were identified. Of the 168 designed SSR primer pairs, 153 were validated and 16 pairs were utilized for the polymorphic analysis of 25 rattan accessions. These findings can be used to interpret the molecular mechanisms of cirrus development, and the developed microsatellites markers provide valuable data for assisting rattan taxonomy and expanding the understanding of genomic study in rattan.

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 bamboos adaptation to different light environment.

The chromosome-level genome assemblies of two rattans (Calamus simplicifolius and Daemonorops jenkinsiana)

Abstract Background Calamus simplicifolius and Daemonorops jenkinsiana are two representative rattans, the most significant material sources for the rattan industry. However, the lack of reference genome sequences is a major obstacle for basic and applied biology on rattan. Findings We produced two chromosome-level genome assemblies of C. simplicifolius and D. jenkinsiana using Illumina, Pacific Biosciences, and Hi-C sequencing data. A total of ∼730 Gb and ∼682 Gb of raw data covered the predicted genome lengths (∼1.98 Gb of C. simplicifolius and ∼1.61 Gb of D. jenkinsiana) to ∼372 × and ∼426 × read depths, respectively. The two de novo genome assemblies, ∼1.94 Gb and ∼1.58 Gb, were generated with scaffold N50s of ∼160 Mb and ∼119 Mb in C. simplicifolius and D. jenkinsiana, respectively. The C. simplicifolius and D. jenkinsiana genomes were predicted to harbor  51,235 and  53,342 intact protein-coding gene models, respectively. Benchmarking Universal Single-Copy Orthologs evaluation demonstrated that genome completeness reached 96.4% and 91.3% in the C. simplicifolius and D. jenkinsiana genomes, respectively. Genome evolution showed that four Arecaceae plants clustered together, and the divergence time between the two rattans was ∼19.3 million years ago. Additionally, we identified 193 and 172 genes involved in the lignin biosynthesis pathway in the C. simplicifolius and D. jenkinsiana genomes, respectively. Conclusions We present the first de novo assemblies of two rattan genomes (C. simplicifolius and D. jenkinsiana). These data will not only provide a fundamental resource for functional genomics, particularly in promoting germplasm utilization for breeding, but also serve as reference genomes for comparative studies between and among different species.

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

Abstract 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.