Identification of the causal mutation in early heading mutant of bread wheat (Triticum aestivum L.) using MutMap approach.

In bread wheat (Triticum aestivum L.), fine-tuning the heading time is essential to maximize grain yield. Photoperiod-1 (Ppd-1) and VERNALIZATION 1 (Vrn-1) are major genes affecting photoperiod sensitivity and vernalization requirements, respectively. These genes have predominantly governed heading timing. However, Ppd-1 and Vrn-1 significantly impact heading dates, necessitating another gene that can slightly modify heading dates for fine-tuning. In this study, we developed an early heading mutant from the ethyl methanesulfonate-mutagenized population of the Japanese winter wheat cultivar "Kitahonami." MutMap analysis identified a nonsense mutation in the clock component gene Wheat PHYTOCLOCK 1/LUX ARRHYTHMO (WPCL-D1) as the probable SNP responsible for the early heading mutant on chromosome 3D. Segregation analysis using F2 and F3 populations confirmed that plants carrying the wpcl-D1 allele headed significantly earlier than those with the functional WPCL-D1. The early heading mutant exhibited increased expression levels of Ppd-1 and circadian clock genes, such as WPCL1 and LATE ELONGATED HYPOCOTYL (LHY). Notably, the transcript accumulation levels of Ppd-A1 and Ppd-D1 were influenced by the copy number of the functional WPCL1 gene. These results suggest that a loss-of-function mutation in WPCL-D1 is the causal mutation for the early heading phenotype. Adjusting the functional copy number of WPCL1 will be beneficial in fine-tuning of heading dates. Supplementary Information: The online version contains supplementary material available at 10.1007/s11032-024-01478-5.

Series of Resistance Genes in Barley (Hordeum vulgare) that Control Barley Yellow Mosaic Virus Multiplication and the Root-to-Leaf Systemic Movement.

The infection of young winter barley (Hordeum vulgare L.) root system in winter by barley yellow mosaic virus (BaYMV) can lead to high yield losses. Resistance breeding is critical for managing this virus, but there are only a few reports on resistance genes that describe how the genes control BaYMV propagation and the systemic movement from the roots to the leaves. Here we report a real-time quantitative PCR analysis of the virus in barley roots and leaves carrying BaYMV resistance genes (rym1 to rym15 and an unknown gene) to elucidate the molecular mechanisms underlying the barley response to BaYMV. The resistance mechanism directly targets the virus. Moreover, the resistance genes/cultivars were classified into the following three groups according to their BaYMV titer: (i) immune (BaYMV was undetectable in the roots or leaves), (ii) partially immune (BaYMV was detected in the roots but not in the leaves), and (iii) susceptible (BaYMV was detected in the roots and leaves). Our results clarified the functions of the resistance genes in barley roots and leaves following a BaYMV infection. We anticipate our analysis to be a starting point for more understanding of the correspondence between resistance genes of Triticeae and the soil-borne viruses.

Mutations within the miR172 target site of wheat AP2 homoeologs regulate lodicule size and rachis internode length.

Closed fertilization in flowers, or cleistogamy, reduces the risk of fungal infection in Triticeae crops. In barley (Hordeum vulgare), cleistogamy is determined by a single recessive gene, cly1, which results from a single nucleotide polymorphism within the microRNA172 target site of the Apetala2 (AP2) transcription factor gene. The recessive cly1 allele negatively regulates the development of lodicules, keeping florets closed at anthesis. However, cleistogamy is not evident in hexaploid wheat (Triticum aestivum) cultivars. This study aimed at identifying mutations in wheat AP2 orthologs by ethyl methane sulfonate-induced mutagenesis and high-resolution melt analysis. Although flowers of AP2 mutants induced in the A and D genomes opened at anthesis, their lodicule size was significantly smaller, especially in the direction of depth, than that of wild-type plants. One of the mutants that carried a nucleotide replacement in AP2 from the D genome produced a compact spike caused by a substantial decrease in rachis internode length, analogous to the barley dense spike. Cleistogamous hexaploid wheat might be generated by combining effective mutant alleles of AP2-homoeologous genes.

Organ-enriched gene expression during floral morphogenesis in wild barley.

Floral morphology varies considerably between dicots and monocots. The ABCDE model explaining how floral organ development is controlled was formulated using core eudicots and applied to grass crops. Barley (Hordeum. vulgare) has unique floral morphogenesis. Wild barley (H. vulgare ssp. spontaneum), which is the immediate ancestor of cultivated barley (H. vulgare ssp. vulgare), contains a rich reservoir of genetic diversity. However, the wild barley genes involved in floral organ development are still relatively uncharacterized. In this study, we generated an organ-specific transcriptome atlas for wild barley floral organs. Genome-wide transcription profiles indicated that 22 838 protein-coding genes were expressed in at least one organ. These genes were grouped into seven clusters according to the similarities in their expression patterns. Moreover, 5619 genes exhibited organ-enriched expression, 677 of which were members of 47 transcription factor families. Gene ontology analyses suggested that the functions of the genes with organ-enriched expression influence the biological processes in floral organs. The co-expression regulatory network showed that the expression of 690 genes targeted by MADS-box proteins was highly positively correlated with the expression of ABCDE model genes during floral morphogenesis. Furthermore, the expression of 138 genes was specific to the wild barley OUH602 genome and not the Morex genome; most of these genes were highly expressed in the glume, awn, lemma, and palea. This study revealed the global gene expression patterns underlying floral morphogenesis in wild barley. On the basis of the study findings, a molecular mechanism controlling floral morphology in barley was proposed.

