Reaction norm for genomic prediction of plant growth: modeling drought stress response in soybean.

KEY MESSAGE: We proposed models to predict the effects of genomic and environmental factors on daily soybean growth and applied them to soybean growth data obtained with unmanned aerial vehicles. Advances in high-throughput phenotyping technology have made it possible to obtain time-series plant growth data in field trials, enabling genotype-by-environment interaction (G × E) modeling of plant growth. Although the reaction norm is an effective method for quantitatively evaluating G × E and has been implemented in genomic prediction models, no reaction norm models have been applied to plant growth data. Here, we propose a novel reaction norm model for plant growth using spline and random forest models, in which daily growth is explained by environmental factors one day prior. The proposed model was applied to soybean canopy area and height to evaluate the influence of drought stress levels. Changes in the canopy area and height of 198 cultivars were measured by remote sensing using unmanned aerial vehicles. Multiple drought stress levels were set as treatments, and their time-series soil moisture was measured. The models were evaluated using three cross-validation schemes. Although accuracy of the proposed models did not surpass that of single-trait genomic prediction, the results suggest that our model can capture G × E, especially the latter growth period for the random forest model. Also, significant variations in the G × E of the canopy height during the early growth period were visualized using the spline model. This result indicates the effectiveness of the proposed models on plant growth data and the possibility of revealing G × E in various growth stages in plant breeding by applying statistical or machine learning models to time-series phenotype data.

Genetic control of morphological traits useful for improving sorghum.

Global climate change and global warming, coupled with the growing population, have raised concerns about sustainable food supply and bioenergy demand. Sorghum [Sorghum bicolor (L.) Moench] ranks fifth among cereals produced worldwide; it is a C4 crop with a higher stress tolerance than other major cereals and has a wide range of uses, such as grains, forage, and biomass. Therefore, sorghum has attracted attention as a promising crop for achieving sustainable development goals (SDGs). In addition, sorghum is a suitable genetic model for C4 grasses because of its high morphological diversity and relatively small genome size compared to other C4 grasses. Although sorghum breeding and genetic studies have lagged compared to other crops such as rice and maize, recent advances in research have identified several genes and many quantitative trait loci (QTLs) that control important agronomic traits in sorghum. This review outlines traits and genetic information with a focus on morphogenetic aspects that may be useful in sorghum breeding for grain and biomass utilization.

DOMINANT AWN INHIBITOR Encodes the ALOG Protein Originating from Gene Duplication and Inhibits AWN Elongation by Suppressing Cell Proliferation and Elongation in Sorghum.

The awn, a needle-like structure extending from the tip of the lemma in grass species, plays a role in environmental adaptation and fitness. In some crops, awns appear to have been eliminated during domestication. Although numerous genes involved in awn development have been identified, several dominant genes that eliminate awns are also known to exist. For example, in sorghum (Sorghum bicolor), the dominant awn-inhibiting gene has been known since 1921; however, its molecular features remain uncharacterized. In this study, we conducted quantitative trait locus analysis and a genome-wide association study of awn-related traits in sorghum and identified DOMINANT AWN INHIBITOR (DAI), which encodes the ALOG family protein on chromosome 3. DAI appeared to be present in most awnless sorghum cultivars, likely because of its effectiveness. Detailed analysis of the ALOG protein family in cereals revealed that DAI originated from a duplication of its twin paralog (DAIori) on chromosome 10. Observations of immature awns in near-isogenic lines revealed that DAI inhibits awn elongation by suppressing both cell proliferation and elongation. We also found that only DAI gained a novel function to inhibit awn elongation through an awn-specific expression pattern distinct from that of DAIori. Interestingly, heterologous expression of DAI with its own promoter in rice inhibited awn elongation in the awned cultivar Kasalath. We found that DAI originated from gene duplication, providing an interesting example of gain-of-function that occurs only in sorghum but shares its functionality with rice and sorghum.

