Life-Cycle Multiomics of Rice Shoots Reveals Growth Stage-Specific Effects of Drought Stress and Time-Lag Drought Responses.

Field-grown rice plants are exposed to various stresses at different stages of their life cycle, but little is known about the effects of stage-specific stresses on phenomes and transcriptomes. In this study, we performed integrated time-course multiomics on rice at 3-d intervals from seedling to heading stage under six drought conditions in a well-controlled growth chamber. Drought stress at seedling and reproductive stages reduced yield performance by reducing seed number and setting rate, respectively. High temporal resolution analysis revealed that drought response occurred in two steps: a rapid response via the abscisic acid (ABA) signaling pathway and a slightly delayed DEHYDRATION-RESPONSIVE ELEMENT-BINDING PROTEIN (DREB) pathway, allowing plants to respond flexibly to deteriorating soil water conditions. Our long-term time-course multiomics showed that temporary drought stress delayed flowering due to prolonged expression of the flowering repressor gene GRAIN NUMBER, PLANT HEIGHT AND HEADING DATE 7 (Ghd7) and delayed expression of the florigen genes HEADING DATE 3a (Hd3a) and RICE FLOWERING LOCUS T 1 (RFT1). Our life-cycle multiomics dataset on rice shoots under drought conditions provides a valuable resource for further functional genomic studies to improve crop resilience to drought stress.

Development of 12 sets of chromosome segment substitution lines that enhance allele mining in Asian cultivated rice.

Many agronomic traits that are important in rice breeding are controlled by multiple genes. The extensive time and effort devoted so far to identifying and selecting such genes are still not enough to target multiple agronomic traits in practical breeding in Japan because of a lack of suitable plant materials in which to efficiently detect and validate beneficial alleles from diverse genetic resources. To facilitate the comprehensive analysis of genetic variation in agronomic traits among Asian cultivated rice, we developed 12 sets of chromosome segment substitution lines (CSSLs) with the japonica background, 11 of them in the same genetic background, using donors representing the genetic diversity of Asian cultivated rice. Using these materials, we overviewed the chromosomal locations of 1079 putative QTLs for seven agronomic traits and their allelic distribution in Asian cultivated rice through multiple linear regression analysis. The CSSLs will allow the effects of putative QTLs in the highly homogeneous japonica background to be validated.

Transcriptome profiles of rice roots under simulated microgravity conditions and following gravistimulation.

Root system architecture affects the efficient uptake of water and nutrients in plants. The root growth angle, which is a critical component in determining root system architecture, is affected by root gravitropism; however, the mechanism of root gravitropism in rice remains largely unknown. In this study, we conducted a time-course transcriptome analysis of rice roots under conditions of simulated microgravity using a three-dimensional clinostat and following gravistimulation to detect candidate genes associated with the gravitropic response. We found that HEAT SHOCK PROTEIN (HSP) genes, which are involved in the regulation of auxin transport, were preferentially up-regulated during simulated microgravity conditions and rapidly down-regulated by gravistimulation. We also found that the transcription factor HEAT STRESS TRANSCRIPTION FACTOR A2s (HSFA2s) and HSFB2s, showed the similar expression patterns with the HSPs. A co-expression network analysis and an in silico motif search within the upstream regions of the co-expressed genes revealed possible transcriptional control of HSPs by HSFs. Because HSFA2s are transcriptional activators, whereas HSFB2s are transcriptional repressors, the results suggest that the gene regulatory networks governed by HSFs modulate the gravitropic response through transcriptional control of HSPs in rice roots.

Quantification of Soil-surface Roots in Seedlings and Mature Rice Plants.

Soil-surface roots (SORs) in rice are primary roots that elongate over or near the soil surface. SORs help avoid excessive reduction of stress that occurs in paddy, such as in saline conditions. SORs may also be beneficial for rice growth in phosphorus-deficient paddy fields. Thus, SOR is a useful trait for crop adaptation to certain environmental stresses. To identify a promising genetic material showing SOR, we established methods for evaluating SOR under different growth conditions. We introduced procedures to evaluate the genetic diversity of SOR in various growth stages and conditions: the Cup method allowed us to quantify SOR at the seedling stage, and the Basket method, using a basket buried in a pot or field, is useful in quantifying SOR at the adult stage. These protocols are expected to contribute not only to the evaluation of the genetic diversity of SOR, but also the isolation of related genes in rice.

