miR172 downregulates the translation of cleistogamy 1 in barley.

Background and Aims: Floret opening in barley is induced by the swelling of the lodicule, a trait under the control of the cleistogamy1 (cly1) gene. The product of cly1 is a member of the APETALA2 (AP2) transcription factor family, which inhibits lodicule development. A sequence polymorphism at the miR172 target site within cly1 has been associated with variation in lodicule development and hence with the cleistogamous phenotype. It was unclear whether miR172 actually functions in cly1 regulation and, if it does, which miR172 gene contributes to cleistogamy. It was also interesting to explore whether miR172-mediated cly1 regulation occurs at transcriptional level or at translational level. Methods: Deep sequencing of small RNA identified the miR172 sequences expressed in barley immature spikes. miR172 genes were confirmed by computational and expression analysis. miR172 and cly1 expression profiles were determined by in situ hybridization and quantitative expression analysis. Immunoblot analysis provided the CLY1 protein quantifications. Definitive evidence of the role of miR172 in cleistogamy was provided by a transposon Ds-induced mutant of Hv-miR172a. Key Results: A small RNA analysis of the immature barley spike revealed three isomers, miR172a, b and c, of which miR172a was the most abundant. In situ hybridization analysis showed that miR172 and cly1 co-localize in the lodicule primordium, suggesting that these two molecules potentially interact with one another. Immunoblot analysis showed that the sequence polymorphism at the miR172 target site within cly1 reduced the abundance of the CLY1 protein, but not that of its transcript. In a Ds-induced mutant of Hv-miR172a, which generates no mature miR172a, the lodicules fail to grow, resulting in a very small lodicule. Conclusions: Direct evidence is presented to show that miR172a acts to reduce the abundance of the CLY1 protein, which enables open flowering in barley.

Genome-wide analysis of rice cis-natural antisense transcription under cadmium exposure using strand-specific RNA-Seq.

BACKGROUND: The elucidation of novel transcripts and their expression in response to various stress conditions is necessary to understand the transcriptional network of plants as an adaptation to biotic and abiotic stresses. We performed strand-specific RNA-Seq (ssRNA-Seq) on rice exposed to cadmium (Cd) for 24 h and investigated the expression of cis-natural antisense transcripts (cis-NATs), a class of endogenous coding or non-protein-coding RNAs with sequence complementarity to the opposite strands of RAP transcripts. RESULTS: Many RAP transcripts possessed cis-NATs and these cis-NATs were responsive to some extent. Cis-NATs were upregulated from 26, 266 and 409 RAP gene loci, while 2054, 2501 and 2825 RAP transcripts were upregulated from 38,123 RAP loci under high Cd exposure in roots at 1, 12 and 24 h, respectively. In addition, most of the upregulated cis-NATs showed little upregulation under ABA or cold treatment. A number of cis-NATs were upregulated from less than 35 RAP gene loci in different tissue and time-point combinations under low Cd exposure, suggesting that cis-NATs respond to environmental stress. Furthermore, 409 RAP transcripts with upregulated cis-NATs were classified into three groups based on the expression of the RAP transcripts from the opposite DNA strand, including 138 upregulated, 128 invariable, and 143 downregulated transcripts, although the responses of cis-NATs and RAP transcripts were not always correlated. CONCLUSIONS: We have shown that the cis-NATs identified by ssRNA-Seq analysis are novel genes and that some of them are stress-specific and show different responses depending on the degree of stress and tissue. These results improve our understanding of the complete molecular mechanism of plant adaptation to Cd exposure.

Expression of Flavone Synthase II and Flavonoid 3'-Hydroxylase Is Associated with Color Variation in Tan-Colored Injured Leaves of Sorghum.