Sorghum Ionomics Reveals the Functional SbHMA3a Allele that Limits Excess Cadmium Accumulation in Grains.

Understanding uptake and redistribution of essential minerals or sequestering of toxic elements is important for optimized crop production. Although the mechanisms controlling mineral transport have been elucidated in rice and other species, little is understood in sorghum-an important C4 cereal crop. Here, we assessed the genetic factors that govern grain ionome profiles in sorghum using recombinant inbred lines (RILs) derived from a cross between BTx623 and NOG (Takakibi). Pairwise correlation and clustering analysis of 22 elements, measured in sorghum grains harvested under greenhouse conditions, indicated that the parental lines, as well as the RILs, show different ionomes. In particular, BTx623 accumulated significantly higher levels of cadmium (Cd) than NOG, because of differential root-to-shoot translocation factors between the two lines. Quantitative trait locus (QTL) analysis revealed a prominent QTL for grain Cd concentration on chromosome 2. Detailed analysis identified SbHMA3a, encoding a P1B-type ATPase heavy metal transporter, as responsible for low Cd accumulation in grains; the NOG allele encoded a functional HMA3 transporter (SbHMA3a-NOG) whose Cd-transporting activity was confirmed by heterologous expression in yeast. BTx623 possessed a truncated, loss-of-function SbHMA3a allele. The functionality of SbHMA3a in NOG was confirmed by Cd concentrations of F2 grains derived from the reciprocal cross, in which the NOG allele behaved in a dominant manner. We concluded that SbHMA3a-NOG is a Cd transporter that sequesters excess Cd in root tissues, as shown in other HMA3s. Our findings will facilitate the isolation of breeding cultivars with low Cd in grains or in exploiting high-Cd cultivars for phytoremediation.

Targeted gene disruption of ATP synthases 6-1 and 6-2 in the mitochondrial genome of Arabidopsis thaliana by mitoTALENs.

We recently achieved targeted disruptions of cytoplasmic male sterility (CMS)-associated genes in the mitochondrial genomes of rice and rapeseed by using mitochondria-targeted transcription activator-like effector nucleases (mitoTALENs). It was the first report of stable and heritable targeted gene modification of plant mitochondrial genomes. Here, we attempted to use mitoTALENs to disrupt two mitochondrial genes in the model plant Arabidopsis thaliana(Arabidopsis) using three different promoters and two types of TALENs. The targets were the two isoforms of the ATP synthase subunit 6 gene, atp6-1 and atp6-2. Each of these genes was successfully deleted and the mitochondrial genomes were recovered in a homoplasmic state. The nuclear genome also has a copy of atp6-1, and we were able to confirm that it was the mitochondrial gene and not the nuclear pseudogene that was knocked out. Among the three mitoTALEN promoters tried, the RPS5A promoter was the most effective. Conventional mitoTALENs were more effective than single-molecule mito-compactTALENs. Targeted mitochondrial gene deletion was achieved by crossing as well as by floral-dip transformation to introduce the mitoTALEN constructs into the nucleus. The gene disruptions were caused by large (kb-size) deletions. The ends of the remaining sequences were connected to distant loci, mostly by illegitimate homologous recombinations between repeats.

Transcriptional switch for programmed cell death in pith parenchyma of sorghum stems.

Pith parenchyma cells store water in various plant organs. These cells are especially important for producing sugar and ethanol from the sugar juice of grass stems. In many plants, the death of pith parenchyma cells reduces their stem water content. Previous studies proposed that a hypothetical D gene might be responsible for the death of stem pith parenchyma cells in Sorghum bicolor, a promising energy grass, although its identity and molecular function are unknown. Here, we identify the D gene and note that it is located on chromosome 6 in agreement with previous predictions. Sorghum varieties with a functional D allele had stems enriched with dry, dead pith parenchyma cells, whereas those with each of six independent nonfunctional D alleles had stems enriched with juicy, living pith parenchyma cells. D expression was spatiotemporally coupled with the appearance of dead, air-filled pith parenchyma cells in sorghum stems. Among D homologs that are present in flowering plants, Arabidopsis ANAC074 also is required for the death of stem pith parenchyma cells. D and ANAC074 encode previously uncharacterized NAC transcription factors and are sufficient to ectopically induce programmed death of Arabidopsis culture cells via the activation of autolytic enzymes. Taken together, these results indicate that D and its Arabidopsis ortholog, ANAC074, are master transcriptional switches that induce programmed death of stem pith parenchyma cells. Thus, targeting the D gene will provide an approach to breeding crops for sugar and ethanol production.

