Maize Gene Atlas Developed by RNA Sequencing and Comparative Evaluation of Transcriptomes Based on RNA Sequencing and Microarrays

Transcriptome analysis is a valuable tool for identification and characterization of genes and pathways underlying plant growth and development. We previously published a microarray-based maize gene atlas from the analysis of 60 unique spatially and temporally separated tissues from 11 maize organs [1]. To enhance the coverage and resolution of the maize gene atlas, we have analyzed 18 selected tissues representing five organs using RNA sequencing (RNA-Seq). For a direct comparison of the two methodologies, the same RNA samples originally used for our microarray-based atlas were evaluated using RNA-Seq. Both technologies produced similar transcriptome profiles as evident from high Pearson's correlation statistics ranging from 0.70 to 0.83, and from nearly identical clustering of the tissues. RNA-Seq provided enhanced coverage of the transcriptome, with 82.1% of the filtered maize genes detected as expressed in at least one tissue by RNA-Seq compared to only 56.5% detected by microarrays. Further, from the set of 465 maize genes that have been historically well characterized by mutant analysis, 427 show significant expression in at least one tissue by RNA-Seq compared to 390 by microarray analysis. RNA-Seq provided higher resolution for identifying tissue-specific expression as well as for distinguishing the expression profiles of closely related paralogs as compared to microarray-derived profiles. Co-expression analysis derived from the microarray and RNA-Seq data revealed that broadly similar networks result from both platforms, and that co-expression estimates are stable even when constructed from mixed data including both RNA-Seq and microarray expression data. The RNA-Seq information provides a useful complement to the microarray-based maize gene atlas and helps to further understand the dynamics of transcription during maize development.

[1]  D. Ware,et al.  Digital Gene Expression Signatures for Maize Development1[W][OA] , 2010, Plant Physiology.

[2]  A. Mohanty,et al.  Pep1, a Secreted Effector Protein of Ustilago maydis, Is Required for Successful Invasion of Plant Cells , 2009, PLoS pathogens.

[3]  Olga G. Troyanskaya,et al.  A scalable method for integration and functional analysis of multiple microarray datasets , 2006, Bioinform..

[4]  V. Walbot,et al.  Rapid transcriptome responses of maize (Zea mays) to UV-B in irradiated and shielded tissues , 2004, Genome Biology.

[5]  Bronwyn R Frame,et al.  Maize (Zea mays L.). , 2015, Methods in molecular biology.

[6]  M. Gerstein,et al.  RNA-Seq: a revolutionary tool for transcriptomics , 2009, Nature Reviews Genetics.

[7]  S. Moose,et al.  Glossy15, an APETALA2-like gene from maize that regulates leaf epidermal cell identity. , 1996, Genes & development.

[8]  J. Jeddeloh,et al.  Arabidopsis MET1 cytosine methyltransferase mutants. , 2003, Genetics.

[9]  Olga G. Troyanskaya,et al.  The Sleipnir library for computational functional genomics , 2008, Bioinform..

[10]  M. Freeling,et al.  Division and differentiation during normal and liguleless-1 maize leaf development. , 1990, Development.

[11]  Trupti Joshi,et al.  An integrated transcriptome atlas of the crop model Glycine max, and its use in comparative analyses in plants. , 2010, The Plant journal : for cell and molecular biology.

[12]  Matteo Pellegrini,et al.  Whole-Genome Analysis of Histone H3 Lysine 27 Trimethylation in Arabidopsis , 2007, PLoS biology.

[13]  J. Prioul,et al.  Maize genes encoding the small subunit of ADP-glucose pyrophosphorylase. , 2001, Plant physiology.

[14]  M. Giroux,et al.  Cloning and characterization of the brittle-2 gene of maize. , 1990 .

[15]  B. Williams,et al.  Mapping and quantifying mammalian transcriptomes by RNA-Seq , 2008, Nature Methods.

[16]  Lisa C. Harper,et al.  MaizeGDB: curation and outreach go hand-in-hand , 2011, Database J. Biol. Databases Curation.

[17]  N. Stanietsky,et al.  The interaction of TIGIT with PVR and PVRL2 inhibits human NK cell cytotoxicity , 2009, Proceedings of the National Academy of Sciences.

[18]  Xuehui Huang,et al.  Function annotation of the rice transcriptome at single-nucleotide resolution by RNA-seq. , 2010, Genome research.

[19]  Hongyu Zhao,et al.  A transcriptome atlas of rice cell types uncovers cellular, functional and developmental hierarchies , 2009, Nature Genetics.

[20]  Cole Trapnell,et al.  Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. , 2010, Nature biotechnology.

[21]  A. Molina,et al.  A protective role for the embryo surrounding region of the maize endosperm, as evidenced by the characterisation of ZmESR-6, a defensin gene specifically expressed in this region , 2005, Plant Molecular Biology.

[22]  Cole Trapnell,et al.  Ultrafast and memory-efficient alignment of short DNA sequences to the human genome , 2009, Genome Biology.

[23]  Rex T. Nelson,et al.  RNA-Seq Atlas of Glycine max: A guide to the soybean transcriptome , 2010, BMC Plant Biology.

