Positive effects of public breeding on US rice yields under future climate scenarios.

In this study, we model and predict rice yields by integrating molecular marker variation, varietal productivity, and climate, focusing on the Southern U.S. rice-growing region. This region spans the states of Arkansas, Louisiana, Texas, Mississippi, and Missouri and accounts for 85% of total U.S. rice production. By digitizing and combining four decades of county-level variety acreage data (1970 to 2015) with varietal information from genotyping-by-sequencing data, we estimate annual historical county-level allele frequencies. These allele frequencies are used together with county-level weather and yield data to develop ten machine learning models for yield prediction. A two-layer meta-learner ensemble model that combines all ten methods is externally evaluated against observations from historical Uniform Regional Rice Nursery trials (1980 to 2018) conducted in the same states. Finally, the ensemble model is used with forecasted weather from the Coupled Model Intercomparison Project across the 110 rice-growing counties to predict production in the coming decades for Composite Variety Groups assembled based on year of release, breeding program, and several breeding trends. Results indicate positive effects over time of public breeding on rice resilience to future climates, and potential reasons are discussed.

Regulator of Awn Elongation 3, an E3 ubiquitin ligase, is responsible for loss of awns during African rice domestication.

Two species of rice have been independently domesticated from different ancestral wild species in Asia and Africa. Comparison of mutations that underlie phenotypic and physiological alterations associated with domestication traits in these species gives insights into the domestication history of rice in both regions. Asian cultivated rice, Oryza sativa, and African cultivated rice, Oryza glaberrima, have been modified and improved for common traits beneficial for humans, including erect plant architecture, nonshattering seeds, nonpigmented pericarp, and lack of awns. Independent mutations in orthologous genes associated with these traits have been documented in the two cultivated species. Contrary to this prevailing model, selection for awnlessness targeted different genes in O. sativa and O. glaberrima. We identify Regulator of Awn Elongation 3 (RAE3) a gene that encodes an E3 ubiquitin ligase and is responsible for the awnless phenotype only in O. glaberrima. A 48-bp deletion may disrupt the substrate recognition domain in RAE3 and diminish awn elongation. Sequencing analysis demonstrated low nucleotide diversity in a ~600-kb region around the derived rae3 allele on chromosome 6 in O. glaberrima compared with its wild progenitor. Identification of RAE3 sheds light on the molecular mechanism underlying awn development and provides an example of how selection on different genes can confer the same domestication phenotype in Asian and African rice.

Genome wide association analysis of root hair traits in rice reveals novel genomic regions controlling epidermal cell differentiation.

BACKGROUND: Genome wide association (GWA) studies demonstrate linkages between genetic variants and traits of interest. Here, we tested associations between single nucleotide polymorphisms (SNPs) in rice (Oryza sativa) and two root hair traits, root hair length (RHL) and root hair density (RHD). Root hairs are outgrowths of single cells on the root epidermis that aid in nutrient and water acquisition and have also served as a model system to study cell differentiation and tip growth. Using lines from the Rice Diversity Panel-1, we explored the diversity of root hair length and density across four subpopulations of rice (aus, indica, temperate japonica, and tropical japonica). GWA analysis was completed using the high-density rice array (HDRA) and the rice reference panel (RICE-RP) SNP sets. RESULTS: We identified 18 genomic regions related to root hair traits, 14 of which related to RHD and four to RHL. No genomic regions were significantly associated with both traits. Two regions overlapped with previously identified quantitative trait loci (QTL) associated with root hair density in rice. We identified candidate genes in these regions and present those with previously published expression data relevant to root hair development. We re-phenotyped a subset of lines with extreme RHD phenotypes and found that the variation in RHD was due to differences in cell differentiation, not cell size, indicating genes in an associated genomic region may influence root hair cell fate. The candidate genes that we identified showed little overlap with previously characterized genes in rice and Arabidopsis. CONCLUSIONS: Root hair length and density are quantitative traits with complex and independent genetic control in rice. The genomic regions described here could be used as the basis for QTL development and further analysis of the genetic control of root hair length and density. We present a list of candidate genes involved in root hair formation and growth in rice, many of which have not been previously identified as having a relation to root hair growth. Since little is known about root hair growth in grasses, these provide a guide for further research and crop improvement.