Spatial kernel models capturing field heterogeneity for accurate estimation of genetic potential.

According to Fisher's principles, an experimental field is typically divided into multiple blocks for local control. Although homogeneity is supposed within a block, this assumption may not be practical for large blocks, such as those including hundreds of plots. In line evaluation trials, which are essential in plant breeding, field heterogeneity must be carefully treated, because it can cause bias in the estimation of genetic potential. To more accurately estimate genotypic values in a large field trial, we developed spatial kernel models incorporating genome-wide markers, which consider continuous heterogeneity within a block and over the field. In the simulation study, the spatial kernel models were robust under various conditions. Although heritability, spatial autocorrelation range, replication number, and missing plots directly affected the estimation accuracy of genotypic values, the spatial kernel models always showed superior performance over the classical block model. We also employed these spatial kernel models for quantitative trait locus mapping. Finally, using field experimental data of bioenergy sorghum lines, we validated the performance of the spatial kernel models. The results suggested that a spatial kernel model is effective for evaluating the genetic potential of lines in a heterogeneous field.

RAD-seq-Based High-Density Linkage Map Construction and QTL Mapping of Biomass-Related Traits in Sorghum using the Japanese Landrace Takakibi NOG.

Sorghum [Sorghum bicolor (L.) Moench] grown locally by Japanese farmers is generically termed Takakibi, although its genetic diversity compared with geographically distant varieties or even within Takakibi lines remains unclear. To explore the genomic diversity and genetic traits controlling biomass and other physiological traits in Takakibi, we focused on a landrace, NOG, in this study. Admixture analysis of 460 sorghum accessions revealed that NOG belonged to the subgroup that represented Asian sorghums, and it was only distantly related to American/African accessions including BTx623. In an attempt to dissect major traits related to biomass, we generated a recombinant inbred line (RIL) from a cross between BTx623 and NOG, and we constructed a high-density linkage map based on 3,710 single-nucleotide polymorphisms obtained by restriction-site-associated DNA sequencing of 213 RIL individuals. Consequently, 13 fine quantitative trait loci (QTLs) were detected on chromosomes 2, 3, 6, 7, 8 and 9, which included five QTLs for days to heading, three for plant height (PH) and total shoot fresh weight and two for Brix. Furthermore, we identified two dominant loci for PH as being identical to the previously reported dw1 and dw3. Together, these results corroborate the diversified genome of Japanese Takakibi, while the RIL population and high-density linkage map generated in this study will be useful for dissecting other important traits in sorghum.

Impacts of dominance effects on genomic prediction of sorghum hybrid performance.

Non-additive (dominance and epistasis) effects have remarkable influences on hybrid performance, e.g., via heterosis. Nevertheless, only additive effects are often considered in genomic predictions (GP). In this study, we demonstrated the importance of dominance effects in the prediction of hybrid performance in bioenergy sorghum [Sorghum bicolor (L.) Moench]. The dataset contained more than 400 hybrids between 200 inbred lines and two testers. The hybrids exhibited considerable heterosis in culm length and fresh weight, and the degree of heterosis was consistent with the genetic distance from the corresponding tester. The degree of heterosis was further different among subpopulations. Conversely, Brix exhibited limited heterosis. Regarding GP, we examined three statistical models and four training dataset types. In most of the dataset types, genomic best linear unbiased prediction (GBLUP) with additive effects had lower prediction accuracy than GBLUP with additive and dominance effects (GBLUP-AD) and Gaussian kernel regression (GK). The superiority of GBLUP-AD and GK depended on the level of dominance variance, which was high for culm length and fresh weight, and low for Brix. Considering subpopulations, the influence of dominance was more complex. Our findings highlight the importance of considering dominance effects in GP models for sorghum hybrid breeding.