Genetic resistance in barley against Japanese soil-borne wheat mosaic virus functions in the roots.

Infection by the Japanese soil-borne wheat mosaic virus (JSBWMV) can lead to substantial losses in the grain yield of barley and wheat crops. While genetically based resistance to this virus has been documented, its mechanistic basis remains obscure. In this study, the deployment of a quantitative PCR assay showed that the resistance acts directly against the virus rather than by inhibiting the colonization of the roots by the virus' fungal vector Polymyxa graminis. In the susceptible barley cultivar (cv.) Tochinoibuki, the JSBWMV titre was maintained at a high level in the roots during the period December-April, and the virus was translocated from the root to the leaf from January onwards. In contrast, in the roots of both cv. Sukai Golden and cv. Haruna Nijo, the titre was retained at a low level, and translocation of the virus to the shoot was strongly suppressed throughout the host's entire life cycle. The roots of wild barley (Hordeum vulgare ssp. spontaneum) accession H602 responded in the early stages of infection similarly to those of the resistant cultivated forms, but the host was unable to suppress the translocation of the virus to the shoot from March onwards. The virus titre in the root was presumed to have been restricted by the action of the gene product of Jmv1 (on chromosome 2H), while the stochastic nature of the infection was suppressed by the action of that of Jmv2 (on chromosome 3H), a gene harbored by cv. Sukai Golden but not by either cv. Haruna Nijo or accession H602.

Wheat Ym2 originated from Aegilops sharonensis and confers resistance to soil-borne Wheat yellow mosaic virus infection to the roots.

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

Genome sequencing-based coverage analyses facilitate high-resolution detection of deletions linked to phenotypes of gamma-irradiated wheat mutants.

BACKGROUND: Gamma-irradiated mutants of Triticum aestivum L., hexaploid wheat, provide novel and agriculturally important traits and are used as breeding materials. However, the identification of causative genomic regions of mutant phenotypes is challenging because of the large and complicated genome of hexaploid wheat. Recently, the combined use of high-quality reference genome sequences of common wheat and cost-effective resequencing technologies has made it possible to evaluate genome-wide polymorphisms, even in complex genomes. RESULTS: To investigate whether the genome sequencing approach can effectively detect structural variations, such as deletions, frequently caused by gamma irradiation, we selected a grain-hardness mutant from the gamma-irradiated population of Japanese elite wheat cultivar "Kitahonami." The Hardness (Ha) locus, including the puroindoline protein-encoding genes Pina-D1 and Pinb-D1 on the short arm of chromosome 5D, primarily regulates the grain hardness variation in common wheat. We performed short-read genome sequencing of wild-type and grain-hardness mutant plants, and subsequently aligned their short reads to the reference genome of the wheat cultivar "Chinese Spring." Genome-wide comparisons of depth-of-coverage between wild-type and mutant strains detected ~ 130 Mbp deletion on the short arm of chromosome 5D in the mutant genome. Molecular markers for this deletion were applied to the progeny populations generated by a cross between the wild-type and the mutant. A large deletion in the region including the Ha locus was associated with the mutant phenotype, indicating that the genome sequencing is a powerful and efficient approach for detecting a deletion marker of a gamma-irradiated mutant phenotype. In addition, we investigated a pre-harvest sprouting tolerance mutant and identified a 67.8 Mbp deletion on chromosome 3B where Viviparous-B1 and GRAS family transcription factors are located. Co-dominant markers designed to detect the deletion-polymorphism confirmed the association with low germination rate, leading to pre-harvest sprouting tolerance. CONCLUSIONS: Short read-based genome sequencing of gamma-irradiated mutants facilitates the identification of large deletions linked to mutant phenotypes when combined with segregation analyses in progeny populations. This method allows effective application of mutants with agriculturally important traits in breeding using marker-assisted selection.

Genome-Wide Analysis of Snf2 Gene Family Reveals Potential Role in Regulation of Spike Development in Barley.

Sucrose nonfermenting 2 (Snf2) family proteins, as the catalytic core of ATP-dependent chromatin remodeling complexes, play important roles in nuclear processes as diverse as DNA replication, transcriptional regulation, and DNA repair and recombination. The Snf2 gene family has been characterized in several plant species; some of its members regulate flower development in Arabidopsis. However, little is known about the members of the family in barley (Hordeum vulgare). Here, 38 Snf2 genes unevenly distributed among seven chromosomes were identified from the barley (cv. Morex) genome. Phylogenetic analysis categorized them into 18 subfamilies. They contained combinations of 21 domains and consisted of 3 to 34 exons. Evolution analysis revealed that segmental duplication contributed predominantly to the expansion of the family in barley, and the duplicated gene pairs have undergone purifying selection. About eight hundred Snf2 family genes were identified from 20 barley accessions, ranging from 38 to 41 genes in each. Most of these genes were subjected to purification selection during barley domestication. Most were expressed abundantly during spike development. This study provides a comprehensive characterization of barley Snf2 family members, which should help to improve our understanding of their potential regulatory roles in barley spike development.