Can biochemical traits bridge the gap between genomics and plant performance? A study in rice under drought.

The possibility of introducing metabolic/biochemical phenotyping to complement genomics-based predictions in breeding pipelines has been considered for years. Here we examine to what extent and under what environmental conditions metabolic/biochemical traits can effectively contribute to understanding and predicting plant performance. In this study, multivariable statistical models based on flag leaf central metabolism and oxidative stress status were used to predict grain yield (GY) performance for 271 indica rice (Oryza sativa) accessions grown in the field under well-watered and reproductive stage drought conditions. The resulting models displayed significantly higher predictability than multivariable models based on genomic data for the prediction of GY under drought (Q2 = 0.54-0.56 versus 0.35) and for stress-induced GY loss (Q2 = 0.59-0.64 versus 0.03-0.06). Models based on the combined datasets showed predictabilities similar to metabolic/biochemical-based models alone. In contrast to genetic markers, models with enzyme activities and metabolite values also quantitatively integrated the effect of physiological differences such as plant height on GY. The models highlighted antioxidant enzymes of the ascorbate-glutathione cycle and a lipid oxidation stress marker as important predictors of rice GY stability under drought at the reproductive stage, and these stress-related variables were more predictive than leaf central metabolites. These findings provide evidence that metabolic/biochemical traits can integrate dynamic cellular and physiological responses to the environment and can help bridge the gap between the genome and the phenome of crops as predictors of GY performance under drought.

A massively parallel barcoded sequencing pipeline enables generation of the first ORFeome and interactome map for rice.

Systematic mappings of protein interactome networks have provided invaluable functional information for numerous model organisms. Here we develop PCR-mediated Linkage of barcoded Adapters To nucleic acid Elements for sequencing (PLATE-seq) that serves as a general tool to rapidly sequence thousands of DNA elements. We validate its utility by generating the ORFeome for Oryza sativa covering 2,300 genes and constructing a high-quality protein-protein interactome map consisting of 322 interactions between 289 proteins, expanding the known interactions in rice by roughly 50%. Our work paves the way for high-throughput profiling of protein-protein interactions in a wide range of organisms.

Genetic mapping identifies a rice naringenin O-glucosyltransferase that influences insect resistance.

Naringenin, the biochemical precursor for predominant flavonoids in grasses, provides protection against UV damage, pathogen infection and insect feeding. To identify previously unknown loci influencing naringenin accumulation in rice (Oryza sativa), recombinant inbred lines derived from the Nipponbare and IR64 cultivars were used to map a quantitative trait locus (QTL) for naringenin abundance to a region of 50 genes on rice chromosome 7. Examination of candidate genes in the QTL confidence interval identified four predicted uridine diphosphate-dependent glucosyltransferases (Os07g31960, Os07g32010, Os07g32020 and Os07g32060). In vitro assays demonstrated that one of these genes, Os07g32020 (UGT707A3), encodes a glucosyltransferase that converts naringenin and uridine diphosphate-glucose to naringenin-7-O-ß-d-glucoside. The function of Os07g32020 was verified with CRISPR/Cas9 mutant lines, which accumulated more naringenin and less naringenin-7-O-ß-d-glucoside and apigenin-7-O-ß-d-glucoside than wild-type Nipponbare. Expression of Os12g13800, which encodes a naringenin 7-O-methyltransferase that produces sakuranetin, was elevated in the mutant lines after treatment with methyl jasmonate and insect pests, Spodoptera litura (cotton leafworm), Oxya hyla intricata (rice grasshopper) and Nilaparvata lugens (brown planthopper), leading to a higher accumulation of sakuranetin. Feeding damage from O. hyla intricata and N. lugens was reduced on the Os07g32020 mutant lines relative to Nipponbare. Modification of the Os07g32020 gene could be used to increase the production of naringenin and sakuranetin rice flavonoids in a more targeted manner. These findings may open up new opportunities for selective breeding of this important rice metabolic trait.

