Natural antisense transcripts are significantly involved in regulation of drought stress in maize

Abstract Natural antisense transcripts (NATs) are a prominent and complex class of regulatory RNAs. Using strand-specific RNA sequencing, we identified 1769 sense and antisense transcript pairs (NAT pairs) in two maize inbreds with different sensitivity to drought, as well as in two derivative recombination inbred lines (RILs). A significantly higher proportion of NATs relative to non-NATs are specifically expressed under water stress (WS). Surprisingly, expression of sense and antisense transcripts produced by NAT pairs is significantly correlated, particularly under WS. We found an unexpected large proportion of NATs with protein coding potential, as estimated by ribosome release scores. Small RNAs significantly accumulate within NAT pairs, with 21 nt smRNA particularly enriched in overlapping regions of these pairs of genes. The abundance of these smRNAs is significantly altered in the leafbladeless1 mutant, suggesting that these genes may be regulated by the tasiRNA pathway. Further, NATs are significantly hypomethylated and include fewer transposable element sequences relative to non-NAT genes. NAT gene regions also exhibit higher levels of H3K36me3, H3K9ac, and H3K4me3, but lower levels of H3K27me3, indicating that NAT gene pairs generally exhibit an open chromatin configuration. Finally, NAT pairs in 368 diverse maize inbreds and 19 segregating populations were specifically enriched for polymorphisms associated with drought tolerance. Taken together, the data highlight the potential impact of that small RNAs and histone modifications have in regulation of NAT expression, and the significance of NATs in response to WS.

[1]  M. Axtell,et al.  Genome-Wide Characterization of Maize Small RNA Loci and Their Regulation in the required to maintain repression6-1 (rmr6-1) Mutant and Long-Term Abiotic Stresses1[OPEN] , 2016, Plant Physiology.

[2]  I. Makałowska,et al.  Biological functions of natural antisense transcripts. , 2016, Acta biochimica Polonica.

[3]  P. Kumari,et al.  cncRNAs: Bi-functional RNAs with protein coding and non-coding functions , 2015, Seminars in cell & developmental biology.

[4]  Xiaojie Li,et al.  Ribosome profiling reveals dynamic translational landscape in maize seedlings under drought stress. , 2015, The Plant journal : for cell and molecular biology.

[5]  Robert W. Williams,et al.  Overexpression of ARGOS Genes Modifies Plant Sensitivity to Ethylene, Leading to Improved Drought Tolerance in Both Arabidopsis and Maize[OPEN] , 2015, Plant Physiology.

[6]  S. Dhanasekaran,et al.  The landscape of antisense gene expression in human cancers , 2015, Genome research.

[7]  J. Mellor,et al.  Sense and antisense transcription are associated with distinct chromatin architectures across genes , 2015, Nucleic acids research.

[8]  Andrew P. Harrison,et al.  Genome-wide view of natural antisense transcripts in Arabidopsis thaliana , 2015, DNA research : an international journal for rapid publication of reports on genes and genomes.

[9]  S. Salzberg,et al.  StringTie enables improved reconstruction of a transcriptome from RNA-seq reads , 2015, Nature Biotechnology.

[10]  J Crossa,et al.  Genomic prediction in biparental tropical maize populations in water-stressed and well-watered environments using low-density and GBS SNPs , 2014, Heredity.

[11]  M. Hammell,et al.  Genome-Wide Analysis of leafbladeless1-Regulated and Phased Small RNAs Underscores the Importance of the TAS3 ta-siRNA Pathway to Maize Development , 2014, PLoS genetics.

[12]  Alexandre Lomsadze,et al.  Identification of protein coding regions in RNA transcripts , 2014, BCB.

[13]  Lu Jiang,et al.  Genome-wide identification of housekeeping genes in maize , 2014, Plant Molecular Biology.

[14]  F. Ali,et al.  Genome Wide Association Studies Using a New Nonparametric Model Reveal the Genetic Architecture of 17 Agronomic Traits in an Enlarged Maize Association Panel , 2014, PLoS genetics.

[15]  Lukas Burger,et al.  Short sequences can efficiently recruit histone H3 lysine 27 trimethylation in the absence of enhancer activity and DNA methylation , 2014, Proceedings of the National Academy of Sciences.

[16]  Qingli Guo,et al.  Identification of Maize Long Non-Coding RNAs Responsive to Drought Stress , 2014, PloS one.

