De novo Transcriptome Analysis of Miscanthus lutarioriparius Identifies Candidate Genes in Rhizome Development

Comparison of two contrasting Leymus chinensis accessions reveals the roles of the cell wall and auxin in rhizome development

Leymus chinensis, a perennial native forage grass, is widely distributed in the steppes of Inner Mongolia as the dominant species. The main reproductive strategy of this grass is clonal propagation, which occurs via the proliferation of subterranean horizontal stems known as rhizomes. To elucidate the mechanism underlying rhizome development in this grass, we collected 60 accessions of L. chinensis and evaluated their rhizome development. One accession, which we named SR-74 (Strong Rhizomes), had significantly better rhizome development capacity than the accession WR-16 (Weak Rhizomes) in terms of rhizome number, total and primary rhizome length, and number of rhizome seedlings. Rhizome elongation was positively correlated with the number of internodes in the rhizome, which affected plant biomass. Compared to WR-16, SR-74 had higher rhizome tip hardness, higher abundance of transcripts participating in the biosynthesis of cell wall components, and higher levels of the metabolites L-phenylalanine, trans-cinnamic acid, 3-coumaric acid, ferulic acid, and coniferin. These metabolites in the phenylpropanoid biosynthesis pathway are precursors of lignin. In addition, SR-74 rhizomes contained higher amounts of auxin and auxin metabolites, including L-Trp, IPA, IBA, IAA and IAA-Asp, as well as upregulated expression of the auxin biosynthesis and signaling genes YUCCA6, YUCCA8, YUCCA10, YUCCA11, PIN1, PIN2, UGT1, UGT2, UGT4, UGT10, GH3, IAA7, IAA23, and IAA30. We propose a network between auxin signaling and the cell wall underlying rhizome development in L. chinensis.

Comparison between the Gametophyte and the Sporophyte Transcriptomes of the Endangered Fern Vandenboschia speciosa

Genomic resources are essential to understanding the evolution and functional biology of organisms. Nevertheless, generating genomic resources from endangered species may be challenging due to the scarcity of available specimens and sampling difficulties. In this study, we compare the transcriptomes of the sporophyte and the gametophyte of the endangered fern Vandenboschia speciosa. After Illumina sequencing and de novo transcriptome assembly of the gametophyte, annotation proved the existence of cross-species contamination in the gametophyte sample. Thus, we developed an in silico decontamination step for the gametophyte sequences. Once the quality check of the decontaminated reads passed, we produced a de novo assembly with the decontaminated gametophyte reads (with 43,139 contigs) and another combining the sporophyte and in silico decontaminated gametophyte reads (with 42,918 contigs). A comparison of the enriched GO terms from the top 1000 most expressed transcripts from both tissues showed that the gametophyte GO term set was enriched in sequences involved in development, response to stress, and plastid organization, while the sporophyte GO term set had a larger representation of more general metabolic functions. This study complements the available genomic resources on the life cycle of the endangered fern Vandenboschia speciosa.

Molecular Traits Underlying the Growth Promotion and Metabolite Accumulation in Rheum palmatum Inoculated with Endophytic Trichoderma citrinoviride HT-1

Trichoderma spp. are an important plant-growth-promoting fungi. Trichoderma citrinoviride HT-1 was isolated from Rheum palmatum root, which has beneficial effects on growth and metabolite accumulation. However, the improvement mechanisms for growth and metabolite accumulation of T. citrinoviride HT-1 are unclear. In this study, RNA sequencing (RNA-seq) and high-performance liquid chromatography (HPLC) were used to measure the effect of different concentrations of conidial suspension of the HT-1 strain on the growth promotion and metabolite accumulation of R. palmatum seedlings. The results showed that the highest biomass and metabolites of R. palmatum seedlings were obtained through treatment with the HT-1 strain at a final spore concentration of 107 spores/mL. RNA sequencing indicated that 1662 genes were upregulated and 2155 genes were downregulated after inoculation with 107 spores/mL of the HT-1 strain. This strain induced significant upregulation of related genes in the phenylpropanoid biosynthesis pathway, plant hormone signal transduction pathway, biosynthesis of secondary metabolites pathway, and plant–pathogen interaction pathway in R. palmatum. The gene expression trends were revealed through quantitative real-time polymerase chain reaction (qRT-PCR) and were consistent with those determined by RNA-seq. Our results will help us to understand the growth-promoting mechanisms of the HT-1 strain on R. palmatum and provide a theoretical basis for the application of T. citrinoviride HT-1 as a biological fertilizer.