Effect of salt tolerance on biomass production in a large population of sorghum accessions.

Salinity causes major reductions in cultivated land area, crop productivity, and crop quality, and salt-tolerant crops have been required to sustain agriculture in salinized areas. The annual C4 crop plant Sorghum bicolor (L.) Moench is salt tolerant, with large variation among accessions. Sorghum's salt tolerance is often evaluated during early growth, but such evaluations are weakly related to overall performance. Here, we evaluated salt tolerance of 415 sorghum accessions grown in saline soil (0, 50, 100, and 150 mM NaCl) for 3 months. Some accessions produced up to 400 g per plant of biomass and showed no growth inhibition at 50 mM NaCl. Our analysis indicated that the genetic factors that affected biomass production under 100 mM salt stress were more different from those without salt stress, comparing to the differences between those under 50 mM and 100 mM salt stress. A genome-wide association study for salt tolerance identified two single-nucleotide polymorphisms (SNPs) that were significantly associated with biomass production, only at 50 mM NaCl. Additionally, two SNPs were significantly associated with salt tolerance index as an indicator for growth response of each accession to salt stress. Our results offer candidate genetic resources and SNP markers for breeding salt-tolerant sorghum.

miRNAs control HAM1 functions at the single-cell-layer level and are essential for normal embryogenesis in Arabidopsis.

KEY MESSAGE: miR171a controls HAM1 functions within the protodermal cells of the embryo, and these controls are essential for normal embryogenesis in Arabidopsis. Arabidopsis thaliana miR171a is known to bind to and cleave mRNAs of three HAIRY MERISTEM (HAM) genes that encode members of the GRAS family transcriptional regulators. The molecular functions of the HAM genes are still being elucidated in Arabidopsis. However, detailed expression patterns of miR171a and the effects of the failure of miR171a to suppress HAM genes were unknown till now. Here, we show the detailed expression patterns of miR171a and HAM1 using green fluorescent protein and confocal scanning microscopy. Our observations revealed that miR171a was expressed in the surface cell layer of the embryo and shoot apical meristem, and it controlled HAM1 functions. To determine the impact of the failure of miR171a to suppress of HAM1, we introduced seven synonymous mutations into the miR171a target site of the HAM1 gene (modified HAM1, mHAM1) and generated transgenic plants that had mHAM1 driven by HAM1 native promoter. The mHAM1 transgenic plants showed organogenic defects. These results indicate that the control of HAM1 functions at the single-cell-layer level by miR171a is essential for proper organ formation in Arabidopsis.

RCN1/OsABCG5, an ATP-binding cassette (ABC) transporter, is required for hypodermal suberization of roots in rice (Oryza sativa).

Suberin is a complex polymer composed of aliphatic and phenolic compounds. It is a constituent of apoplastic plant interfaces. In many plant species, including rice (Oryza sativa), the hypodermis in the outer part of roots forms a suberized cell wall (the Casparian strip and/or suberin lamellae), which inhibits the flow of water and ions and protects against pathogens. To date, there is no genetic evidence that suberin forms an apoplastic transport barrier in the hypodermis. We discovered that a rice reduced culm number1 (rcn1) mutant could not develop roots longer than 100 mm in waterlogged soil. The mutated gene encoded an ATP-binding cassette (ABC) transporter named RCN1/OsABCG5. RCN1/OsABCG5 gene expression in the wild type was increased in most hypodermal and some endodermal roots cells under stagnant deoxygenated conditions. A GFP-RCN1/OsABCG5 fusion protein localized at the plasma membrane of the wild type. Under stagnant deoxygenated conditions, well suberized hypodermis developed in wild types but not in rcn1 mutants. Under stagnant deoxygenated conditions, apoplastic tracers (periodic acid and berberine) were blocked at the hypodermis in the wild type but not in rcn1, indicating that the apoplastic barrier in the mutant was impaired. The amount of the major aliphatic suberin monomers originating from C(28) and C(30) fatty acids or ω-OH fatty acids was much lower in rcn1 than in the wild type. These findings suggest that RCN1/OsABCG5 has a role in the suberization of the hypodermis of rice roots, which contributes to formation of the apoplastic barrier.

