Diagnosing arsenic-mediated biochemical responses in rice cultivars using Raman spectroscopy

Rice (Oryza sativa) is the primary crop for nearly half of the world's population. Groundwater in many rice-growing parts of the world often has elevated levels of arsenite and arsenate. At the same time, rice can accumulate up to 20 times more arsenic compared to other staple crops. This places an enormous amount of people at risk of chronic arsenic poisoning. In this study, we investigated whether Raman spectroscopy (RS) could be used to diagnose arsenic toxicity in rice based on biochemical changes that were induced by arsenic accumulation. We modeled arsenite and arsenate stresses in four different rice cultivars grown in hydroponics over a nine-day window. Our results demonstrate that Raman spectra acquired from rice leaves, coupled with partial least squares-discriminant analysis, enabled accurate detection and identification of arsenic stress with approximately 89% accuracy. We also performed high-performance liquid chromatography (HPLC)-analysis of rice leaves to identify the key molecular analytes sensed by RS in confirming arsenic poisoning. We found that RS primarily detected a decrease in the concentration of lutein and an increase in the concentration of vanillic and ferulic acids due to the accumulation of arsenite and arsenate in rice. This showed that these molecules are detectable indicators of biochemical response to arsenic accumulation. Finally, a cross-correlation of RS with HPLC and ICP-MS demonstrated RS's potential for a label-free, non-invasive, and non-destructive quantification of arsenic accumulation in rice.

Genetic factors underlying anaerobic germination in rice: Genome-wide association study and transcriptomic analysis

The success of rice (Oryza sativa L.) germination and survival under submerged conditions is mainly determined by the rapid growth of the coleoptile to reach the water surface. Previous reports have shown the presence of genetic variability within rice accessions in the levels of flooding tolerance during germination or anaerobic germination (AG). Although many studies have focused on the physiological mechanisms of oxygen stress, few studies have explored the breadth of natural variation in AG tolerance-related traits in rice. In this study, we evaluated the coleoptile lengths of a geographically diverse rice panel of 241 accessions, including global accessions along with elite breeding lines and released cultivars from the United States, under the normal and flooded conditions in laboratory and greenhouse environments. A genome-wide association study (GWAS) was performed using a 7K single-nucleotide polymorphism (SNP) array and the phenotypic data of normal coleoptile length, flooded coleoptile length, flooding tolerance index, and survival at 14 d after seeding (DAS). Out of the 30 significant GWAS quantitative trait loci (QTL) regions identified, 14 colocalized with previously identified candidate genes of AG tolerance, whereas 16 were potentially novel. Two rice accessions showing contrasting phenotypic responses to AG stress were selected for the transcriptomics study. The combined approach of GWAS and transcriptomics analysis identified 77 potential candidate genes related to AG tolerance. The findings of our study may assist rice improvement programs in developing rice cultivars with robust tolerance under flooding stress during germination and the early seedling stage.

Guidelines for Performing CRISPR/Cas9 Genome Editing for Gene Validation and Trait Improvement in Crops

With the rapid advances in plant genome editing techniques over the past 10 years, more efficient and powerful crop genome editing applications are now possible. Candidate genes for key traits can be validated using CRISPR/Cas9-based knockouts and through the up- and down-regulation of gene expression. Likewise, new trait improvement approaches can take advantage of targeted editing to improve stress tolerance, disease resistance, and nutritional traits. However, several key steps in the process can prove tricky for researchers who might be new to plant genome editing. Here, we present step-by-step guidelines and best practices for a crop genome editing pipeline that should help to improve the rate of success. Important factors in the process include proper target sequence analysis and single guide RNA (sgRNA) design, sequencing of the target site in the genotypes of interest, performing an in vitro CRISPR/Cas9 ribonucleoprotein (RNP) assay to validate the designed sgRNAs, preparing the transformation constructs, considering a protoplast editing step as further validation, and, finally, stable plant transformation and mutation detection by Sanger and/or next-generation sequencing. With these detailed guidelines, a new user should be able to quickly set up a genome editing pipeline in their crop of interest and start making progress with the different CRISPR/Cas-based editing variants for gene validation and trait improvement purposes.

Comparative transcriptomic analysis of germinating rice seedlings to individual and combined anaerobic and cold stress

BACKGROUND

Rice is one of the most important cereals consumed worldwide. Two major abiotic factors affecting rice plants in different growth stages are flooding stress and cold stress. These abiotic stresses can take place independently or simultaneously and significantly affect rice plants during germination and seedling growth. Fortunately, a wide array of phenotypic responses conferring flooding stress and chilling stress tolerance exist within the rice germplasm, indicating the presence of different molecular mechanisms underlying tolerance to these stresses. Understanding these differences may assist in developing improved rice cultivars having higher tolerance to both stresses. In this study, we conducted a comparative global gene expression analysis of two rice genotypes with contrasting phenotypes under cold stress, anaerobic stress, and combined cold and anaerobic stress during germination.

RESULTS

The differential gene expression analysis revealed that 5571 differentially expressed genes (DEGs), 7206 DEGs, and 13279 DEGs were identified under anaerobic stress, cold stress, and combined stress, respectively. Genes involved in the carbohydrate metabolic process, glucosyltransferase activity, regulation of nitrogen compound metabolic process, protein metabolic process, lipid metabolic process, cellular nitrogen compound biosynthetic process, lipid biosynthetic process, and a microtubule-based process were enriched across all stresses. Notably, the common Gene Ontology (GO) analysis identified three hub genes, namely Os08g0176800 (similar to mRNA-associated protein mrnp 41), Os11g0454200 (dehydrin), and OS10g0505900 (expressed protein).

CONCLUSION

A large number of differentially expressed genes were identified under anaerobic, cold conditions during germination and the combination of the two stress conditions in rice. These results will assist in the identification of promising candidate genes for possible manipulation toward rice crops that are more tolerant under flooding and cold during germination, both independently and concurrently.

Optimization of gene editing in cowpea through protoplast transformation and agroinfiltration by targeting the phytoene desaturase gene

Cowpea (Vigna unguiculata) is a legume staple widely grown across Sub-Saharan Africa and other tropical and sub-tropical regions. Considering projected climate change and global population increases, cowpea's adaptation to hot climates, resistance to drought, and nitrogen-fixing capabilities make it an especially attractive crop for facing future challenges. Despite these beneficial traits, efficient varietal improvement is challenging in cowpea due to its recalcitrance to transformation and long regeneration times. Transient gene expression assays can provide solutions to alleviate these issues as they allow researchers to test gene editing constructs before investing in the time and resource- intensive process of transformation. In this study, we developed an improved cowpea protoplast isolation protocol, a transient protoplast assay, and an agroinfiltration assay to be used for initial testing and validation of gene editing constructs and for gene expression studies. To test these protocols, we assessed the efficacy of a CRISPR-Cas9 construct containing four multiplexed single-guide RNA (sgRNA) sequences using polyethylene glycol (PEG)-mediated transformation and agroinfiltration with phytoene desaturase (PDS) as the target gene. Sanger sequencing of DNA from transformed protoplasts and agroinfiltrated cowpea leaves revealed several large deletions in the target sequences. The protoplast system and agroinfiltration protocol developed in this study provide versatile tools to test gene editing components before initiating plant transformation, thus improving the chance of using active sgRNAs and attaining the desired edits and target phenotype.

Optimization of Prime Editing in Rice, Peanut, Chickpea, and Cowpea Protoplasts by Restoration of GFP Activity

Precise editing of the plant genome has long been desired for functional genomic research and crop breeding. Prime editing is a newly developed precise editing technology based on CRISPR-Cas9, which uses an engineered reverse transcriptase (RT), a catalytically impaired Cas9 endonuclease (nCas9), and a prime editing guide RNA (pegRNA). In addition, prime editing has a wider range of editing types than base editing and can produce nearly all types of edits. Although prime editing was first established in human cells, it has recently been applied to plants. As a relatively new technique, optimization will be needed to increase the editing efficiency in different crops. In this study, we successfully edited a mutant GFP in rice, peanut, chickpea, and cowpea protoplasts. In rice, up to 16 times higher editing efficiency was achieved with a dual pegRNA than the single pegRNA containing vectors. Edited-mutant GFP protoplasts have also been obtained in peanut, chickpea, and cowpea after transformation with the dual pegRNA vectors, albeit with much lower editing efficiency than in rice, ranging from 0.2% to 0.5%. These initial results promise to expedite the application of prime editing in legume breeding programs to accelerate crop improvement.

Increasing the level of resistant starch in 'Presidio' rice through multiplex CRISPR-Cas9 gene editing of starch branching enzyme genes