The transcriptomic landscapes of rice cultivars with diverse root system architectures grown in upland field conditions.

Root system architecture affects plant drought resistance and other key agronomic traits such as lodging. However, although phenotypic and genomic variation has been extensively analyzed, few field studies have integrated phenotypic and transcriptomic information, particularly for below-ground traits such as root system architecture. Here, we report the phenotypic and transcriptomic landscape of 61 rice (Oryza sativa) accessions with highly diverse below-ground traits grown in an upland field. We found that four principal components explained the phenotypic variation and that accessions could be classified into four subpopulations (indica, aus, japonica and admixed) based on their tiller numbers and crown root diameters. Transcriptome analysis revealed that differentially expressed genes associated with specific subpopulations were enriched with stress response-related genes, suggesting that subpopulations have distinct stress response mechanisms. Root growth was negatively correlated with auxin-inducible genes, suggesting an association between auxin signaling and upland field conditions. A negative correlation between crown root diameter and stress response-related genes suggested that thicker crown root diameter is associated with resistance to mild drought stress. Finally, co-expression network analysis implemented with DNA affinity purification followed by sequencing analysis identified phytohormone signaling networks and key transcription factors negatively regulating crown root diameter. Our datasets provide a useful resource for understanding the genomic and transcriptomic basis of phenotypic variation under upland field conditions.

Root angle modifications by the DRO1 homolog improve rice yields in saline paddy fields.

The root system architecture (RSA) of crops can affect their production, particularly in abiotic stress conditions, such as with drought, waterlogging, and salinity. Salinity is a growing problem worldwide that negatively impacts on crop productivity, and it is believed that yields could be improved if RSAs that enabled plants to avoid saline conditions were identified. Here, we have demonstrated, through the cloning and characterization of qSOR1 (quantitative trait locus for SOIL SURFACE ROOTING 1), that a shallower root growth angle (RGA) could enhance rice yields in saline paddies. qSOR1 is negatively regulated by auxin, predominantly expressed in root columella cells, and involved in the gravitropic responses of roots. qSOR1 was found to be a homolog of DRO1 (DEEPER ROOTING 1), which is known to control RGA. CRISPR-Cas9 assays revealed that other DRO1 homologs were also involved in RGA. Introgression lines with combinations of gain-of-function and loss-of-function alleles in qSOR1 and DRO1 demonstrated four different RSAs (ultra-shallow, shallow, intermediate, and deep rooting), suggesting that natural alleles of the DRO1 homologs could be utilized to control RSA variations in rice. In saline paddies, near-isogenic lines carrying the qSOR1 loss-of-function allele had soil-surface roots (SOR) that enabled rice to avoid the reducing stresses of saline soils, resulting in increased yields compared to the parental cultivars without SOR. Our findings suggest that DRO1 homologs are valuable targets for RSA breeding and could lead to improved rice production in environments characterized by abiotic stress.

High-throughput three-dimensional visualization of root system architecture of rice using X-ray computed tomography.

BACKGROUND: X-ray computed tomography (CT) allows us to visualize root system architecture (RSA) beneath the soil, non-destructively and in a three-dimensional (3-D) form. However, CT scanning, reconstruction processes, and root isolation from X-ray CT volumes, take considerable time. For genetic analyses, such as quantitative trait locus mapping, which require a large population size, a high-throughput RSA visualization method is required. RESULTS: We have developed a high-throughput process flow for the 3-D visualization of rice (Oryza sativa) RSA (consisting of radicle and crown roots), using X-ray CT. The process flow includes use of a uniform particle size, calcined clay to reduce the possibility of visualizing non-root segments, use of a higher tube voltage and current in the X-ray CT scanning to increase root-to-soil contrast, and use of a 3-D median filter and edge detection algorithm to isolate root segments. Using high-performance computing technology, this analysis flow requires only 10 min (33 s, if a rough image is acceptable) for CT scanning and reconstruction, and 2 min for image processing, to visualize rice RSA. This reduced time allowed us to conduct the genetic analysis associated with 3-D RSA phenotyping. In 2-week-old seedlings, 85% and 100% of radicle and crown roots were detected, when 16 cm and 20 cm diameter pots were used, respectively. The X-ray dose per scan was estimated at < 0.09 Gy, which did not impede rice growth. Using the developed process flow, we were able to follow daily RSA development, i.e., 4-D RSA development, of an upland rice variety, over 3 weeks. CONCLUSIONS: We developed a high-throughput process flow for 3-D rice RSA visualization by X-ray CT. The X-ray dose assay on plant growth has shown that this methodology could be applicable for 4-D RSA phenotyping. We named the RSA visualization method 'RSAvis3D' and are confident that it represents a potentially efficient application for 3-D RSA phenotyping of various plant species.

