A Sorghum bicolorexpression atlas reveals dynamic genotype-specific expression profiles for vegetative tissues of grain, sweet and bioenergy sorghums
© Shakoor et al.; licensee BioMed Central Ltd. 2014
Received: 4 November 2013
Accepted: 15 January 2014
Published: 23 January 2014
Effective improvement in sorghum crop development necessitates a genomics-based approach to identify functional genes and QTLs. Sequenced in 2009, a comprehensive annotation of the sorghum genome and the development of functional genomics resources is key to enable the discovery and deployment of regulatory and metabolic genes and gene networks for crop improvement.
This study utilizes the first commercially available whole-transcriptome sorghum microarray (Sorgh-WTa520972F) to identify tissue and genotype-specific expression patterns for all identified Sorghum bicolor exons and UTRs. The genechip contains 1,026,373 probes covering 149,182 exons (27,577 genes) across the Sorghum bicolor nuclear, chloroplast, and mitochondrial genomes. Specific probesets were also included for putative non-coding RNAs that may play a role in gene regulation (e.g., microRNAs), and confirmed functional small RNAs in related species (maize and sugarcane) were also included in our array design. We generated expression data for 78 samples with a combination of four different tissue types (shoot, root, leaf and stem), two dissected stem tissues (pith and rind) and six diverse genotypes, which included 6 public sorghum lines (R159, Atlas, Fremont, PI152611, AR2400 and PI455230) representing grain, sweet, forage, and high biomass ideotypes.
Here we present a summary of the microarray dataset, including analysis of tissue-specific gene expression profiles and associated expression profiles of relevant metabolic pathways. With an aim to enable identification and functional characterization of genes in sorghum, this expression atlas presents a new and valuable resource to the research community.
KeywordsSorghum bicolor Gene atlas Transcriptome Gene expression Functional genomics Microarray
Sorghum [Sorghum bicolor (L.) Moench] is a staple cereal crop for millions of people in the marginal, semi-arid environments of Africa and South Asia. Its unique and advanced ability to grow in regions of low and variable rainfall highlight its potential to impact agricultural productivity in widespread water-limited environments [1, 2]. Originating and evolving across the diverse environmental landscape of Africa, morphological and physiological adaptation strategies has advanced sorghum as a naturally heat and drought-tolerant warm season C4 grass that is more efficient at utilizing water, nitrogen and energy resources with respect to other major crops, including maize and wheat [1, 3, 4]. Occupying seven million hectares of farmland, the United States is currently the world’s top sorghum producer (8.8 million annual metric tons), followed by India (7.0), Mexico (6.9), and Nigeria (4.8) (http://cgiar.org/sorghum). Cultivated in diverse climates and environmental conditions, the challenges of increasing performance and yield on marginal lands and cooler climates remains at the forefront of sorghum improvement efforts worldwide [5, 6].
Sorghum is globally established as an important source of food, feed, sugar and fiber, and recent interest in bioenergy feedstocks also spotlights sorghum as an attractive alternative for sustainable biofuel production . Framed upon the 2009 sorghum reference genome , translational genomic resources have been developed that directly impact research in other closely related C4 feedstock grasses, including switchgrass and Miscanthus[8, 9]. Comprehensive understanding of the genetic and molecular mechanisms that regulate metabolite biosynthesis, transport and storage in these species is essential for the efficient development of biofuel feedstocks.
Global transcriptome profiling further provides a means to access gene networks for the discovery of functional connections between genes, mRNAs and their regulatory proteins, and complex traits expressed through coordinated and dynamic gene networks across different tissues and developmental stages . Over the last decade, microarray-based expression profiling has provided a reliable high-throughput platform for genome-wide analysis of gene expression in many organisms. Microarrays offer substantial advantages for functional genomics, as they are increasingly cost-effective, provide a comparable accuracy of expression profiling to RNA-sequencing, and have been shown to provide comprehensive expression data (up to 90% of the transcriptome) in a given tissue . Well-established microarray data analysis tools are also available for querying, visualizing and analyzing the genomes and predicted genes [12, 13], as well as for analyzing the transcriptome profiling data and integrating with other public datasets [14–17].
To provide insight into the sorghum transcriptome, we generated a record of gene expression in a set of seven tissues and six diverse sorghum genotypes. The choice of samples reflects our aim to develop and enrich the current sorghum transcriptome literature. Previous studies have predominantly focused on reproductive tissues, and the majority of these reports do not represent the complete sorghum transcriptome. Several of these studies have also been limited to the reference genome (BTx623) or Keller, a recently resequenced sweet sorghum variety [18–22].
Comparable whole plant transcriptome maps are available for a number of other model species, including Arabidopsis thaliana, maize (Zea mays), barley (Hordeum vulgare) , rice (Oryza sativa) [26, 27], and soybean (Glycine max) . These recent transcriptome surveys were constructed with only one genotype or line/accession for their respective species of interest, whereas the present study aims to highlight the practical importance of examining expression profiles across diverse tissue types, developmental stages, as well as genotypes in order to accurately target genes and metabolic pathways for the efficient development of improved feedstocks.
Fundamental understanding of sorghum genomics is necessary for improving sorghum for agronomic and compositional traits. Specifically, genotypes with high biomass and increased levels of fermentable stem sugars are ideal for developing feedstocks for the biofuel industry. We developed this genomic resource, the whole-transcriptome array as well as the vegetative transcriptome in diverse genotypes and tissues, in order to facilitate the characterization of molecular networks and regulatory mechanisms governing important metabolic pathways including, but not limited to, cell wall biosynthesis for lignocellulosic biomass as well as synthesis, translocation, and storage of fermentable photosynthates for energy content. The relevance of our dataset is demonstrated by genotype and tissue-specific expression of the phenylpropanoid and lignin biosynthetic pathway genes.
Intended as readily available public resource for functional gene characterization, the transcriptome data presented here is available through NCBI's Gene Expression Omnibus (GEO) under accession number GSE49879, and the Sorghum Genome Array is available through Affymetrix (http://affymetrix.com).
Results and discussion
Generation and quality assessment of data
A whole-transcriptome exon array for Sorghum bicolor was custom-designed by Chromatin, Inc. (http://chromatininc.com) and Affymetrix: Sorgh-WTa520972F. This genechip contains 1,026,373 probes covering 149,182 exons (27,577 genes) across the Sorghum bicolor genome (10 chromosomes), chloroplast and mitochondria. The sequences used to construct the probesets included all identified Sorghum bicolor exons from the Sbi1 assembly (http://www.phytozome.net). Multiple probes were chosen for each exon, with a minimum of one probe per exon and 25 probes per gene. In addition to standard Affymetrix controls, positive controls in the microarray design included probes for constitutively expressed Sorghum bicolor genes (actin, ubiquitin and eIF4a1). Probes for intronic regions of actin and ubiquitin were also utilized to determine background expression levels.
To study the sorghum transcriptome and build a gene expression atlas, we collected 78 diverse samples from various developmental stages and tissue types (Additional file 1). In order to broadly capture sorghum genetic diversity, we included genotypes representing three major ideotypes, including grain, sweet, and bioenergy sorghums. Our study includes R159, an elite grain sorghum characterized by the valuable agronomic traits of uniform growth and disease resistance . Grain sorghum is cultivated primarily for its high starch content, applications in human/animal health and nutrition, and as biofuel feedstock for ethanol production . We also included two sweet sorghums, Fremont and Atlas, that produce increased biomass and accumulate high levels of fermentable carbohydrates in the stem. Additionally, Fremont is drought resistant and flowers early, while Atlas is less susceptible to lodging (due to a stiff stalk phenotype) and flowers later . We also selected three bioenergy or high biomass lines, PI455230, PI152611, and AR2400 that produce increased levels of cellulosic material and are photoperiod sensitive, which allows the plant to produce higher amounts of vegetative matter under long day conditions (Additional file 2). PI152611 is specifically a forage line, a fast-growing, highly digestible grass utilized for livestock feed [5, 29].
The primary goal of this study was to obtain relevant and applicable data for the research community developing sorghum as a global feedstock; this research interest guided our sample selection towards vegetative tissues, with a strong bias for stem tissues. A comprehensive trancriptomic profile of sorghum inflorescence and leaf data was recently made available to the community . We compared the leaf RNA sequencing dataset with the present leaf dataset to demonstrate and confirm that our microarray analysis approach towards transcriptome profiling was appropriate. The Spearman correlation of the transcriptome across technologies is 0.61 (Additional file 3), which is consistent with several studies comparing RNA-seq and microarray methods for genome-wide transcriptome profiling [31–33]. The present comparison corroborates these studies and demonstrates that the microarray platform for expression profiling correlates well with current sequencing methods. With a common goal of crop improvement, complementary datasets such as these generate a core of information that can be explored for the functional characterization of genes and genetic pathways.
We assessed data quality for hybridization by comparing normalized signals of all probe sets between biological replicates using Pearson’s correlation analysis. The biological replicates were highly correlated, with an average Pearson’s correlation coefficient of 0.99 (Additional file 4). The highly reproducible results of the replicate data further validate the quality of the microarray platform and present dataset. Previous studies have consistently established strong correlations between qRT-PCR data and microarray data processed using robust multi-array analysis (RMA) [34, 35]. However, we also tested a small subset of these genes via qRT-PCR to validate the array-generated expression data and expression patterns across multiple tissue types (Additional file 5A and 5B).
To further assure data quality, we also examined the general expression patterns of well-characterized genes that have been highlighted for tissue-specific expression in previous studies. In microarray experiments with RNA isolated from shoot tips, we observed high expression levels for homologs of SPATULA, a shoot tip transcription factor that is strongly expressed in shoot tips and young leaf primordia . Similarly, the sorghum homolog for TIP2-3, a root-specific aquaporin gene , was also expressed at higher levels in our study using root-isolated RNA (Additional file 6).
Global gene expression patterns
We detected the expression of 19,354 genes in at least one of the 78 samples, representing 70.2% of all genes on the array (27,577 genes). The number of expressed transcripts detected in the various tissues ranged between 10,850 and 11,587 (representing 56 to 60% of all expressed genes on the array). Expressed genes were determined following established methods , and with a conservative and arbitrary expression threshold cutoff of 320 (five times the mean normalized signal from intronic gene probes used as controls), we found that 15.4% of genes on the array were detected in all tissues (4256/27,577) (Additional file 7). Gene ontology (GO) annotation analysis of these constitutively expressed genes reveals that most are involved in basic biological processes including development, protein synthesis/modification, and signal transduction (Additional file 8). Similar to published work in maize, expression of constitutive genes varied among the samples, with the coefficient of variation (CV) ranging from 5% to 129%. With a CV of 10.4%, we identified a ubiquitin-conjugating enzyme, Sb09g023560, as one of the most stably expressed genes (Additional file 9). This class of genes was also identified in the maize atlas as the most stably expressed among variable tissues .
A diverse range of plant tissues was sampled in this study; however, 29.8% of the probesets were not detected above our designated expression threshold level. Several plausible explanations can account for this incomplete expression coverage, including gene expression from specific tissues and/or developmental stages not included in this study, false positive gene models, and levels of expression below detection threshold limits. Further, the arrays were developed utilizing the BTx623 reference sequence and do not capture polymorphisms, copy number variation and presence-absence variation across all the sampled genotypes.
Transcriptome-based classification of sorghum tissues
Interestingly, three out of 46 stem samples clustered near the group of meristematic tissues (roots and shoot tips). All three of these 'outlier’ samples were collected at the top internode, 61 days after planting (DAP) in three of the six sampled genotypes (PI455230, PI152611, and AR2400). At 70DAP, the stem samples from same genotypes clustered with the other stem samples. These lines are characterized as high biomass genotypes, whereas the remaining three genotypes can be characterized as either grain or sweet lines (R159, Atlas, and Fremont). The PCA indicates that at 61DAP, the patterns of gene expression in the stem of the high biomass lines are more related to meristematic regions, or regions of active growth. While it is possible that these three stem samples were collected too close to the meristematic shoot tip region, further study may indicate that the differential transcriptome in the stems of these lines capture a transition zone of gene expression in which sorghum commits to post-reproductive pathways of sugar production and grain fill versus continued biomass production. This result further demonstrates the importance of examining genotype, tissue type, as well as temporal expression patterns when targeting transcriptional programs of interest.
Tissue and genotype-specific patterns of gene expression
We next attempted to determine whether functional gene classes were over-represented in specific genotypes. GO analysis did not reveal statistical differences in the enrichment of GO slim terms using agriGO (Fisher’s exact test and the Yekutieli (false-discovery rate under dependency) multi-test adjustment method) . However, this can partially be attributed to the incomplete annotation of the sorghum genome, as well as stage and tissue-specific expression not captured in our sample collection.
Tissue-specific expression of genes involved in the phenylpropanoid-monolignol pathway
To exemplify the functional utility of this data, we highlighted the expression data of 10 key enzymes associated with the phenylpropanoid-monolignol biosynthesis pathway (Additional file 13). Currently, one of the primary strategies for bioenergy feedstock improvement is through lignin modification. Alterations in lignin content and composition aim to improve the digestibility of forage and saccharification efficiency of lignocellulosic biofuels [41, 42]. Thus, modifying the expression of genes in the lignin biosynthesis pathway is an attractive approach to achieving this goal.
Comprehensive transcriptome profiling provides a global overview of gene networks and allows for the discovery of functional connections between genes, mRNAs and their regulatory proteins. In the present study, we constructed a gene expression atlas covering an array of tissues, developmental stages and genotypes using the first commercially available sorghum microarray (Sorgh-WTa520972F). We observed tissue and genotype-specific expression patterns of relevant metabolic pathways that highlight the significance of intra-species variation in sorghum.
Developed as a new resource for crop breeding and genomic discovery, Sorgh-WTa520972F is produced by Affymetrix and is available to the public research community. We are currently utilizing this microarray to identify differential gene expression related to key metabolic processes (e.g., starch/lignin biosynthesis) for the identification of regulatory regions. Additional avenues for future study with this array are wide-ranging and can include gene expression profiling during abiotic/biotic stress, plant infection and disease establishment to investigate genetic mechanisms and applications to plant breeding and crop improvement. Detailed expression analysis of small RNAs included in the array design may also reveal key insights in diverse biological processes, including RNA-guided gene regulation. Sorgh-WTa520972F can also be utilized in quantitative trait locus (QTL) mapping and validation methods (e.g., identify differentially expressed genes from 'tolerant’ versus 'sensitive’ varieties). Minimal costs associated with microarray analysis allow for the generation of high-throughput expression profiles or combinations of profiles of elite breeding lines for accelerated crop-breeding efforts. Applications of this resource can target numerous agronomic traits in sorghum as well as provide insight in closely related grasses (e.g., sugarcane, switchgrass, Miscanthus x giganteus) for improved feedstock development.
To study the sorghum transcriptome and build the present gene expression atlas, we collected 78 samples from various developmental stages and tissue types (Additional file 1). Six diverse sorghum genotypes were grown in Chromatin’s greenhouse and field sites (Champaign, IL). These six genotypes were chosen to represent ideotypes of sorghum cultivation, including sweet, grain and high biomass sorghum varieties. Greenhouse grown seedling shoot and root samples were collected at 10DAP, which is roughly five days after plant emergence. Whole leaf and meristematic shoot tip samples were collected at 38DAP. This time-point captures the active growth phase of vegetative structures, including leaves, shoots and tillers. The stem tissue samples were collected at two time points: 61 and 70DAP. At 61DAP, the stem is fully formed in both flowering and non-flowering types. In flowering types, the head is also fully formed, and the period between 61 and 70DAP is a stage of active metabolism, capturing the transition between flowering (61DAP) and active grain filling (70DAP) . The stem tissue was further dissected into the pith and the rind. As a bioenergy crop, the majority of fermentable sugar available in sorghum is present in the pith. The majority of lignin, however, is found in the rind . Two tissue types (shoot and root) were represented by two biological replicates.
A whole-transcriptome exon array for Sorghum bicolor: Sorgh-WTa520972F was designed and utilized for the present expression study. The array contains 1,026,373 probes covering 149,182 exons (27,577 genes) across the Sorghum bicolor nuclear, chloroplast and mitochondrial genome. The sequences used to construct the probesets included all identified Sorghum bicolor exons from the Sbi1 assembly and Sbi1.4 annotation (http://phytozome.net). We also added sequences for putative non-coding RNAs in Sorghum bicolor that may play a role in gene regulation (e.g., rRNAs, tRNAs, snoRNAs and microRNAs). Confirmed functional small RNAs in closely related species (maize, sugarcane) were also included in our array design (http://bioinformatics.cau.edu.cn/PMRD, http://www.ncrna.org/frnadb) (Additional file 10).
RNA Isolation and hybridization
Total RNA from all tissue types was extracted using a NucleoSpin RNA Plant Kit (Maxherey-Nagel, Germany). RNA integrity, as indicated by the detection of discrete ribosomal subunits, was verified electrophoretically. The RNA quality and quantity was further validated with a NanoDrop spectrophotometer (NanoDrop Technologies, Wilmington, DE). Prior to hybridization, the total RNA profile was also analyzed with Agilent 2100 Bioanalyzer (Agilent technologies, Waldbronn, Germany). Synthesis of cDNA, probe labeling and hybridization was performed by Precision Biomarker (Precision Biomarker Resources, Inc. Evanston, Illinois)
Data extraction and evaluation of gene expression
Background correction and normalization were performed using a robust multi-chip average (RMA) algorithm in the Bioconductor Affy package . Present calls for expressed genes were determined following established methods . In brief, an expressed gene was identified by a RMA-normalized linear expression of >/= 320 in at least one of the 78 samples. The expression cut-off was five times the mean RMA-normalized signal from 576 negative-control oligos selected from the intronic regions of known constitutive genes (e.g., actin, ubiquitin, and eIF4a1). A mean signal intensity of 64 was determined for the negative control oligos analyzed across all 78 slides. Constitutively expressed genes were identified by a RMA-normalized linear expression value of >/= 320 in all 78 samples.
Principal component analysis, hierarchical clustering and z-scores
To study the biological relatedness and identify expression trends among the samples, we utilized the cmdscale function and then plotted using R. We used RMA-normalized log2 normalized expression values in the PCA analysis. Hierarchical clustering was performed using RMA-normalized log2 normalized expression values and clustered using Pearson’s correlation analysis. The Z scores were calculated as follows: Z = (X-Xmean)/SD, where X is the average expression of a given gene in a tissue, and Xmean and SD are the mean expression and standard deviation respectively of that gene across all the selected tissues.
GO Slim enrichment analysis
We evaluated enrichment of GO slim terms of biological process category (http://geneontology.org/GO.slims) in agriGO (http://bioinfo.cau.edu.cn/agriGO/) by Fisher’s exact test (p-value ≤0.05) and the Yekutieli (false-discovery rate under dependency) multi-test adjustment method .
The relative mRNA expression was measured using Peltier Thermal Cycler PTC-200 PCR machine (MJ Research, Waltham, MA, USA) and the SuperScript III Platinum SYBR Green One-Step qRT- PCR kit (Invitrogen, Carlsbad, CA). Three independent reverse transcription reactions were performed for each RNA sample, and qRT-PCR was carried out under the following conditions: 100 nanograms of each RNA sample was reverse transcribed at 60°C for 3 minutes, and reverse transcription was followed by initial activation at 95°C for 5 minutes, and 40 amplification cycles at 95°C for 15 s and 50°C for 30s. Results were analysed using MJ Opticon Monitor 3.1.32 software, and relative expression of mRNA was calculated by the comparative Ct method (2-[Δ][Δ]Ct) . Gene expression values across tissue types were normalized to ubiquitin expression.
Availability of supporting data
The transcriptome dataset supporting the results of this article is available through NCBI's Gene Expression Omnibus (GEO) under accession number GSE49879, and the Sorghum Genome Array is available through Affymetrix (http://affymetrix.com).
We would like to thank Otto Folkerts and Jennifer Mach for critically reviewing and editing this manuscript. We would also like to thank Shu-wen Huang for her contribution to data procurement. For their contributions to sample collection and RNA preparation, we also thank Zhanmin Fan, Piero Barone, Soumya Bhattacharya, Michael Root, David Sobczak, Kaushal Patel, Duyen Tran, and Shivani Verma. This project is partially funded through ARPA-E Cooperative Agreement DE-AR0470-1541.
- Rosenow D, Quisenberry J, Wendt C, Clark L: Drought tolerant sorghum and cotton germplasm. Agric Water Manage. 1983, 7 (1): 207-222.View ArticleGoogle Scholar
- Mann J, Kimber C, Miller F: The origin and early cultivation of sorghums in Africa. Texas Agric Exp Station. 1983, 1: 128074-Available electronically from http://hdl.handle.net/1969.1/1280Google Scholar
- de Vries SC, van de Ven GW, van Ittersum MK, Giller KE: Resource use efficiency and environmental performance of nine major biofuel crops, processed by first-generation conversion techniques. Biomass and Bioenergy. 2010, 34 (5): 588-601. 10.1016/j.biombioe.2010.01.001.View ArticleGoogle Scholar
- Kimber CT, Dahlberg JA, Kresovich S: The gene pool of Sorghum bicolorand its improvement. In Genomics of the Saccharinae. New York: Springer;2013:23–41.View ArticleGoogle Scholar
- Shoemaker C, Bransby DI: The role of sorghum as a bioenergy feedstock.In Sustainable alternative fuel feedstock opportunities, challenges androadmaps for six US regions. Ankeny, IA: Soil and Water ConservationSociety; 2010:149–159.Google Scholar
- Burow G, Burke JJ, Xin Z, Franks CD: Genetic dissection of early-season cold tolerance in sorghum (Sorghum bicolor (L) Moench). Mol Breeding. 2011, 28 (3): 391-402. 10.1007/s11032-010-9491-4.View ArticleGoogle Scholar
- Paterson AH, Bowers JE, Bruggmann R, Dubchak I, Grimwood J, Gundlach H, Haberer G, Hellsten U, Mitros T, Poliakov A: The Sorghum bicolor genome and the diversification of grasses. Nature. 2009, 457 (7229): 551-556. 10.1038/nature07723.View ArticlePubMedGoogle Scholar
- Zhang J-Y, Lee Y-C, Torres-Jerez I, Wang M, Yin Y, Chou W-C, He J, Shen H, Srivastava AC, Pennacchio C: Development of an integrated transcript sequence database and a gene expression atlas for gene discovery and analysis in switchgrass (Panicum virgatum L). Plant J. 2013, 74 (1): 160-173. 10.1111/tpj.12104.View ArticlePubMedGoogle Scholar
- Swaminathan K, Alabady MS, Varala K, De Paoli E, Ho I, Rokhsar DS, Arumuganathan AK, Ming R, Green PJ, Meyers BC: Genomic and small RNA sequencing of Miscanthus x giganteus shows the utility of sorghum as a reference genome sequence for Andropogoneae grasses. Genome Biol. 2010, 11 (2): R12-10.1186/gb-2010-11-2-r12.PubMed CentralView ArticlePubMedGoogle Scholar
- Komili S, Silver PA: Coupling and coordination in gene expression processes: a systems biology view. Nat Rev Genet. 2008, 9 (1): 38-48. 10.1038/nrg2223.View ArticlePubMedGoogle Scholar
- Mane S, Evans C, Cooper K, Crasta O, Folkerts O, Hutchison S, Harkins T, Thierry-Mieg D, Thierry-Mieg J, Jensen R: Transcriptome sequencing of the Microarray Quality Control (MAQC) RNA reference samples using next generation sequencing. BMC Genomics. 2009, 10 (1): 264-10.1186/1471-2164-10-264.PubMed CentralView ArticlePubMedGoogle Scholar
- Goodstein DM, Shu S, Howson R, Neupane R, Hayes RD, Fazo J, Mitros T, Dirks W, Hellsten U, Putnam N: Phytozome: a comparative platform for green plant genomics. Nucleic Acids Res. 2012, 40 (D1): D1178-D1186. 10.1093/nar/gkr944.PubMed CentralView ArticlePubMedGoogle Scholar
- Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J: Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004, 5 (10): R80-10.1186/gb-2004-5-10-r80.PubMed CentralView ArticlePubMedGoogle Scholar
- Barrett T, Troup DB, Wilhite SE, Ledoux P, Rudnev D, Evangelista C, Kim IF, Soboleva A, Tomashevsky M, Edgar R: NCBI GEO: mining tens of millions of expression profiles–database and tools update. Nucleic Acids Res. 2007, 35 (suppl 1): D760-D765.PubMed CentralView ArticlePubMedGoogle Scholar
- Zimmermann P, Laule O, Schmitz J, Hruz T, Bleuler S, Gruissem W: Genevestigator transcriptome meta-analysis and biomarker search using rice and barley gene expression databases. Mol Plant. 2008, 1 (5): 851-857. 10.1093/mp/ssn048.View ArticlePubMedGoogle Scholar
- Joung J-G, Corbett AM, Fellman SM, Tieman DM, Klee HJ, Giovannoni JJ, Fei Z: Plant MetGenMAP: an integrative analysis system for plant systems biology. Plant Physiol. 2009, 151 (4): 1758-1768. 10.1104/pp.109.145169.PubMed CentralView ArticlePubMedGoogle Scholar
- Katari MS, Nowicki SD, Aceituno FF, Nero D, Kelfer J, Thompson LP, Cabello JM, Davidson RS, Goldberg AP, Shasha DE: VirtualPlant: a software platform to support systems biology research. Plant Physiol. 2010, 152 (2): 500-515. 10.1104/pp.109.147025.PubMed CentralView ArticlePubMedGoogle Scholar
- Buchanan CD, Lim S, Salzman RA, Kagiampakis I, Morishige DT, Weers BD, Klein RR, Pratt LH, Cordonnier-Pratt M-M, Klein PE: Sorghum bicolor's transcriptome response to dehydration, high salinity and ABA. Plant Mol Biol. 2005, 58 (5): 699-720. 10.1007/s11103-005-7876-2.View ArticlePubMedGoogle Scholar
- Davidson RM, Gowda M, Moghe G, Lin H, Vaillancourt B, Shiu SÄ, Jiang N, Robin Buell C: Comparative transcriptomics of three Poaceae species reveals patterns of gene expression evolution. Plant J. 2012, 71 (3): 492-502.PubMedGoogle Scholar
- Salzman RA, Brady JA, Finlayson SA, Buchanan CD, Summer EJ, Sun F, Klein PE, Klein RR, Pratt LH, Cordonnier-Pratt M-M: Transcriptional profiling of sorghum induced by methyl jasmonate, salicylic acid, and aminocyclopropane carboxylic acid reveals cooperative regulation and novel gene responses. Plant Physiol. 2005, 138 (1): 352-368. 10.1104/pp.104.058206.PubMed CentralView ArticlePubMedGoogle Scholar
- Jiang S-Y, Ma Z, Vanitha J, Ramachandran S: Genetic variation and expression diversity between grain and sweet sorghum lines. BMC Genomics. 2013, 14 (1): 18-10.1186/1471-2164-14-18.PubMed CentralView ArticlePubMedGoogle Scholar
- Dugas DV, Monaco MK, Olson A, Klein RR, Kumari S, Ware D, Klein PE: Functional annotation of the transcriptome of Sorghum bicolor in response to osmotic stress and abscisic acid. BMC Genomics. 2011, 12 (1): 514-10.1186/1471-2164-12-514.PubMed CentralView ArticlePubMedGoogle Scholar
- Schmid M, Davison TS, Henz SR, Pape UJ, Demar M, Vingron M, Schölkopf B, Weigel D, Lohmann JU: A gene expression map of Arabidopsis thaliana development. Nat Genet. 2005, 37 (5): 501-506. 10.1038/ng1543.View ArticlePubMedGoogle Scholar
- Sekhon RS, Lin H, Childs KL, Hansey CN, Buell CR, de Leon N, Kaeppler SM: Genome-wide atlas of transcription during maize development. Plant J. 2011, 66 (4): 553-563. 10.1111/j.1365-313X.2011.04527.x.View ArticlePubMedGoogle Scholar
- Druka A, Muehlbauer G, Druka I, Caldo R, Baumann U, Rostoks N, Schreiber A, Wise R, Close T, Kleinhofs A: An atlas of gene expression from seed to seed through barley development. Funct Integr Genomics. 2006, 6 (3): 202-211. 10.1007/s10142-006-0025-4.View ArticlePubMedGoogle Scholar
- Jiao Y, Tausta SL, Gandotra N, Sun N, Liu T, Clay NK, Ceserani T, Chen M, Ma L, Holford M: A transcriptome atlas of rice cell types uncovers cellular, functional and developmental hierarchies. Nat Genet. 2009, 41 (2): 258-263. 10.1038/ng.282.View ArticlePubMedGoogle Scholar
- Wang L, Xie W, Chen Y, Tang W, Yang J, Ye R, Liu L, Lin Y, Xu C, Xiao J: A dynamic gene expression atlas covering the entire life cycle of rice. Plant J. 2010, 61 (5): 752-766. 10.1111/j.1365-313X.2009.04100.x.View ArticlePubMedGoogle Scholar
- Libault M, Farmer A, Joshi T, Takahashi K, Langley RJ, Franklin LD, He J, Xu D, May G, Stacey G: An integrated transcriptome atlas of the crop model Glycine max, and its use in comparative analyses in plants. Plant J. 2010, 63 (1): 86-99.PubMedGoogle Scholar
- USDA, ARS, National Genetic Resources Program: Germplasm ResourcesInformation Network - (GRIN). Beltsville, Maryland: National GermplasmResources Laboratory; 2014.Google Scholar
- Martin JH, Stephens JC: The culture and use of sorghums for forage.Washington D.C: US Department of Agriculture; 1955.Google Scholar
- Agarwal A, Koppstein D, Rozowsky J, Sboner A, Habegger L, Hillier LW, Sasidharan R, Reinke V, Waterston RH, Gerstein M: Comparison and calibration of transcriptome data from RNA-Seq and tiling arrays. BMC Genomics. 2010, 11 (1): 383-10.1186/1471-2164-11-383.PubMed CentralView ArticlePubMedGoogle Scholar
- Kogenaru S, Yan Q, Guo Y, Wang N: RNA-seq and microarray complement each other in transcriptome profiling. BMC Genomics. 2012, 13 (1): 629-10.1186/1471-2164-13-629.PubMed CentralView ArticlePubMedGoogle Scholar
- Marioni JC, Mason CE, Mane SM, Stephens M, Gilad Y: RNA-seq: an assessment of technical reproducibility and comparison with gene expression arrays. Genome Res. 2008, 18 (9): 1509-1517. 10.1101/gr.079558.108.PubMed CentralView ArticlePubMedGoogle Scholar
- Dallas PB, Gottardo NG, Firth MJ, Beesley AH, Hoffmann K, Terry PA, Freitas JR, Boag JM, Cummings AJ, Kees UR: Gene expression levels assessed by oligonucleotide microarray analysis and quantitative real-time RT-PCR - how well do they correlate?. BMC Genomics. 2005, 6 (1): 59-10.1186/1471-2164-6-59.PubMed CentralView ArticlePubMedGoogle Scholar
- Morey JS, Ryan JC, Van Dolah FM: Microarray validation: factors influencing correlation between oligonucleotide microarrays and real-time PCR. Biol Proced Online. 2006, 8 (1): 175-193. 10.1251/bpo126.PubMed CentralView ArticlePubMedGoogle Scholar
- Ichihashi Y, Horiguchi G, Gleissberg S, Tsukaya H: The bHLH transcription factor SPATULA controls final leaf size in Arabidopsis thaliana. Plant Cell Physiol. 2010, 51 (2): 252-261. 10.1093/pcp/pcp184.View ArticlePubMedGoogle Scholar
- Lopez F, Bousser A, Sissoéff I, Hoarau J, Mahé A: Characterization in maize of ZmTIP2-3, a root-specific tonoplast intrinsic protein exhibiting aquaporin activity. J Exp Bot. 2004, 55 (396): 539-541. 10.1093/jxb/052.View ArticlePubMedGoogle Scholar
- Meyer E, Logan TL, Juenger TE: Transcriptome analysis and gene expression atlas for Panicum hallii var. filipes, a diploid model for biofuel research. Plant J. 2012, 70 (5): 879-890. 10.1111/j.1365-313X.2012.04938.x.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 (suppl 2): W64-W70.PubMed CentralView ArticlePubMedGoogle Scholar
- Qazi HA, Paranjpe S, Bhargava S: Stem sugar accumulation in sweet sorghum - activity and expression of sucrose metabolizing enzymes and sucrose transporters. J Plant Physiol. 2012, 169 (6): 605-613. 10.1016/j.jplph.2012.01.005.View ArticlePubMedGoogle Scholar
- Zhao Q, Tobimatsu Y, Zhou R, Pattathil S, Gallego-Giraldo L, Fu C, Jackson LA, Hahn MG, Kim H, Chen F: Loss of function of cinnamyl alcohol dehydrogenase 1 leads to unconventional lignin and a temperature-sensitive growth defect in Medicago truncatula. Proc Natl Acad Sci USA. 2013, 110 (33): 13660-13665. 10.1073/pnas.1312234110.PubMed CentralView ArticlePubMedGoogle Scholar
- Saballos A, Sattler SE, Sanchez E, Foster TP, Xin Z, Kang C, Pedersen JF, Vermerris W: Brown midrib2 (Bmr2) encodes the major 4 coumarate: coenzyme A ligase involved in lignin biosynthesis in sorghum (Sorghum bicolor (L.) Moench). Plant J. 2012, 70 (5): 818-830. 10.1111/j.1365-313X.2012.04933.x.View ArticlePubMedGoogle Scholar
- Jouanin L, Lapierre C: Lignins: biosynthesis. Biodegradation andBioengineering: Elsevier Science; 2012.Google Scholar
- Sattler SE, Funnell-Harris DL: Modifying lignin to improve bioenergy feedstocks: strengthening the barrier against pathogens?. Front Plant Sci. 2013, 4: 70.PubMed CentralView ArticlePubMedGoogle Scholar
- Vanderlip R, Reeves H: Growth stages of sorghum [Sorghum bicolor (L.) Moench]. Agronomy J. 1972, 64 (1): 13-16. 10.2134/agronj1972.00021962006400010005x.View ArticleGoogle Scholar
- Whitfield MB, Chinn MS, Veal MW: Processing of materials derived from sweet sorghum for biobased products. Ind Crops Products. 2012, 37 (1): 362-375. 10.1016/j.indcrop.2011.12.011.View ArticleGoogle 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. 2001, 25 (4): 402-408. 10.1006/meth.2001.1262.View ArticlePubMedGoogle 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. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.