Youko Oono

Monitoring expression profiles of Arabidopsis genes during cold acclimation and deacclimation using DNA microarrays

A comparative analysis of gene expression profiles during cold acclimation and deacclimation is necessary to elucidate the molecular mechanisms of cold stress responses in higher plants. We analyzed gene expression profiles in the process of cold acclimation and deacclimation (recovery from cold stress) using two microarray systems, the 7K RAFL cDNA microarray and the Agilent 22K oligonucleotide array. By both microarray analyses, we identified 292 genes up-regulated and 320 genes down-regulated during deacclimation, and 445 cold up-regulated genes and 341 cold down-regulated genes during cold acclimation. Many genes up-regulated during deacclimation were found to be down-regulated during cold acclimation, and vice versa. The genes up-regulated during deacclimation were classified into (1) regulatory proteins involved in further regulation of signal transduction and gene expression and (2) functional proteins involved in the recovery process from cold-stress-induced damages and plant growth. We also applied expression profiling studies to identify the key genes involved in the biosynthesis of carbohydrates and amino acids that are known to play important roles in cold acclimation. We compared genes that are regulated during deacclimation with those regulated during rehydration after dehydration to discuss the similarity and difference of each recovery process.

Expression of ER quality control-related genes in response to changes in BiP1 levels in developing rice endosperm.

Binding protein (BiP) is the key chaperone involved in folding of secretory proteins such as seed storage proteins in the ER lumen. To obtain functional information about BiP1, a gene that is predominantly expressed during rice seed maturation, we generated several transgenic rice plants in which various levels of BiP1 protein accumulated in an endosperm-specific manner. Severe suppression (BiP1 KD) or significant over-expression (BiP1 OEmax) of BiP1 not only altered seed phenotype and the intracellular structure of endosperm cells, but also reduced seed storage protein content, starch accumulation and grain weight. Microarray and RT-PCR analyses indicated that expression of many chaperone and co-chaperone genes was induced in transgenic plants, with more prominent expression in the BiP1 KD line than in the BiP1 OEmax line. Transcriptional induction of most chaperones was observed in calli treated with dithiothreitol or tunicamycin, treatments that trigger ER stress, indicating that induction of the chaperone genes in transgenic rice was caused by an ER stress response. In transient assays using rice protoplasts, the ortholog (Os06g0622700) of the AtbZIP60 transcription factor was shown to be involved in activation of some chaperone genes. Slight increases in the BiP1 level compared with wild-type, accompanied by increased levels of calnexin and protein disulfide isomerase-like proteins, resulted in significant enhancement of seed storage protein content, without any change in intracellular structure or seed phenotype. Judicious modification of BiP1 levels in transgenic rice can provide suitable conditions for the production of secretory proteins by alleviating ER stress.

Analysis of ER stress in developing rice endosperm accumulating β‐amyloid peptide

The common neurodegenerative disorder known as Alzheimers disease is characterized by cerebral neuritic plaques of amyloid beta (Abeta) peptide. Plaque formation is related to the highly aggregative property of this peptide, because it polymerizes to form insoluble plaques or fibrils causing neurotoxicity. Here, we expressed Abeta peptide as a new causing agent to endoplasmic reticulum (ER) stress to study ER stress occurred in plant. When the dimer of Abeta(1-42) peptide was expressed in maturing seed under the control of the 2.3-kb glutelin GluB-1 promoter containing its signal peptide, a maximum of about 8 mug peptide per grain accumulated and was deposited at the periphery of distorted ER-derived PB-I protein bodies. Synthesis of Abeta peptide in the ER lumen severely inhibited the synthesis and deposition of seed storage proteins, resulting in the generation of many small and abnormally appearing PB bodies. This ultrastructural change was accounted for by ER stress leading to the accumulation of aggregated Abeta peptide in the ER lumen and a coordinated increase in ER-resident molecular chaperones such as BiPs and PDIs in Abeta-expressing plants. Microarray analysis also confirmed that expression of several BiPs, PDIs and OsbZIP60 containing putative transmembrane domains was affected by the ER stress response. Abeta-expressing transgenic rice kernels exhibited an opaque and shrunken phenotype. When grain phenotype and expression levels were compared among transgenic rice grains expressing several different recombinant peptides, such detrimental effects on grain phenotype were correlated with the expressed peptide causing ER stress rather than expression levels.

Massive parallel sequencing of mRNA in identification of unannotated salinity stress-inducible transcripts in rice (Oryza sativa L.)