Unveiling the transcriptomic complexity of Miscanthus sinensis using a combination of PacBio long read- and Illumina short read sequencing platforms

Background

Miscanthus sinensis Andersson is a perennial grass that exhibits remarkable lignocellulose characteristics suitable for sustainable bioenergy production. However, knowledge of the genetic resources of this species is relatively limited, which considerably hampers further work on its biology and genetic improvement.

Results

In this study, through analyzing the transcriptome of mixed samples of leaves and stems using the latest PacBio Iso-Seq sequencing technology combined with Illumina HiSeq, we report the first full-length transcriptome dataset of M. sinensis with a total of 58.21 Gb clean data. An average of 15.75 Gb clean reads of each sample were obtained from the PacBio Iso-Seq system, which doubled the data size (6.68 Gb) obtained from the Illumina HiSeq platform. The integrated analyses of PacBio- and Illumina-based transcriptomic data uncovered 408,801 non-redundant transcripts with an average length of 1,685 bp. Of those, 189,406 transcripts were commonly identified by both methods, 169,149 transcripts with an average length of 619 bp were uniquely identified by Illumina HiSeq, and 51,246 transcripts with an average length of 2,535 bp were uniquely identified by PacBio Iso-Seq. Approximately 96 % of the final combined transcripts were mapped back to the Miscanthus genome, reflecting the high quality and coverage of our sequencing results. When comparing our data with genomes of four species of Andropogoneae, M. sinensis showed the closest relationship with sugarcane with up to 93 % mapping ratios, followed by sorghum with up to 80 % mapping ratios, indicating a high conservation of orthologs in these three genomes. Furthermore, 306,228 transcripts were successfully annotated against public databases including cell wall related genes and transcript factor families, thus providing many new insights into gene functions. The PacBio Iso-Seq data also helped identify 3,898 alternative splicing events and 2,963 annotated AS isoforms within 10 function categories.

Conclusions

Taken together, the present study provides a rich data set of full-length transcripts that greatly enriches our understanding of M. sinensis transcriptomic resources, thus facilitating further genetic improvement and molecular studies of the Miscanthus species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07971-x.

De Novo Sporophyte Transcriptome Assembly and Functional Annotation in the Endangered Fern Species Vandenboschia speciosa (Willd.) G. Kunkel

We sequenced the sporophyte transcriptome of Killarney fern (Vandenboschia speciosa (Willd.) G. Kunkel). In addition to being a rare endangered Macaronesian-European endemism, this species has a huge genome (10.52 Gb) as well as particular biological features and extreme ecological requirements. These characteristics, together with the systematic position of ferns among vascular plants, make it of high interest for evolutionary, conservation and functional genomics studies. The transcriptome was constructed de novo and contained 36,430 transcripts, of which 17,706 had valid BLAST hits. A total of 19,539 transcripts showed at least one of the 7362 GO terms assigned to the transcriptome, whereas 6547 transcripts showed at least one of the 1359 KEGG assigned terms. A prospective analysis of functional annotation results provided relevant insights on genes involved in important functions such as growth and development as well as physiological adaptations. In this context, a catalogue of genes involved in the genetic control of plant development, during the vegetative to reproductive transition, in stress response as well as genes coding for transcription factors is given. Altogether, this study provides a first step towards understanding the gene expression of a significant fern species and the in silico functional and comparative analyses reported here provide important data and insights for further comparative evolutionary studies in ferns and land plants in general.

