Genome skimming-based STMS marker discovery and its validation in temperate hill bamboo Drepanostachyum falcatum.

Drepanostachyum falcatum (Nees) Keng f. is one of the most widely distributed shrubby bamboo species in the temperate region of northwest (NW) Himalayas. Along with the other three temperate bamboo species, namely Yushania anceps, Thamnocalamus spathiflorus and Himalayacalamus falconeri, commonly called as 'ringal', and utilized for making various articles of household and commercial purpose by local artisans. Despite huge ecological and socio-economic importance, they are least studied and lacks baseline genetic information. In this study, ~10 Gb genome sequence data with 70.68 million reads were generated for D. falcatum, through genome skimming approach based on high throughput next-generation sequencing technology with Illumina protocol. The high-quality reads were de novo assembled into 31,997 contigs, which comprised 1943 microsatellite repeats. The dinucleotide and trinucleotide repeats were most abundantly distributed in the genome with 52.95 and 41.17%, respectively. Depending on the sufficient flanking sequence, only 1123 repeats were successfully tagged with primer pairs and these sites were designated as sequence-tagged microsatellite (STMS) markers. Further, a subset of 106 STMS markers were validated through PCR amplification; 77 marker loci were successfully amplified, and 48 of these showed polymorphism. Same set of marker loci were also tested for their cross-amplification in other three temperate bamboo species of the NW Himalayas, which revealed good level of transferability (27-48%) but lesser polymorphism (4-12%). In addition, the genomewide in silico cross-amplification revealed poor cross-transferability in other bamboo taxa representing four different phylogenetic lineages, namely Phyllostachys edulis (10.2%), Bonia amplexicaulis (3.03%), Guadua angustifolia (1.60%), Olyra latifolia (0.89%) and Raddia guianensis (0.36%). Ten polymorphic markers were further used to estimate the measures of genetic diversity in two natural populations, which revealed high genetic diversity (polymorphic information content, PIC = 0.889; expected heterozygosity, He = 0.756) and low genetic differentiation (FST=0.061; Nm = 5.445). To the best of our knowledge, this is one of the pioneer studies carried out for the development of genomic STMS markers through genome skimming approach in Indian bamboo species. The marker information generated here is novel and of paramount importance for future genetic studies in D. falcatum as well as other temperate bamboo species through cross-transferability.

Bamboo Transposon Research: Current Status and Perspectives.

Bamboo, a fast-growing non-timber forest plant with many uses, is a valuable species for green development. However, bamboo flowering is very infrequent, extending, in general, for up to 120 years. Ecologically, bamboo species are generally better adapted to various environments than other grasses. Therefore, the species deserves a special status in what could be called Ecological Bioeconomy. An understanding of the genetic processes of bamboo can help us sustainably develop and manage bamboo forests. Transposable elements (TEs), jumping genes or transposons, are major genetic elements in plant genomes. The rapid development of the bamboo reference genome, at the chromosome level, reveals that TEs occupy over 63.24% of the genome. This is higher than found in rice, Brachypodium, and sorghum. The bamboo genome contains diverse families of TEs, which play a significant role in bamboo's biological processes including growth and development. TEs provide important clues for understanding the evolution of the bamboo genome. In this chapter, we briefly describe the current status of research on TEs in the bamboo genome, their regulation, and transposition mechanisms. Perspectives for future research are also provided.

Genome-Wide Characterization and Gene Expression Analyses of GATA Transcription Factors in Moso Bamboo (Phyllostachys edulis).