De novo Genome Assembly of the indica Rice Variety IR64 Using Linked-Read Sequencing and Nanopore Sequencing.

IR64 is a rice variety with high-yield that has been widely cultivated around the world. IR64 has been replaced by modern varieties in most growing areas. Given that modern varieties are mostly progenies or relatives of IR64, genetic analysis of IR64 is valuable for rice functional genomics. However, chromosome-level genome sequences of IR64 have not been available previously. Here, we sequenced the IR64 genome using synthetic long reads obtained by linked-read sequencing and ultra-long reads obtained by nanopore sequencing. We integrated these data and generated the de novo assembly of the IR64 genome of 367 Mb, equivalent to 99% of the estimated size. Continuity of the IR64 genome assembly was improved compared with that of a publicly available IR64 genome assembly generated by short reads only. We annotated 41,458 protein-coding genes, including 657 IR64-specific genes, that are missing in other high-quality rice genome assemblies IRGSP-1.0 of japonica cultivar Nipponbare or R498 of indica cultivar Shuhui498. The IR64 genome assembly will serve as a genome resource for rice functional genomics as well as genomics-driven and/or molecular breeding.

Backhoe-assisted monolith method for plant root phenotyping under upland conditions.

Root system architecture (RSA) is one of the most important traits determining water and nutrient availability for plants. Modification of RSA is known to be a useful approach for improving root performance of crops. However, for conducting root phenotyping, there are few alternatives for the rapid collection of root samples from a constant soil volume. In this report, we propose a rapid root-sampling method, which uses a steel cylinder known as round monolith and backhoes to reduce the physical effort. The monolith was set on the ground surrounding individual rice plants and vertically driven back by a backhoe. Soil samples with 20 cm width and 25 cm depth were excavated by the monolith, from which root samples were then isolated. This backhoe-assisted monolith method requires at most five minutes to collect root samples from one plant. Using this method, we quantified the root traits of three rice lines, reported to form different types of root system such as shallow-, intermediate-, and deep-roots, using a root image analysis software. The data obtained through this method, which showed the same trend as previously reported, clearly demonstrated that this method is useful for quantitative evaluation of roots in the soil.

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.

Fine Mapping of QUICK ROOTING 1 and 2, Quantitative Trait Loci Increasing Root Length in Rice.

The volume that the root system can occupy is associated with the efficiency of water and nutrient uptake from soil. Genetic improvement of root length, which is a limiting factor for root distribution, is necessary for increasing crop production. In this report, we describe identification of two quantitative trait loci (QTLs) for maximal root length, QUICK ROOTING 1 (QRO1) on chromosome 2 and QRO2 on chromosome 6, in cultivated rice (Oryza sativa L.). We measured the maximal root length in 26 lines carrying chromosome segments from the long-rooted upland rice cultivar Kinandang Patong in the genetic background of the short-rooted lowland cultivar IR64. Five lines had longer roots than IR64. By rough mapping of the target regions in BC4F2 populations, we detected putative QTLs for maximal root length on chromosomes 2, 6, and 8. To fine-map these QTLs, we used BC4F3 recombinant homozygous lines. QRO1 was mapped between markers RM5651 and RM6107, which delimit a 1.7-Mb interval on chromosome 2, and QRO2 was mapped between markers RM20495 and RM3430-1, which delimit an 884-kb interval on chromosome 6. Both QTLs may be promising gene resources for improving root system architecture in rice.

Jasmonate regulates juvenile-to-adult phase transition in rice.