Organ-Specific Transcriptome Analysis Identifies Candidate Genes Involved in the Stem Specialization of Bermudagrass (Cynodon dactylon L.)

As an important warm-season turfgrass and forage grass species with wide applications, bermudagrass (Cynodon dactylon L.) simultaneously has shoot, stolon and rhizome, three types of stems with different physiological functions. To better understand how the three types of stems differentiate and specialize, we generated an organ-specific transcriptome dataset of bermudagrass encompassing 114,169 unigenes, among which 100,878 and 65,901 could be assigned to the Kyoto Encyclopedia of Genes and Genomes (KEGG) and the Gene Ontology (GO) terms, respectively. Using the dataset, we comprehensively analyzed the gene expression of different organs, especially the shoot, stolon and rhizome. The results indicated that six organs of bermudagrass all contained more than 52,000 significantly expressed unigenes, however, only 3,028 unigenes were enrich-expressed in different organs. Paired comparison analyses further indicated that 11,762 unigenes were differentially expressed in the three types of stems. Gene enrichment analysis revealed that 39 KEGG pathways were enriched with the differentially expressed unigenes (DEGs). Specifically, 401 DEGs were involved in plant hormone signal transduction, whereas 1,978 DEGs were transcription factors involved in gene expression regulation. Furthermore, in agreement with the starch content and starch synthase assay results, DEGs encoding starch synthesis-related enzymes all showed the highest expression level in the rhizome. These results not only provided new insights into the specialization of stems in bermudagrass but also made solid foundation for future gene functional studies in this important grass species and other stoloniferous/rhizomatous plants.

Proteome analysis of bermudagrass stolons and rhizomes provides new insights into the adaptation of plant stems to aboveground and underground growth

As an important perennial warm-season turfgrass species, bermudagrass (Cynodon dactylon L.) forms underground-growing rhizomes and aboveground-growing stolons simultaneously, making it a fast propagating clonal plant with strong regeneration ability. In the current study, we compared the internode proteomes of rhizomes and stolons at the same developmental stage in the bermudagrass cultivar Yangjiang using iTRAQ. The results indicated that 228 protein species were differentially accumulated in the two specialized stems. In agreement with the different contents of starch, chlorophyll, anthocyanin and H₂O₂ in the two types of stems, photosynthesis and flavonoid biosynthesis were enriched with differentially accumulated protein species (DAPs) in stolons, whereas starch and sucrose metabolism, glycolysis, and H₂O₂ metabolism were enriched with DAPs in rhizomes. Burying stolons in the soil resulted in the gradual degradation of chlorophyll and anthocyanin, accumulation of starch, and increment of H₂O₂, which is similar to the physiological characteristics of rhizomes. These results collectively revealed that stolons and rhizomes of bermudagrass have significant differences at the proteome level and light might play important regulatory roles in the discrepancy of the proteome profiles and specialization of the two stems, providing new insights into the adaptation of plant stems to aboveground and underground growth. BIOLOGICAL SIGNIFICANCE: As two types of specialized stems that grow underground and aboveground respectively, rhizomes and stolons play important roles in overwintering and ecological invasion of many perennial and clonal plants. However, because rhizomes and stolons rarely coexist in single plant species, the differences between the two stems remain unclear at the molecular level. In this study, through an iTRAQ comparative proteomic analysis, we reported the identification of 228 differentially accumulated protein species (DAPs) in rhizomes and stolons of bermudagrass for the first time. We found that the 228 DAPs were interconnected to form protein networks in regulating diverse cellular activities and biochemical reactions. We also observed that stolons growing underground showed similar physiological activities and DAP expression as those of underground-growing rhizomes, suggesting that light might play important regulatory roles in the specialization of stolons and rhizomes. These results expanded our understanding of the mysterious adaption of plant stems to different growth conditions.

Revealing the full-length transcriptome of caucasian clover rhizome development

BACKGROUND