Genetic architecture of root and shoot ionomes in rice (Oryza sativa L.).

KEY MESSAGE: Association analysis for ionomic concentrations of 20 elements identified independent genetic factors underlying the root and shoot ionomes of rice, providing a platform for selecting and dissecting causal genetic variants. Understanding the genetic basis of mineral nutrient acquisition is key to fully describing how terrestrial organisms interact with the non-living environment. Rice (Oryza sativa L.) serves both as a model organism for genetic studies and as an important component of the global food system. Studies in rice ionomics have primarily focused on above ground tissues evaluated from field-grown plants. Here, we describe a comprehensive study of the genetic basis of the rice ionome in both roots and shoots of 6-week-old rice plants for 20 elements using a controlled hydroponics growth system. Building on the wealth of publicly available rice genomic resources, including a panel of 373 diverse rice lines, 4.8 M genome-wide single-nucleotide polymorphisms, single- and multi-marker analysis pipelines, an extensive tome of 321 candidate genes and legacy QTLs from across 15 years of rice genetics literature, we used genome-wide association analysis and biparental QTL analysis to identify 114 genomic regions associated with ionomic variation. The genetic basis for root and shoot ionomes was highly distinct; 78 loci were associated with roots and 36 loci with shoots, with no overlapping genomic regions for the same element across tissues. We further describe the distribution of phenotypic variation across haplotypes and identify candidate genes within highly significant regions associated with sulfur, manganese, cadmium, and molybdenum. Our analysis provides critical insight into the genetic basis of natural phenotypic variation for both root and shoot ionomes in rice and provides a comprehensive resource for dissecting and testing causal genetic variants.

Bridging old and new: diversity and evaluation of high iron-associated stress response of rice cultivated in West Africa.

Adoption of rice varieties that perform well under high iron-associated (HIA) stress environments can enhance rice production in West Africa. This study reports the genetic characterization of 323 rice accessions and breeding lines cultivated in West Africa using genotyping-by-sequencing and their phenotypic response to HIA treatments in hydroponic solution (1500 mg l-1 FeSO4·7H2O) and hot-spot fields. The germplasm consisted of four genetic subpopulations: Oryza glaberrima (14%), O. sativa-japonica (7%), O. sativa-indica Group 1 (45%), and O. sativa-indica Group 2 (25%). Severe versus mild stress in the field was associated with a reduced SPAD value (12%), biomass (56%), and grain yield (57%), with leaf bronzing explaining 30% and 21% of the variation for biomass and grain yield, respectively. Association mapping using 175 indica genotypes identified 23 significant single nucleotide polymorphism (SNP) markers that mapped to 14 genomic regions. Genome-wide association study (GWAS) signals associated with leaf bronzing, a routinely used indicator of HIA stress, differed in hydroponic compared with field conditions. Contrastingly, six significant SNPs on chromosomes 8 and 9 were associated with the SPAD value under HIA stress in both field and hydroponic experiments, and a candidate potassium transporter gene mapped under the peak on chromosome 8. This study helps define criteria for assessing rice performance under HIA environments.

Association mapping and genetic dissection of drought-induced canopy temperature differences in rice.

Drought-stressed plants display reduced stomatal conductance, which results in increased leaf temperature by limiting transpiration. In this study, thermal imaging was used to quantify the differences in canopy temperature under drought in a rice diversity panel consisting of 293 indica accessions. The population was grown under paddy field conditions and drought stress was imposed for 2 weeks at flowering. The canopy temperature of the accessions during stress negatively correlated with grain yield (r= -0.48) and positively with plant height (r=0.56). Temperature values were used to perform a genome-wide association (GWA) analysis using a 45K single nucleotide polynmorphism (SNP) map. A quantitative trait locus (QTL) for canopy temperature under drought was detected on chromosome 3 and fine-mapped using a high-density imputed SNP map. The candidate genes underlying the QTL point towards differences in the regulation of guard cell solute intake for stomatal opening as the possible source of temperature variation. Genetic variation for the significant markers of the QTL was present only within the tall, low-yielding landraces adapted to drought-prone environments. The absence of variation in the shorter genotypes, which showed lower leaf temperature and higher grain yield, suggests that breeding for high grain yield in rice under paddy conditions has reduced genetic variation for stomatal response under drought.