Juvenile-to-adult phase transition is an important shift for the acquisition of adult vegetative characteristics and subsequent reproductive competence. We identified a recessive precocious (pre) mutant exhibiting a long leaf phenotype in rice. The long leaf phenotype is conspicuous in the second to the fourth leaves, which are juvenile and juvenile-to-adult transition leaves. We found that morphological and physiological traits, such as midrib formation, shoot meristem size, photosynthetic rate and plastochron, in juvenile and juvenile-to-adult transition stages of the pre mutant have precociously acquired adult characteristics. In agreement with these results, expression patterns of miR156 and miR172, which are microRNAs regulating phase change, support the accelerated juvenile-to-adult phase change in the pre mutant. The mutated gene encodes an allene oxide synthase (OsAOS1), which is a key enzyme for the biosynthesis of jasmonic acid (JA). The pre mutant showed a low level of JA and enhanced sensitivity to gibberellic acid, which promotes the phase change in some plant species. We also show that prolonged plastochron in the pre mutant is caused by accelerated PLASTOCHRON1 (PLA1) function. The present study reveals a substantial role of JA as a negative regulator of vegetative phase change.

Genetic improvement for root growth angle to enhance crop production.

The root system is an essential organ for taking up water and nutrients and anchoring shoots to the ground. On the other hand, the root system has rarely been regarded as breeding target, possibly because it is more laborious and time-consuming to evaluate roots (which require excavation) in a large number of plants than aboveground tissues. The root growth angle (RGA), which determines the direction of root elongation in the soil, affects the area in which roots capture water and nutrients. In this review, we describe the significance of RGA as a potential trait to improve crop production, and the physiological and molecular mechanisms that regulate RGA. We discuss the prospects for breeding to improve RGA based on current knowledge of quantitative trait loci for RGA in rice.

A QTL for root growth angle on rice chromosome 7 is involved in the genetic pathway of DEEPER ROOTING 1.

BACKGROUND: Root growth angle (RGA) is an important trait that influences the ability of rice to avoid drought stress. DEEPER ROOTING 1 (DRO1), which is a major quantitative trait locus (QTL) for RGA, is responsible for the difference in RGA between the shallow-rooting cultivar IR64 and the deep-rooting cultivar Kinandang Patong. However, the RGA differences between these cultivars cannot be fully explained by DRO1. The objective of this study was to identify new QTLs for RGA explaining the difference in RGA between these cultivars. RESULTS: By crossing IR64 (which has a non-functional allele of DRO1) with Kinandang Patong (which has a functional allele of DRO1), we developed 26 chromosome segment substitution lines (CSSLs) that carried a particular chromosome segment from Kinandang Patong in the IR64 genetic background. Using these CSSLs, we found only one chromosomal region that was related to RGA: on chromosome 9, which includes DRO1. Using an F2 population derived from a cross between Kinandang Patong and the Dro1-NIL (near isogenic line), which had a functional DRO1 allele in the IR64 genetic background, we identified a new QTL for RGA (DRO3) on the long arm of chromosome 7. CONCLUSIONS: DRO3 may only affect RGA in plants with a functional DRO1 allele, suggesting that DRO3 is involved in the DRO1 genetic pathway.

QTLs underlying natural variation of root growth angle among rice cultivars with the same functional allele of DEEPER ROOTING 1.

BACKGROUND: The functional allele of the rice gene DEEPER ROOTING 1 (DRO1) increases the root growth angle (RGA). However, wide natural variation in RGA is observed among rice cultivars with the functional DRO1 allele. To elucidate genetic factors related to such variation, we quantitatively measured RGA using the basket method and analyzed quantitative trait loci (QTLs) for RGA in three F2 mapping populations derived from crosses between the large RGA-type cultivar Kinandang Patong and each of three accessions with varying RGA: Momiroman has small RGA and was used to produce the MoK-F2 population; Yumeaoba has intermediate RGA (YuK-F2 population); Tachisugata has large RGA (TaK-F2 population). All four accessions belong to the same haplotype group of functional DRO1 allele. RESULTS: We detected the following statistically significant QTLs: one QTL on chromosome 4 in MoK-F2, three QTLs on chromosomes 2, 4, and 6 in YuK-F2, and one QTL on chromosome 2 in TaK-F2. Among them, the two QTLs on chromosome 4 were located near DRO2, which has been previously reported as a major QTL for RGA, whereas the two major QTLs for RGA on chromosomes 2 (DRO4) and 6 (DRO5) were novel. With the LOD threshold reduced to 3.0, several minor QTLs for RGA were also detected in each population. CONCLUSION: Natural variation in RGA in rice cultivars carrying functional DRO1 alleles may be controlled by a few major QTLs and by several additional minor QTLs.