Caucasian clover (Trifolium ambiguum M. Bieb.) is a strongly rhizomatous, low-crowned perennial leguminous and ground-covering grass. The species may be used as an ornamental plant and is resistant to cold, arid temperatures and grazing due to a well-developed underground rhizome system and a strong clonal reproduction capacity. However, the posttranscriptional mechanism of the development of the rhizome system in caucasian clover has not been comprehensively studied. Additionally, a reference genome for this species has not yet been published, which limits further exploration of many important biological processes in this plant.

RESULT

We adopted PacBio sequencing and Illumina sequencing to identify differentially expressed genes (DEGs) in five tissues, including taproot (T1), horizontal rhizome (T2), swelling of taproot (T3), rhizome bud (T4) and rhizome bud tip (T5) tissues, in the caucasian clover rhizome. In total, we obtained 19.82 GB clean data and 80,654 nonredundant transcripts were analysed. Additionally, we identified 78,209 open reading frames (ORFs), 65,227 coding sequences (CDSs), 58,276 simple sequence repeats (SSRs), 6821 alternative splicing (AS) events, 2429 long noncoding RNAs (lncRNAs) and 4501 putative transcription factors (TFs) from 64 different families. Compared with other tissues, T5 exhibited more DEGs, and co-upregulated genes in T5 are mainly annotated as involved in phenylpropanoid biosynthesis. We also identified betaine aldehyde dehydrogenase (BADH) as a highly expressed gene-specific to T5. A weighted gene co-expression network analysis (WGCNA) of transcription factors and physiological indicators were combined to reveal 11 hub genes (MEgreen-GA3), three of which belong to the HB-KNOX family, that are up-regulated in T3. We analysed 276 DEGs involved in hormone signalling and transduction, and the largest number of genes are associated with the auxin (IAA) signalling pathway, with significant up-regulation in T2 and T5.

CONCLUSIONS

This study contributes to our understanding of gene expression across five different tissues and provides preliminary insight into rhizome growth and development in caucasian clover.

Differential regulatory pathways associated with drought-inhibition and post-drought recuperation of rhizome development in perennial grass

BACKGROUND AND AIMS

Rhizomes are key organs for the establishment of perennial grass stands and adaptation to environmental stress. However, mechanisms regulating rhizome initiation and elongation under drought stress and during post-drought recovery remain unclear. The objective of this study is to investigate molecular factors and metabolic processes involved in drought effects and post-drought recovery in rhizome growth in perennial grass species by comparative transcriptomic and proteomic profiling.

METHODS

Tall fescue (Festuca arundinacea) (B-type rhizome genotype, 'BR') plants were exposed to drought stress and re-watering in growth chambers. The number and length of rhizomes were measured following drought stress and re-watering. Hormone and sugar contents were analysed, and transcriptomic and proteomic analyses were performed to identify metabolic factors, genes and proteins associated with rhizome development.

KEY RESULTS

Rhizome initiation and elongation were inhibited by drought stress, and were associated with increases in the contents of abscisic acid (ABA) and soluble sugars, but declines in the contents of indoleacetic acid (IAA), zeatin riboside (ZR) and gibberellin (GA4). Genes involved in multiple metabolic processes and stress defence systems related to rhizome initiation exhibited different responses to drought stress, including ABA signalling, energy metabolism and stress protection. Drought-inhibition of rhizome elongation could be mainly associated with the alteration of GA4 and antioxidants contents, energy metabolism and stress response proteins. Upon re-watering, new rhizomes were regenerated from rhizome nodes previously exposed to drought stress, which was accompanied by the decline in ABA content and increases in IAA, ZR and GA4, as well as genes and proteins for auxin, lipids, lignin and nitrogen metabolism.

CONCLUSIONS

Drought-inhibition of rhizome initiation and elongation in tall fescue was mainly associated with adjustments in hormone metabolism, carbohydrate metabolism and stress-defence systems. Rhizome regeneration in response to re-watering involved reactivation of hormone and lipid metabolism, secondary cell-wall development, and nitrogen remobilization and cycling.

Characterization of rhizome transcriptome and identification of a rhizomatous ER body in the clonal plant Cardamine leucantha