Moso bamboo is well-known for its rapid-growth shoots and widespread rhizomes. However, the regulatory genes of these two processes are largely unexplored. GATA transcription factors regulate many developmental processes, but their roles in moso bamboo height control and rhizome development remains unexplored. Here, thirty-one bamboo GATA factors (PeGATAs) were identified, which are evolutionarily closer to rice than Arabidopsis, and their gene expression patterns were analyzed in bamboo development and phytohormone response with bioinformatics and molecular methods. Interestingly, PeGATAs could only be classified into three groups. Phytohormone responsive cis-elements were found in PeGATA promoters and the expression profiles showed that PeGATA genes might respond to gibberellin acid and abscisic acid but not to auxin at the transcriptional level. Furthermore, PeGATA genes have a tissue-specific expression pattern in bamboo rhizomes. Interestingly, most PeGATA genes were down-regulated during the rapid-growth of bamboo shoots. In addition, over-expressing one of the PeGATA genes, PeGATA26, significantly repressed the primary root length and plant height of transgenic Arabidopsis plants, which may be achieved by promoting the gibberellin acid turnover. Overall, our results provide insight into the function of GATA transcription factors in bamboo, and into genetic resources for engineering plant height.

Genome-Wide Identification and Expression Analyses of the bZIP Transcription Factor Genes in moso bamboo (Phyllostachys edulis).

The basic leucine zipper (bZIP) transcription factor (TF) family is one of the largest gene families, and play crucial roles in many processes, including stress responses, hormone effects. The TF family also participates in plant growth and development. However, limited information is available for these genes in moso bamboo (Phyllostachys edulis), one of the most important non-timber forest products in the world. In the present study, 154 putative PhebZIP genes were identified in the moso bamboo genome. The phylogenetic analyses indicate that the PhebZIP gene proteins classify into 9 subfamilies and the gene structures and conserved motifs that analyses identified among all PhebZIP proteins suggested a high group-specificity. Microsynteny and evolutionary patterns analyses of the non-synonymous (Ka) and synonymous (Ks) substitution rates and their ratios indicated that paralogous pairs of PhebZIP genes in moso bamboo underwent a large-scale genome duplication event that occurred 7-15 million years ago (MYA). According to promoter sequence analysis, we further selected 18 genes which contain the higher number of cis-regulatory elements for expression analysis. The result showed that these genes are extensively involved in GA-, ABA- and MeJA-responses, with possibly different mechanisms. The tissue-specific expression profiles of PhebZIP genes in five plant tissues/organs/developmental stages suggested that these genes are involved in moso bamboo organ development, especially seed development. Subcellular localization and transactivation activity analysis showed that PhebZIP47 and PhebZIP126 were localized in the nucleus and PhebZIP47 with no transcriptional activation in yeast. Our research provides a comprehensive understanding of PhebZIP genes and may aid in the selection of appropriate candidate genes for further cloning and functional analysis in moso bamboo growth and development, and improve their resistance to stress during their life.

Molecular phylogeny of 21 tropical bamboo species reconstructed by integrating non-coding internal transcribed spacer (ITS1 and 2) sequences and their consensus secondary structure.

The unavailability of the reproductive structure and unpredictability of vegetative characters for the identification and phylogenetic study of bamboo prompted the application of molecular techniques for greater resolution and consensus. We first employed internal transcribed spacer (ITS1, 5.8S rRNA and ITS2) sequences to construct the phylogenetic tree of 21 tropical bamboo species. While the sequence alone could grossly reconstruct the traditional phylogeny amongst the 21-tropical species studied, some anomalies were encountered that prompted a further refinement of the phylogenetic analyses. Therefore, we integrated the secondary structure of the ITS sequences to derive individual sequence-structure matrix to gain more resolution on the phylogenetic reconstruction. The results showed that ITS sequence-structure is the reliable alternative to the conventional phenotypic method for the identification of bamboo species. The best-fit topology obtained by the sequence-structure based phylogeny over the sole sequence based one underscores closer clustering of all the studied Bambusa species (Sub-tribe Bambusinae), while Melocanna baccifera, which belongs to Sub-Tribe Melocanneae, disjointedly clustered as an out-group within the consensus phylogenetic tree. In this study, we demonstrated the dependability of the combined (ITS sequence+structure-based) approach over the only sequence-based analysis for phylogenetic relationship assessment of bamboo.

Transcriptome Sequencing and Analysis for Culm Elongation of the World's Largest Bamboo (Dendrocalamus sinicus).