[17]  M. Albà,et al.  Long non-coding RNAs as a source of new peptides , 2014, eLife.

[18]  M. Axtell,et al.  Identification and annotation of small RNA genes using ShortStack. , 2014, Methods.

[19]  Huan Wang,et al.  Genome-wide identification of long noncoding natural antisense transcripts and their responses to light in Arabidopsis , 2014, Genome research.

[20]  M. A. Pedraza,et al.  Insights into the Maize Pan-Genome and Pan-Transcriptome[W][OPEN] , 2014, Plant Cell.

[21]  Wei Shi,et al.  featureCounts: an efficient general purpose program for assigning sequence reads to genomic features , 2013, Bioinform..

[22]  Xiaosen Guo,et al.  Identification of candidate genes for drought tolerance by whole-genome resequencing in maize , 2014, BMC Plant Biology.

[23]  Steven R. Eichten,et al.  Genome-wide discovery and characterization of maize long non-coding RNAs , 2014, Genome Biology.

[24]  Yuke He,et al.  Global analysis of cis-natural antisense transcripts and their heat-responsive nat-siRNAs in Brassica rapa , 2013, BMC Plant Biology.

[25]  L. Steinmetz,et al.  Gene regulation by antisense transcription , 2013, Nature Reviews Genetics.

[26]  Y. Poirier,et al.  A Rice cis-Natural Antisense RNA Acts as a Translational Enhancer for Its Cognate mRNA and Contributes to Phosphate Homeostasis and Plant Fitness[C][W] , 2013, Plant Cell.

[27]  K. Shinozaki,et al.  Genome-Wide Analysis of ZmDREB Genes and Their Association with Natural Variation in Drought Tolerance at Seedling Stage of Zea mays L , 2013, PLoS genetics.

[28]  Colin N. Dewey,et al.  De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis , 2013, Nature Protocols.

[29]  P. Benfey,et al.  Integrated detection of natural antisense transcripts using strand-specific RNA sequencing data , 2013, Genome research.

[30]  J. Rinn,et al.  Ribosome profiling reveals resemblance between long non-coding RNAs and 5′ leaders of coding RNAs , 2013, Development.

[31]  J. Kendall,et al.  The maize methylome influences mRNA splice sites and reveals widespread paramutation-like switches guided by small RNA , 2013, Genome research.

[32]  Nicholas T. Ingolia,et al.  Ribosome Profiling Provides Evidence that Large Noncoding RNAs Do Not Encode Proteins , 2013, Cell.

[33]  S. Puranik,et al.  Epigenetic mechanisms of plant stress responses and adaptation , 2013, Plant Cell Reports.

[34]  Stefano Lonardi,et al.  Mechanisms of small RNA generation from cis-NATs in response to environmental and developmental cues. , 2013, Molecular plant.

[35]  Nathan M. Springer,et al.  Genomic Distribution of Maize Facultative Heterochromatin Marked by Trimethylation of H3K27[W] , 2013, Plant Cell.

[36]  Stephen H. Hughes,et al.  H3K4me3 Interactions with TAF3 Regulate Preinitiation Complex Assembly and Selective Gene Activation , 2013, Cell.

[37]  Xiaohong Yang,et al.  Genome-wide association study dissects the genetic architecture of oil biosynthesis in maize kernels , 2012, Nature Genetics.

[38]  Jina,et al.  Mechanisms of small RNA generation from cis-NATs in response to environmental and developmental cues. , 2013, Molecular plant.

[39]  Thomas R. Gingeras,et al.  STAR: ultrafast universal RNA-seq aligner , 2013, Bioinform..

[40]  Ying Lu,et al.  Strand-specific RNA-seq reveals widespread occurrence of novel cis- natural antisense transcripts in rice , 2012, BMC Genomics.

[41]  Manel Esteller,et al.  Cis-acting noncoding RNAs: friends and foes , 2012, Nature Structural &Molecular Biology.

[42]  Zhe Wu,et al.  Flowering time control: another window to the connection between antisense RNA and chromatin. , 2012, Trends in genetics : TIG.

[43]  Wei Li,et al.  RSeQC: quality control of RNA-seq experiments , 2012, Bioinform..

[44]  Anna M. McGeachy,et al.  The ribosome profiling strategy for monitoring translation in vivo by deep sequencing of ribosome-protected mRNA fragments , 2012, Nature Protocols.