Rice (Oryza sativa L.) is an excellent source of starch, which is composed of amylopectin and amylose. Resistant starch (RS) is a starch product that is not easily digestible and absorbed in the stomach or small intestine and instead is passed on directly to the large intestine. Cereals high in RS may be beneficial to improve human health and reduce the risk of diet-related chronic diseases. It has been reported through chemical mutagenesis and RNA interference studies that starch branching enzymes (SBEs) play a major role in contributing to higher levels of RS in cereal crops. In this study, we used multiplex clustered regularly interspaced short palindromic repeat (CRISPR)-CRISPR associated protein 9 (Cas9) genome editing to simultaneously target all four SBE genes in rice using the endogenous transfer RNA (tRNA)-processing system for expressing the single-guide RNAs (sgRNAs) targeting these genes. The CRISPR-Cas9 vector construct with four SBE gene sgRNAs was transformed into the U.S. rice cultivar Presidio using Agrobacterium-mediated transformation. Knockout mutations were identified at all four SBE genes across eight transgene-positive T₀ plants. Transgene-free T₁ lines with different combinations of disrupted SBE genes were identified, with several SBE-edited lines showing significantly increased RS content up to 15% higher than the wild-type (WT) cultivar Presidio. Although further efforts are needed to fix all of the mutant alleles as homozygous, our study demonstrated the potential of multiplex genome editing to develop high-RS lines.

Functional Allele Validation by Gene Editing to Leverage the Wealth of Genetic Resources for Crop Improvement

Advances in molecular technologies over the past few decades, such as high-throughput DNA marker genotyping, have provided more powerful plant breeding approaches, including marker-assisted selection and genomic selection. At the same time, massive investments in plant genetics and genomics, led by whole genome sequencing, have led to greater knowledge of genes and genetic pathways across plant genomes. However, there remains a gap between approaches focused on forward genetics, which start with a phenotype to map a mutant locus or QTL with the goal of cloning the causal gene, and approaches using reverse genetics, which start with large-scale sequence data and work back to the gene function. The recent establishment of efficient CRISPR-Cas-based gene editing promises to bridge this gap and provide a rapid method to functionally validate genes and alleles identified through studies of natural variation. CRISPR-Cas techniques can be used to knock out single or multiple genes, precisely modify genes through base and prime editing, and replace alleles. Moreover, technologies such as protoplast isolation, in planta transformation, and the use of developmental regulatory genes promise to enable high-throughput gene editing to accelerate crop improvement.

Phenotypic variation and genome-wide association studies of main culm panicle node number, maximum node production rate, and degree-days to heading in rice

Background

Grain yield is a complex trait that results from interaction between underlying phenotypic traits and climatic, edaphic, and biotic variables. In rice, main culm panicle node number (MCPNN; the node number on which the panicle is borne) and maximum node production rate (MNPR; the number of leaves that emerge per degree-day > 10°C) are primary phenotypic plant traits that have significant positive direct effects on yield-related traits. Degree-days to heading (DDTH), which has a significant positive effect on grain yield, is influenced by the interaction between MCPNN and MNPR. The objective of this research is to assess the phenotypic variation of MCPNN, MNPR, and DDTH in a panel of diverse rice accessions, determine regions in the rice genome associated with these traits using genome-wide association studies (GWAS), and identify putative candidate genes that control these traits.

Results

Considerable variation was observed for the three traits in a 220-genotype diverse rice population. MCPNN ranged from 8.1 to 20.9 nodes in 2018 and from 9.9 to 21.0 nodes in 2019. MNPR ranged from 0.0097 to 0.0214 nodes/degree day > 10°C in 2018 and from 0.0108 to 0.0193 nodes/degree-day > 10°C in 2019. DDTH ranged from 713 to 2,345 degree-days > 10°C in 2018 and from 778 to 2,404 degree-days > 10°C in 2019.

Thirteen significant (P < 2.91 x 10^-7) trait-single nucleotide polymorphism (SNP) associations were identified using the multilocus mixed linear model for GWAS. Significant associations between MCPNN and three SNPs in chromosome 2 (S02_12032235, S02_11971745, and S02_12030176) were detected with both the 2018 and best linear unbiased prediction (BLUP) datasets. Nine SNPs in chromosome 6 (S06_1970442, S06_2310856, S06_2550351, S06_1968653, S06_2296852, S06_1968680, S06_1968681, S06_1970597, and S06_1970602) were significantly associated with MNPR in the 2019 dataset. One SNP in chromosome 11 (S11_29358169) was significantly associated with the DDTH in the BLUP dataset.

Conclusions

This study identifies SNP markers that are putatively associated with MCPNN, MNPR, and DDTH. Some of these SNPs were located within or near gene models, which identify possible candidate genes involved in these traits. Validation of the putative candidate genes through expression and gene editing analyses are necessary to confirm their roles in regulating MCPNN, MNPR, and DDTH. Identifying the underlying genetic basis for primary phenotypic traits MCPNN and MNPR could lead to the development of fast and efficient approaches for their estimation, such as marker-assisted selection and gene editing, which is essential in increasing breeding efficiency and enhancing grain yield in rice. On the other hand, DDTH is a resultant variable that is highly affected by nitrogen and water management, plant density, and several other factors.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08629-y.

Raman Spectroscopy Enables Non-invasive and Confirmatory Diagnostics of Aluminum and Iron Toxicities in Rice

Metal toxicities can be detrimental to a plant health, as well as to the health of animals and humans that consume such plants. Metal content of plants can be analyzed using colorimetric, atomic absorption- or mass spectroscopy-based methods. However, these techniques are destructive, costly and laborious. In the current study, we investigate the potential of Raman spectroscopy (RS), a modern spectroscopic technique, for detection and identification of metal toxicities in rice. We modeled medium and high levels of iron and aluminum toxicities in hydroponically grown plants. Spectroscopic analyses of their leaves showed that both iron and aluminum toxicities can be detected and identified with ∼100% accuracy as early as day 2 after the stress initiation. We also showed that diagnostics accuracy was very high not only on early, but also on middle (day 4-day 8) and late (day 10-day 14) stages of the stress development. Importantly this approach only requires an acquisition time of 1 s; it is non-invasive and non-destructive to plants. Our findings suggest that if implemented in farming, RS can enable pre-symptomatic detection and identification of metallic toxins that would lead to faster recovery of crops and prevent further damage.

Carbon Nanotube-Mediated Plasmid DNA Delivery in Rice Leaves and Seeds

CRISPR-Cas gene editing technologies offer the potential to modify crops precisely; however, in vitro plant transformation and regeneration techniques present a bottleneck due to the lengthy and genotype-specific tissue culture process. Ideally, in planta transformation can bypass tissue culture and directly lead to transformed plants, but efficient in planta delivery and transformation remains a challenge. This study investigates transformation methods that have the potential to directly alter germline cells, eliminating the challenge of in vitro plant regeneration. Recent studies have demonstrated that carbon nanotubes (CNTs) loaded with plasmid DNA can diffuse through plant cell walls, facilitating transient expression of foreign genetic elements in plant tissues. To test if this approach is a viable technique for in planta transformation, CNT-mediated plasmid DNA delivery into rice tissues was performed using leaf and excised-embryo infiltration with reporter genes. Quantitative and qualitative data indicate that CNTs facilitate plasmid DNA delivery in rice leaf and embryo tissues, resulting in transient GFP, YFP, and GUS expression. Experiments were also initiated with CRISPR-Cas vectors targeting the phytoene desaturase (PDS) gene for CNT delivery into mature embryos to create heritable genetic edits. Overall, the results suggest that CNT-based delivery of plasmid DNA appears promising for in planta transformation, and further optimization can enable high-throughput gene editing to accelerate functional genomics and crop improvement activities.

Optimization of Protoplast Isolation and Transformation for a Pilot Study of Genome Editing in Peanut by Targeting the Allergen Gene Ara h 2

The cultivated peanut (Arachis hypogaea L.) is a legume consumed worldwide in the form of oil, nuts, peanut butter, and candy. Improving peanut production and nutrition will require new technologies to enable novel trait development. Clustered regularly interspaced short palindromic repeats and CRISPR-associated protein 9 (CRISPR-Cas9) is a powerful and versatile genome-editing tool for introducing genetic changes for studying gene expression and improving crops, including peanuts. An efficient in vivo transient CRISPR-Cas9- editing system using protoplasts as a testbed could be a versatile platform to optimize this technology. In this study, multiplex CRISPR-Cas9 genome editing was performed in peanut protoplasts to disrupt a major allergen gene with the help of an endogenous tRNA-processing system. In this process, we successfully optimized protoplast isolation and transformation with green fluorescent protein (GFP) plasmid, designed two sgRNAs for an allergen gene, Ara h 2, and tested their efficiency by in vitro digestion with Cas9. Finally, through deep-sequencing analysis, several edits were identified in our target gene after PEG-mediated transformation in protoplasts with a Cas9 and sgRNA-containing vector. These findings demonstrated that a polyethylene glycol (PEG)-mediated protoplast transformation system can serve as a rapid and effective tool for transient expression assays and sgRNA validation in peanut.

Reference-Guided De Novo Genome Assembly to Dissect a QTL Region for Submergence Tolerance Derived from Ciherang-Sub1