BackgroundMicroarray technology is limited to monitoring the expression of previously annotated genes that have corresponding probes on the array. Computationally annotated genes have not fully been validated, because ESTs and full-length cDNAs cannot cover entire transcribed regions. Here, mRNA-Seq (an Illumina cDNA sequencing application) was used to monitor whole mRNAs of salinity stress-treated rice tissues.ResultsThirty-six-base-pair reads from whole mRNAs were mapped to the rice genomic sequence: 72.0% to 75.2% were mapped uniquely to the genome, and 5.0% to 5.7% bridged exons. From the piling up of short reads mapped on the genome, a series of programs (Bowtie, TopHat, and Cufflinks) comprehensively predicted 51,301 (shoot) and 54,491 (root) transcripts, including 2,795 (shoot) and 3,082 (root) currently unannotated in the Rice Annotation Project database. Of these unannotated transcripts, 995 (shoot) and 1,052 (root) had ORFs similar to those encoding the amino acid sequences of functional proteins in a BLASTX search against UniProt and RefSeq databases. Among the unannotated genes, 213 (shoot) and 436 (root) were differentially expressed in response to salinity stress. Sequence-based and array-based measurements of the expression ratios of previously annotated genes were highly correlated.ConclusionUnannotated transcripts were identified on the basis of the piling up of mapped reads derived from mRNAs in rice. Some of these unannotated transcripts encoding putative functional proteins were expressed differentially in response to salinity stress.

Genome-wide transcriptome analysis reveals that cadmium stress signaling controls the expression of genes in drought stress signal pathways in rice.

Plant growth is severely affected by toxic concentrations of the non-essential heavy metal cadmium (Cd). Comprehensive transcriptome analysis by RNA-Seq following cadmium exposure is required to further understand plant responses to Cd and facilitate future systems-based analyses of the underlying regulatory networks. In this study, rice plants were hydroponically treated with 50 µM Cd for 24 hours and ∼60,000 expressed transcripts, including transcripts that could not be characterized by microarray-based approaches, were evaluated. Upregulation of various ROS-scavenging enzymes, chelators and metal transporters demonstrated the appropriate expression profiles to Cd exposure. Gene Ontology enrichment analysis of the responsive transcripts indicated the upregulation of many drought stress-related genes under Cd exposure. Further investigation into the expression of drought stress marker genes such as DREB suggested that expression of genes in several drought stress signal pathways was activated under Cd exposure. Furthermore, qRT-PCR analyses of randomly selected Cd-responsive metal transporter transcripts under various metal ion stresses suggested that the expression of Cd-responsive transcripts might be easily affected by other ions. Our transcriptome analysis demonstrated a new transcriptional network linking Cd and drought stresses in rice. Considering our data and that Cd is a non-essential metal, the network underlying Cd stress responses and tolerance, which plants have developed to adapt to other stresses, could help to acclimate to Cd exposure. Our examination of this transcriptional network provides useful information for further studies of the molecular mechanisms of plant adaptation to Cd exposure and the improvement of tolerance in crop species.

Construction of Pseudomolecule Sequences of the aus Rice Cultivar Kasalath for Comparative Genomics of Asian Cultivated Rice

Having a deep genetic structure evolved during its domestication and adaptation, the Asian cultivated rice (Oryza sativa) displays considerable physiological and morphological variations. Here, we describe deep whole-genome sequencing of the aus rice cultivar Kasalath by using the advanced next-generation sequencing (NGS) technologies to gain a better understanding of the sequence and structural changes among highly differentiated cultivars. The de novo assembled Kasalath sequences represented 91.1% (330.55 Mb) of the genome and contained 35 139 expressed loci annotated by RNA-Seq analysis. We detected 2 787 250 single-nucleotide polymorphisms (SNPs) and 7393 large insertion/deletion (indel) sites (>100 bp) between Kasalath and Nipponbare, and 2 216 251 SNPs and 3780 large indels between Kasalath and 93-11. Extensive comparison of the gene contents among these cultivars revealed similar rates of gene gain and loss. We detected at least 7.39 Mb of inserted sequences and 40.75 Mb of unmapped sequences in the Kasalath genome in comparison with the Nipponbare reference genome. Mapping of the publicly available NGS short reads from 50 rice accessions proved the necessity and the value of using the Kasalath whole-genome sequence as an additional reference to capture the sequence polymorphisms that cannot be discovered by using the Nipponbare sequence alone.

mRNA-Seq Reveals a Comprehensive Transcriptome Profile of Rice under Phosphate Stress

