A detailed gene expression study of the Miscanthusgenus reveals changes in the transcriptome associated with the rejuvenation of spring rhizomes
- Adam Barling†1, 2,
- Kankshita Swaminathan†1,
- Therese Mitros3, 4,
- Brandon T James1, 2,
- Juliette Morris1,
- Ornella Ngamboma1,
- Megan C Hall3, 7,
- Jessica Kirkpatrick1, 2,
- Magdy Alabady1, 8,
- Ashley K Spence5, 9,
- Matthew E Hudson1, 2,
- Daniel S Rokhsar3, 4, 6 and
- Stephen P Moose1, 2Email author
© Barling et al.; licensee BioMed Central Ltd. 2013
Received: 25 January 2013
Accepted: 4 December 2013
Published: 9 December 2013
The Miscanthus genus of perennial C4 grasses contains promising biofuel crops for temperate climates. However, few genomic resources exist for Miscanthus, which limits understanding of its interesting biology and future genetic improvement. A comprehensive catalog of expressed sequences were generated from a variety of Miscanthus species and tissue types, with an emphasis on characterizing gene expression changes in spring compared to fall rhizomes.
Illumina short read sequencing technology was used to produce transcriptome sequences from different tissues and organs during distinct developmental stages for multiple Miscanthus species, including Miscanthus sinensis, Miscanthus sacchariflorus, and their interspecific hybrid Miscanthus × giganteus. More than fifty billion base-pairs of Miscanthus transcript sequence were produced. Overall, 26,230 Sorghum gene models (i.e., ~ 96% of predicted Sorghum genes) had at least five Miscanthus reads mapped to them, suggesting that a large portion of the Miscanthus transcriptome is represented in this dataset. The Miscanthus × giganteus data was used to identify genes preferentially expressed in a single tissue, such as the spring rhizome, using Sorghum bicolor as a reference. Quantitative real-time PCR was used to verify examples of preferential expression predicted via RNA-Seq. Contiguous consensus transcript sequences were assembled for each species and annotated using InterProScan. Sequences from the assembled transcriptome were used to amplify genomic segments from a doubled haploid Miscanthus sinensis and from Miscanthus × giganteus to further disentangle the allelic and paralogous variations in genes.
This large expressed sequence tag collection creates a valuable resource for the study of Miscanthus biology by providing detailed gene sequence information and tissue preferred expression patterns. We have successfully generated a database of transcriptome assemblies and demonstrated its use in the study of genes of interest. Analysis of gene expression profiles revealed biological pathways that exhibit altered regulation in spring compared to fall rhizomes, which are consistent with their different physiological functions. The expression profiles of the subterranean rhizome provides a better understanding of the biological activities of the underground stem structures that are essentials for perenniality and the storage or remobilization of carbon and nutrient resources.
Miscanthus is a perennial C4 grass that belongs to the Andropogoneae tribe within the Poaceae family, which includes important agricultural crops for food and fuel such as sugarcane, sorghum, and maize. Following their introduction into the Western world in the 1930s , members of the Miscanthus genus are now grown as ornamental crops in many regions of the United States due to their characteristically robust growth and attractive late-season inflorescence.
The Miscanthus genus consists of approximately fifteen species, most of which are either diploids or tetraploids . The grass is an obligate outcrosser with a large, highly repetitive 2.5 Gbp (giga base pairs) genome that is distributed across nineteen chromosomes [3, 4]. Natural hybridization events between the two most predominant Miscanthus species, M. sinensis and M. sacchariflorus, have been reported [5, 6]. Ribosomal DNA evidence suggests that the large statured, cold tolerant, sterile triploid hybrid M. × giganteus (3n = 57) is the result of a natural hybridization event between a diploid M. sinensis (2n = 38) and a tetraploid M. sacchariflorus (4n = 76) [2, 4, 7].
Plants of the Miscanthus genus, especially Miscanthus × giganteus, have generated interest as a source of lignocellulosic biomass for the bioenergy industry. Although Miscanthus has been of horticultural interest for some time, it essentially remains a genus of wild species. Genetic selections for the genus have largely concentrated on traits desirable to the horticultural and landscaping industry; there have been few focused breeding efforts targeting traits that would enhance the potential of Miscanthus as a perennial bioenergy feedstock. The availability of molecular tools for Miscanthus will accelerate improvement of biofuel-centric traits in Miscanthus. Recent advances in Miscanthus genomics have enabled the construction of complete genetic maps for M. sinensis[8–10]. These genetic maps revealed a recent allotetraploidization event in Miscanthus in which pairs of homeologous chromosomes show extensive synteny to the Sorghum bicolor genome, with a single chromosome fusion accounting for the nineteen linkage groups.
Deep sequencing technologies applied to gene discovery through transcriptome sequencing has efficiently increased genetic information for many non-model plant organisms such as barley, grape, wheat, and lodgepole pine [11–15]. Importantly, the high degree of sequence similarity and genome organization between Miscanthus and Sorghum make Sorghum bicolor a suitable reference genome sequence for the analysis of the Miscanthus transcriptome [4, 9, 10]. A preliminary study of dormant Miscanthus × giganteus rhizomes was used to assess variation among available Miscanthus × giganteus accessions , but a comprehensive catalog of expressed sequences in the Miscanthus genus is not yet available. We report here high-depth sequencing of expressed mRNAs from a variety of M. × giganteus tissues as well as multiple accessions of M. sinensis and one accession of M. sacchariflorus. The data generated enable a robust assembly of the Miscanthus transcriptome with demonstrated utility in the analysis of changes in gene expression and evolution of genic sequences within the genus.
Results and discussion
Sequencing the Miscanthustranscriptome
Miscanthus RNA-Seq libraries sequenced for this study
(Billion base pairs)
Miscanthus x giganteus ‘Illinois clone’
RO2, RZ1 and 3, RB2, ES1 and 2, VA1, SA1, ST2, PA2, II1, MI2, ML1, FR3
Miscanthus sacchariflorus ‘Golf Course’
Miscanthus sinensis ‘White Kaskade’
Miscanthus sinensis ‘Goliath’
Miscanthus sinensis ‘Amur Silvergrass’
Miscanthus sinensis ‘Grosse Fontaine’
Miscanthus sinensis ‘Undine’
Miscanthus sinensis ‘Zebrinus’
Tissue specific expression profile of the Miscanthus × giganteus transcriptome using the Sorghumgenome as a reference
The M. × giganteus tissues were sequenced in two separate Illumina short-read sequencing runs, both to assemble the Miscanthus transcriptome (Table 1, Figure 1A) and to identify genes preferentially expressed in a single M. × giganteus tissue-type. Approximately ten million reads were obtained for each tissue. Although Miscanthus does not currently have a completed genome the high nucleotide identity of Miscanthus to Sorghum suggests that the Sorghum genome can be used as a suitable reference for profiling tissue specific transcript expression in Miscanthus.
When expression profiles for each library are subjected to hierarchical clustering, the libraries tend to group primarily by organ type (Figure 1B). However, because some libraries were sequenced at different read lengths (36 versus 76-bp, Table 1), relative mapping efficiencies to the Sorghum reference could contribute to apparent relationships among libraries. We assessed this directly by three analyses. First, Figure 2A shows that libraries sequenced to 36-bp produced approximately half the proportion of reads mapping to the Sorghum reference compared to the 76-bp libraries. Second, when the number of Sorghum gene models with a minimum of five matching reads are compared among libraries from similar tissues, a substantial number of gene models appear to be uniquely represented in only one library (Figure 2C and D). This observation is particularly noteworthy for the comparison of Emerging Shoot (1, 36-bp reads) and Emerging Shoot (2, 76-bp reads), where the same RNA sample was used to independently construct two libraries. Although many Sorghum gene models were sampled at read depths greater than 10, a substantial number show lesser depth (Figure 2E). It is predominately these low-coverage gene models that account for the apparent differences among closely related (e.g. Vegetative Shoot Apex and Sub-Apex Shoot) or identical (e.g. Emerging shoots) RNA samples (Figure 2F).
While most transcripts are ubiquitously expressed in all tissues, transcripts that are differentially expressed yet abundant in at least one tissue are interesting as markers for developmental programs or tissue-specific biology. The Rank Products (RP) method [17, 18] is a useful non-parametric test to evaluate the significance of differential expression by a series of fold change comparisons. Rankings arise from consistencies in fold change differences between samples; as such, a series of pairwise comparisons for each individual tissue against the rest of the sequenced tissues in our study identifies high-ranking transcripts that are preferentially expressed in a tissue compared to the rest. The RP method has been used recently to help develop expression profiles for plants such as soybean , aspen trees , and the study of hormonal responses in Arabidopsis. We employed RP to identify genes preferentially expressed in one particular tissue compared to the other sampled tissues, i.e. the “rest of the plant” (Additional file 1).
The highly ranked genes from this analysis included many whose expression is known to be associated with biological processes that occur primarily in one of the sampled tissues. Examples include photosynthetic genes like phosphoenolpyruvate carboxylase (PEPC) and pyruvate orthophosphate dikinase (PPDK) in Mature Leaf, genes involved in floral organ development like APETALA3 and PISTILLATA in the Inflorescence samples, and regulators of flowering like APETALA1 in the Pre-Flowering Apex [22–27] (Additional file 1). Overall, we believe that we have generated a good repertoire of gene expression in Miscanthus for a number of stages and tissues. The primary appeal of this information is its potential use in the future investigation of the Miscanthus genus’ unique traits and characteristics. The high rankings of genes known to be highly expressed in certain tissue types in other plant species strengthens confidence in our approach to identify genes preferentially expressed in lesser-studied organs such as the subterranean rhizome; thus, we choose to focus our validation experiments on genes preferentially expressed in the Spring Rhizome and associated organs (Rhizome Buds, Emerging Shoot and Root).
As no housekeeping genes have been tested or verified for use in M. × giganteus, five potential control candidates were deduced from the Rank Product data. These potential control candidates contained Sorghum gene models with near-equal RPKM (Reads per Kilo-base per Million) values in each of the five tested tissues used in this verification. From these five candidates, the two best-performing gene models (in terms of amplification efficiency via RT-qPCR and closest-to-equivalent expression) were chosen as control genes for this study.
The RT-qPCR results correlated well with the expression patterns estimated by the RNA-Seq analysis (Figure 3), confirming that the expression variation observed from RNA-Seq provides a good representation of changes in transcript profiles among samples. Occasionally, gene expression for the root tissue appeared higher in the RT-qPCR. We attribute this discrepancy to the differences in the growth conditions for the root tissues sampled for RNA-Seq and RT-qPCR. The RNA-Seq library was prepared from roots of greenhouse plants grown in Turface, whereas the RT-qPCR analysis was performed with root tissue harvested from the same long-standing M. × giganteus field plot from which the majority of other tissue samples were obtained. In addition to the aforementioned tests, two additional leaf specific genes were assayed and both methods showed consistent results (Figure 3).
Seasonal transcription responses in Miscanthus × giganteusrhizomes
We noticed that a number of the genes identified by the Rank Product analysis as preferentially expressed in Spring Rhizomes were annotated with functions associated with the biosynthesis or signaling of plant hormones. Such pathways might be expected to be highly active in rejuvenating rhizomes. To assess this hypothesis more directly, we obtained biological replicate samples from rhizomes harvested in both Spring (May 5) and Fall (October 29) during the 2012 growing season and used RNA-Seq for transcript profiling (Additional file 2). A Gene Ontology analysis of these samples shows an enrichment in Spring Rhizomes of transcripts associated with cell wall biogenesis, root development, and both the biogenesis and signaling of jasmonic acid (Additional file 3). These findings confirm observations from the initial Rank Product analysis. In contrast, rhizomes collected in late fall show an enrichment of transcripts associated with seed maturation and dormancy. Overall, the upregulation of hormonal signaling in the spring and dormancy in the fall is consistent with seasonal changes in the physiological functions of rhizomes.
Quantitative trait analysis in a Sorghum bicolor by Sorghum propinqum population identified a 15 MB interval on Sorghum chromosome 1 associated with rhizomatousness and cold tolerance . It is interesting to note that many genes in this interval are highly expressed in the M. × giganteus Rhizome and also are differentially expressed between Spring and Fall Rhizomes. Noteworthy among these genes are three predicted ZIM domain proteins (Sb01g033020, Sb01g045190, Sb01g045180) with homology to Arabidopsis JAZ/TIFY transcription factors associated with jasmonic acid biosynthesis and signaling (Additional file 2). Conversely the M. × giganteus homolog of Sb01g038670 is highly expressed in Fall Rhizomes (Additional file 2). Sb01g038670 encodes a putative small hydrophobic membrane protein that belongs to a low temperature and salt responsive protein family and shows similarity to Arabidopsis RCI2s and Maize PMP3s [29–33].
De-novo assembly of the short read data
The M. × giganteus genotype was formed via hybridization of M. sinensis with M. sacchariflorus. Thus, we expect that the detailed assembly produced for M. × giganteus should also be broadly useful for investigating expression variation in other Miscanthus accessions. We evaluated this in two ways. First we generated libraries from a single tissue (expanding leaves containing both mature and immature portions) for four M. sinensis accessions and then mapped the reads to either an assembly produced only from that accession or to the M. × giganteus assembly. Leaf samples clearly have a reduced representation of the full transcriptome of M. × giganteus, as evidenced by the fewer number of contigs produced and their shorter N50 (Figure 4A). This is not unexpected, as most leaf tissue reads likely come from a small number of very highly expressed genes; as a result, less abundant transcripts will be more poorly represented. Importantly, when leaf-only libraries are mapped to M. × giganteus, the proportion of mapped reads rises to the level observed for M. × giganteus onto itself (Figure 4B), suggesting that nearly all reads in the leaf libraries are in fact represented within the M. × giganteus assembly. We reasoned there might be two approaches to improve accession-specific assemblies, greater read depth of the same tissue, or the inclusion of more tissues. Figure 4B shows that more than doubling the read depth of the leaf libraries had no impact on the proportion of mapped reads (those within the green circles); however, even a single library containing a mixture of tissues (those within purple circle) sequenced at moderate depth yields accession-specific assemblies comparable to M. × giganteus. Having established that moderate depth sequencing of mixed tissues offers the best assembly, we generated such a library from M. sacchariflorus accession ‘Golf Course’ and confirmed that the M. × giganteus assembly is of sufficient quality to obtain high proportions of read-mapping for both M. sacchariflorus and M. sinensis accessions.
A practical challenge of having many closely related para-alleles in Miscanthus spp. is the propensity with which chimeric products can be generated during PCR amplification due to the aberrant pairing of incompletely amplified fragments from the para-alleles during successive PCR cycles (Additional file 4). Whereas such PCR chimeras are easy to identify with Sanger sequencing of multiple clones from PCR amplicons, less rigorous methods of genotyping polyploids based on sizing of PCR-amplified fragments (e.g., SSRs) are likely to have a high error rate due to the incidence of such artifacts.
Annotation of the Miscanthusassemblies
The Miscanthus contigs were annotated using InterProScan version 4.8 [38, 39]. Eighty-eight percent of contigs were assigned at least one annotation (ftp://ftp.jgi-psf.org/pub/JGI_data/Miscanthus/transcriptome/). The top twenty most common Gene Ontology (GO) assignments in the three main categories (Cellular Component, Molecular Function, and Biological Process) in the assembled Miscanthus transcriptome are available in Additional file 5 and provide additional evidence that we have a comprehensive collection of transcripts.
Although most repeats in the genome are silenced, it is not uncommon for some repetitive elements to show expression, particularly in actively developing tissues. Of the 269,530 Miscanthus contigs, 1,693 were annotated by InterProScan to contain one or more elements found in retrotransposons: Integrase, RNase H, Reverse Transcriptase and the gag structural protein (ftp://ftp.jgi-psf.org/pub/JGI_data/Miscanthus/transcriptome/). Three of these contigs (GrosseFontaine_TContig13633, Mxg_TContig47918 and Undine_TContig8294) contained all four polypeptides, suggesting they could potentially represent intact functional retrotransposons. To further investigate the presence of putative repetitive elements in the assembly, we compared the assembly to the Plant Repeat Database, which provides a comprehensive well-characterized list of the most common plant repeats . Less that 2% of the contigs matched the repeat database (Additional file 6A), and more than half of these contigs were residual ribosomal RNA, likely due to incomplete removal of non-poly-adenylated RNAs during the library preparation. Aside from ribosomal RNAs, the most common matches were typically unclassified retrotransposons, transposons, and MITES of the Tourist type (Additional file 6B).
The grasses of the Andropogoneae tribe—maize, Sorghum, sugarcane, and Miscanthus—are among the world’s most economically important crops. An abundance of genomic resources exist for the two annual crops in this group, maize and Sorghum. In contrast, the perennials sugarcane and Miscanthus have lagged behind, in part because of the size and complexity of their genomes. The Miscanthus transcriptome reported in this study represents a major new genomic resource for the perennial Andropogoneae and will enable comparative genomic studies that advance our understanding of perenniality in grasses.
This Miscanthus expression study provides a first glance at the transcriptome of active subterranean tissues collected during an annual seasonal cycle. It is interesting to note that these tissues show preferred expression of genes involved in jasmonic acid signaling, indole biosynthesis, auxin responses, abscisic acid pathways, and osmo-sensing. The transcripts preferentially expressed in the tissues underground suggest that changes in plant hormone pathways are associated with nutrient remobilization and growth in spring. Jasmonate synthesis and signalling appears to be particularly active in the Spring Rhizomes. Exogenous jasmonate has been shown to induce underground tubers in rhubarb, yams and potatoes, and to promote shoot and bulb formation in garlic grown via tissue culture [41–44]. It is also interesting that three ZIM/tify domain containing proteins located in the Sorghum rhizomatousness interval  are highly expressed in Spring Rhizomes while the homolog of low temperature and salt responsive protein, RCI2 [30–33, 45], in the interval is expressed in Fall Rhizomes. ZIM domain proteins are transcription factors in the jasmonic acid signaling pathway, which usually function as transcriptional repressors [46–49]. The role of jasmonate and other plant hormones in rhizome biology and nutrient cycling in Miscanthus deserves further investigation. In general, while hormones appear to rage in Spring Rhizomes, genes involved in amino acid metabolism and seed maturation are high in the Fall Rhizomes (Additional files 2 and 3).
As the transcriptome assembly presented here is based solely on short-read sequencing, there are situations where the paralogous transcripts are collapsed in regions of high similarity and are represented as separate contigs in regions of greater variation. It is apparent that longer read sequencing is required to produce transcript assemblies that consistently separate alleles from paralogs. Nevertheless, the information on gene expression in Miscanthus reported here will be valuable in exploring Miscanthus biology and aid in the further sequencing and annotation of the Miscanthus genomes.
Sample collection and processing
Tissue samples used in this study were collected either from a M. × giganteus test plot that was established in 1980 in Urbana, Illinois at the University of Illinois Turf Farm or from individual plants grown in the Plant Science Laboratory greenhouse at the University of Illinois. Specific collection information, including sampling location, tissue type, sampling time, and application are shown in Additional file 7. Root samples used in the M. × giganteus sequencing project were collected from rhizomes grown in the greenhouse in calcinated clay (Turface) in order to increase the efficiency of root-tissue sampling. The samples were flash frozen in liquid nitrogen immediately following their excision. Total RNA was extracted from a pool of ten biological replicates per tissue to curb the possible bias from one sample, using an RNA extraction protocol developed for pine . Following the manufacturer’s protocol, Dynabeads (Invitrogen catalog number 61005) were used to purify the mRNA . The yield of the mRNA was quantified with a NanoDrop Spectrophotometer ND-1000 and the quality verified on an Agilent 2100 Bioanalyzer. To ensure the highest quality possible mRNA would be used for sequencing, only samples with a 260/280 of 2 ± 0.1 and a minimum RNA integrity number of 8 were used. The libraries were made and sequenced on an Illumina Genome Analyzer IIx by the W. M. Keck Center at the University of Illinois.
For Miscanthus × giganteus, RNA from the various tissues was extracted and sequenced separately, with a minimum of one lane of short read data obtained for each tissue type. All samples were sequenced on an Illumina Genome Analyzer IIx. For Rhizome, Emerging Shoot 1, Vegetative Shoot Apex, Sub-Apex Shoot, Immature Inflorescence, and Mature Leaf, 36 bp paired end reads were obtained, whereas 76 bp paired end reads were obtained for the rest of the tissues. In the case of M. sacchariflorus ‘Golf Course,’ M. sinensis ‘White Kaskade,’ and M. sinensis ‘Goliath,’ tissues were pooled before the RNA extraction (Table 1, Additional file 7). For the rest of the M. sinensis accessions, expanding leaves containing both mature and immature tissues were sampled for RNA extraction and sequencing.
Transcriptome assembly and annotation
A total of 106 billion base pairs of sequence distributed in 767 million Illumina reads were generated (Table 1, Additional file 7, SRP023501, SRP023470, SRP017791). De novo assemblies of the raw reads were performed separately for each accession using ABySS  and Phrap version 1.080721 (Phil Green, http://www.phrap.org/) as previously described in Swaminathan, 2012 . Each contig was translated in all six open reading frames (ORFs) and re-oriented based on homology to a Sorghum gene model using BLAST, with a minimum e-value of 1E-10. If the contig showed no homology to Sorghum, the contig was reoriented based on the longest ORF. A FASTA file of the reoriented assembly is provided. The contigs were annotated using InterProScan version 4.8 [38, 39] Both the assembly and annotation files are available for download from ftp://ftp.jgi-psf.org/pub/JGI_data/Miscanthus/transcriptome/. The number of putative expressed repeats was identified based on homology to a repeat in the Plant Repeat Databases (ftp://ftp.plantbiology.msu.edu/pub/data/TIGR_Plant_Repeats/) using blastn with an E-value cutoff of 1E-6.
Clustering of the contigs with the Sorghumannotated transcriptome and sugarcane ESTs
Single linkage  was used to cluster Miscanthus sequences with S. bicolor gene models  (ftp://ftp.jgi-psf.org/pub/compgen/phytozome/v9.0/Sbicolor_v1.4/) and sugarcane gene index (http://compbio.dfci.harvard.edu/tgi/cgi-bin/tgi/gimain.pl?gudb=s_officinarum). An all by all BLAT  alignment was used to find contigs that were 95% identical for over 90% of the length of the smaller of the two contigs from the same species or were 90% identical for over 90% of the length between species were assigned to the same cluster. Clusters with more than 300 members were discarded, as they are more likely to be an artifact caused by repetitive or low-complexity sequences. Clusters (701) that only contained Miscanthus and sugarcane sequences were re-matched to the Sorghum bicolor genome using Blat . Clusters (449) that did not align to the Sorghum genome at 90% identity over 90% of the length were classified as clusters with no match to Sorghum.
Cloning and sequencing genic loci from Miscanthusspp
Eleven genes present in a single copy within Sorghum were matched to the Miscanthus transcriptome assemblies using nucleotide-nucleotide BLAST (blastn). The best match for each gene from each Miscanthus assembly was aligned using Sequencher [Gene Codes Corporation version 5.0.1] with a minimum identity cutoff of 90%. Splice junctions were identified by aligning the Miscanthus contigs to the Sorghum genome using BLAT  with minIdentity set to 98. Thirteen primer pairs were then designed using IDT’s PrimerQuest (http://www.idtdna.com/Scitools/Applications/Primerquest/), taking care to minimize SNPs and avoid splice junctions. To confirm the primers were unique, the Novoalign program (Novoalign 2.05.13 http://www.novocraft.com/main/index.php) was used to map each primer pair to the Sorghum genome. The primer sequences are available in Additional file 8.
Genomic amplified PCR products were cleaned using the QIAprep Spin Miniprep kit (Qiagen catalog # 27106) and transformed using the pGem T easy Vector System II kit (Promega catalog # A1380). A minimum of eight colonies was chosen per accession for each primer; plasmids were extracted using the QIAprep 96 Turbo Miniprep Kit (Qiagen catalog # 27191). Plasmids were Sanger sequenced from both ends by the Roy J. Carver Biotechnology Center at the University of Illinois. Sequences were trimmed and aligned to the contig from which their primers were designed using Sequencher. All sequences have been deposited in Genbank (Accession numbers KF299554 - KF299740). For the two genes shown in Figure 5, genetic diversity was increased by including additional Miscanthus species and accessions.
Sequence ends were truncated so that every sequence was the same length; where two or more sequences from the same accession shared 100% identity, they were collapsed. Contigs were then exported in FASTA format and MEGA5 (http://www.megasoftware.net/)  was used for the evolutionary analyses. The evolutionary history was inferred by using the Maximum Likelihood method based on the Hasegawa-Kishino-Yano model , with the number of bootstrap replications set to 1,000, the number of discrete gamma categories set to five, the site coverage cutoff set at 20%, and the Close-Neighbor-Interchange set as the heuristic method.
Expression Analysis of Miscanthus × giganteus
Reads were adapter-trimmed and quality controlled with Perl scripts prior to import to the CLC Genomics Workbench Version 3.7. (CLC bio 2010). Low-quality bases and bad reads were discarded from input files through the use of Trim.pl (http://wiki.bioinformatics.ucdavis.edu/index.php/Trim.pl), trimming bases with quality below 10 (phred) using windowed adaptive trimming. Reads were aligned to the unmasked Sorghum bicolor genome, with exon subfeatures included, downloaded from phytozome (ftp://ftp.jgi-psf.org/pub/compgen/phytozome/v9.0/Sbicolor_v1.4/), using the following settings: 94.4% identity, extend annotated gene regions 300 flanking residues both upstream and downstream, and only use reads with a maximum of five hits. Exon discovery was enabled with a required relative expression level of 0.2 with a minimum of ten reads of at least 50 nucleotides in length. Unique gene map counts were exported from CLC for each tissue file.
For the M. × giganteus tissue preferred expression, RPKM values were calculated based on these unique counts and subsequently used in a differential expression analysis performed via the non-parametric rank products (RP) methodology  using the Perl script provided by the authors. With the RP method, genes in each individual sample are ranked based on the gene-length normalized expression consistencies and differences observed when juxtaposed against the normalized expression of the other samplings by means of a series of pairwise comparisons. As a result, the final rankings for each sample identifies, preferentially expressed genes within a single tissue by comparing each tissue to all other tissue types with the exception of Emerging Shoot 1 and 2, which were treated as a single sample with expression values averaged between the two. Listings of RP results are provided in Additional file 1.
Three biological replicates of M. × giganteus were used for the Spring versus Fall Rhizome comparison. Reads were again mapped with CLC Genomics Workbench using identical parameters to those outlined above. In total, 23,015 out of the 27,609 S. bicolor gene models had at least one read that would map in a sample. Of these, 9,264 genes had twenty or more counts per million in at least 3 samples and were considered for the differential expression analysis using two Bioconductor packages: LIMMA and edgeR (Robinson, et al.). The LIMMA (Smyth, et al.) package was used with both FPKM (fragments per kilobase of transcript per million mapped reads) and VOOM (Law, et al.) normalization methods. A total of 3,381 genes were differentially expressed in all three methods under a false discovery rate of 0.05 and a fold change value of at least two (Additional file 2). A GO analysis was performed on the 9,264 genes using the Parametric Analysis of Gene Set Enrichment (PAGE) tool in agriGO  (Additional files 2 and 3).
RT-qPCR on genes preferentially expressed in the rhizome
Total RNA was extracted from newly collected tissue-stock of M. × giganteus Emerging Shoot, Mature Leaf, Rhizome Bud, Root, and Spring Rhizome, all of which were sampled in April and May of 2011 from three dissimilar locations at the University of Illinois Turf Farm. Primers were designed for nine genes preferentially expressed in the rhizome according to the rank product analysis (Additional file 1, Additional file 8). For controls, five genes with near-equal RPKM expression values in each of the five sampled tissues were chosen. In addition, two primer sets for genes with known preferential leaf expression were added to this study (Additional file 1, Additional file 8). The primers were evaluated for amplification efficiency using the LightCycler Software package (ver. 220.127.116.11) on a Roche LightCycler 480. Five of the nine primer pairs designed to rhizome-preferred genes (Sb07g004190, Sb01g005150, Sb04g025430, Sb10g022200 and Sb03g043280), both the leaf genes (Sb09g028720 and Sb10g028120), and two of the controls genes (Sb09g019750 and Sb02g041180) had an amplification efficiency of 2 ± 0.1 and were chosen for RT-qPCR. As the other four of the nine primer pairs designed to rhizome-preferred genes did not possess adequate amplification efficiency, likely due to non-specific amplification, they could not be used effectively in RT-qPCR and were therefore discarded.
RT-qPCR was performed using four technical replicates and three biological replicates for every sampled tissue on a Roche LightCycler 480. Gene expression was determined by exporting data from the LightCycler Software package (ver. 18.104.22.168) into Microsoft Excel and performing a relative gene expression analysis using the ΔΔCt method .
Data access and visualization
Funding was provided by the Energy Biosciences Institute to SPM, MEH, DSR, AB, KS, TM, BTJ, JM, ON, MCH, JK, and MA. Some of the RNA-Seq data used was funded by Award DE-SC0005433 from the joint Department of Energy and USDA Feedstock Genomics program. We would like to thank Dr. Thomas Voigt and members of the EBI Agronomy Program for growing and maintaining the accessions used in this study. We would also like to thank Won Byoung Chae and John A. Juvik for the greenhouse-grown M. sinensis leaf tissue samples and Katarzyna Glowacka and Stanislaw Jezowski, from the Institute of Plant Genetics, Polish Academy of Sciences, Strzeszyńska 34, 60–479 Poznań, Poland, for leaf tissue from the M. sinensis double haploid lines (DH1 and DH2) and their parents (DH1P and DH2P). We thank Alvaro Hernandez and the UIUC Keck Center for Illumina RNA sequencing, Kranthi Varala for sharing scripts to aid the assembly, and David Goodstein and the JGI Phytozome team for database access and numerous pipeline scripts.
- Lewandowski I, Clifton-Brown JC, Scurlock JMO, Huisman W: Miscanthus: European experience with a novel energy crop. Biomass Bioenergy. 2000, 19: 209-227.View ArticleGoogle Scholar
- Hodkinson TR, Chase MW, Renvoize SA: Characterization of a genetic resource collection for Miscanthus (saccharinae, andropogoneae, poaceae) using AFLP and ISSR PCR. Ann Bot. 2002, 89: 627-636.PubMed CentralView ArticlePubMedGoogle Scholar
- Rayburn A, Crawford J, Rayburn C, Juvik J: Genome size of three Miscanthus species. Plant Mol Biol Report. 2009, 27: 184-188.View ArticleGoogle Scholar
- Swaminathan K, Alabady MS, Varala K, De Paoli E, Ho I, Rokhsar DS, Arumuganathan AK, Ming R, Green PJ, Meyers BC, Moose SP, Hudson ME, others: Genomic and small RNA sequencing of Miscanthus\times giganteus shows the utility of sorghum as a reference genome sequence for Andropogoneae grasses. Genome Biol. 2010, 11: R12-PubMed CentralView ArticlePubMedGoogle Scholar
- Nishiwaki A, Mizuguti A, Kuwabara S, Toma Y, Ishigaki G, Miyashita T, Yamada T, Matuura H, Yamaguchi S, Rayburn AL, Akashi R, Stewart JR: Discovery of natural Miscanthus (Poaceae) triploid plants in sympatric populations of Miscanthus sacchariflorus and Miscanthus sinensis in southern Japan. Am J Botany. 2011, 98: 154-159.View ArticleGoogle Scholar
- Dwiyanti MS, Rudolph A, Swaminathan K, Nishiwaki A, Shimono Y, Kuwabara S, Matuura H, Nadir M, Moose S, Stewart JR, Yamada T: Genetic analysis of putative triploid Miscanthus hybrids and tetraploid M. Sacchariflorus collected from sympatric populations of kushima, Japan. BioEnergy Res. 2012, 6: 1-8.Google Scholar
- Greef J, Deuter M, Jung C, Schondelmaier J: Genetic diversity of European Miscanthus species revealed by AFLP fingerprinting. Genet Resour Crop Evol. 1997, 44: 185-195.View ArticleGoogle Scholar
- Kim C, Zhang D, Auckland SA, Rainville LK, Jakob K, Kronmiller B, Sacks EJ, Deuter M, Paterson AH: SSR-based genetic maps of Miscanthus sinensis and M. sacchariflorus, and their comparison to sorghum. Theoretical and Applied Genetics. 2012, 124: 1325-1338.View ArticlePubMedGoogle Scholar
- Ma X-F, Jensen E, Alexandrov N, Troukhan M, Zhang L, Thomas-Jones S, Farrar K, Clifton-Brown J, Donnison I, Swaller T, Flavell R: High resolution genetic mapping by genome sequencing reveals genome duplication and tetraploid genetic structure of the diploid Miscanthus sinensis. PloS One. 2012, 7: e33821-PubMed CentralView ArticlePubMedGoogle Scholar
- Swaminathan K, Chae WB, Mitros T, Varala K, Xie L, Barling A, Glowacka K, Hall M, Jezowski S, Ming R, Hudson M, Juvik JA, Rokhsar DS, Moose SP: A framework genetic map for Miscanthus sinensis from RNA-Seq-based markers shows recent tetraploidy. BMC Genomics. 2012, 13: 142-PubMed CentralView ArticlePubMedGoogle Scholar
- Mayer KFX, Waugh R, Brown JWS, Schulman A, Langridge P, Platzer M, Fincher GB, Muehlbauer GJ, Sato K, Close TJ, Wise RP, Stein N: A physical, genetic and functional sequence assembly of the barley genome. Nature. 2012, 491: 711-716.PubMedGoogle Scholar
- Venturini L, Ferrarini A, Zenoni S, Tornielli GB, Fasoli M, Dal Santo S, Minio A, Buson G, Tononi P, Zago ED, Zamperin G, Bellin D, Pezzotti M, Delledonne M: De novo transcriptome characterization of Vitis vinifera cv. Corvina unveils varietal diversity. BMC Genomics. 2013, 14: 41-PubMed CentralView ArticlePubMedGoogle Scholar
- Schreiber AW, Sutton T, Caldo RA, Kalashyan E, Lovell B, Mayo G, Muehlbauer GJ, Druka A, Waugh R, Wise RP, Langridge P, Baumann U: Comparative transcriptomics in the Triticeae. BMC Genomics. 2009, 10: 285-PubMed CentralView ArticlePubMedGoogle Scholar
- Parchman TL, Geist KS, Grahnen JA, Benkman CW, Buerkle CA: Transcriptome sequencing in an ecologically important tree species: assembly, annotation, and marker discovery. BMC Genomics. 2010, 11: 180-PubMed CentralView ArticlePubMedGoogle Scholar
- Sekhon RS, Briskine R, Hirsch CN, Myers CL, Springer NM, Buell CR, De Leon N, Kaeppler SM: Maize gene atlas developed by RNA sequencing and comparative evaluation of transcriptomes based on RNA sequencing and microarrays. PloS One. 2013, 8: e61005-PubMed CentralView ArticlePubMedGoogle Scholar
- Chouvarine P, Cooksey AM, McCarthy FM, Ray DA, Baldwin BS, Burgess SC, Peterson DG: Transcriptome-based differentiation of closely-related Miscanthus lines. PloS One. 2012, 7: e29850-PubMed CentralView ArticlePubMedGoogle Scholar
- Breitling R, Armengaud P, Amtmann A, Herzyk P: Rank products: a simple, yet powerful, new method to detect differentially regulated genes in replicated microarray experiments. FEBS Letters. 2004, 573: 83-92.View ArticlePubMedGoogle Scholar
- Eisinga R, Breitling R, Heskes T: The exact probability distribution of the rank product statistics for replicated experiments. FEBS Letters. 2013, 587: 677-682.View ArticlePubMedGoogle Scholar
- Asakura T, Tamura T, Terauchi K, Narikawa T, Yagasaki K, Ishimaru Y, Abe K: Global gene expression profiles in developing soybean seeds. Plant Physiol Biochem. 2012, 52: 147-153.View ArticlePubMedGoogle Scholar
- Wei H, Gou J, Yordanov Y, Zhang H, Thakur R, Jones W, Burton A: Global transcriptomic profiling of aspen trees under elevated [CO2] to identify potential molecular mechanisms responsible for enhanced radial growth. J Plant Res. 2013, 126: 305-320.View ArticlePubMedGoogle Scholar
- Nemhauser JL, Hong F, Chory J: Different plant hormones regulate similar processes through largely nonoverlapping transcriptional responses. Cell. 2006, 126: 467-475.View ArticlePubMedGoogle Scholar
- Sims TL, Hague DR: Light-stimulated increase of translatable mRNA for phosphoenolpyruvate carboxylase in leaves of maize. J Biolog Chem. 1981, 256: 8252-8255.Google Scholar
- Hague DR, Uhler M, Collins PD: Cloning of cDNA for pyruvate, Pi dikinase from maize leaves. Nucleic Acids Res. 1983, 11: 4853-4865.PubMed CentralView ArticlePubMedGoogle Scholar
- Sheen J: C4 Gene expression. Ann Rev Plant Physiol Plant Mol Biol. 1999, 50: 187-217.View ArticleGoogle Scholar
- Goto K, Meyerowitz EM: Function and regulation of the Arabidopsis floral homeotic gene PISTILLATA. Genes Dev. 1994, 8: 1548-1560.View ArticlePubMedGoogle Scholar
- Jack T, Fox GL, Meyerowitz EM: Arabidopsis homeotic gene APETALA3 ectopic expression: transcriptional and posttranscriptional regulation determine floral organ identity. Cell. 1994, 76: 703-716.View ArticlePubMedGoogle Scholar
- Mandel MA, Gustafson-Brown C, Savidge B, Yanofsky MF: Molecular characterization of the Arabidopsis floral homeotic gene APETALA1. Nature. 1992, 360: 273-277.View ArticlePubMedGoogle Scholar
- Washburn JD, Murray SC, Burson BL, Klein RR, Jessup RW: Targeted mapping of quantitative trait locus regions for rhizomatousness in chromosome SBI-01 and analysis of overwintering in a Sorghum bicolor × S. propinquum population. Mol Breeding. 2013, 31: 153-162.View ArticleGoogle Scholar
- Fu J, Zhang D-F, Liu Y-H, Ying S, Shi Y-S, Song Y-C, Li Y, Wang T-Y: Isolation and characterization of maize PMP3 genes involved in salt stress tolerance. PloS One. 2012, 7: e31101-PubMed CentralView ArticlePubMedGoogle Scholar
- Capel J, Jarillo JA, Salinas J, Martinez-Zapater JM: Two homologous low-temperature-inducible genes from Arabidopsis encode highly hydrophobic proteins. Plant Physiology. 1997, 115: 569-576.PubMed CentralView ArticlePubMedGoogle Scholar
- Mitsuya S, Taniguchi M, Miyake H, Takabe T: Disruption of RCI2A leads to over-accumulation of Na + and increased salt sensitivity in Arabidopsis thaliana plants. Planta. 2005, 222: 1001-1009.View ArticlePubMedGoogle Scholar
- Mitsuya S, Taniguchi M, Miyake H, Takabe T: Overexpression of RCI2A decreases Na + uptake and mitigates salinity-induced damages in Arabidopsis thaliana plants. Physiol Plant. 2006, 128: 95-102.View ArticleGoogle Scholar
- Medina J, Catalá R, Salinas J: Developmental and stress regulation of RCI2A and RCI2B, two cold-inducible genes of arabidopsis encoding highly conserved hydrophobic proteins. Plant Physiology. 2001, 125: 1655-1666.PubMed CentralView ArticlePubMedGoogle Scholar
- Birol I, Jackman SD, Nielsen CB, Qian JQ, Varhol R, Stazyk G, Morin RD, Zhao Y, Hirst M, Schein JE, Horsman DE, Connors JM, Gascoyne RD, Marra M a, Jones SJM: De novo transcriptome assembly with ABySS. Bioinformatics. 2009, 25: 2872-2877.View ArticlePubMedGoogle Scholar
- Surget-Groba Y, Montoya-Burgos JI: Optimization of de novo transcriptome assembly from next-generation sequencing data. Genome Research. 2010, 20: 1432-1440.PubMed CentralView ArticlePubMedGoogle Scholar
- Głowacka K, Kaczmarek Z, Jeżowski S: Androgenesis in the bioenergy plant: from Calli induction to plant regeneration. Crop Sci. 2012, 52: 2659-View ArticleGoogle Scholar
- Camacho C, Coulouris G, Avagyan V, Ma N, Papadopoulos J, Bealer K, Madden TL: BLAST+: architecture and applications. BMC Bioinformatics. 2009, 10: 421-PubMed CentralView ArticlePubMedGoogle Scholar
- Quevillon E, Silventoinen V, Pillai S, Harte N, Mulder N, Apweiler R, Lopez R: InterProScan: protein domains identifier. Nucleic Acids Res. 2005, 33: W116-W120.PubMed CentralView ArticlePubMedGoogle Scholar
- Zdobnov EM, Apweiler R: InterProScan--an integration platform for the signature-recognition methods in InterPro. Bioinform (Oxford, England). 2001, 17: 847-848.View ArticleGoogle Scholar
- Ouyang S, Buell CR: The TIGR Plant Repeat Databases: a collective resource for the identification of repetitive sequences in plants. Nucleic Acids Res. 2004, 32: D360-D363.PubMed CentralView ArticlePubMedGoogle Scholar
- Koda Y, Kikuta Y: Possible involvement of jasmonic acid in tuberization of Yam plants. Plant Cell Physiol. 1991, 32: 629-633.Google Scholar
- Koda Y, Kikuta Y, Tazaki H, Tsujino Y, Sakamura S, Yoshihara T: Potato tuber-inducing activities of jasmonic acid and related compounds. Phytochemistry. 1991, 30: 1435-1438.View ArticleGoogle Scholar
- Ravnikar M, Žel J, Plaper I, Špacapan A: Jasmonic acid stimulates shoot and bulb formation of garlic in vitro. J Plant Growth Regul. 1993, 12: 73-77.View ArticleGoogle Scholar
- Rayirath UP, Lada RR, Caldwell CD, Asiedu SK, Sibley KJ: Role of ethylene and jasmonic acid on rhizome induction and growth in rhubarb (Rheum rhabarbarum L.). Plant Cell Tissue Organ Cult. 2010, 105: 253-263.View ArticleGoogle Scholar
- Medina J, Ballesteros ML, Salinas J: Phylogenetic and functional analysis of Arabidopsis RCI2 genes. J Exper Botany. 2007, 58: 4333-4346.View ArticleGoogle Scholar
- Chini A, Fonseca S, Fernández G, Adie B, Chico JM, Lorenzo O, García-Casado G, López-Vidriero I, Lozano FM, Ponce MR, Micol JL, Solano R: The JAZ family of repressors is the missing link in jasmonate signalling. Nature. 2007, 448: 666-671.View ArticlePubMedGoogle Scholar
- Fonseca S, Chico JM, Solano R: The jasmonate pathway: the ligand, the receptor and the core signalling module. Curr Opin Plant Biol. 2009, 12: 539-547.View ArticlePubMedGoogle Scholar
- Chini A, Fonseca S, Chico JM, Fernández-Calvo P, Solano R: The ZIM domain mediates homo- and heteromeric interactions between Arabidopsis JAZ proteins. Plant J. 2009, 59: 77-87.View ArticlePubMedGoogle Scholar
- Katsir L, Chung HS, Koo AJK, Howe GA: Jasmonate signaling: a conserved mechanism of hormone sensing. Curr Opin Plant Biol. 2008, 11: 428-435.PubMed CentralView ArticlePubMedGoogle Scholar
- Chang S, Puryear J, Cairney J: A simple and efficient method for isolating RNA from pine trees. Plant Mol Biol Reporter. 1993, 11: 113-116.View ArticleGoogle Scholar
- Jakobsen KS, Breivold E, Hornes E: Purification of mRNA directly from crude plant tissues in 15 minutes using magnetic oligo dT microspheres. Nucleic Acids Res. 1990, 18: 3669-PubMed CentralView ArticlePubMedGoogle Scholar
- Johnson SC: Hierarchical clustering schemes. Psychometrika. 1967, 32: 241-254.View ArticlePubMedGoogle Scholar
- Paterson AH, Bowers JE, Bruggmann R, Dubchak I, Grimwood J, Gundlach H, Haberer G, Hellsten U, Mitros T, Poliakov A, Schmutz J, Spannagl M, Tang H, Wang X, Wicker T, Bharti AK, Chapman J, Feltus FA, Gowik U, Grigoriev IV, Lyons E, Maher CA, Martis M, Narechania A, Otillar RP, Penning BW, Salamov AA, Wang Y, Zhang L, Carpita NC, et al: The Sorghum bicolor genome and the diversification of grasses. Nature. 2009, 457: 551-556.View ArticlePubMedGoogle Scholar
- Kent WJ: BLAT–the BLAST-like alignment tool. Genome Res. 2002, 12: 656-664.PubMed CentralView ArticlePubMedGoogle Scholar
- Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011, 28: 2731-2739.PubMed CentralView ArticlePubMedGoogle Scholar
- Hasegawa M, Kishino H, Yano T: Dating of the human-ape splitting by a molecular clock of mitochondrial DNA. J Mol Evol. 1985, 22: 160-174.View ArticlePubMedGoogle Scholar
- Du Z, Zhou X, Ling Y, Zhang Z, Su Z: agriGO: a GO analysis toolkit for the agricultural community. Nucleic Acids Res. 2010, 38: W64-W70.PubMed CentralView ArticlePubMedGoogle Scholar
- Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods (San Diego, Calif.). 2001, 25: 402-408.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.