[24]  S. Tingey,et al.  DNA array profiling of gene expression changes during maize embryo development , 2002, Functional & Integrative Genomics.

[25]  James C. Schnable,et al.  Genes Identified by Visible Mutant Phenotypes Show Increased Bias toward One of Two Subgenomes of Maize , 2011, PloS one.

[26]  B. Burr,et al.  Maize anthocyanin regulatory gene pl is a duplicate of c1 that functions in the plant. , 1993, The Plant cell.

[27]  Ute Baumann,et al.  An atlas of gene expression from seed to seed through barley development , 2006, Functional & Integrative Genomics.

[28]  Kay Denyer,et al.  Two Paralogous Genes Encoding Small Subunits of ADP-glucose Pyrophosphorylase in Maize, Bt2 and L2, Replace the Single Alternatively Spliced Gene Found in Other Cereal Species , 2007, Journal of Molecular Evolution.

[29]  M. Bucan,et al.  Promoter features related to tissue specificity as measured by Shannon entropy , 2005, Genome Biology.

[30]  J. Prioul,et al.  Expression of ADP-Glucose Pyrophosphorylase in Maize (Zea mays L.) Grain and Source Leaf during Grain Filling , 1994, Plant physiology.

[31]  Stefan R. Henz,et al.  A gene expression map of Arabidopsis thaliana development , 2005, Nature Genetics.

[32]  Jun Zheng,et al.  Genome-wide transcriptome analysis of two maize inbred lines under drought stress , 2010, Plant Molecular Biology.

[33]  R. Sekhon,et al.  Genome-wide atlas of transcription during maize development. , 2011, The Plant journal : for cell and molecular biology.

[34]  Huanming Yang,et al.  Deep RNA sequencing at single base-pair resolution reveals high complexity of the rice transcriptome. , 2010, Genome research.

[35]  V. Walbot,et al.  Gene Expression Profiling in Response to Ultraviolet Radiation in Maize Genotypes with Varying Flavonoid Content1[w] , 2003, Plant Physiology.

[36]  G. Weiller,et al.  A gene expression atlas of the model legume Medicago truncatula. , 2008, The Plant journal : for cell and molecular biology.

[37]  James C. Schnable,et al.  Differentiation of the maize subgenomes by genome dominance and both ancient and ongoing gene loss , 2011, Proceedings of the National Academy of Sciences.

[38]  Henry D. Priest,et al.  Genome-wide mapping of alternative splicing in Arabidopsis thaliana. , 2010, Genome research.

[39]  F. Vignols,et al.  The brown midrib3 (bm3) mutation in maize occurs in the gene encoding caffeic acid O-methyltransferase. , 1995, The Plant cell.

[40]  James C. Schnable,et al.  Genome-Wide Analysis of Syntenic Gene Deletion in the Grasses , 2012, Genome biology and evolution.

[41]  R. Sekhon,et al.  Utility of RNA Sequencing for Analysis of Maize Reproductive Transcriptomes , 2011 .

[42]  J. Wall,et al.  Comparison of the protein compositions of selected corns and their wild relatives, teosinte and Tripsacum , 1977 .

[43]  C. Robin Buell,et al.  Maize (Zea mays L.) Genome Diversity as Revealed by RNA-Sequencing , 2012, PloS one.

[44]  B. Larkins,et al.  The Development of Endosperm in Grasses1 , 2009, Plant Physiology.

[45]  Robert Turgeon,et al.  The developmental dynamics of the maize leaf transcriptome , 2010, Nature Genetics.

[46]  Dawn H. Nagel,et al.  The B73 Maize Genome: Complexity, Diversity, and Dynamics , 2009, Science.

[47]  M. Giroux,et al.  ADP-glucose pyrophosphorylase in shrunken-2 and brittle-2 mutants of maize , 1994, Molecular and General Genetics MGG.

[48]  Virginia Walbot,et al.  Genetic diversity contribution to errors in short oligonucleotide microarray analysis. , 2006, Plant biotechnology journal.

[49]  R. Fischer,et al.  Imprinting and Seed Development , 2004, The Plant Cell Online.

[50]  M. Stephens,et al.  RNA-seq: an assessment of technical reproducibility and comparison with gene expression arrays. , 2008, Genome research.

[51]  Lior Pachter,et al.  Sequence Analysis , 2020, Definitions.

[52]  B. Snel,et al.  A global definition of expression context is conserved between orthologs, but does not correlate with sequence conservation , 2006, BMC Genomics.

[53]  V. Solovyev,et al.  Ab initio gene finding in Drosophila genomic DNA. , 2000, Genome research.

[54]  T. Setter,et al.  Comparative Transcriptional Profiling of Placenta and Endosperm in Developing Maize Kernels in Response to Water Deficit1 , 2003, Plant Physiology.

[55]  D. Galbraith,et al.  DESIGN, PRODUCTION, AND UTILIZATION OF LONG OLIGONUCLEOTIDE MICROARRAYS FOR EXPRESSION ANALYSIS IN MAIZE , 2005 .