Global climate change has increased the number of severe flooding events that affect agriculture, including rice production in the U.S. and internationally. Heavy rainfall can cause rice plants to be completely submerged, which can significantly affect grain yield or completely destroy the plants. Recently, a major effect submergence tolerance QTL during the vegetative stage, qSub8.1, which originated from Ciherang-Sub1, was identified in a mapping population derived from a cross between Ciherang-Sub1 and IR10F365. Ciherang-Sub1 was, in turn, derived from a cross between Ciherang and IR64-Sub1. Here, we characterize the qSub8.1 region by analyzing the sequence information of Ciherang-Sub1 and its two parents (Ciherang and IR64-Sub1) and compare the whole genome profile of these varieties with the Nipponbare and Minghui 63 (MH63) reference genomes. The three rice varieties were sequenced with 150 bp pair-end whole-genome shotgun sequencing (Illumina HiSeq4000), followed by performing the Trimmomatic-SOAPdenovo2-MUMmer3 pipeline for genome assembly, resulting in approximate genome sizes of 354.4, 343.7, and 344.7 Mb, with N50 values of 25.1, 25.4, and 26.1 kb, respectively. The results showed that the Ciherang-Sub1 genome is composed of 59–63% Ciherang, 22–24% of IR64-Sub1, and 15–17% of unknown sources. The genome profile revealed a more detailed genomic composition than previous marker-assisted breeding and showed that the qSub8.1 region is mostly from Ciherang, with some introgressed segments from IR64-Sub1 and currently unknown source(s).

Optimizing Agrobacterium-Mediated Transformation and CRISPR-Cas9 Gene Editing in the tropical japonica Rice Variety Presidio

Bottlenecks in plant transformation and regeneration have slowed progress in applying CRISPR/Cas-based genome editing for crop improvement. Rice (Oryza sativa L.) has highly efficient temperate japonica transformation protocols, along with reasonably efficient indica protocols using immature embryos. However, rapid and efficient protocols are not available for transformation and regeneration in tropical japonica varieties, even though they represent the majority of rice production in the U.S. and South America. The current study has optimized a protocol using callus induction from mature seeds with both Agrobacterium-mediated and biolistic transformation of the high-yielding U.S. tropical japonica cultivar Presidio. Gene editing efficiency was tested by evaluating knockout mutations in the phytoene desaturase (PDS) and young seedling albino (YSA) genes, which provide a visible phenotype at the seedling stage for successful knockouts. Using the optimized protocol, transformation of 648 explants with particle bombardment and 532 explants with Agrobacterium led to a 33% regeneration efficiency. The YSA targets had ambiguous phenotypes, but 60% of regenerated plants for PDS showed an albino phenotype. Sanger sequencing of edited progeny showed a number of insertions, deletions, and substitutions at the gRNA target sites. These results pave the way for more efficient gene editing of tropical japonica rice varieties.

OsCOP1 regulates embryo development and flavonoid biosynthesis in rice (Oryza sativa L.)

Novel mutations of OsCOP1 were identified to be responsible for yellowish pericarp and embryo lethal phenotype, which revealed that OsCOP1 plays a crucial role in flavonoid biosynthesis and embryogenesis in rice seed. Successful production of viable seeds is a major component of plant life cycles, and seed development is a complex, highly regulated process that affects characteristics such as seed viability and color. In this study, three yellowish-pericarp embryo lethal (yel) mutants, yel-hc, yel-sk, and yel-cc, were produced from three different japonica cultivars of rice (Oryza sativa L). Mutant seeds had yellowish pericarps and exhibited embryonic lethality, with significantly reduced grain size and weight. Morphological aberrations were apparent by 5 days after pollination, with abnormal embryo development and increased flavonoid accumulation observed in the yel mutants. Genetic analysis and mapping revealed that the phenotype of the three yel mutants was controlled by a single recessive gene, LOC_Os02g53140, an ortholog of Arabidopsis thaliana CONSTITUTIVE PHOTOMORPHOGENIC 1 (COP1). The yel-hc, yel-sk, and yel-cc mutants carried mutations in the RING finger, coiled-coil, and WD40 repeat domains, respectively, of OsCOP1. CRISPR/Cas9-targeted mutagenesis was used to knock out OsCOP1 by targeting its functional domains, and transgenic seed displayed the yel mutant phenotype. Overexpression of OsCOP1 in a homozygous yel-hc mutant background restored pericarp color, and the aberrant flavonoid accumulation observed in yel-hc mutant was significantly reduced in the embryo and endosperm. These results demonstrate that OsCOP1 is associated with embryo development and flavonoid biosynthesis in rice grains. This study will facilitate a better understanding of the functional roles of OsCOP1 involved in early embryogenesis and flavonoid biosynthesis in rice seeds.

Transcriptome profiling of two rice genotypes under mild field drought stress during grain-filling stage

Drought is one of the most critical abiotic stresses that threaten crop production worldwide. This stress affects the rice crop in all stages of rice development; however, the occurrence during reproductive and grain-filling stages has the most impact on grain yield. Although many global transcriptomic studies have been performed during the reproductive stage in rice, very limited information is available for the grain-filling stage. Hence, we intend to investigate how the rice plant responds to drought stress during the grain-filling stage and how the responses change over time under field conditions. Two rice genotypes were selected for RNA-seq analysis: '4610', previously reported as a moderately tolerant breeding line, and Rondo, an elite indica rice cultivar susceptible to drought conditions. Additionally, 10 agronomic traits were evaluated under normal irrigated and drought conditions. Leaf tissues were collected during grain-filling stages at two time points, 14 and 21 days after the drought treatment, from both the drought field and normal irrigated field conditions. Based on agronomic performances, '4610' was less negatively affected than Rondo under mild drought conditions, and expression profiling largely aligned with the phenotypic data. The transcriptomic data indicated that, in general, '4610' had much earlier responses than its counterpart in mitigating the impact of drought stress. Several key genes and gene families related to drought stress or stress-related conditions were found differentially expressed in this study, including transcription factors, drought tolerance genes and reactive oxygen species scavengers. Furthermore, this study identified novel differentially expressed genes (DEGs) without function annotations that may play roles in drought tolerance-related functions. Some of the important DEGs detected in this study can be targeted for future research.

Improved Transformation and Regeneration of Indica Rice: Disruption of SUB1A as a Test Case via CRISPR-Cas9

Gene editing by use of clustered regularly interspaced short palindromic repeats (CRISPR) has become a powerful tool for crop improvement. However, a common bottleneck in the application of this approach to grain crops, including rice (Oryza sativa), is efficient vector delivery and calli regeneration, which can be hampered by genotype-dependent requirements for plant regeneration. Here, methods for Agrobacterium-mediated and biolistic transformation and regeneration of indica rice were optimized using CRISPR-Cas9 gene-editing of the submergence tolerance regulator SUBMERGENCE 1A-1 gene of the cultivar Ciherang-Sub1. Callus induction and plantlet regeneration methods were optimized for embryogenic calli derived from immature embryos and mature seed-derived calli. Optimized regeneration (95%) and maximal editing efficiency (100%) were obtained from the immature embryo-derived calli. Phenotyping of T₁ seeds derived from the edited T₀ plants under submergence stress demonstrated inferior phenotype compared to their controls, which phenotypically validates the disruption of SUB1A-1 function. The methods pave the way for rapid CRISPR-Cas9 gene editing of recalcitrant indica rice cultivars.

Non-Invasive Identification of Nutrient Components in Grain

Digital farming is a modern agricultural concept that aims to maximize the crop yield while simultaneously minimizing the environmental impact of farming. Successful implementation of digital farming requires development of sensors to detect and identify diseases and abiotic stresses in plants, as well as to probe the nutrient content of seeds and identify plant varieties. Experimental evidence of the suitability of Raman spectroscopy (RS) for confirmatory diagnostics of plant diseases was previously provided by our team and other research groups. In this study, we investigate the potential use of RS as a label-free, non-invasive and non-destructive analytical technique for the fast and accurate identification of nutrient components in the grains from 15 different rice genotypes. We demonstrate that spectroscopic analysis of intact rice seeds provides the accurate rice variety identification in ~86% of samples. These results suggest that RS can be used for fully automated, fast and accurate identification of seeds nutrient components.

Utilizing PacBio Iso-Seq for Novel Transcript and Gene Discovery of Abiotic Stress Responses in Oryza sativa L

The wide natural variation present in rice is an important source of genes to facilitate stress tolerance breeding. However, identification of candidate genes from RNA-Seq studies is hampered by the lack of high-quality genome assemblies for the most stress tolerant cultivars. A more targeted solution is the reconstruction of transcriptomes to provide templates to map RNA-seq reads. Here, we sequenced transcriptomes of ten rice cultivars of three subspecies on the PacBio Sequel platform. RNA was isolated from different organs of plants grown under control and abiotic stress conditions in different environments. Reconstructed de novo reference transcriptomes resulted in 37,500 to 54,600 plant-specific high-quality isoforms per cultivar. Isoforms were collapsed to reduce sequence redundancy and evaluated, e.g., for protein completeness (BUSCO). About 40% of all identified transcripts were novel isoforms compared to the Nipponbare reference transcriptome. For the drought/heat tolerant aus cultivar N22, 56 differentially expressed genes in developing seeds were identified at combined heat and drought in the field. The newly generated rice transcriptomes are useful to identify candidate genes for stress tolerance breeding not present in the reference transcriptomes/genomes. In addition, our approach provides a cost-effective alternative to genome sequencing for identification of candidate genes in highly stress tolerant genotypes.

Raman Spectroscopy Enables Non-invasive and Confirmatory Diagnostics of Salinity Stresses, Nitrogen, Phosphorus, and Potassium Deficiencies in Rice