Arabidopsis dynamin-related proteins DRP2B and DRP1A participate together in clathrin-coated vesicle formation during endocytosis.

Endocytosis performs a wide range of functions in animals and plants. Clathrin-coated vesicle (CCV) formation is an initial step of endocytosis, and in animal cells is largely achieved by dynamins. However, little is known of its molecular mechanisms in plant cells. To identify dynamin-related proteins (DRPs) involved in endocytic CCV formation in plant cells, we compared the behaviors of two structurally different Arabidopsis DRPs, DRP2B and DRP1A, with those of the clathrin light chain (CLC), a marker of CCVs, at the plasma membrane by variable incidence angle fluorescent microscopy (VIAFM). DRP2B shares domain organization with animal dynamins whereas DRP1A is plant-specific. We show that green fluorescent protein (GFP)-tagged DRP2B and DRP1A colocalized with CLC tagged with monomeric Kusabira Orange (mKO) in Arabidopsis cultured cells. Time-lapse VIAFM observations suggested that both GFP-DRP2B and GFP-DRP1A appeared and accumulated on the existing mKO-CLC foci and disappeared at the same time as or immediately after the disappearance of mKO-CLC. Moreover, DRP2B and DRP1A colocalized and assembled/disassembled together at the plasma membrane in Arabidopsis cells. A yeast two-hybrid assay showed that DRP2B and DRP1A interacted with each other. An inhibitor of clathrin-mediated endocytosis, tyrphostin A23, disturbed the localization of DRP1A, but had little effect on the localization of DRP2B, indicating that DRP1A and DRP2B have different molecular properties. These results suggest that DRP2B and DRP1A participate together in endocytic CCV formation in Arabidopsis cells despite the difference of their molecular properties.

Multispectral Phenotyping and Genetic Analyses of Spring Appearance in Greening Plant, Phedimus spp.

The change in appearance during the seasonal transitions in ornamental greening plants is an important characteristic. In particular, the early onset of green leaf color is a desirable trait for a cultivar. In this study, we established a method for phenotyping leaf color change by multispectral imaging and performed genetic analysis based on the phenotypes to clarify the potential of the approach in breeding greening plants. We performed multispectral phenotyping and quantitative trait locus (QTL) analysis of an F1 population derived from 2 parental lines of Phedimus takesimensis, known to be a drought and heat-tolerant rooftop plant species. The imaging was conducted in April of 2019 and 2020 when dormancy breakage occurs and growth extension begins. Principal component analysis of 9 different wavelength values showed a high contribution from the first principal component (PC1), which captured variation in the visible light range. The high interannual correlation in PC1 and in the intensity of visible light indicated that the multispectral phenotyping captured genetic variation in the color of leaves. We also performed restriction site-associated DNA sequencing and obtained the first genetic linkage map of Phedimus spp. QTL analysis revealed 2 QTLs related to early dormancy breakage. Based on the genotypes of the markers underlying these 2 QTLs, the F1 phenotypes with early (late) dormancy break, green (red or brown) leaves, and a high (low) degree of vegetative growth were classified. The results suggest the potential of multispectral phenotyping in the genetic dissection of seasonal leaf color changes in greening plants.

Genetic Dissection of Seasonal Changes in a Greening Plant Based on Time-Series Multispectral Imaging.