Loss of function at RAE2, a previously unidentified EPFL, is required for awnlessness in cultivated Asian rice.

Domestication of crops based on artificial selection has contributed numerous beneficial traits for agriculture. Wild characteristics such as red pericarp and seed shattering were lost in both Asian (Oryza sativa) and African (Oryza glaberrima) cultivated rice species as a result of human selection on common genes. Awnedness, in contrast, is a trait that has been lost in both cultivated species due to selection on different sets of genes. In a previous report, we revealed that at least three loci regulate awn development in rice; however, the molecular mechanism underlying awnlessness remains unknown. Here we isolate and characterize a previously unidentified EPIDERMAL PATTERNING FACTOR-LIKE (EPFL) family member named REGULATOR OF AWN ELONGATION 2 (RAE2) and identify one of its requisite processing enzymes, SUBTILISIN-LIKE PROTEASE 1 (SLP1). The RAE2 precursor is specifically cleaved by SLP1 in the rice spikelet, where the mature RAE2 peptide subsequently induces awn elongation. Analysis of RAE2 sequence diversity identified a highly variable GC-rich region harboring multiple independent mutations underlying protein-length variation that disrupt the function of the RAE2 protein and condition the awnless phenotype in Asian rice. Cultivated African rice, on the other hand, retained the functional RAE2 allele despite its awnless phenotype. Our findings illuminate the molecular function of RAE2 in awn development and shed light on the independent domestication histories of Asian and African cultivated rice.

Genomic regions responsible for seminal and crown root lengths identified by 2D & 3D root system image analysis.

BACKGROUND: Genetic improvement of root system architecture is a promising approach for improved uptake of water and mineral nutrients distributed unevenly in the soil. To identify genomic regions associated with the length of different root types in rice, we quantified root system architecture in a set of 26 chromosome segment substitution lines derived from a cross between lowland indica rice, IR64, and upland tropical japonica rice, Kinandang Patong, (IK-CSSLs), using 2D & 3D root phenotyping platforms. RESULTS: Lengths of seminal and crown roots in the IK-CSSLs grown under hydroponic conditions were measured by 2D image analysis (RootReader2D). Twelve CSSLs showed significantly longer seminal root length than the recurrent parent IR64. Of these, 8 CSSLs also exhibited longer total length of the three longest crown roots compared to IR64. Three-dimensional image analysis (RootReader3D) for these CSSLs grown in gellan gum revealed that only one CSSL, SL1003, showed significantly longer total root length than IR64. To characterize the root morphology of SL1003 under soil conditions, SL1003 was grown in Turface, a soil-like growth media, and roots were quantified using RootReader3D. SL1003 had larger total root length and increased total crown root length than did IR64, although its seminal root length was similar to that of IR64. The larger TRL in SL1003 may be due to increased crown root length. CONCLUSIONS: SL1003 carries an introgression from Kinandang Patong on the long arm of chromosome 1 in the genetic background of IR64. We conclude that this region harbors a QTL controlling crown root elongation.

LABA1, a Domestication Gene Associated with Long, Barbed Awns in Wild Rice.