Proper management of nutrients in agricultural systems is critically important for maximizing crop yields while simultaneously minimizing the health and environmental impacts of pollution from fertilizers. These goals can be achieved by timely confirmatory diagnostics of nutrient deficiencies in plants, which enable precise administration of fertilizers and other supplementation in fields. Traditionally, nutrient diagnostics are performed by wet-laboratory analyses, which are both time- and labor-consuming. Unmanned aerial vehicle (UAV) and satellite imaging have offered a non-invasive alternative. However, these imaging approaches do not have sufficient specificity, and they are only capable of detecting symptomatic stages of nutrient deficiencies. Raman spectroscopy (RS) is a non-invasive and non-destructive technique that can be used for confirmatory detection and identification of both biotic and abiotic stresses on plants. Herein, we show the use of a hand-held Raman spectrometer for highly accurate pre-symptomatic diagnostics of nitrogen, phosphorus, and potassium deficiencies in rice (Oryza sativa). Moreover, we demonstrate that RS can also be used for pre symptomatic diagnostics of medium and high salinity stresses. A Raman-based analysis is fast (1 s required for spectral acquisition), portable (measurements can be taken directly in the field), and label-free (no chemicals are needed). These advantages will allow RS to transform agricultural practices, enabling precision agriculture in the near future.

Flood resilience loci SUBMERGENCE 1 and ANAEROBIC GERMINATION 1 interact in seedlings established underwater

Crops with resilience to multiple climatic stresses are essential for increased yield stability. Here, we evaluate the interaction between two loci associated with flooding survival in rice (Oryza sativa L.). ANAEROBIC GERMINATION 1 (AG1), encoding trehalose 6-phosphate phosphatase 7 (TPP7), promotes mobilization of endosperm reserves to enhance the elongation of a hollow coleoptile in seeds that are seeded directly into shallow paddies. SUBMERGENCE 1 (SUB1), encoding the ethylene-responsive transcription factor SUB1A-1, confers tolerance to complete submergence by dampening carbohydrate catabolism, to enhance recovery upon desubmergence. Interactions between AG1/TPP7 and SUB1/SUB1A-1 were investigated under three flooding scenarios using four near-isogenic lines by surveying growth and survival. Pyramiding of the two loci does not negatively affect anaerobic germination or vegetative-stage submergence tolerance. However, the pyramided AG1 SUB1 genotype displays reduced survival when seeds are planted underwater and maintained under submergence for 16 d. To better understand the roles of TPP7 and SUB1A-1 and their interaction, temporal changes in carbohydrates and shoot transcriptomes were monitored in the four genotypes varying at the two loci at four developmental timeponts, from day 2 after seeding through day 14 of complete submergence. TPP7 enhances early coleoptile elongation, whereas SUB1A-1 promotes precocious photoautotrophy and then restricts underwater elongation. By contrast, pyramiding of the AG1 and SUB1 slows elongation growth, the transition to photoautotrophy, and survival. mRNA-sequencing highlights time-dependent and genotype-specific regulation of mRNAs associated with DNA repair, cell cycle, chromatin modification, plastid biogenesis, carbohydrate catabolism and transport, elongation growth, and other processes. These results suggest that interactions between AG1/TPP7 and SUB1/SUB1A-1 could impact seedling establishment if paddy depth is not effectively managed after direct seeding.

Responses of AG1 and AG2 QTL introgression lines and seed pre-treatment on growth and physiological processes during anaerobic germination of rice under flooding

Rice seeds germinating in flooded soils encounter hypoxia or even anoxia leading to poor seed germination and crop establishment. Introgression of AG1 and AG2 QTLs associated with tolerance of flooding during germination, together with seed pre-treatment via hydro-priming or presoaking can enhance germination and seedling growth in anaerobic soils. This study assessed the performance of elite lines incorporating AG1, AG2 and their combination when directly seeded in flooded soils using dry seeds. The QTLs were in the background of two popular varieties PSB Rc82 and Ciherang-Sub1, evaluated along with the donors Kho Hlan On (AG1) and Ma-Zhan Red (AG2) and recipient parents PSB Rc82 and Ciherang-Sub1. In one set of experiments conducted in the greenhouse, seedling emergence, growth, and carbohydrate mobilization from seeds were assessed. Metabolites associated with reactive oxygen species (ROS) scavenging including malondialdehyde (MDA) as a measure of lipid peroxidation, ascorbate, total phenolic concentration (TPC), and activities of ROS scavenging enzymes were quantified in seeds germinating under control (saturated) and flooded (10 cm) soils. In another set of experiments conducted in a natural field with 3-5 cm flooding depths, control and pretreated seeds of Ciherang-Sub1 introgression lines and checks were used. Flooding reduced seedling emergence of all genotypes, though emergence of AG1 + AG2 introgression lines was greater than the other AG lines. Soluble sugars increased, while starch concentration decreased gradually under flooding especially in the tolerant checks and in AG1 + AG2 introgression lines. Less lipid peroxidation and higher α-amylase activity, higher ascorbate (RAsA) and TPC were observed in the tolerant checks and in the AG1 + AG2 introgression lines. Lipid peroxidation correlated negatively with ascorbate, TPC, and with ROS scavengers. Seed hydro-priming or pre-soaking increased emergence by 7-10% over that of dry seeds. Introgression of AG2 and AG1 + AG2 QTLs with seed pretreatment showed 101-153% higher emergence over dry seeds of intolerant genotypes in the field. Lines carrying AG1 + AG2 QTLs showed higher α-amylase activity, leading to rapid starch degradation and increase in soluble sugars, ascorbate, and TPC, together leading to higher germination and seedling growth in flooded soils. Seed hydro-priming or pre-soaking for 24 h also improved traits associated with flooding tolerance. Combining tolerance with seed management could therefore, improve crop establishment in flooded soils and encourage large-scale adoption of direct seeded rice system.

Allelic sequence variation in the Sub1A, Sub1B and Sub1C genes among diverse rice cultivars and its association with submergence tolerance

Erratic rainfall leading to flash flooding causes huge yield losses in lowland rice. The traditional varieties and landraces of rice possess variable levels of tolerance to submergence stress, but gene discovery and utilization of these resources has been limited to the Sub1A-1 allele from variety FR13A. Therefore, we analysed the allelic sequence variation in three Sub1 genes in a panel of 179 rice genotypes and its association with submergence tolerance. Population structure and diversity analysis based on a 36-plex genome wide genic-SNP assay grouped these genotypes into two major categories representing Indica and Japonica cultivar groups with further sub-groupings into Indica, Aus, Deepwater and Aromatic-Japonica cultivars. Targetted re-sequencing of the Sub1A, Sub1B and Sub1C genes identfied 7, 7 and 38 SNPs making 8, 9 and 67 SNP haplotypes, respectively. Haplotype networks and phylogenic analysis revealed evolution of Sub1B and Sub1A genes by tandem duplication and divergence of the ancestral Sub1C gene in that order. The alleles of Sub1 genes in tolerant reference variety FR13A seem to have evolved most recently. However, no consistent association could be found between the Sub1 allelic variation and submergence tolerance probably due to low minor allele frequencies and presence of exceptions to the known Sub1A-1 association in the genotype panel. We identified 18 cultivars with non-Sub1A-1 source of submergence tolerance which after further mapping and validation in bi-parental populations will be useful for development of superior flood tolerant rice cultivars.

An improved 7K SNP array, the C7AIR, provides a wealth of validated SNP markers for rice breeding and genetics studies

Single nucleotide polymorphisms (SNPs) are highly abundant, amendable to high-throughput genotyping, and useful for a number of breeding and genetics applications in crops. SNP frequencies vary depending on the species and populations under study, and therefore target SNPs need to be carefully selected to be informative for each application. While multiple SNP genotyping systems are available for rice (Oryza sativa L. and its relatives), they vary in their informativeness, cost, marker density, speed, flexibility, and data quality. In this study, we report the development and performance of the Cornell-IR LD Rice Array (C7AIR), a second-generation SNP array containing 7,098 markers that improves upon the previously released C6AIR. The C7AIR is designed to detect genome-wide polymorphisms within and between subpopulations of O. sativa, as well as O. glaberrima, O. rufipogon and O. nivara. The C7AIR combines top-performing SNPs from several previous rice arrays, including 4,007 SNPs from the C6AIR, 2,056 SNPs from the High Density Rice Array (HDRA), 910 SNPs from the 384-SNP GoldenGate sets, 189 SNPs from the 44K array selected to add information content for elite U.S. tropical japonica rice varieties, and 8 trait-specific SNPs. To demonstrate its utility, we carried out a genome-wide association analysis for plant height, employing the C7AIR across a diversity panel of 189 rice accessions and identified 20 QTLs contributing to plant height. The C7AIR SNP chip has so far been used for genotyping >10,000 rice samples. It successfully differentiates the five subpopulations of Oryza sativa, identifies introgressions from wild and exotic relatives, and is useful for quantitative trait loci (QTL) and association mapping in diverse materials. Moreover, data from the C7AIR provides valuable information that can be used to select informative and reliable SNP markers for conversion to lower-cost genotyping platforms for genomic selection and other downstream applications in breeding.

Genome-Wide Association Mapping to Identify Genetic Loci for Cold Tolerance and Cold Recovery During Germination in Rice

