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).

Background: Microarray 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. Results: Thirty-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. Conclusion: Unannotated 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.

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.

The high toxicity of cadmium (Cd) and its ready uptake by plants has become a major agricultural problem. To investigate the genetic architecture and genetic regulation of Cd tolerance in barley, we conducted quantitative trait loci (QTL) analysis in the phenotypically polymorphic Oregon Wolfe Barley (OWB) mapping population, derived from a cross between Rec and Dom parental genotypes. Through evaluating the Cd tolerance of 87 available doubled haploid lines of the OWB mapping population at the seedling stage, one minor and one major QTL were detected on chromosomes 2H and 6H, respectively. For chlorosis and necrosis traits, the major QTL explained 47.24% and 38.59% of the phenotypic variance, respectively. RNA-Seq analysis of the parental seedlings under Cd treatment revealed 542 differentially expressed genes between Cd-tolerant Rec and Cd-susceptible Dom genotypes. By analyzing sequence variations in transcribed sequences of the parental genotypes, 155,654 SNPs and 1,525 InDels were identified between the two contrasting genotypes and may contribute to Cd tolerance. Finally, by integrating the data from the identified QTLs and RNA-Seq analysis, 16 Cd tolerance-related candidate genes were detected, nine of which were metal ion transporters. These results provide promising candidate genes for further gene cloning and improving Cd tolerance in barley.

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.

Wheat yellow mosaic virus (WYMV) is a pathogen transmitted into its host’s roots by the soil-borne vector Polymyxa graminis . Ym1 and Ym2 genes protect the host from the significant yield losses caused by the virus, but the mechanistic basis of these resistance genes remains poorly understood. Here, it has been shown that Ym1 and Ym2 act within the root either by hindering the initial movement of WYMV from the vector into the root and/or by suppressing viral multiplication. A mechanical inoculation experiment on the leaf revealed that the presence of Ym1 reduced viral infection incidence, rather than viral titer, while that of Ym2 was ineffective in the leaf. To understand the basis of the root specificity of the Ym2 product, the gene was isolated from bread wheat using a positional cloning approach. The candidate gene encodes a CC-NBS-LRR protein and it correlated allelic variation with respect to its sequence with the host’s disease response. Ym2 (B37500) and its paralog (B35800) are found in the near-relatives, respectively, Aegilops sharonensi s and Aegilops speltoides (a close relative of the donor of bread wheat’s B genome), while both sequences, in a concatenated state, are present in several accessions of the latter species . Structural diversity in Ym2 has been generated via translocation and recombination between the two genes and enhanced by the formation of a chimeric gene resulting from an intralocus recombination event. The analysis has revealed how the Ym2 region has evolved during the polyploidization events leading to the creation of cultivated wheat.