Plants have developed several morphological and physiological strategies to adapt to phosphate stress. We analyzed the inducible transcripts associated with phosphate starvation and over-abundant phosphate supply to characterize the transcriptome in rice seedlings using the mRNA-Seq strategy. Fifty-three million reads obtained from 16 libraries under various phosphate stress and recovery treatments were uniquely mapped to the rice genome. Transcripts identified specifically tagged to 40,574 (root) and 39,748 (shoot) Rice Annotation Project (RAP) transcripts. Additionally, we detected uniquely 10,388 transcripts with no match to any RAP transcript. These transcripts that showed specific response to Pi stress include those without ORFs that may act as non-protein coding transcripts. With an accompanying browser of the transcriptome under Pi stress, a deeper understanding of the structural and functional features of both annotated and unannotated Pi stress-responsive transcripts can provide useful information in improving Pi acquisition and utilization in rice and other cereal crops.

TENOR: Database for Comprehensive mRNA-Seq Experiments in Rice

Here we present TENOR (Transcriptome ENcyclopedia Of Rice, http://tenor.dna.affrc.go.jp), a database that encompasses large-scale mRNA sequencing (mRNA-Seq) data obtained from rice under a wide variety of conditions. Since the elucidation of the ability of plants to adapt to various growing conditions is a key issue in plant sciences, it is of great interest to understand the regulatory networks of genes responsible for environmental changes. We used mRNA-Seq and performed a time-course transcriptome analysis of rice, Oryza sativa L. (cv. Nipponbare), under 10 abiotic stress conditions (high salinity; high and low phosphate; high, low and extremely low cadmium; drought; osmotic; cold; and flood) and two plant hormone treatment conditions (ABA and jasmonic acid). A large number of genes that were responsive to abiotic stresses and plant hormones were detected by differential expression analysis. Furthermore, several responsive genes were found to encode transcription factors that could control the transcriptional network of stress responses, but the timing of the induction of these genes was not uniform across conditions. A significant number of cis-regulatory elements were enriched in the promoter regions of the responsive genes and were shared among conditions. These data suggest that some key components of gene regulation networks are shared between different stress signaling pathways. All the resources (novel genes identified from mRNA-Seq data, expression profiles, co-expressed genes and cis-regulatory elements) can be searched for and are available in TENOR.

RNA sequencing-mediated transcriptome analysis of rice plants in endoplasmic reticulum stress conditions

BackgroundThe endoplasmic reticulum (ER) stress response is widely known to function in eukaryotes to maintain the homeostasis of the ER when unfolded or misfolded proteins are overloaded in the ER. To understand the molecular mechanisms of the ER stress response in rice (Oryza sativa L.), we previously analyzed the expression profile of stably transformed rice in which an ER stress sensor/transducer OsIRE1 was knocked-down, using the combination of preliminary microarray and quantitative RT-PCR. In this study, to obtain more detailed expression profiles of genes involved in the initial stages of the ER stress response in rice, we performed RNA sequencing of wild-type and transgenic rice plants produced by homologous recombination in which endogenous genomic OsIRE1 was replaced by missense alleles defective in ribonuclease activity.ResultsAt least 38,076 transcripts were investigated by RNA sequencing, 380 of which responded to ER stress at a statistically significant level (195 were upregulated and 185 were downregulated). Furthermore, we successfully identified 17 genes from the set of 380 ER stress-responsive genes that were not included in the probe set of the currently available microarray chip in rice. Notably, three of these 17 genes were non-annotated genes, even in the latest version of the Rice Annotation Project Data Base (RAP-DB, version IRGSP-1.0).ConclusionsTherefore, RNA sequencing-mediated expression profiling provided valuable information about the ER stress response in rice plants and led to the discovery of new genes related to ER stress.

Diversity in the complexity of phosphate starvation transcriptomes among rice cultivars based on RNA-Seq profiles

Rice has developed several morphological and physiological strategies to adapt to phosphate starvation in the soil. In order to elucidate the molecular basis of response to phosphate starvation, we performed mRNA sequencing of 4 rice cultivars with variation in growth response to Pi starvation as indicated by the shoot/root dry weight ratio. Approximately 254 million sequence reads were mapped onto the IRGSP-1.0 reference rice genome sequence and an average of about 5,000 transcripts from each cultivar were found to be responsive under phosphate starvation. Comparative analysis of the RNA-Seq profiles of the 4 cultivars revealed similarities as well as distinct differences in expression of these responsive transcripts. We elucidated a set of core responsive transcripts including annotated and unannotated transcripts commonly expressed in the 4 cultivars but with different levels of expression. De novo assembly of unmapped reads to the Nipponbare genome generated a set of sequence contigs representing potential new transcripts that may be involved in tolerance to phosphate starvation. This study can be used for identification of genes and gene networks associated with environmental stress and the development of novel strategies for improving tolerance to phosphate starvation in rice and other cereal crops.