Search

[This corrects the article DOI: 10.1371/journal.pgen.1003202.].

Salinity and drought stress are the primary cause of crop losses worldwide. In sodic saline soils sodium chloride (NaCl) disrupts normal plant growth and development. The complex interactions of plant systems with abiotic stress have made RNA sequencing a more holistic and appealing approach to study transcriptome level responses in a single cell and/or tissue. In this work, we determined the Petunia transcriptome response to NaCl stress by sequencing leaf samples and assembling 196 million Illumina reads with Trinity software. Using our reference transcriptome we identified more than 7,000 genes that were differentially expressed within 24 h of acute NaCl stress. The proposed transcriptome can also be used as an excellent tool for biological and bioinformatics in the absence of an available Petunia genome and it is available at the SOL Genomics Network (SGN) http://solgenomics.net. Genes related to regulation of reactive oxygen species, transport, and signal transductions as well as novel and undescribed transcripts were among those differentially expressed in response to salt stress. The candidate genes identified in this study can be applied as markers for breeding or to genetically engineer plants to enhance salt tolerance. Gene Ontology analyses indicated that most of the NaCl damage happened at 24 h inducing genotoxicity, affecting transport and organelles due to the high concentration of Na+ ions. Finally, we report a modification to the library preparation protocol whereby cDNA samples were bar-coded with non-HPLC purified primers, without affecting the quality and quantity of the RNA-seq data. The methodological improvement presented here could substantially reduce the cost of sample preparation for future high-throughput RNA sequencing experiments.

Transcriptome variation plays an important role in affecting the phenotype of an organism. However, an understanding of the underlying mechanisms regulating transcriptome variation in segregating populations is still largely unknown. We sought to assess and map variation in transcript abundance in maize shoot apices in the intermated B73 × Mo17 recombinant inbred line population. RNA-based sequencing (RNA-seq) allowed for the detection and quantification of the transcript abundance derived from 28,603 genes. For a majority of these genes, the population mean, coefficient of variation, and segregation patterns could be predicted by the parental expression levels. Expression quantitative trait loci (eQTL) mapping identified 30,774 eQTL including 96 trans-eQTL "hotspots," each of which regulates the expression of a large number of genes. Interestingly, genes regulated by a trans-eQTL hotspot tend to be enriched for a specific function or act in the same genetic pathway. Also, genomic structural variation appeared to contribute to cis-regulation of gene expression. Besides genes showing Mendelian inheritance in the RIL population, we also found genes whose expression level and variation in the progeny could not be predicted based on parental difference, indicating that non-Mendelian factors also contribute to expression variation. Specifically, we found 145 genes that show patterns of expression reminiscent of paramutation such that all the progeny had expression levels similar to one of the two parents. Furthermore, we identified another 210 genes that exhibited unexpected patterns of transcript presence/absence. Many of these genes are likely to be gene fragments resulting from transposition, and the presence/absence of their transcripts could influence expression levels of their ancestral syntenic genes. Overall, our results contribute to the identification of novel expression patterns and broaden the understanding of transcriptional variation in plants.

The shoot apical meristem (SAM) maintains a pool of indeterminate cells within the SAM proper, while lateral organs are initiated from the SAM periphery. Laser microdissection-microarray technology was used to compare transcriptional profiles within these SAM domains to identify novel maize genes that function during leaf development. Nine hundred and sixty-two differentially expressed maize genes were detected; control genes known to be upregulated in the initiating leaf (P0/P1) or in the SAM proper verified the precision of the microdissections. Genes involved in cell division/growth, cell wall biosynthesis, chromatin remodeling, RNA binding, and translation are especially upregulated in initiating leaves, whereas genes functioning during protein fate and DNA repair are more abundant in the SAM proper. In situ hybridization analyses confirmed the expression patterns of six previously uncharacterized maize genes upregulated in the P0/P1. P0/P1-upregulated genes that were also shown to be downregulated in leaf-arrested shoots treated with an auxin transport inhibitor are especially implicated to function during early events in maize leaf initiation. Reverse genetic analyses of asceapen1 (asc1), a maize D4-cyclin gene upregulated in the P0/P1, revealed novel leaf phenotypes, less genetic redundancy, and expanded D4-CYCLIN function during maize shoot development as compared to Arabidopsis. These analyses generated a unique SAM domain-specific database that provides new insight into SAM function and a useful platform for reverse genetic analyses of shoot development in maize.

Microarrays enable comparative analyses of gene expression on a genomic scale, however these experiments frequently identify an abundance of differentially expressed genes such that it may be difficult to identify discrete functional networks that are hidden within large microarray datasets. Microarray analyses in which mutant organisms are compared to nonmutant siblings can be especially problematic when the gene of interest is expressed in relatively few cells. Here, we describe the use of laser microdissection microarray to perform transcriptional profiling of the maize shoot apical meristem (SAM), a ~100-microm pillar of organogenic cells that is required for leaf initiation. Microarray analyses compared differential gene expression within the SAM and incipient leaf primordium of nonmutant and narrow sheath mutant plants, which harbored mutations in the duplicate genes narrow sheath1 (ns1) and narrow sheath2 (ns2). Expressed in eight to ten cells within the SAM, ns1 and ns2 encode paralogous WUSCHEL1-like homeobox (WOX) transcription factors required for recruitment of leaf initials that give rise to a large lateral domain within maize leaves. The data illustrate the utility of laser microdissection-microarray analyses to identify a relatively small number of genes that are differentially expressed within the SAM. Moreover, these analyses reveal potentially conserved WOX gene functions and implicate specific hormonal and signaling pathways during early events in maize leaf development.