Sorghum (Sorghum bicolor L. Moench) exhibits various color changes in injured leaves in response to cutting stress. Here, we aimed to identify key genes for the light brown and dark brown color variations in tan-colored injured leaves of sorghum. For this purpose, sorghum M36001 (light brown injured leaves), Nakei-MS3B (purple), and a progeny, #7 (dark brown), from Nakei-MS3B × M36001, were used. Accumulated pigments were detected by using high-performance liquid chromatography: M36001 accumulated only apigenin in its light brown leaves; #7 accumulated both luteolin and a small amount of apigenin in its dark brown leaves, and Nakei-MS3B accumulated 3-deoxyanthocyanidins (apigeninidin and luteolinidin) in its purple leaves. Apigenin or luteolin glucoside derivatives were also accumulated, in different proportions. Differentially expressed genes before and after cutting stress were identified by using RNA sequencing (RNA-seq). Integration of our metabolic and RNA-seq analyses suggested that expression of only flavone synthase II (FNSII) led to the synthesis of apigenin in M36001, expression of both FNSII and flavonoid 3'-hydroxylase (F3'H) led to the synthesis of apigenin and luteolin in #7, and expression of both flavanone 4-reductase and F3'H led to the synthesis of 3-deoxyanthocyanidins in Nakei-MS3B. These results suggest that expression of FNSII is related to the synthesis of flavones (apigenin and luteolin) and the expression level of F3'H is related to the balance of apigenin and luteolin. Expression of FNSII and F3'H is thus associated with dark or light brown coloration in tan-colored injured leaves of sorghum.

Genome-Wide Transcriptome Analysis of Cadmium Stress in Rice.

Rice growth is severely affected by toxic concentrations of the nonessential heavy metal cadmium (Cd). To elucidate the molecular basis of the response to Cd stress, we performed mRNA sequencing of rice following our previous study on exposure to high concentrations of Cd (Oono et al., 2014). In this study, rice plants were hydroponically treated with low concentrations of Cd and approximately 211 million sequence reads were mapped onto the IRGSP-1.0 reference rice genome sequence. Many genes, including some identified under high Cd concentration exposure in our previous study, were found to be responsive to low Cd exposure, with an average of about 11,000 transcripts from each condition. However, genes expressed constitutively across the developmental course responded only slightly to low Cd concentrations, in contrast to their clear response to high Cd concentration, which causes fatal damage to rice seedlings according to phenotypic changes. The expression of metal ion transporter genes tended to correlate with Cd concentration, suggesting the potential of the RNA-Seq strategy to reveal novel Cd-responsive transporters by analyzing gene expression under different Cd concentrations. This study could help to develop novel strategies for improving tolerance to Cd exposure in rice and other cereal crops.

Expression level of a flavonoid 3'-hydroxylase gene determines pathogen-induced color variation in sorghum.

BACKGROUND: Sorghum (Sorghum bicolor L. Moench) accumulates 3-deoxyanthocyanidins and exhibits orange to purple coloration on parts of the leaf in response to infection with the fungus Bipolaris sorghicola. We aimed to identify the key genes determining this color variation. RESULTS: Sorghum populations derived from Nakei-MS3B and M36001 accumulated apigeninidin, or both apigeninidin and luteolinidin, in different proportions in lesions caused by B. sorghicola infection, suggesting that the relative proportions of the two 3-deoxyanthocyanidins determine color variation. QTL analysis and genomic sequencing indicated that two closely linked loci on chromosome 4, containing the flavonoid 3'-hydroxylase (F3'H) and Tannin1 (Tan1) genes, were responsible for the lesion color variation. The F3'H locus in Nakei-MS3B had a genomic deletion resulting in the fusion of two tandemly arrayed F3'H genes. The recessive allele at the Tan1 locus derived from M36001 had a genomic insertion and encoded a non-functional WD40 repeat transcription factor. Whole-mRNA sequencing revealed that expression of the fused F3'H gene was conspicuously induced in purple sorghum lines. The levels of expression of F3'H matched the relative proportions of apigeninidin and luteolinidin. CONCLUSIONS: Expression of F3'H is responsible for the synthesis of luteolinidin; the expression level of this gene is therefore critical in determining color variation in sorghum leaves infected with B. sorghicola.

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.

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

BACKGROUND: The 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. RESULTS: At 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). CONCLUSIONS: Therefore, 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.

Genome-wide expression analysis of reactive oxygen species gene network in Mizuna plants grown in long-term spaceflight.