Two dominant genes in barley (Hordeum vulgare L.) complementarily encode perfect resistance to Japanese soil-borne wheat mosaic virus.

Japanese soil-borne wheat mosaic virus (Furovirus) is a damaging pathogen of wheat and barley. This virus can survive in the soil for several decades, so the deployment of resistant cultivars represents the only practical control measure. Here, a genetic analysis has identified two regions of the barley genome-one on chromosome 2H and the other on chromosome 3H-as harboring gene(s) encoding resistance to this virus. The joint presence of both loci, termed Jmv1 and Jmv2, made the plants essentially immune, with resistance being dominant over susceptibility at each locus. Phylogenetic analysis showed that the virus is not closely related to the type Furovirus species Soil-borne wheat mosaic virus. There was a difference between the RNA1- and RNA2-based phylogenies of the virus species in Furovirus implying the independent segregation of the virus subgenomes.

Comparative Analysis of Root Transcriptome Reveals Candidate Genes and Expression Divergence of Homoeologous Genes in Response to Water Stress in Wheat.

Crop cultivars with larger root systems have an increased ability to absorb water and nutrients under conditions of water deficit. To unravel the molecular mechanism of water-stress tolerance in wheat, we performed RNA-seq analysis on the two genotypes, Colotana 296-52 (Colotana) and Tincurrin, contrasting the root growth under polyethylene-glycol-induced water-stress treatment. Out of a total of 35,047 differentially expressed genes, 3692 were specifically upregulated in drought-tolerant Colotana under water stress. Transcription factors, pyrroline-5-carboxylate reductase and late-embryogenesis-abundant proteins were among upregulated genes in Colotana. Variant calling between Colotana and Tincurrin detected 15,207 SNPs and Indels, which may affect protein function and mediate the contrasting root length phenotype. Finally, the expression patterns of five triads in response to water, high-salinity, heat, and cold stresses were analyzed using qRT-PCR to see if there were differences in homoeologous gene expression in response to those conditions. The five examined triads showed variation in the contribution of homoeologous genes to water, high-salinity, heat, and cold stresses in the two genotypes. The variation of homoeologous gene expression in response to environmental stresses may enable plants to better cope with stresses in their natural environments.

Combined QTL mapping and RNA-Seq profiling reveals candidate genes associated with cadmium tolerance in barley.

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.

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.

Gene expression of rice seeds surviving 13- and 20-month exposure to space environment.

Rice seeds were exposed outside of the international space station to assess the risk of space environment exposure on gene expression associated with seed germination. The germination percentages of the space-stored and ground-stored seeds exposed for 13 months were 48 and 96% respectively. Those for 20 months were 7 and 76%, respectively. Germination was defined 3 days after imbibition, except for the space-stored seeds exposed for 20 months, which germinated 5 days after imbibition. Subsequent RNA-seq analyses of the dry seeds, germinated seeds, and roots and shoots of seedlings revealed that the mutation rates of mRNA sequences were not significantly different between space-stored and ground-stored samples exposed for 13 months and 20 months. In all, 4 and 16 transcripts of glycolysis-related genes were increased in the germinated seeds after 13-month and 20-month exposure, respectively. Also, 2 and 39 transcripts of long-lived mRNA required for germination were decreased more than 2-fold in the dry seeds after 13-month and 20-month exposure, respectively. These results suggest that damage to long-lived mRNA in seeds by a space environment delays and reduces germination.

Variation of DNA methylation patterns associated with gene expression in rice (Oryza sativa) exposed to cadmium.

We report genome-wide single-base resolution maps of methylated cytosines and transcriptome change in Cd-exposed rice. Widespread differences were identified in CG and non-CG methylation marks between Cd-exposed and Cd-free rice genomes. There are 2320 non-redundant differentially methylated regions detected in the genome. RNA sequencing revealed 2092 DNA methylation-modified genes differentially expressed under Cd exposure. More genes were found hypermethylated than those hypomethylated in CG, CHH and CHG (where H is A, C or T) contexts in upstream, gene body and downstream regions. Many of the genes were involved in stress response, metal transport and transcription factors. Most of the DNA methylation-modified genes were transcriptionally altered under Cd stress. A subset of loss of function mutants defective in DNA methylation and histone modification activities was used to identify transcript abundance of selected genes. Compared with wide type, mutation of MET1 and DRM2 resulted in general lower transcript levels of the genes under Cd stress. Transcripts of OsIRO2, OsPR1b and Os09g02214 in drm2 were significantly reduced. A commonly used DNA methylation inhibitor 5-azacytidine was employed to investigate whether DNA demethylation affected physiological consequences. 5-azacytidine provision decreased general DNA methylation levels of selected genes, but promoted growth of rice seedlings and Cd accumulation in rice plant.

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.

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.

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.

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.

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.