The rhizome is a plant organ that develops from a shoot apical meristem but penetrates into belowground environments. To characterize the gene expression profile of rhizomes, we compared the rhizome transcriptome with those of the leaves, shoots and roots of a rhizomatous Brassicaceae plant, Cardamine leucantha. Overall, rhizome transcriptomes were characterized by the absence of genes that show rhizome-specific expression and expression profiles intermediate between those of shoots and roots. Our results suggest that both endogenous developmental factors and external environmental factors are important for controlling the rhizome transcriptome. Genes that showed relatively high expression in the rhizome compared to shoots and roots included those related to belowground defense, control of reactive oxygen species and cell elongation under dark conditions. A comparison of transcriptomes further allowed us to identify the presence of an ER body, a defense-related belowground organelle, in epidermal cells of the C. leucantha rhizome, which is the first report of ER bodies in rhizome tissue.

Comparative analysis of sugarcane root transcriptome in response to the plant growth-promoting Burkholderia anthina MYSP113

The diazotrophic Burkholderia anthina MYSP113 is a vital plant growth-promoting bacteria and sugarcane root association. The present study based on a detailed analysis of sugarcane root transcriptome by using the HiSeq-Illumina platform in response to the strain MYSP113. The bacterium was initially isolated from the rhizosphere of sugarcane. To better understand biological, cellular, and molecular mechanisms, a de novo transcriptomic assembly of sugarcane root was performed. HiSeq-Illumina platformwas employed for the sequencing of an overall of 16 libraries at a 2×100 bp configuration. Differentially expressed genes (DEGs) analysis identified altered gene expression in 370 genes (total of 199 up-regulated genes and 171 down-regulated genes). Deciphering the huge datasets, concerning the functioning and production of biological systems, a high throughput genome sequencing analysis was attempted with Gene ontology functional analyses and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. The report revealed a total of 148930 unigenes. 70414 (47.28%) of them were annotated successfully to Gene Ontology (GO) terms. 774 at 45 days, 4985 of 30 days and 15 days of 6846 terms were significantly regulated. GO analysis revealed that many genes involved in the metabolic, oxidation-reduction process and biological regulatory processes in response to strain MYSP113 and significantly enriched as compare to the control. Moreover, KEGG enriched results show that differentially expressed genes were classified into different pathway categories involved in various processes, such as nitrogen metabolism, plant hormone signal transduction, etc. The sample correlation analyses could help examine the similarity at the gene expression level. The reliability of the observed differential gene expression patterns was validated with quantitative real-time PCR (qRT-PCR). Additionally, plant enzymes activities such as peroxidase and superoxide dismutase were significantly increased in plant roots after the inoculation of strain MYSP113. The results of the research may help in understanding the plant growth-promoting rhizobacteria and plant interaction.

Identification of candidate reference genes for quantitative RT-PCR in Miscanthus sinensis subjected to various abiotic stresses

Quantitative real-time PCR (qRT-PCR) has been widely used for studying gene expression at the transcript level. Its accuracy usually relies on the reference genes that are utilized for data normalization. Miscanthus sinensis, a perennial C4 grass with high biomass and strong resistance to adversities, is often utilized as a high value energy crop. However, no reliable reference genes have been investigated for normalizing gene expression for this species. In this study, 12 candidate reference genes were selected to identify their stability under five different abiotic stress treatments (drought, salt, cadmium, chromium and arsenic) by using geNorm, NormFinder, BestKeeper and RefFinder softwares. The results showed that 18S rRNA and Unigene33312 were the best reference genes under drought treatments. Unigene33312 and Unigene33024 were found to be the most stably expressed genes under salt stress and Cd stress. Moreover, Unigene33024 and PP2A were the most suitable reference genes under Cr stress and Unigene33024 and Sb09g019750 were deemed more suitable reference genes under As stress. In total, considering all the samples, Unigene33024 and PP2A were the most stable genes while ACTIN and Unigene26576 were the least stable reference genes for internal control. The expression patterns of two target genes (Cu/Zn SOD and CAT) were used to further verify those selected reference genes under different conditions. The results showed that the most and the least stable reference genes had clearly different expression patterns. This work comprehensively estimated the stability of reference genes in M. sinensis which may give insight to the reference genes selection in other tissues as well as other related varieties. These suggested reference genes would assist in further putative gene expression validation in M. sinensis.