BACKGROUND: Spaceflight environment have been shown to generate reactive oxygen species (ROS) and induce oxidative stress in plants, but little is known about the gene expression of the ROS gene network in plants grown in long-term spaceflight. The molecular response and adaptation to the spaceflight environment of Mizuna plants harvested after 27 days of cultivation onboard the International Space Station (ISS) were measured using genome-wide mRNA expression analysis (mRNA-Seq). RESULTS: Total reads of transcripts from the Mizuna grown in the ISS as well as on the ground by mRNA-Seq showed 8,258 and 14,170 transcripts up-regulated and down-regulated, respectively, in the space-grown Mizuna when compared with those from the ground-grown Mizuna. A total of 20 in 32 ROS oxidative marker genes were up-regulated, including high expression of four hallmarks, and preferentially expressed genes associated with ROS-scavenging including thioredoxin, glutaredoxin, and alternative oxidase genes. In the transcription factors of the ROS gene network, MEKK1-MKK4-MPK3, OXI1-MKK4-MPK3, and OXI1-MPK3 of MAP cascades, induction of WRKY22 by MEKK1-MKK4-MPK3 cascade, induction of WRKY25 and repression of Zat7 by Zat12 were suggested. RbohD and RbohF genes were up-regulated preferentially in NADPH oxidase genes, which produce ROS. CONCLUSIONS: This large-scale transcriptome analysis revealed that the spaceflight environment induced oxidative stress and the ROS gene network activation in the space-grown Mizuna. Among transcripts altered in expression by space conditions, some were common genes response to abiotic and biotic stress. Furthermore, certain genes were exclusively up-regulated in Mizuna grown on the ISS. Surprisingly, Mizuna grew in space normally, as well as on the ground, demonstrating that plants can acclimate to long-term exposure in the spaceflight environment by reprogramming the expression of the ROS gene network.

A BAC physical map of aus rice cultivar 'Kasalath', and the map-based genomic sequence of 'Kasalath' chromosome 1.

Comparative analysis using available genomic resources within closely related species is an effective way to investigate genomic sequence and structural diversity. Rice (Oryza sativa L.) has undergone significant physiological and morphological changes during its domestication and local adaptation. We present a complete bacterial artificial chromosome (BAC) physical map for the aus rice cultivar 'Kasalath', which covers 90% of the sequence of temperate japonica rice cultivar 'Nipponbare'. Examination of physical distances between computational and experimental measurements of 'Kasalath' BAC insert size revealed the presence of more than 500 genomic regions that appear to have significant chromosome structural changes between the two cultivars. In particular, a genomic region on the long arm of 'Kasalath' chromosome 11 carrying a disease-resistance gene cluster was greatly expanded relative to the 'Nipponbare' genome. We also decoded 41.37 Mb of high-quality genomic sequence from 'Kasalath' chromosome 1. Extensive comparisons of chromosome 1 between 'Kasalath' and 'Nipponbare' led to the discovery of 317,843 single-nucleotide polymorphisms (SNPs) and 66,331 insertion/deletion (indel) sites. Nearly two-thirds of the expressed genes on rice chromosome 1 carried natural variations involving SNPs and/or indels that resulted in substitutions, insertions or deletions of amino acids in one cultivar relative to the other. We also observed gain and loss of genes caused by large indels. This study provides an important framework and an invaluable dataset for further understanding of the molecular mechanisms underlying the evolution and functions of the rice genome.

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.

Simultaneous transcriptome analysis of Sorghum and Bipolaris sorghicola by using RNA-seq in combination with de novo transcriptome assembly.