[45]  Arjun Krishnan,et al.  Effects of Drought on Gene Expression in Maize Reproductive and Leaf Meristem Tissue Revealed by RNA-Seq1[W][OA] , 2012, Plant Physiology.

[46]  A. Berr,et al.  Chromatin modification and remodelling: a regulatory landscape for the control of Arabidopsis defence responses upon pathogen attack , 2012, Cellular microbiology.

[47]  M. Nishiyama,et al.  Identification of sense and antisense transcripts regulated by drought in sugarcane , 2012, Plant Molecular Biology.

[48]  Steven L Salzberg,et al.  Fast gapped-read alignment with Bowtie 2 , 2012, Nature Methods.

[49]  O. Mittelsten Scheid,et al.  Stress-Induced Chromatin Changes: A Critical View on Their Heritability , 2012, Plant & cell physiology.

[50]  Shang Gao,et al.  Genome-wide analysis of plant nat-siRNAs reveals insights into their distribution, biogenesis and function , 2012, Genome Biology.

[51]  David R. Kelley,et al.  Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks , 2012, Nature Protocols.

[52]  Yu Zhang,et al.  Miniature Inverted–Repeat Transposable Elements (MITEs) Have Been Accumulated through Amplification Bursts and Play Important Roles in Gene Expression and Species Diversity in Oryza sativa , 2011, Molecular biology and evolution.

[53]  Damon Lisch,et al.  How important are transposons for plant evolution? , 2012, Nature Reviews Genetics.

[54]  Lin Fang,et al.  Resequencing 50 accessions of cultivated and wild rice yields markers for identifying agronomically important genes , 2011, Nature Biotechnology.

[55]  F. Ariel,et al.  Dual RNAs in plants. , 2011, Biochimie.

[56]  C. K. Vanderpool,et al.  Dual-function RNA regulators in bacteria. , 2011, Biochimie.

[57]  J. Araus,et al.  Large-scale screening for maize drought resistance using multiple selection criteria evaluated under water-stressed and well-watered environments , 2011 .

[58]  Ning Li,et al.  The panorama of physiological responses and gene expression of whole plant of maize inbred line YQ7-96 at the three-leaf stage under water deficit and re-watering , 2011, Theoretical and Applied Genetics.

[59]  Marcel Martin.  Cutadapt removes adapter sequences from high-throughput sequencing reads , 2011 .

[60]  Felix Krueger,et al.  Bismark: a flexible aligner and methylation caller for Bisulfite-Seq applications , 2011, Bioinform..

[61]  Stéphane Robin,et al.  Integrative epigenomic mapping defines four main chromatin states in Arabidopsis , 2011, The EMBO journal.

[62]  Howard Y. Chang,et al.  A long noncoding RNA maintains active chromatin to coordinate homeotic gene expression , 2011, Nature.

[63]  O. Clarenz,et al.  Dynamic Regulation of H3K27 Trimethylation during Arabidopsis Differentiation , 2011, PLoS genetics.

[64]  Wolfgang Huber,et al.  Antisense expression increases gene expression variability and locus interdependency , 2011, Molecular systems biology.

[65]  Lovelace J. Luquette,et al.  Comprehensive analysis of the chromatin landscape in Drosophila , 2010, Nature.

[66]  Hideaki Sugawara,et al.  The Sequence Read Archive , 2010, Nucleic Acids Res..

[67]  Dynamics of phytohormone and DNA methylation patterns changes during dormancy induction in strawberry (Fragaria × ananassa Duch.) , 2011, Plant Cell Reports.

[68]  P. Kapranov,et al.  Comprehensive Polyadenylation Site Maps in Yeast and Human Reveal Pervasive Alternative Polyadenylation , 2010, Cell.

[69]  Jean-Marcel Ribaut,et al.  Joint linkage–linkage disequilibrium mapping is a powerful approach to detecting quantitative trait loci underlying drought tolerance in maize , 2010, Proceedings of the National Academy of Sciences.

[70]  N. Friedman,et al.  Strand-specific RNA sequencing reveals extensive regulated long antisense transcripts that are conserved across yeast species , 2010, Genome Biology.

[71]  H. Hakonarson,et al.  ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data , 2010, Nucleic acids research.