Low temperature significantly affects rice growth and yield. Temperatures lower than 15°C are generally detrimental for germination and uniform seedling stand. To investigate the genetic architecture underlying cold tolerance during germination in rice, we conducted a genome-wide association study using a novel diversity panel of 257 rice accessions from around the world and the 7K SNP marker array. Phenotyping was conducted in controlled growth chambers under dark conditions at 13°C. The rice accessions were measured for low-temperature germinability, germination index, coleoptile length under cold stress, plumule length at 4-day recovery, and plumule length recovery rate. A total of 51 QTLs were identified at p < 0.001 and 17 QTLs were identified using an FDR < 0.05 across the different chilling indices with the whole panel of accessions. At the threshold of p < 0.001, a total of 20 QTLs were identified in the subset of japonica accessions, while 9 QTLs were identified in the subset of indica accessions. Considering the recurring SNPs and linked SNPs across different chilling indices, we identified 31 distinct QTL regions in the whole panel, 13 QTL regions in the japonica subset, and 7 distinct QTL regions in the indica subset. Among these QTL regions, three regions were common between the whole panel and japonica, three regions were common between the whole panel and indica, and one region was common between indica and japonica. A subset of QTL regions was potentially colocalized with previously identified genes and QTLs, including 10 from the japonica subset, 4 from the indica subset, and 6 from the whole panel. On the other hand, a total of 21 potentially novel QTL regions from the whole panel, 10 from the japonica subset, and 1 from the indica subset were identified. The results of our study provide useful information on the genetic architecture underlying cold tolerance during germination in rice, which in turn can be used for further molecular study and crop improvement for low-temperature stressed environments.

Trait-based mapping to identify the genetic factors underlying anaerobic germination of rice: Phenotyping, GXE, and QTL mapping

BACKGROUND

Anaerobic germination is one of the most important traits for rice under direct-seeded conditions. The trait reduces risk of crop failure due to waterlogged conditions after seeding and allows water to be used as a means of weed control. The identification of QTLs and causal genes for anaerobic germination will facilitate breeding for improved direct-seeded rice varieties. In this study, we explored a BC₁F_2:3 population developed from a cross between BJ1, an indica landrace, and NSIC Rc222, a high-yielding recurrent parent. The population was phenotyped under different screening methods (anaerobic screenhouse, anaerobic tray, and aerobic screenhouse) to establish the relationship among the methods and to identify the most suitable screening method, followed by bulk segregant analysis (BSA) to identify large-effect QTLs.

RESULTS

The study showed high heritability for survival (SUR) under all three phenotyping conditions. Although high correlation was observed within screening environments between survival at 14 and 21 days after seeding, the correlation across environments was low. Germination under aerobic and anaerobic conditions showed very low correlation, indicating the independence of their genetic control. The results were further confirmed through AMMI analysis. Four significant markers with an effect on anaerobic germination were identified through BSA. CIM analysis revealed qAG1-2, qAG6-2, qAG7-4, and qAG10-1 having significant effects on the trait. qAG6-2 and qAG10-1 were consistent across screening conditions and seedling age while qAG1-2 and qAG7-4 were specific to screening methods. All QTLs showed an effect when survival across all screening methods was analyzed. Together, the QTLs explained 39 to 55% of the phenotypic variation for survival under anaerobic conditions. No QTL effects were observed under aerobic conditions.

CONCLUSIONS

The study helped us understand the effect of phenotyping method on anaerobic germination, which will lead to better phenotyping for this trait in future studies. The QTLs identified through this study will allow the improvement of breeding lines for the trait through marker-assisted selection or through forward breeding approaches such as genomic selection. The high frequency of the BJ1 allele of these QTLs will enhance the robustness of germination under anaerobic conditions in inbred and hybrid rice varieties.

Increasing flooding tolerance in rice: combining tolerance of submergence and of stagnant flooding

BACKGROUND AND AIMS

Rice ecosystems in the tropical coastal areas are subject to two types of flooding stress: transient complete submergence and long-term water stagnation (stagnant flooding). Here, we aimed to dissect the mechanisms for stagnant flooding tolerance of rice genotypes carrying SUB1, a quantitative trait locus for submergence tolerance.

METHODS

We screened 80 elite genotypes under stagnant flooding stress in the lowland rice fields in the wet and dry seasons, and examined the tolerance mechanisms of promising genotypes for the two following seasons.

KEY RESULTS

Yield reduction under stagnant flooding averaged 48 % in the dry season and 89 % in the wet season. Elite genotypes carrying SUB1 showed 49 % lower yield than those without SUB1 under stagnant flooding, with no differences under shallow water conditions. However, we identified a few high-yielding Sub1 genotypes that were as tolerant of stagnant flooding as a reference genotype that lacked SUB1. These genotypes had intermediate stature with more shoot elongation in response to rising water than a moderately tolerant Sub1 reference variety, resulting in greater canopy expansion and higher yield. It was important to increase lodging resistance, since plant height >140 cm increased lodging under stagnant flooding. The culm diameter was closely associated with culm strength; reduced aerenchyma formation and increased lignin accumulation in the culm should increase lodging resistance.

CONCLUSIONS

The study demonstrated a successful combination of submergence and stagnant flooding tolerance in a rice breeding programme, and identified elite Sub1 genotypes that also tolerate stagnant flooding. Our results will support genetic improvement of Sub1 varieties for stagnant flooding tolerance.

Deciphering Genetics Underlying Stable Anaerobic Germination in Rice: Phenotyping, QTL Identification, and Interaction Analysis

Anaerobic germination (AG) is an important trait for direct-seeded rice (DSR) to be successful. Rice usually has low germination under anaerobic conditions, which leads to a poor crop stand in DSR when rain occurs after seeding. The ability of rice to germinate under water reduces the risk of poor crop stand. Further, this allows the use of water as a method of weed control. The identification of the genetic factors leading to high anaerobic germination is required to develop improved DSR varieties. In the present study, two BC₁F_2:3 mapping families involving a common parent with anaerobic germination potential, Kalarata, an indica landrace, and two recurrent parents, NSIC Rc222 and NSIC Rc238, were used. Phenotyping was done under two environmental conditions and genotyping was carried out through the KASP SNP genotyping platform. A total of 185 and 189 individuals genotyped with 170 and 179 polymorphic SNPs were used for QTL analysis for the two populations, Kalarata/NSIC Rc238 and Kalarata/NSIC Rc222, respectively. A total of five QTLs on chromosomes 3, 5, 6, 7, and 8 for survival (SUR) and four QTLs on chromosomes 1, 3 (two locations), and 7 for the trait seedling height (SH) across the populations and over the screening conditions were identified. Except for the QTLs on chromosomes 5 and 8, the parent with AG potential, Kalarata, contributed all the other QTLs. Among the five QTLs for SUR, the second-largest QTL (qSUR6-1) was novel for AG potential in rice, showing a stable expression in terms of genetic background and screening conditions explaining 11.96% to 16.01% of the phenotypic variation. The QTL for SH (qSH1-1) was also novel. Considering different genetic backgrounds and different screening conditions, the QTLs identified for the trait SUR explained phenotypic variation in the range of 57.60% to 73.09% while that for the trait SH ranged from 13.53% to 34.30%.

QTL mapping for tolerance to anaerobic germination in rice from IR64 and the aus landrace Kharsu 80A

Direct seeding of rice often results in poor crop establishment due to unlevelled fields, unpredicted heavy rains after sowing, and weed and pest invasion. Thus, it is important to develop varieties able to tolerate flooding during germination, also known as anaerobic germination (AG), to address these constraints. A study was conducted to identify QTLs associated with AG tolerance from an IR64/Kharsu 80A F2:3 mapping population using 190 lines phenotyped for seedling survival under the stress. Genotyping was performed using a genomewide 384-plex Indica/Indica SNP set. Four QTLs derived from Kharsu 80A providing increased tolerance to anaerobic germination were identified: three on chromosome 7 (qAG7.1, qAG7.2 and qAG7.3) and one on chromosome 3 (qAG3), with LOD values ranging from 5.7 to 7.7, and phenotypic variance explained (R2) from 8.1% to 12.6%. The QTLs identified in this study can be further investigated to better understand the genetic bases of AG tolerance in rice, and used for marker-assisted selection to develop more robust direct-seeded rice varieties.

Optimizing Sowing and Flooding Depth for Anaerobic Germination-Tolerant Genotypes to Enhance Crop Establishment, Early Growth, and Weed Management in Dry-Seeded Rice (Oryza sativa L.)