The recent development of RNA sequencing (RNA-seq) technology has enabled us to analyze the transcriptomes of plants and their pathogens simultaneously. However, RNA-seq often relies on aligning reads to a reference genome and is thus unsuitable for analyzing most plant pathogens, as their genomes have not been fully sequenced. Here, we analyzed the transcriptomes of Sorghum bicolor (L.) Moench and its pathogen Bipolaris sorghicola simultaneously by using RNA-seq in combination with de novo transcriptome assembly. We sequenced the mixed transcriptome of the disease-resistant sorghum cultivar SIL-05 and B. sorghicola in infected leaves in the early stages of infection (12 and 24 h post-inoculation) by using Illumina mRNA-Seq technology. Sorghum gene expression was quantified by aligning reads to the sorghum reference genome. For B. sorghicola, reads that could not be aligned to the sorghum reference genome were subjected to de novo transcriptome assembly. We identified genes of B. sorghicola for growth of this fungus in sorghum, as well as genes in sorghum for the defense response. The genes of B. sorghicola included those encoding Woronin body major protein, LysM domain-containing intracellular hyphae protein, transcriptional factors CpcA and HacA, and plant cell-wall degrading enzymes. The sorghum genes included those encoding two receptors of the simple eLRR domain protein family, transcription factors that are putative orthologs of OsWRKY45 and OsWRKY28 in rice, and a class III peroxidase that is a homolog involved in disease resistance in the Poaceae. These defense-related genes were particularly strongly induced among paralogs annotated in the sorghum genome. Thus, in the absence of genome sequences for the pathogen, simultaneous transcriptome analysis of plant and pathogen by using de novo assembly was useful for identifying putative key genes in the plant-pathogen interaction.

Characterisation of the wheat (Triticum aestivum L.) transcriptome by de novo assembly for the discovery of phosphate starvation-responsive genes: gene expression in Pi-stressed wheat.

BACKGROUND: Phosphorus (P) is an essential macronutrient for plant growth and development. To modulate their P homeostasis, plants must balance P uptake, mobilisation, and partitioning to various organs. Despite the worldwide importance of wheat as a cultivated food crop, molecular mechanisms associated with phosphate (Pi) starvation in wheat remain unclear. To elucidate these mechanisms, we used RNA-Seq methods to generate transcriptome profiles of the wheat variety 'Chinese Spring' responding to 10 days of Pi starvation. RESULTS: We carried out de novo assembly on 73.8 million high-quality reads generated from RNA-Seq libraries. We then constructed a transcript dataset containing 29,617 non-redundant wheat transcripts, comprising 15,047 contigs and 14,570 non-redundant full-length cDNAs from the TriFLDB database. When compared with barley full-length cDNAs, 10,656 of the 15,047 contigs were unalignable, suggesting that many might be distinct from barley transcripts. The average expression level of the contigs was lower than that of the known cDNAs, implying that these contigs included transcripts that were rarely represented in the full-length cDNA library. Within the non-redundant transcript set, we identified 892-2,833 responsive transcripts in roots and shoots, corresponding on average to 23.4% of the contigs not covered by cDNAs in TriFLDB under Pi starvation. The relative expression level of the wheat IPS1 (Induced by Phosphate Starvation 1) homologue, TaIPS1, was 341-fold higher in roots and 13-fold higher in shoots; this finding was further confirmed by qRT-PCR analysis. A comparative analysis of the wheat- and rice-responsive transcripts for orthologous genes under Pi-starvation revealed commonly upregulated transcripts, most of which appeared to be involved in a general response to Pi starvation, namely, an IPS1-mediated signalling cascade and its downstream functions such as Pi remobilisation, Pi uptake, and changes in Pi metabolism. CONCLUSIONS: Our transcriptome profiles demonstrated the impact of Pi starvation on global gene expression in wheat. This study revealed that enhancement of the Pi-mediated signalling cascade using IPS1 is a potent adaptation mechanism to Pi starvation that is conserved in both wheat and rice and validated the effectiveness of using short-read next-generation sequencing data for wheat transcriptome analysis in the absence of reference genome information.

Two-wavelength laser diode interferometer with time-sharing sinusoidal phase modulation.

We describe an interferometer system that uses two separate wavelengths to measure step height. The overlapping interference images detected by a CCD camera are easily separated by an ordinary integrating-bucket method and time-sharing sinusoidal phase modulation, in which two laser diodes are alternately modulated with a sinusoidal signal. A phase map is obtained only for the laser diode into which the modulation signal is injected. In this instance, a 1-microm step height was accurately detected.