Common wild rice (Oryza rufipogon), the wild relative of Asian cultivated rice (Oryza sativa), flaunts long, barbed awns, which are necessary for efficient propagation and dissemination of seeds. By contrast, O. sativa cultivars have been selected to be awnless or to harbor short, barbless awns, which facilitate seed processing and storage. The transition from long, barbed awns to short, barbless awns was a crucial event in rice domestication. Here, we show that the presence of long, barbed awns in wild rice is controlled by a major gene on chromosome 4, LONG AND BARBED AWN1 (LABA1), which encodes a cytokinin-activating enzyme. A frame-shift deletion in LABA1 of cultivated rice reduces the cytokinin concentration in awn primordia, disrupting barb formation and awn elongation. Sequencing analysis demonstrated low nucleotide diversity and a selective sweep encompassing an ∼800-kb region around the derived laba1 allele in cultivated rice. Haplotype analysis revealed that the laba1 allele originated in the japonica subspecies and moved into the indica gene pool via introgression, suggesting that humans selected for this locus in early rice domestication. Identification of LABA1 provides new insights into rice domestication and also sheds light on the molecular mechanism underlying awn development.

The Tyrosine Aminomutase TAM1 Is Required for ß-Tyrosine Biosynthesis in Rice.

Non-protein amino acids, often isomers of the standard 20 protein amino acids, have defense-related functions in many plant species. A targeted search for jasmonate-induced metabolites in cultivated rice (Oryza sativa) identified (R)-ß-tyrosine, an isomer of the common amino acid (S)-α-tyrosine in the seeds, leaves, roots, and root exudates of the Nipponbare cultivar. Assays with 119 diverse cultivars showed a distinct presence/absence polymorphism, with ß-tyrosine being most prevalent in temperate japonica cultivars. Genetic mapping identified a candidate gene on chromosome 12, which was confirmed to encode a tyrosine aminomutase (TAM1) by transient expression in Nicotiana benthamiana and in vitro enzyme assays. A point mutation in TAM1 eliminated ß-tyrosine production in Nipponbare. Rice cultivars that do not produce ß-tyrosine have a chromosome 12 deletion that encompasses TAM1. Although ß-tyrosine accumulation was induced by the plant defense signaling molecule jasmonic acid, bioassays with hemipteran and lepidopteran herbivores showed no negative effects at physiologically relevant ß-tyrosine concentrations. In contrast, root growth of Arabidopsis thaliana and other tested dicot plants was inhibited by concentrations as low as 1 µM. As ß-tyrosine is exuded into hydroponic medium at higher concentrations, it may contribute to the allelopathic potential of rice.

Genomic selection and association mapping in rice (Oryza sativa): effect of trait genetic architecture, training population composition, marker number and statistical model on accuracy of rice genomic selection in elite, tropical rice breeding lines.

Genomic Selection (GS) is a new breeding method in which genome-wide markers are used to predict the breeding value of individuals in a breeding population. GS has been shown to improve breeding efficiency in dairy cattle and several crop plant species, and here we evaluate for the first time its efficacy for breeding inbred lines of rice. We performed a genome-wide association study (GWAS) in conjunction with five-fold GS cross-validation on a population of 363 elite breeding lines from the International Rice Research Institute's (IRRI) irrigated rice breeding program and herein report the GS results. The population was genotyped with 73,147 markers using genotyping-by-sequencing. The training population, statistical method used to build the GS model, number of markers, and trait were varied to determine their effect on prediction accuracy. For all three traits, genomic prediction models outperformed prediction based on pedigree records alone. Prediction accuracies ranged from 0.31 and 0.34 for grain yield and plant height to 0.63 for flowering time. Analyses using subsets of the full marker set suggest that using one marker every 0.2 cM is sufficient for genomic selection in this collection of rice breeding materials. RR-BLUP was the best performing statistical method for grain yield where no large effect QTL were detected by GWAS, while for flowering time, where a single very large effect QTL was detected, the non-GS multiple linear regression method outperformed GS models. For plant height, in which four mid-sized QTL were identified by GWAS, random forest produced the most consistently accurate GS models. Our results suggest that GS, informed by GWAS interpretations of genetic architecture and population structure, could become an effective tool for increasing the efficiency of rice breeding as the costs of genotyping continue to decline.

Hybrid assembly with long and short reads improves discovery of gene family expansions.