Poor crop establishment, high weed infestation, and consequent yield loss are major concerns for dry-seeded rice (DSR). Flooding after seeding helps in managing weeds but reduces seed germination and crop stand. Anaerobic germination (AG)-tolerant rice genotypes could overcome these problems in DSR. Screenhouse experiments were established to evaluate the effect of seed sowing depth (SD) (0.5 cm, 1 cm, and 2 cm) and flooding depth (FD) (saturated, 2 cm, and 5 cm) on crop establishment, early growth, and weed competitiveness in DSR using AG-tolerant genotypes (Khao Hlan On, Ma-Zhan Red, IR64+AG1, and IR64). Echinochloa crus-galli, Ludwigia hyssopifolia, and Cyperus difformis were used in the weedy treatment. Rice plants reached maximum emergence 9-13 days later under flooding compared with saturated conditions. Crop emergence decreased by 12-22% at 0.5 and 1 cm SD and by 48-60% at 2 cm SD, when combined with 2 or 5 cm FD compared with saturated conditions. The 2 cm SD reduced seedling emergence by 23-42% in Khao Hlan On and Ma-Zhan Red, by 62-70% in IR64+AG1, and by 90-92% in IR64 under flooding. Initial growth in rice plant height was slow under flooding but increased progressively after the seedlings emerged from water and the final height was not affected by FD. Leaf area, total shoot biomass, tiller density, and leaf number per pot of rice were higher at 1 cm SD (P< 0.05), but decreased drastically at 2 cm SD under flooding. The emergence of E. crus-galli and L. hyssopifolia decreased by 53-65% and 89-95%, respectively, but increased by 49-68% in C. difformis under 2 and 5 cm FD, respectively, compared with saturated conditions. The shoot biomass of the weeds followed the same trend. Khao Hlan On showed the highest weed-competitive ability under all FD while the biomass of IR64+AG1 and IR64 decreased by 10-14% due to weed competition under 2 cm FD. The 1 cm SD showed better growth for all genotypes under different FD. The 2 cm FD is sufficient to have a significant control of problematic weed species. The tolerance of AG of rice genotypes should be further enhanced to increase their weed-competitive ability.

Identification of QTLs for yield and agronomic traits in rice under stagnant flooding conditions

BACKGROUND

Stagnant flooding, where water of 25-50 cm remains until harvest time, is a major problem in rainfed lowland areas. Most of the Sub1 varieties, which can withstand around 2 weeks of complete submergence, perform poorly in these conditions. Hence, varieties tolerant of stagnant flooding are essential.

RESULTS

This paper presents the first study to map QTLs associated with tolerance to stagnant flooding, along with a parallel study under normal irrigation, using an F₇ mapping population consisting of 148 RILs derived from a cross of Ciherang-Sub1 and the stagnant-flooding tolerant line IR10F365. Phenotypic data was collected for 15 key traits under both environments. Additionally, survival rate was measured under stress conditions. Genotyping was performed using the Illumina Infinium genotyping platform with a 6 K SNP chip, resulting in 469 polymorphic SNPs. Under stress and irrigated conditions, 38 and 46 QTLs were identified, respectively. Clusters of QTLs were detected in both stress and normal conditions, especially on chromosomes 3 and 5.

CONCLUSIONS

Unique and common QTLs were identified and their physiological consequences are discussed. These beneficial QTLs can be used as targets for molecular breeding and can be further investigated to understand the underlying molecular mechanisms involved in stagnant flooding tolerance in rice.

Genetics, Physiological Mechanisms and Breeding of Flood-Tolerant Rice (Oryza sativa L.)

Flooding of rice fields is a serious problem in the river basins of South and South-East Asia where about 15 Mha of lowland rice cultivation is regularly affected. Flooding creates hypoxic conditions resulting in poor germination and seedling establishment. Flash flooding, where rice plants are completely submerged for 10-15 d during their vegetative stage, causes huge losses. Water stagnation for weeks to months also leads to substantial yield losses when large parts of rice aerial tissues are inundated. The low-yielding traditional varieties and landraces of rice adapted to these flooding conditions have been replaced by flood-sensitive high-yielding rice varieties. The 'FR13A' rice variety and the Submergence 1A (SUB1A) gene were identified for flash flooding and subsequently introgressed to high-yielding rice varieties. The challenge is to find superior alleles of the SUB1A gene, or even new genes that may confer greater tolerance to submergence. Similarly, genes have been identified in tolerant landraces of rice for their ability to survive by rapid stem elongation (SNORKEL1 and SNORKEL2) during deep-water flooding, and for anaerobic germination ability (TPP7). Research on rice genotypes and novel genes that are tolerant to prolonged water stagnation is in progress. These studies will greatly assist in devising more efficient and precise molecular breeding strategies for developing climate-resilient high-yielding rice varieties for flood-prone regions. Here we review the state of our knowledge of flooding tolerance in rice and its application in varietal improvement.

EGRINs (Environmental Gene Regulatory Influence Networks) in Rice That Function in the Response to Water Deficit, High Temperature, and Agricultural Environments

Environmental gene regulatory influence networks (EGRINs) coordinate the timing and rate of gene expression in response to environmental signals. EGRINs encompass many layers of regulation, which culminate in changes in accumulated transcript levels. Here, we inferred EGRINs for the response of five tropical Asian rice (Oryza sativa) cultivars to high temperatures, water deficit, and agricultural field conditions by systematically integrating time-series transcriptome data, patterns of nucleosome-free chromatin, and the occurrence of known cis-regulatory elements. First, we identified 5447 putative target genes for 445 transcription factors (TFs) by connecting TFs with genes harboring known cis-regulatory motifs in nucleosome-free regions proximal to their transcriptional start sites. We then used network component analysis to estimate the regulatory activity for each TF based on the expression of its putative target genes. Finally, we inferred an EGRIN using the estimated transcription factor activity (TFA) as the regulator. The EGRINs include regulatory interactions between 4052 target genes regulated by 113 TFs. We resolved distinct regulatory roles for members of the heat shock factor family, including a putative regulatory connection between abiotic stress and the circadian clock. TFA estimation using network component analysis is an effective way of incorporating multiple genome-scale measurements into network inference.

Multiple abiotic stimuli are integrated in the regulation of rice gene expression under field conditions

Plants rely on transcriptional dynamics to respond to multiple climatic fluctuations and contexts in nature. We analyzed the genome-wide gene expression patterns of rice (Oryza sativa) growing in rainfed and irrigated fields during two distinct tropical seasons and determined simple linear models that relate transcriptomic variation to climatic fluctuations. These models combine multiple environmental parameters to account for patterns of expression in the field of co-expressed gene clusters. We examined the similarities of our environmental models between tropical and temperate field conditions, using previously published data. We found that field type and macroclimate had broad impacts on transcriptional responses to environmental fluctuations, especially for genes involved in photosynthesis and development. Nevertheless, variation in solar radiation and temperature at the timescale of hours had reproducible effects across environmental contexts. These results provide a basis for broad-based predictive modeling of plant gene expression in the field.

DOI:http://dx.doi.org/10.7554/eLife.08411.001

Plants need to be able to sense and respond to changes in temperature, light levels and other aspects of their environment. One way in which plants can rapidly respond to these changes is to modify how genes involved in growth and other processes are expressed. Therefore, understanding how this happens may help us to improve the ability of crops to grow when exposed to drought or other extreme environmental conditions.

Most previous studies into the effect of the environment on plant gene expression have been carried out under controlled conditions in a laboratory. These findings cannot reflect the full range of gene expression patterns that occur in the natural environment, where multiple factors (e.g. sunlight, water, nutrients) may vary at the same time. Therefore, it is important to also analyze the effect of fluctuations in multiple environmental factors in more complex field experiments.

Plessis et al. developed mathematical models to analyze the gene expression patterns of rice plants grown in the tropical environment of the Philippines using two different farming practices. One field of rice was flooded and constantly supplied with fresh water (referred to as the irrigated field), while the other field was dry and only received water from rainfall (the rainfed field). The experiments show that temperature and levels of sunlight (including UV radiation) have a strong impact on gene expression in the rice plants. Short-term variations in temperature and sunlight levels also have the most consistent effect across the different fields and seasons tested. However, for many genes, the plants grown in the irrigated field responded to the changes in environmental conditions in a different way to the plants grown in the rainfed field.

Further analysis identified groups of genes whose expression combined responses to several environmental factors at the same time. For example, certain genes that responded to increases in sunlight in the absence of drought responded to both sunlight levels and the shortage of water when a drought occurred. The next step is to test more types of environments and climates to be able to predict gene expression responses under future climatic conditions.

DOI:http://dx.doi.org/10.7554/eLife.08411.002

A trehalose-6-phosphate phosphatase enhances anaerobic germination tolerance in rice

Global socioeconomic developments create strong incentives for farmers to shift from transplanted to direct-seeded rice (DSR) as a means of intensification and economization(1). Rice production must increase to ensure food security(2) and the bulk of this increase will have to be achieved through intensification of cultivation, because expansion of cultivated areas is reaching sustainable limits(3). Anaerobic germination tolerance, which enables uniform germination and seedling establishment under submergence(4), is a key trait for the development of tropical DSR varieties(5,6). Here, we identify a trehalose-6-phosphate phosphatase gene, OsTPP7, as the genetic determinant in qAG-9-2, a major quantitative trait locus (QTL) for anaerobic germination tolerance(7). OsTPP7 is involved in trehalose-6-phosphate (T6P) metabolism, central to an energy sensor that determines anabolism or catabolism depending on local sucrose availability(8,9). OsTPP7 activity may increase sink strength in proliferating heterotrophic tissues by indicating low sugar availability through increased T6P turnover, thus enhancing starch mobilization to drive growth kinetics of the germinating embryo and elongating coleoptile, which consequently enhances anaerobic germination tolerance.