Good appearance throughout the year is important for perennial ornamental plants used for rooftop greenery. However, the methods for evaluating appearance throughout the year, such as plant color and growth activity, are not well understood. In this study, evergreen and winter-dormant parents of Phedimus takesimensis and 94 F1 plants were used for multispectral imaging. We took 16 multispectral image measurements from March 2019 to April 2020 and used them to calculate 15 vegetation indices and the area of plant cover. QTL analysis was also performed. Traits such as the area of plant cover and vegetation indices related to biomass were high during spring and summer (growth period), whereas vegetation indices related to anthocyanins were high in winter (dormancy period). According to the PCA, changes in the intensity of light reflected from the plants at different wavelengths over the course of a year were consistent with the changes in plant color and growth activity. Seven QTLs were found to be associated with major seasonal growth changes. This approach, which monitors not only at a single point in time but also over time, can reveal morphological changes during growth, senescence, and dormancy throughout the year.

Novel QTL for Lateral Root Density and Length Improve Phosphorus Uptake in Rice (Oryza sativa L.).

The rice root system consists of two types of lateral roots, indeterminate larger L-types capable of further branching, and determinate, short, unbranched S-types. L-type laterals correspond to the typical lateral roots of cereals whereas S-type laterals are unique to rice. Both types contribute to nutrient and water uptake and genotypic variation for density and length of these laterals could be exploited in rice improvement to enhance adaptations to nutrient and water-limited environments. Our objectives were to determine how best to screen for lateral root density and length and to identify markers linked to genotypic variation for these traits. Using different growing media showed that screening in nutrient solution exposed genotypic variation for S-type and L-type density, but only the lateral roots of soil-grown plants varied for their lengths. A QTL mapping population developed from parents contrasting for lateral root traits was grown in a low-P field, roots were sampled, scanned and density and length of lateral roots measured. One QTL each was detected for L-type density (LDC), S-type density on crown root (SDC), S-type density on L-type (SDL), S-type length on L-type (SLL), and crown root number (RNO). The QTL for LDC on chromosome 5 had a major effect, accounting for 46% of the phenotypic variation. This strong positive effect was confirmed in additional field experiments, showing that lines with the donor parent allele at qLDC5 had 50% higher LDC. Investigating the contribution of lateral root traits to P uptake using stepwise regressions indicated LDC and RNO were most influential, followed by SDL. Simulating effects of root trait differences conferred by the main QTL in a P uptake model confirmed that qLDC5 was most effective in improving P uptake followed by qRNO9 for RNO and qSDL9 for S-type lateral density on L-type laterals. Pyramiding qLDC5 with qRNO9 and qSDL9 would be possible given that trade-offs between traits were not detected. Phenotypic selection for the RNO trait during variety development would be feasible, however, the costs of doing so reliably for lateral root density traits is prohibitive and markers identified here therefore provide the first opportunity to incorporate such traits into a breeding program.

Random regression for modeling soybean plant response to irrigation changes using time-series multispectral data.

Plant response to drought is an important yield-related trait under abiotic stress, but the method for measuring and modeling plant responses in a time series has not been fully established. The objective of this study was to develop a method to measure and model plant response to irrigation changes using time-series multispectral (MS) data. We evaluated 178 soybean (Glycine max (L.) Merr.) accessions under three irrigation treatments at the Arid Land Research Center, Tottori University, Japan in 2019, 2020 and 2021. The irrigation treatments included W5: watering for 5 d followed by no watering 5 d, W10: watering for 10 d followed by no watering 10 d, D10: no watering for 10 d followed by watering 10 d, and D: no watering. To capture the plant responses to irrigation changes, time-series MS data were collected by unmanned aerial vehicle during the irrigation/non-irrigation switch of each irrigation treatment. We built a random regression model (RRM) for each of combination of treatment by year using the time-series MS data. To test the accuracy of the information captured by RRM, we evaluated the coefficient of variation (CV) of fresh shoot weight of all accessions under a total of nine different drought conditions as an indicator of plant's stability under drought stresses. We built a genomic prediction model (MTRRM model) using the genetic random regression coefficients of RRM as secondary traits and evaluated the accuracy of each model for predicting CV. In 2020 and 2021,the mean prediction accuracies of MTRRM models built in the changing irrigation treatments (r = 0.44 and 0.49, respectively) were higher than that in the continuous drought treatment (r = 0.34 and 0.44, respectively) in the same year. When the CV was predicted using the MTRRM model across 2020 and 2021 in the changing irrigation treatment, the mean prediction accuracy (r = 0.46) was 42% higher than that of the simple genomic prediction model (r =0.32). The results suggest that this RRM method using the time-series MS data can effectively capture the genetic variation of plant response to drought.