BACKGROUND: Long-read and short-read sequencing technologies offer competing advantages for eukaryotic genome sequencing projects. Combinations of both may be appropriate for surveys of within-species genomic variation. METHODS: We developed a hybrid assembly pipeline called "Alpaca" that can operate on 20X long-read coverage plus about 50X short-insert and 50X long-insert short-read coverage. To preclude collapse of tandem repeats, Alpaca relies on base-call-corrected long reads for contig formation. RESULTS: Compared to two other assembly protocols, Alpaca demonstrated the most reference agreement and repeat capture on the rice genome. On three accessions of the model legume Medicago truncatula, Alpaca generated the most agreement to a conspecific reference and predicted tandemly repeated genes absent from the other assemblies. CONCLUSION: Our results suggest Alpaca is a useful tool for investigating structural and copy number variation within de novo assemblies of sampled populations.

The buffering capacity of stems: genetic architecture of nonstructural carbohydrates in cultivated Asian rice, Oryza sativa.

Harnessing stem carbohydrate dynamics in grasses offers an opportunity to help meet future demands for plant-based food, fiber and fuel production, but requires a greater understanding of the genetic controls that govern the synthesis, interconversion and transport of such energy reserves. We map out a blueprint of the genetic architecture of rice (Oryza sativa) stem nonstructural carbohydrates (NSC) at two critical developmental time-points using a subpopulation-specific genome-wide association approach on two diverse germplasm panels followed by quantitative trait loci (QTL) mapping in a biparental population. Overall, 26 QTL are identified; three are detected in multiple panels and are associated with starch-at-maturity, sucrose-at-maturity and NSC-at-heading. They tag OsHXK6 (rice hexokinase), ISA2 (rice isoamylase) and a tandem array of sugar transporters. This study provides the foundation for more in-depth molecular investigation to validate candidate genes underlying rice stem NSC and informs future comparative studies in other agronomically vital grass species.

Harvesting the promising fruits of genomics: applying genome sequencing technologies to crop breeding.

Next generation sequencing (NGS) technologies are being used to generate whole genome sequences for a wide range of crop species. When combined with precise phenotyping methods, these technologies provide a powerful and rapid tool for identifying the genetic basis of agriculturally important traits and for predicting the breeding value of individuals in a plant breeding population. Here we summarize current trends and future prospects for utilizing NGS-based technologies to develop crops with improved trait performance and increase the efficiency of modern plant breeding. It is our hope that the application of NGS technologies to plant breeding will help us to meet the challenge of feeding a growing world population.

Natural variation underlies alterations in Nramp aluminum transporter (NRAT1) expression and function that play a key role in rice aluminum tolerance.

Aluminum (Al) toxicity is a major constraint for crop production on acid soils which compose ∼ 40% of arable land in the tropics and subtropics. Rice is the most Al-tolerant cereal crop and offers a good model for identifying Al tolerance genes and mechanisms. Here we investigated natural variation in the rice Nramp aluminum transporter (NRAT1) gene encoding a root plasma membrane Al uptake transporter previously hypothesized to underlie a unique Al tolerance mechanism. DNA sequence variation in the NRAT1 coding and regulatory regions was associated with changes in NRAT1 expression and NRAT1 Al transport properties. These sequence changes resulted in significant differences in Al tolerance that were found to be associated with changes in the Al content of root cell wall and cell sap in 24 representative rice lines from a rice association panel. Expression of the tolerant OsNRAT1 allele in yeast resulted in higher Al uptake than did the sensitive allele and conferred greater Al tolerance when expressed in transgenic Arabidopsis. These findings indicate that NRAT1 plays an important role in rice Al tolerance by reducing the level of toxic Al in the root cell wall and transporting Al into the root cell, where it is ultimately sequestered in the vacuole. Given its ability to enhance Al tolerance in rice and Arabidopsis, this work suggests that the NRAT1 gene or its orthologs may be useful tools for enhancing Al tolerance in a wide range of plant species.