Integrated mRNA and miRNA transcriptome analysis reveals a regulatory network for tuber expansion in Chinese yam (Dioscorea opposita)

BACKGROUND

Yam tuber is a storage organ, derived from the modified stem. Tuber expansion is a complex process, and depends on the expressions of genes that can be influenced by environmental and endogenous factors. However, little is known about the regulatory mechanism of tuber expansion. In order to identify the genes and miRNAs involved in tuber expansion, we examined the mRNAs and small RNAs in Dioscorea opposita (Chinese yam) cv. Guihuai 16 tuber during its initiation and expansion stages.

RESULTS

A total of 14,238 differentially expressed genes in yam tuber at its expansion stage were identified by using RNA sequencing technology. Among them, 5723 genes were up-regulated, and 8515 genes were down-regulated. Functional analysis revealed the coordination of tuber plant involved in processes of cell events, metabolism, biosynthesis, and signal transduction pathways at transcriptional level, suggesting that these differentially expressed genes are somehow involved in response to tuber expansion, including CDPK, CaM, CDL, SAUR, DELLA, SuSy, and expansin. In addition, 541 transcription factor genes showed differential expression during the expansion stage at transcriptional level. MADS, bHLH, and GRAS were involved in cell differentiation, division, and expansion, which may relate to tuber expansion. Noteworthy, data analysis revealed that 22 known tuber miRNAs belong to 10 miRNA families, and 50 novel miRNAs were identified. The integrated analysis of miRNA-mRNA showed that 4 known miRNAs and 11 genes formed 14 miRNA-target mRNA pairs were co-expressed in expansion stage. miRNA160, miRNA396, miRNA535 and miRNA5021 may be involved in complex network to regulate cell division and differentiation in yam during its expansion stage.

CONCLUSION

The mRNA and miRNA datasets presented here identified a subset of candidate genes and miRNAs that are putatively associated with tuber expansion in yam, a hypothetical model of genetic regulatory network associated with tuber expansion in yam was put forward, which may provide a foundation for molecular regulatory mechanism researching on tuber expansion in Dioscorea species.

Integrated mRNA and miRNA transcriptome analysis reveals a regulatory network for tuber expansion in Chinese yam (Dioscorea opposita)

Background

Yam tuber is a storage organ, derived from the modified stem. Tuber expansion is a complex process, and depends on the expressions of genes that can be influenced by environmental and endogenous factors. However, little is known about the regulatory mechanism of tuber expansion. In order to identify the genes and miRNAs involved in tuber expansion, we examined the mRNAs and small RNAs in Dioscorea opposita (Chinese yam) cv. Guihuai 16 tuber during its initiation and expansion stages.

Results

A total of 14,238 differentially expressed genes in yam tuber at its expansion stage were identified by using RNA sequencing technology. Among them, 5723 genes were up-regulated, and 8515 genes were down-regulated. Functional analysis revealed the coordination of tuber plant involved in processes of cell events, metabolism, biosynthesis, and signal transduction pathways at transcriptional level, suggesting that these differentially expressed genes are somehow involved in response to tuber expansion, including CDPK, CaM, CDL, SAUR, DELLA, SuSy, and expansin. In addition, 541 transcription factor genes showed differential expression during the expansion stage at transcriptional level. MADS, bHLH, and GRAS were involved in cell differentiation, division, and expansion, which may relate to tuber expansion. Noteworthy, data analysis revealed that 22 known tuber miRNAs belong to 10 miRNA families, and 50 novel miRNAs were identified. The integrated analysis of miRNA-mRNA showed that 4 known miRNAs and 11 genes formed 14 miRNA-target mRNA pairs were co-expressed in expansion stage. miRNA160, miRNA396, miRNA535 and miRNA5021 may be involved in complex network to regulate cell division and differentiation in yam during its expansion stage.