Dendrocalamus sinicus is the world's largest bamboo species with strong woody culms, and known for its fast-growing culms. As an economic bamboo species, it was popularized for multi-functional applications including furniture, construction, and industrial paper pulp. To comprehensively elucidate the molecular processes involved in its culm elongation, Illumina paired-end sequencing was conducted. About 65.08 million high-quality reads were produced, and assembled into 81,744 unigenes with an average length of 723 bp. A total of 64,338 (79%) unigenes were annotated for their functions, of which, 56,587 were annotated in the NCBI non-redundant protein database and 35,262 were annotated in the Swiss-Prot database. Also, 42,508 and 21,009 annotated unigenes were allocated to gene ontology (GO) categories and clusters of orthologous groups (COG), respectively. By searching against the Kyoto Encyclopedia of Genes and Genomes Pathway database (KEGG), 33,920 unigenes were assigned to 128 KEGG pathways. Meanwhile, 8,553 simple sequence repeats (SSRs) and 81,534 single-nucleotide polymorphism (SNPs) were identified, respectively. Additionally, 388 transcripts encoding lignin biosynthesis were detected, among which, 27 transcripts encoding Shikimate O-hydroxycinnamoyltransferase (HCT) specifically expressed in D. sinicus when compared to other bamboo species and rice. The phylogenetic relationship between D. sinicus and other plants was analyzed, suggesting functional diversity of HCT unigenes in D. sinicus. We conjectured that HCT might lead to the high lignin content and giant culm. Given that the leaves are not yet formed and culm is covered with sheaths during culm elongation, the existence of photosynthesis of bamboo culm is usually neglected. Surprisedly, 109 transcripts encoding photosynthesis were identified, including photosystem I and II, cytochrome b6/f complex, photosynthetic electron transport and F-type ATPase, and 24 transcripts were characterized as antenna proteins that regarded as the main tool for capturing light of plants, implying stem photosynthesis plays a key role during culm elongation due to the unavailability of its leaf. By real-time quantitative PCR, the expression level of 6 unigenes was detected. The results showed the expression level of all genes accorded with the transcriptome data, which confirm the reliability of the transcriptome data. As we know, this is the first study underline the D. sinicus transcriptome, which will deepen the understanding of the molecular mechanisms of culm development. The results may help variety improvement and resource utilization of bamboos.

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.

[Retrieval of leaf area index of moso bamboo forest with Landsat Thematic Mapper image based on PROSAIL canopy radiative transfer model].

The PROSAIL canopy radiative transfer model was used to establish leaf area index (LAI) and canopy reflectance lookup-table for Moso bamboo forest. The combination of Landsat Thematic Mapper (TM) image and this model was then used to retrieve LAI. The results demonstrated that the sensitivity of the input parameters in the PROSAIL model decreased in order of LAI >chlorophyll content (C(ab)) > leaf structure parameters (N) > mean leaf angle (ALA) > equivalent water thickness (C(w)) > dry matter content (C(m)). The most sensitive factors LAI and C(ab) were then used to construct the LAI-canopy reflectance lookup-table. The LAI estimates from the PROSAIL model had good agreement with the reference data, with the coefficient of determination (R2) reached 0.90. The root mean square error (RMSE) and relative RMSE were 0.58 and 13.0%, respectively. However, the mean LAI estimate was higher than the observed value.

The impact of aridification and vegetation type on changes in the community structure of methane-cycling microorganisms in Japanese wetland soils.

Over the years, the wetlands covered by Sphagnum in Bibai, Japan have been turning into areas of aridity, resulting in an invasion of Sasa into the bogs. Yet little is known about the methane-cycling microorganisms in such environments. In this study, the methanotrophic, methanogenic, and archaeal community structures within these two types of wetland vegetation were studied by phylogenetic analysis targeting particulate methane monooxygenase (pmoA), methyl coenzyme M reductase (mcrA), and the archaeal 16S rRNA gene. The pmoA library indicated that Methylomonas and Methylocystis predominated in the Sphagnum-covered and Sasa-invaded areas, respectively. The mcrA and 16S rRNA libraries indicated that Methanoregula were abundant methanogens in the Sphagnum-covered area. In the Sasa-invaded area, by contrast, mcrA genes were not detected, and no 16S rRNA clones were affiliated with previously known methanogens. Because the Sasa-invaded area still produced methane, of the various uncultured populations detected, novel euryarchaeotal lineages are candidate methane producers.