Distinct gene expression profiles in egg and synergid cells of rice as revealed by cell type-specific microarrays.

Double fertilization in flowering plants refers to a process in which two sperm cells, carried by the pollen tube, fertilize both the egg and the central cell after their release into a synergid cell of the female gametophyte. The molecular processes by which the female gametophytic cells express their unique functions during fertilization are not well understood. Genes expressed in egg and synergid cells might be important for multiple stages of the plant reproductive process. Here, we profiled genome-wide gene expression in egg and synergid cells in rice (Oryza sativa), a model monocot, using a nonenzymatic cell isolation technique. We found that the expression profiles of the egg and synergid cells were already specified at the micropylar end of the female gametophyte during the short developmental period that comprises the three consecutive mitotic nuclear divisions after megaspore generation. In addition, we identified a large number of genes expressed in the rice egg and synergid cells and characterized these genes using Gene Ontology analysis. The analysis suggested that epigenetic and posttranscriptional regulatory mechanisms are involved in the specification and/or maintenance of these cells. Comparisons between the rice profiles and reported Arabidopsis (Arabidopsis thaliana) profiles revealed that genes enriched in the egg/synergid cell of rice were distinct from those in Arabidopsis.

QTL mapping for early root and shoot vigor of upland rice (Oryza sativa L.) under P deficient field conditions in Japan and Madagascar.

Upland rice production is limited by the low phosphorus (P) availability of many highly weathered tropical soils and P deficiency is likely to become increasingly limiting in future drier climates because P mobility decreases sharply with soil moisture. Good seedling root development will be crucial to cope with the combined effects of low P and water availability. Upland rice genebank accession DJ123 was used as a donor for P efficiency and root vigor traits in a cross with inefficient local variety Nerica4 and a set of backcross lines were used to characterize the seedling stage response of upland rice to low P availability and to identify associated QTL in field trials in Japan and Madagascar. Ten QTL were detected for crown root number, root, shoot and total dry weight per plant in a highly P deficient field in Japan using the BC1F3 generation. Of these, qPef9 on chromosome 9 affected multiple traits, increasing root number, root weight and total biomass, whereas a neighboring QTL on chromosome 9 (qPef9-2) increased shoot biomass. Field trials with derived BC1F5 lines in a low-P field in Madagascar confirmed a highly influential region on chromosome 9. However, qPef9-2 appeared more influential than qPef9, as the shoot and root biomass contrast between lines carrying DJ123 or Nerica4 alleles at qPef9-2 was +23.8% and +13.5% compared to +19.2% and +14.4% at qPef9. This advantage increased further during the growing season, leading to 46% higher shoot biomass at the late vegetative stage. Results suggest an introgression between 8.0 and 12.9 Mb on chromosome 9 from P efficient donor DJ123 can improve plant performance under P-limited conditions. The QTL identified here have practical relevance because they were confirmed in the target genetic background of the local variety Nerica4 and can therefore be applied directly to improve its performance.

Corrigendum: QTL mapping for early root and shoot vigor of upland rice (Oryza sativa L.) under P deficient field conditions in Japan and Madagascar.

[This corrects the article DOI: 10.3389/fpls.2022.1017419.].