Control of root system architecture by DEEPER ROOTING 1 increases rice yield under drought conditions.

The genetic improvement of drought resistance is essential for stable and adequate crop production in drought-prone areas. Here we demonstrate that alteration of root system architecture improves drought avoidance through the cloning and characterization of DEEPER ROOTING 1 (DRO1), a rice quantitative trait locus controlling root growth angle. DRO1 is negatively regulated by auxin and is involved in cell elongation in the root tip that causes asymmetric root growth and downward bending of the root in response to gravity. Higher expression of DRO1 increases the root growth angle, whereby roots grow in a more downward direction. Introducing DRO1 into a shallow-rooting rice cultivar by backcrossing enabled the resulting line to avoid drought by increasing deep rooting, which maintained high yield performance under drought conditions relative to the recipient cultivar. Our experiments suggest that control of root system architecture will contribute to drought avoidance in crops.

OsIAA13-mediated auxin signaling is involved in lateral root initiation in rice.

The plant hormone auxin is essential for root formation. After auxin perception, transmission of the auxin signal progresses through the degradation of Aux/IAA proteins. In this study, we newly isolated and characterized a rice gain-of-function mutant, Osiaa13, containing a single amino acid substitution in the core sequence required for the degradation of the OsIAA13 protein. The Osiaa13 mutant displayed typical auxin-related phenotypes: the number of lateral roots was significantly reduced and the root gravitropic response was defective. Osiaa13 mutants also exhibited altered GUS staining controlled by the DR5 promoter in lateral root initiation sites. Furthermore, expression levels of several genes that might be associated with lateral root initiation were altered in Osiaa13. Taken together, our results indicate that OsIAA13 is involved in auxin signaling and controls the expression of genes that are required for lateral root initiation in rice.

Molecular mechanism of crown root initiation and the different mechanisms between crown root and radicle in rice.

Monocot plants produce numerous adventitious (crown) roots. The plant hormone auxin has positive effects on crown root formation, while cytokinin suppresses it. We have demonstrated that auxin-induced CROWN ROOTLESS5 (CRL5) regulates crown root initiation in rice through the induction of OsRR1, a negative regulator of cytokinin signaling. CRL5 overexpressing calli formed adventitious roots, although CRL5 overexpressing plants did not induce ectopic roots, suggesting that CRL5, which promotes de novo root initiation, might function only in de-differentiated cells. A radicle initiated normally in a crl5 mutant, in spite of the defect in crown root initiation, whereas crown roots, but not a radicle, were produced in a radicleless1 (ral1) mutant. A crl5 ral1 double mutant displayed an additive phenotype, showing that the formation of each root is regulated by different genetic mechanisms in rice.

The auxin responsive AP2/ERF transcription factor CROWN ROOTLESS5 is involved in crown root initiation in rice through the induction of OsRR1, a type-A response regulator of cytokinin signaling.

Cytokinin is known to have negative effects on de novo auxin-induced root formation. However, the regulatory mechanisms of root initiation by both cytokinin and auxin are poorly understood. In this study, we characterized a rice mutant, termed crown rootless5 (crl5), which produced fewer crown roots and displayed impaired initiation of crown root primordia. The expression of CRL5, which encodes a member of the large AP2/ERF transcription factor family protein, was observed in the stem region where crown root initiation occurs. Exogenous auxin treatment induced CRL5 expression without de novo protein biosynthesis, which also required the degradation of AUX/IAA proteins. A putative auxin response element in the CRL5 promoter region specifically interacted with a rice ARF, demonstrating that CRL5 may be a direct target of an ARF, similar to CRL1/ADVENTITIOUS ROOTLESS1 (ARL1) that also regulates crown root initiation. A crl1 crl5 double mutant displayed an additive phenotype, indicating that these two genes function in different genetic pathways for crown root initiation. In addition, ProACT:CRL5/WT showed a cytokinin-resistant phenotype for crown root initiation, and also up-regulated the expression of two negative regulators of cytokinin signaling, OsRR1 and OsRR2, which were downregulated in crl5. Transgenic plants that over-expressed OsRR1 under the control of the CRL5 promoter in a crl5 mutant background produced a higher number of crown roots than the crl5 plant. Taken together, these results indicate that auxin-induced CRL5 promotes crown root initiation through repression of cytokinin signaling by positively regulating type-A RR, OsRR1.