Evaluation of rice and sugarcane SSR markers for phylogenetic and genetic diversity analyses in bamboo.

Simple sequence repeat (SSR) markers are valuable tools for many purposes such as phylogenetic, fingerprinting, and molecular breeding studies. However, only a few SSR markers are known and available in bamboo species of the tropics (Bambusa spp.). Considering that grass genomes have co-evolved and share large-scale synteny, theoretically it should be possible to use the genome sequence based SSR markers of field crops such as rice (Oryza sativa) and sugarcane (Saccharum spp.) for genome analysis in bamboo. To test this, 98 mapped SSR primers representing 12 linkage groups of rice and 20 EST-derived sugarcane SSR primers were evaluated for transferability to 23 bamboo species. Of the tested markers, 44 (44.9%) rice and 15 (75%) sugarcane SSR primers showed repeatable amplification in at least one species of bamboo and thus were successfully utilized for phylogenetic and genetic diversity analyses. Transferred SSR primers revealed complex amplification patterns in bamboo, with an average of 9.62 fragments per primer, indicating a high level of polyploidy and genetic variability in bamboo. Forty-two of these primers (34 rice and 8 sugarcane SSR primers) detected an average of 2.12 unique fragments per primer and thus could be exploited for species identification. Six bamboo SSR primers exhibited cross transferability, to varying degrees, to different bamboo species. The genetic similarity coefficient indicated a high level of divergence at the species level (73%). However, a relatively low level of diversity was observed within species (25% in 20 accessions of Dendrocalamus hamiltonii). Further, cluster analysis revealed that the major grouping was in accordance with the taxonomical classification of bamboo. Thus, the rice and sugarcane SSRs can be utilized for phylogenetic and genetic diversity studies in bamboo.

Cloning and characterization of a bamboo LEAFY HULL STERILE1 homologous gene.

A cDNA named DlMADS8 was isolated from the young spikelets of the sweet bamboo, Dendrocalamus latiflorus by rapid amplification of cDNA end (RACE). DNA sequence analysis showed that DlMADS8 was composed of full ORF and 3'UTR, but without 5'UTR. The cDNA contained 1059 nucleotides and encoded a putative protein of 244 amino acid residues. The gene displayed the structure of a typical plant MADS-box gene, which consisted of a MADS domain, K domain, a short I region, and the C-terminal region. Phylogenetic analysis of plant MADS-box genes based on amino acid sequences revealed that DlMADS8 was grouped into the AGAMOUS-LIKE 2 (AGL2)-like subfamily. It was homologous to the LEAFY HULL STERILE1 (LHS1) genes of grasses. To study the functions of it, DlMADS8 cDNA clone driven by the CaMV 35S promoter was transformed into Arabidopsis thaliana. Transgenic plants of DlMADS8 exhibited the phenotypes of curled leaves and early flowering. After bolting, three novel phenotypes related to inflorescence development were observed in different transgenic plants. No obvious homeotic conversions of floral organs were observed in all of the 35S::DllMADS8 transgenic Arabidopsis plants. These results indicated that DlMADS8 probably plays a role in floral meristem determinacy and is involved in controlling the flowering time of D. latiflorus.

[Effect of stand condition improvement on surface soil nutrient losses of hilly land Dendrocalamus latiflorus plantations].

A comparative study on the soil nutrient losses of hilly land shoot-used Dendrocalamus lati florus plantations at the Nanjing stated-owned forest farm, Fujian Province showed that the losses of soil N, P and K under three treatments of stand condition improvement varied with season, and 90% or more of these losses appeared from April to September, with the order of overall ploughing > strip ploughing > control. There was no significant difference of N and P losses among the three treatments, while the K loss was markedly higher under overall ploughing than the control.