[72]  Aaron R. Quinlan,et al.  Bioinformatics Applications Note Genome Analysis Bedtools: a Flexible Suite of Utilities for Comparing Genomic Features , 2022 .

[73]  C. Lister,et al.  Targeted 3′ Processing of Antisense Transcripts Triggers Arabidopsis FLC Chromatin Silencing , 2010, Science.

[74]  Mark D. Robinson,et al.  edgeR: a Bioconductor package for differential expression analysis of digital gene expression data , 2009, Bioinform..

[75]  C. Dean,et al.  Cold-induced silencing by long antisense transcripts of an Arabidopsis Polycomb target , 2009, Nature.

[76]  C. Wahlestedt,et al.  Regulatory roles of natural antisense transcripts , 2009, Nature Reviews Molecular Cell Biology.

[77]  D. Spector,et al.  Long noncoding RNAs: functional surprises from the RNA world. , 2009, Genes & development.

[78]  Gonçalo R. Abecasis,et al.  The Sequence Alignment/Map format and SAMtools , 2009, Bioinform..

[79]  X. Liu,et al.  Genome-Wide and Organ-Specific Landscapes of Epigenetic Modifications and Their Relationships to mRNA and Small RNA Transcriptomes in Maize[W] , 2009, The Plant Cell Online.

[80]  David W. Galbraith,et al.  Evidence for Antisense Transcription Associated with MicroRNA Target mRNAs in Arabidopsis , 2009, PLoS genetics.

[81]  M. Settles,et al.  Large-scale analysis of antisense transcription in wheat using the Affymetrix GeneChip Wheat Genome Array , 2009, BMC Genomics.

[82]  Sarah A Kessans,et al.  Small-interfering RNAs from natural antisense transcripts derived from a cellulose synthase gene modulate cell wall biosynthesis in barley , 2008, Proceedings of the National Academy of Sciences.

[83]  Leighton J. Core,et al.  Nascent RNA Sequencing Reveals Widespread Pausing and Divergent Initiation at Human Promoters , 2008, Science.

[84]  Jeannie T. Lee,et al.  Antagonism between DNA and H3K27 Methylation at the Imprinted Rasgrf1 Locus , 2008, PLoS genetics.

[85]  Vladimir Vacic,et al.  Small RNAs and the regulation of cis-natural antisense transcripts in Arabidopsis , 2008, BMC Molecular Biology.

[86]  Edward S. Buckler,et al.  TASSEL: software for association mapping of complex traits in diverse samples , 2007, Bioinform..

[87]  Bing Li,et al.  The Role of Chromatin during Transcription , 2007, Cell.

[88]  Hailing Jin,et al.  A pathogen-inducible endogenous siRNA in plant immunity , 2006, Proceedings of the National Academy of Sciences.

[89]  L. Xiong,et al.  Identification of Drought Tolerance Determinants by Genetic Analysis of Root Response to Drought Stress and Abscisic Acid1 , 2006, Plant Physiology.

[90]  X. Shirley Liu,et al.  Genome-wide in silico identification and analysis of cis natural antisense transcripts (cis-NATs) in ten species , 2006, Nucleic acids research.

[91]  O. Borsani,et al.  Endogenous siRNAs Derived from a Pair of Natural cis-Antisense Transcripts Regulate Salt Tolerance in Arabidopsis , 2005, Cell.

[92]  L. Hurst,et al.  Genome-wide analysis of coordinate expression and evolution of human cis-encoded sense-antisense transcripts. , 2005, Trends in genetics : TIG.

[93]  Thomas D. Wu,et al.  GMAP: a genomic mapping and alignment program for mRNA and EST sequence , 2005, Bioinform..

[94]  Sean R. Eddy,et al.  Rfam: annotating non-coding RNAs in complete genomes , 2004, Nucleic Acids Res..

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

[96]  Joseph M. Dale,et al.  Empirical Analysis of Transcriptional Activity in the Arabidopsis Genome , 2003, Science.

[97]  D. Higgs,et al.  Transcription of antisense RNA leading to gene silencing and methylation as a novel cause of human genetic disease , 2003, Nature Genetics.

[98]  John M. Hancock,et al.  PlantProm: a database of plant promoter sequences , 2003, Nucleic Acids Res..

[99]  Peer Bork,et al.  Systematic identification of novel protein domain families associated with nuclear functions. , 2002, Genome research.