A trehalose-6-phosphate phosphatase enhances anaerobic germination tolerance in rice

Global socioeconomic developments create strong incentives for farmers to shift from transplanted to direct-seeded rice (DSR) as a means of intensification and economization(1). Rice production must increase to ensure food security(2) and the bulk of this increase will have to be achieved through intensification of cultivation, because expansion of cultivated areas is reaching sustainable limits(3). Anaerobic germination tolerance, which enables uniform germination and seedling establishment under submergence(4), is a key trait for the development of tropical DSR varieties(5,6). Here, we identify a trehalose-6-phosphate phosphatase gene, OsTPP7, as the genetic determinant in qAG-9-2, a major quantitative trait locus (QTL) for anaerobic germination tolerance(7). OsTPP7 is involved in trehalose-6-phosphate (T6P) metabolism, central to an energy sensor that determines anabolism or catabolism depending on local sucrose availability(8,9). OsTPP7 activity may increase sink strength in proliferating heterotrophic tissues by indicating low sugar availability through increased T6P turnover, thus enhancing starch mobilization to drive growth kinetics of the germinating embryo and elongating coleoptile, which consequently enhances anaerobic germination tolerance.

A trehalose-6-phosphate phosphatase enhances anaerobic germination tolerance in rice

Global socioeconomic developments create strong incentives for farmers to shift from transplanted to direct-seeded rice (DSR) as a means of intensification and economization(1). Rice production must increase to ensure food security(2) and the bulk of this increase will have to be achieved through intensification of cultivation, because expansion of cultivated areas is reaching sustainable limits(3). Anaerobic germination tolerance, which enables uniform germination and seedling establishment under submergence(4), is a key trait for the development of tropical DSR varieties(5,6). Here, we identify a trehalose-6-phosphate phosphatase gene, OsTPP7, as the genetic determinant in qAG-9-2, a major quantitative trait locus (QTL) for anaerobic germination tolerance(7). OsTPP7 is involved in trehalose-6-phosphate (T6P) metabolism, central to an energy sensor that determines anabolism or catabolism depending on local sucrose availability(8,9). OsTPP7 activity may increase sink strength in proliferating heterotrophic tissues by indicating low sugar availability through increased T6P turnover, thus enhancing starch mobilization to drive growth kinetics of the germinating embryo and elongating coleoptile, which consequently enhances anaerobic germination tolerance.

Physiological analyses of traits associated with tolerance of long-term partial submergence in rice

Floods are major constraints to crop production worldwide. In low-lying, flood-prone areas of the tropics, longer-term partial submergence (stagnant flooding [SF]) greatly reduces rice yield. This study assesses shoot growth and several physiological mechanisms associated with SF tolerance in rice. Five rice genotypes with contrasting responses to SF were evaluated in field ponds. Following transplanting, floodwater was gradually increased at a rate of ∼2 cm day(-1) to reach a final depth of 50 cm and then maintained until maturity. Although plants were not fully submerged, the yield was reduced by 47 % across genotypes compared with those grown under control conditions (6.1 vs. 3.3 t ha(-1)). This reduction was mainly attributed to the reduction in biomass caused by reduced light interception and leaf growth above the water. Stagnant flooding also reduced panicle number per unit area by 52 % because of reduced tillering. Shoot elongation rate kept pace with rising floodwater and correlated positively with leaf growth and biomass production. Conversely, stem non-structural carbohydrate (NSC) concentration correlated negatively with shoot elongation rate, suggesting that fast-elongating genotypes actively consume NSCs to avoid complete submergence. Moderate shoot elongation rate strongly and positively correlated with grain yield under SF; however, elongation at rates >2.0 cm day(-1) was associated with reduced harvest index due to a smaller panicle size and increased lodging. Tolerant varieties were found to be either inherently tall or elongate moderately with rising floodwater. Our studies suggest that to improve tolerance of SF an appropriate phenotype should combine both of these traits. Fine-tuning for optimum shoot elongation with rising floodwater is, therefore, a priority for future work.

QTL mapping and confirmation for tolerance of anaerobic conditions during germination derived from the rice landrace Ma-Zhan Red

Wide adoption of direct-seeded rice practices has been hindered by poorly leveled fields, heavy rainfall and poor drainage, which cause accumulation of water in the fields shortly after sowing, leading to poor crop establishment. This is due to the inability of most rice varieties to germinate and reach the water surface under complete submergence. Hence, tolerance of anaerobic conditions during germination is an essential trait for direct-seeded rice cultivation in both rainfed and irrigated ecosystems. A QTL study was conducted to unravel the genetic basis of tolerance of anaerobic conditions during germination using a population derived from a cross between IR42, a susceptible variety, and Ma-Zhan Red, a tolerant landrace from China. Phenotypic data was collected based on the survival rates of the seedlings at 21 days after sowing of dry seeds under 10 cm of water. QTL analysis of the mapping population consisting of 175 F2:3 families genotyped with 118 SSR markers identified six significant QTLs on chromosomes 2, 5, 6, and 7, and in all cases the tolerant alleles were contributed by Ma-Zhan Red. The largest QTL on chromosome 7, having a LOD score of 14.5 and an R (2) of 31.7 %, was confirmed using a BC2F3 population. The QTLs detected in this study provide promising targets for further genetic characterization and for use in marker-assisted selection to rapidly develop varieties with improved tolerance to anaerobic condition during germination. Ultimately, this trait can be combined with other abiotic stress tolerance QTLs to provide resilient varieties for direct-seeded systems.

Comparison of phenotypic versus marker-assisted background selection for the SUB1 QTL during backcrossing in rice

Marker assisted backcrossing has been used effectively to transfer the submergence tolerance gene SUB1 into popular rice varieties, but the approach can be costly. The selection strategy comprising foreground marker and phenotypic selection was investigated as an alternative. The non-significant correlation coefficients between ranking of phenotypic selection and ranking of background marker selection in BC(2)F(1), BC(3)F(1) and BC(3)F(2) generations indicated inefficiency of phenotypic selection compared to marker-assisted background selection with respect to recovery of the recipient genome. In addition, the introgression size of the chromosome fragment containing SUB1 was approximately 17 Mb, showing the effects of linkage drag. The significant correlation coefficient between rankings of phenotypic selection with the percentage of recipient alleles in the BC(1)F(1) generation suggested that background selection could be avoided in this generation to minimize the genotyping cost. The phenotypically selected best plant of the BC(3)F(1) generation was selfed and backcross recombinant lines were selected in the resulting BC(3)F(4) generation. The selection strategy could be appropriate for the introgression of SUB1 QTL in countries that lack access to high-throughput genotyping facilities.

Identifying novel QTLs for submergence tolerance in rice cultivars IR72 and Madabaru

Short-term submergence is a recurring problem in many rice production areas. The SUB1 gene, derived from the tolerant variety FR13A, has been transferred to a number of widely grown varieties, allowing them to withstand complete submergence for up to 2 weeks. However, in areas where longer-term submergence occurs, improved varieties having higher tolerance levels are needed. To search for novel quantitative trait loci (QTLs) from other donors, an F(2:3) population between IR72 and Madabaru, both moderately tolerant varieties, was investigated. After a repeated phenotyping of 466 families under submergence stress, a subset of 80 families selected from the two extreme phenotypic tails was used for the QTL analysis. Phenotypic data showed transgressive segregation, with several families having an even higher survival rate than the FR13A-derived tolerant check (IR40931). Four QTLs were identified on chromosomes 1, 2, 9, and 12; the largest QTL on chromosome 1 had a LOD score of 11.2 and R (2) of 52.3%. A QTL mapping to the SUB1 region on chromosome 9, with a LOD score of 3.6 and R (2) of 18.6%, had the tolerant allele from Madabaru, while the other three QTLs had tolerant alleles from IR72. The identification of three non-SUB1 QTLs from IR72 suggests that an alternative pathway may be present in this variety that is independent of the ethylene-dependent pathway mediated by the SUB1A gene. These novel QTLs can be combined with SUB1 using marker assisted backcrossing in an effort to enhance the level of submergence tolerance for flood-prone areas.

Molecular marker survey and expression analyses of the rice submergence-tolerance gene SUB1A

The major rice quantitative-trait locus Submergence 1 (Sub1) confers tolerance of submergence for about 2 weeks. To identify novel sources of tolerance, we have conducted a germplasm survey with allele-specific markers targeting SUB1A and SUB1C, two of the three transcription-factor genes within the Sub1 locus. The objective was to identify tolerant genotypes without the SUB1A gene or with the intolerant SUB1A-2 allele. The survey revealed that all tolerant genotypes possessed the tolerant Sub1 haplotype (SUB1A-1/SUB1C-1), whereas all accessions without the SUB1A gene were intolerant. Only the variety James Wee with the SUB1A-2 allele was moderately tolerant. However, some intolerant genotypes with the SUB1A-1 allele were identified and RT-PCR analyses were conducted to compare gene expression in tolerant and intolerant accessions. Initial analyses of leaf samples failed to reveal a clear association of SUB1A transcript abundance and tolerance. Temporal and spatial gene expression analyses subsequently showed that SUB1A expression in nodes and internodes associated best with tolerance across representative genotypes. In James Wee, transcript abundance was high in all tissues, suggesting that some level of tolerance might be conferred by high expression of the SUB1A-2 allele. To further assess tissue-specific expression, we have expressed the GUS reporter gene under the control of the SUB1A-1 promoter. The data revealed highly specific GUS expression at the base of the leaf sheath and in the leaf collar region. Specific expression in the growing part of rice leaves is well in agreement with the role of SUB1A in suppressing leaf elongation under submergence.