Conclusion

The mRNA and miRNA datasets presented here identified a subset of candidate genes and miRNAs that are putatively associated with tuber expansion in yam, a hypothetical model of genetic regulatory network associated with tuber expansion in yam was put forward, which may provide a foundation for molecular regulatory mechanism researching on tuber expansion in Dioscorea species.

Integrated transcriptome and miRNA analysis uncovers molecular regulators of aerial stem-to-rhizome transition in the medical herb Gynostemma pentaphyllum

BACKGROUND

Gynostemma pentaphyllum is an important perennial medicinal herb belonging to the family Cucurbitaceae. Aerial stem-to-rhizome transition before entering the winter is an adaptive regenerative strategy in G. pentaphyllum that enables it to survive during winter. However, the molecular regulation of aerial stem-to-rhizome transition is unknown in plants. Here, integrated transcriptome and miRNA analysis was conducted to investigate the regulatory network of stem-to-rhizome transition.

RESULTS

Nine transcriptome libraries prepared from stem/rhizome samples collected at three stages of developmental stem-to-rhizome transition were sequenced and a total of 5428 differentially expressed genes (DEGs) were identified. DEGs associated with gravitropism, cell wall biosynthesis, photoperiod, hormone signaling, and carbohydrate metabolism were found to regulate stem-to-rhizome transition. Nine small RNA libraries were parallelly sequenced, and seven significantly differentially expressed miRNAs (DEMs) were identified, including four known and three novel miRNAs. The seven DEMs targeted 123 mRNAs, and six pairs of miRNA-target showed significantly opposite expression trends. The GpmiR166b-GpECH2 module involved in stem-to-rhizome transition probably promotes cell expansion by IBA-to-IAA conversion, and the GpmiR166e-GpSGT-like module probably protects IAA from degradation, thereby promoting rhizome formation. GpmiR156a was found to be involved in stem-to-rhizome transition by inhibiting the expression of GpSPL13A/GpSPL6, which are believed to negatively regulate vegetative phase transition. GpmiR156a and a novel miRNA Co.47071 co-repressed the expression of growth inhibitor GpRAV-like during stem-to-rhizome transition. These miRNAs and their targets were first reported to be involved in the formation of rhizomes. In this study, the expression patterns of DEGs, DEMs and their targets were further validated by quantitative real-time PCR, supporting the reliability of sequencing data.

CONCLUSIONS

Our study revealed a comprehensive molecular network regulating the transition of aerial stem to rhizome in G. pentaphyllum. These results broaden our understanding of developmental phase transitions in plants.

Transcriptomics and proteomics reveal genetic and biological basis of superior biomass crop Miscanthus

Miscanthus is a rhizomatous C4 grass which is considered as potential high-yielding energy crop with the low-nutrient requirements, high water-use efficiency, and capability of C mitigation. To better understand the genetic basis, an integrative analysis of the transcriptome and proteome was performed to identify important genes and pathways involved in Miscanthus leaves. At the transcript level, 64,663 transcripts in M. lutarioriparius, 97,043 in M. sacchariflorus, 97,043 in M. sinensis, 67,323 in M. floridulus and 70,021 in M. × giganteus were detected by an RNA sequencing approach. At the protein level, 1964 peptide-represented proteins were identified and 1933 proteins differed by 1.5-fold or more in their relative abundance, as indicated by iTRAQ (isobaric tags for relative and absolute quantitation) analysis. Phylogenies were constructed from the nearly taxa of Miscanthus. A large number of genes closely related to biomass production were found. And SSR markers and their corresponding primers were derived from Miscanthus transcripts and 90% of them were successfully detected by PCR amplification among Miacanthus species. These similarities and variations on the transcriptional and proteomic level between Miscanthus species will serve as a resource for research in Miscanthus and other lignocellulose crops.