Development of submergence-tolerant rice cultivars: the Sub1 locus and beyond

BACKGROUND AND AIMS

Submergence is a recurring problem in the rice-producing rainfed lowlands of south and south-east Asia. Developing rice cultivars with tolerance of submergence and with agronomic and quality traits acceptable to farmers is a feasible approach to address this problem. The objectives of this study were to (a) develop mega varieties with Sub1 introgression that are submergence tolerant, (b) assess the performance of Sub1 in different genetic backgrounds, (c) determine the roles of the Sub1A and Sub1C genes in conferring tolerance, and (d) assess the level of tolerance in F(1) hybrids heterozygous for the Sub1A-1-tolerant allele.

METHODS

Tolerant varieties were developed by marker-assisted backcrossing through two or three backcrosses, and their performance was evaluated to determine the effect of Sub1 in different genetic backgrounds. The roles of Sub1A and Sub1C in conferring the tolerant phenotype were further investigated using recombinants identified within the Sub1 gene cluster based on survival and gene expression data.

KEY RESULTS

All mega varieties with Sub1 introgression had a significantly higher survival rate than the original parents. An intolerant Sub1C allele combined with the tolerant Sub1A-1 allele did not significantly reduce the level of tolerance, and the Sub1C-1 expression appeared to be independent of the Sub1A allele; however, even when Sub1C-1 expression is completely turned off in the presence of Sub1A-2, plants remained intolerant. Survival rates and Sub1A expression were significantly lower in heterozygotes compared with the homozygous tolerant parent.

CONCLUSIONS

Sub1 provided a substantial enhancement in the level of tolerance of all the sensitive mega varieties. Sub1A is confirmed as the primary contributor to tolerance, while Sub1C alleles do not seem important. Lack of dominance of Sub1 suggests that the Sub1A-1 allele should be carried by both parents for developing tolerant rice hybrids.

A marker-assisted backcross approach for developing submergence-tolerant rice cultivars

Submergence stress regularly affects 15 million hectares or more of rainfed lowland rice areas in South and Southeast Asia. A major QTL on chromosome 9, Sub1, has provided the opportunity to apply marker assisted backcrossing (MAB) to develop submergence tolerant versions of rice cultivars that are widely grown in the region. In the present study, molecular markers that were tightly linked with Sub1, flanking Sub1, and unlinked to Sub1 were used to apply foreground, recombinant, and background selection, respectively, in backcrosses between a submergence-tolerant donor and the widely grown recurrent parent Swarna. By the BC(2)F(2) generation a submergence tolerant plant was identified that possessed Swarna type simple sequence repeat (SSR) alleles on all fragments analyzed except the tip segment of rice chromosome 9 that possessed the Sub1 locus. A BC(3)F(2) double recombinant plant was identified that was homozygous for all Swarna type alleles except for an approximately 2.3-3.4 Mb region surrounding the Sub1 locus. The results showed that the mega variety Swarna could be efficiently converted to a submergence tolerant variety in three backcross generations, involving a time of two to three years. Polymorphic markers for foreground and recombinant selection were identified for four other mega varieties to develop a wider range of submergence tolerant varieties to meet the needs of farmers in the flood-prone regions. This approach demonstrates the effective use of marker assisted selection for a major QTL in a molecular breeding program.

Genetic diversity analysis of traditional and improved Indonesian rice (Oryza sativa L.) germplasm using microsatellite markers

The archipelago of Indonesia has a long history of rice production across a broad range of rice-growing environments resulting in a diverse array of local Indonesian rice varieties. Although some have been incorporated into modern breeding programs, the vast majority of these landraces remain untapped. To better understand this rich source of genetic diversity we have characterized 330 rice accessions, including 246 Indonesian landraces and 63 Indonesian improved cultivars, using 30 fluorescently-labeled microsatellite markers. The landraces were selected across 21 provinces and include representatives of the classical subpopulations of cere, bulu, and gundil rices. A total of 394 alleles were detected at the 30 simple sequence repeat loci, with an average number of 13 alleles per locus across all accessions, and an average polymorphism information content value of 0.66. Genetic diversity analysis characterized the Indonesian landraces as 68% indica and 32% tropical japonica, with an indica gene diversity of 0.53 and a tropical japonica gene diversity of 0.56, and a Fst of 0.38 between the two groups. All of the improved varieties sampled were indica, and had an average gene diversity of 0.46. A set of high quality Indonesian varieties, including Rojolele, formed a separate cluster within the tropical japonicas. This germplasm presents a valuable source of diversity for future breeding and association mapping efforts.

Substitution mapping of dth1.1, a flowering-time quantitative trait locus (QTL) associated with transgressive variation in rice, reveals multiple sub-QTL

A quantitative trait locus (QTL), dth1.1, was associated with transgressive variation for days to heading in an advanced backcross population derived from the Oryza sativa variety Jefferson and an accession of the wild rice relative Oryza rufipogon. A series of near-isogenic lines (NILs) containing different O. rufipogon introgressions across the target region were constructed to dissect dth1.1 using substitution mapping. In contrast to the late-flowering O. rufipogon parent, O. rufipogon alleles in the substitution lines caused early flowering under both short- and long-day lengths and provided evidence for at least two distinct sub-QTL: dth1.1a and dth1.1b. Potential candidate genes underlying these sub-QTL include genes with sequence similarity to Arabidopsis GI, FT, SOC1, and EMF1, and Pharbitis nil PNZIP. Evidence from families with nontarget O. rufipogon introgressions in combination with dth1.1 alleles also detected an early flowering QTL on chromosome 4 and a late-flowering QTL on chromosome 6 and provided evidence for additional sub-QTL in the dth1.1 region. The availability of a series of near-isogenic lines with alleles introgressed from a wild relative of rice provides an opportunity to better understand the molecular basis of transgressive variation in a quantitative trait.

Identification of quantitative trait loci for yield and yield components in an advanced backcross population derived from the Oryza sativa variety IR64 and the wild relative O. rufipogon

A BC(2)F(2) population developed from an interspecific cross between Oryza sativa (cv IR64) and O. rufipogon (IRGC 105491) was used in an advanced backcross QTL analysis to identify and introduce agronomically useful genes from this wild relative into the cultivated gene pool. The objectives of this study were: (1) to identify putative yield and yield component QTLs that can be useful to improve the elite cultivar IR64; (2) to compare the QTLs within this study with previously reported QTLs in rice as the basis for identifying QTLs that are stable across different environments and genetic backgrounds; and (3) to compare the identified QTLs with previously reported QTLs from maize to examine the degree of QTL conservation across the grass family. Two hundred eighty-five families were evaluated in two field environments in Indonesia, with two replications each, for 12 agronomic traits. A total of 165 markers consisting of 131 SSRs and 34 RFLPs were used to construct the genetic linkage map. By employing interval mapping and composite interval mapping, 42 QTLs were identified. Despite its inferior performance, 33% of the QTL alleles originating from O. rufipogon had a beneficial effect for yield and yield components in the IR64 background. Twenty-two QTLs (53.4%) were located in similar regions as previously reported rice QTLs, suggesting the existence of stable QTLs across genetic backgrounds and environments. Twenty QTLs (47.6%) were exclusively detected in this study, uncovering potentially novel alleles from the wild, some of which might improve the performance of the tropical indica variety IR64. Additionally, several QTLs for plant height, grain weight, and flowering time detected in this study corresponded to homeologous regions in maize containing previously detected maize QTLs for these traits.

Identification of quantitative trait loci for grain quality in an advanced backcross population derived from the Oryza sativa variety IR64 and the wild relative O. rufipogon

The objective of this study was to identify quantitative trait loci (QTLs) associated with grain quality in rice. Two hundred eighty-five BC(2)F(2 )families developed from an interspecific cross between cv IR64 and Oryza rufipogon (IRGC 105491) were evaluated for 14 seed quality traits. A total of 165 markers consisting of 131 single sequence repeats and 34 restriction fragment length polymorphism markers were used to create a genetic linkage map spanning the 12 rice chromosomes. Twenty-three independent QTLs were identified using single point analysis, interval mapping, and composite interval mapping. These loci consisted of one QTL for filled rough/total rough rice ratio, two for grain density, one for percentage of de-husked rice grains, two for percentage of green rice grains, three for percentage of damaged-yellow rice grains, two for percentage of red rice grains, one for milled rice recovery, three for head rice recovery, four for broken rice grains, two for crushed rice grains, one for amylose content, and one for gel consistency. For most of the QTLs identified in this study, the O. rufipogon-derived allele contributed an undesirable effect. For amylose content and gel consistency, the O. rufipogon allele may be useful in an IR64 background, depending on the cultural preferences of the consumer. Careful selection against the regions associated with negative effects will be required to avoid unwanted grain quality characteristics during the development of improved varieties for yield and yield components using introgressions from O. rufipogon.