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.

Transcriptome Analysis of Rice Root Tips Reveals Auxin, Gibberellin and Ethylene Signaling Underlying Nutritropism.

Nutritropism is a positive tropism toward nutrients in plant roots. An NH4+ gradient is a nutritropic stimulus in rice (Oryza sativa L.). When rice roots are exposed to an NH4+ gradient generated around nutrient sources, root tips bend toward and coil around the sources. The molecular mechanisms are largely unknown. Here, we analyzed the transcriptomes of the inside and outside of bending root tips exhibiting nutritropism to reveal nutritropic signal transduction. Tissues facing the nutrient sources (inside) and away (outside) were separately collected by laser microdissection. Principal component analysis revealed distinct transcriptome patterns between the two tissues. Annotations of 153 differentially expressed genes implied that auxin, gibberellin and ethylene signaling were activated differentially between the sides of the root tips under nutritropism. Exogenous application of transport and/or biosynthesis inhibitors of these phytohormones largely inhibited the nutritropism. Thus, signaling and de novo biosynthesis of the three phytohormones are necessary for nutritropism. Expression patterns of IAA genes implied that auxins accumulated more in the inside tissues, meaning that ammonium stimulus is transduced to auxin signaling in nutritropism similar to gravity stimulus in gravitropism. SAUR and expansin genes, which are known to control cell wall modification and to promote cell elongation in shoot gravitropism, were highly expressed in the inside tissues rather than the outside tissues, and our transcriptome data are unexplainable for differential elongation in root nutritropism.

Single-cell transcriptomic analysis of pea shoot development and cell-type-specific responses to boron deficiency.

Understanding how nutrient stress impacts plant growth is fundamentally important to the development of approaches to improve crop production under nutrient limitation. Here we applied single-cell RNA sequencing to shoot apices of Pisum sativum grown under boron (B) deficiency. We identified up to 15 cell clusters based on the clustering of gene expression profiles and verified cell identity with cell-type-specific marker gene expression. Different cell types responded differently to B deficiency. Specifically, the expression of photosynthetic genes in mesophyll cells (MCs) was down-regulated by B deficiency, consistent with impaired photosynthetic rate. Furthermore, the down-regulation of stomatal development genes in guard cells, including homologs of MUTE and TOO MANY MOUTHS, correlated with a decrease in stomatal density under B deficiency. We also constructed the developmental trajectory of the shoot apical meristem (SAM) cells and a transcription factor interaction network. The developmental progression of SAM to MC was characterized by up-regulation of genes encoding histones and chromatin assembly and remodeling proteins including homologs of FASCIATA1 (FAS1) and SWITCH DEFECTIVE/SUCROSE NON-FERMENTABLE (SWI/SNF) complex. However, B deficiency suppressed their expression, which helps to explain impaired SAM development under B deficiency. These results represent a major advance over bulk-tissue RNA-seq analysis in which cell-type-specific responses are lost and hence important physiological responses to B deficiency are missed. The reported findings reveal strategies by which plants adapt to B deficiency thus offering breeders a set of specific targets for genetic improvement. The reported approach and resources have potential applications well beyond P. sativum species and could be applied to various legumes to improve their adaptability to multiple nutrient or abiotic stresses.

Age-dependent analysis dissects the stepwise control of auxin-mediated lateral root development in rice.

As root elongation rates are different among each individual root, the distance from the root apices does not always reflect the age of root cells. Thus, methods for correcting variations in elongation rates are needed to accurately evaluate the root developmental process. Here, we show that modeling-based age-dependent analysis is effective for dissecting stepwise lateral root (LR) development in rice (Oryza sativa). First, we measured the increases in LR and LR primordium (LRP) numbers, diameters, and lengths in wild type and an auxin-signaling-defective mutant, which has a faster main (crown) root elongation rate caused by the mutation in the gene encoding AUXIN/INDOLE-3-ACETIC ACID protein 13 (IAA13). The longitudinal patterns of these parameters were fitted by the appropriate models and the age-dependent patterns were identified using the root elongation rates. As a result, we found that LR and LRP numbers and lengths were reduced in iaa13. We also found that the duration of the increases in LR and LRP diameters were prolonged in iaa13. Subsequent age-dependent comparisons with gene expression patterns suggest that AUXIN RESPONSE FACTOR11 (ARF11), the homolog of MONOPTEROS (MP)/ARF5 in Arabidopsis (Arabidopsis thaliana), is involved in the initiation and growth of LR(P). Indeed, the arf11 mutant showed a reduction of LR and LRP numbers and lengths. Our results also suggest that PINOID-dependent rootward-to-shootward shift of auxin flux contributes to the increase in LR and LRP diameters. Together, we propose that modeling-based age-dependent analysis is useful for root developmental studies by enabling accurate evaluation of root traits' expression.

Outer apoplastic barriers in roots: prospects for abiotic stress tolerance.

Floods and droughts are becoming more frequent as a result of climate change and it is imperative to find ways to enhance the resilience of staple crops to abiotic stresses. This is crucial to sustain food production during unfavourable conditions. Here, we analyse the current knowledge about suberised and lignified outer apoplastic barriers, focusing on the functional roles of the barrier to radial O2 loss formed as a response to soil flooding and we discuss whether this trait also provides resilience to multiple abiotic stresses. The barrier is composed of suberin and lignin depositions in the exodermal and/or sclerenchyma cell walls. In addition to the important role during soil flooding, the barrier can also restrict radial water loss, prevent phytotoxin intrusion, salt intrusion and the main components of the barrier can impede invasion of pathogens in the root. However, more research is needed to fully unravel the induction pathway of the outer apoplastic barriers and to address potential trade-offs such as reduced nutrient or water uptake. Nevertheless, we suggest that the outer apoplastic barriers might act as a jack of all trades providing tolerance to multiple abiotic and/or biotic stressors.

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.

Triterpenoids in aerenchymatous phellem contribute to internal root aeration and waterlogging adaptability in soybeans.

Soybeans (Glycine max) develop newly differentiated aerenchymatous phellem (AP) in response to waterlogging stress. AP is formed in the hypocotyl and root, thus contributing to internal aeration and adaptation to waterlogging for several legumes. Extensive accumulation of triterpenoids - lupeol and betulinic acid - has been identified in AP. However, their physiological roles in plants remain unclarified. Lupeol is converted from 2,3-oxidosqualene by lupeol synthase (LUS) and oxidized to betulinic acid. Notably, soybeans have two LUS genes (GmLUS1 and GmLUS2). Functional analysis was performed to reveal the biological and physiological functions of triterpenoids in AP using lus mutants. The AP cells of lus1 mutant lacked triterpenoid accumulation and epicuticular wax. Lupeol and betulinic acid were the major components of epicuticular wax and contributed to tissue hydrophobicity and oxygen transport to the roots. Tissue porosity in AP was lower in the lus1 mutant than in the wild-type, which resulted in reduced oxygen transport to the roots via AP. This reduction in oxygen transport resulted in shallow root systems under waterlogged conditions. Triterpenoid accumulation in AP contributes to effective internal aeration and root development for adaptation to waterlogging, suggesting the significance of triterpenoids in improving waterlogging tolerance.

Asymmetric auxin distribution establishes a contrasting pattern of aerenchyma formation in the nodal roots of Zea nicaraguensis during gravistimulation.

Auxin distribution is essential for determining root developmental patterns. The formation of lateral roots and constitutive aerenchyma, which is a gas space developed through cell death, is regulated by auxin in rice (Oryza sativa). However, it is unclear whether the involvement of auxin in constitutive aerenchyma formation is conserved in other species. In this study, we found that constitutive aerenchyma formation was regulated by auxin in the nodal roots of Zea nicaraguensis, a wild relative of maize (Zea mays ssp. mays) grown naturally on frequently flooded coastal plains. Subsequent gravistimulation (root rotation) experiments showed opposite patterns of aerenchyma and lateral root formation. Lateral root formation on the convex side of rotated roots is known to be stimulated by a transient increase in auxin level in the pericycle. We found that aerenchyma formation was accelerated in the cortex on the concave side of the rotated nodal roots of Z. nicaraguensis. A cortex-specific expression analysis of auxin-responsive genes suggested that the auxin level was higher on the concave side than on the convex side. These results suggest that asymmetric auxin distribution underlies the regulation of aerenchyma and lateral root formation in the nodal roots of Z. nicaraguensis. As aerenchyma reduces the respiratory cost of the roots, constitutive aerenchyma on the concave side of the nodal root may balance resource allocation, thereby contributing to the uptake of water and nutrients by newly formed lateral roots. Our study provides insights into auxin-dependent asymmetric root patterning such as that of gravistimulation and hydropatterning response.

The barrier to radial oxygen loss protects roots against hydrogen sulphide intrusion and its toxic effect.

The root barrier to radial O2 loss (ROL) is a key root trait preventing O2 loss from roots to anoxic soils, thereby enabling root growth into anoxic, flooded soils. We hypothesized that the ROL barrier can also prevent intrusion of hydrogen sulphide (H2 S), a potent phytotoxin in flooded soils. Using H2 S- and O2 -sensitive microsensors, we measured the apparent permeance to H2 S of rice roots, tested whether restricted H2 S intrusion reduced its adverse effects on root respiration, and whether H2 S could induce the formation of a ROL barrier. The ROL barrier reduced apparent permeance to H2 S by almost 99%, greatly restricting H2 S intrusion. The ROL barrier acted as a shield towards H2 S; O2 consumption in roots with a ROL barrier remained unaffected at high H2 S concentration (500 µM), compared to a 67% decline in roots without a barrier. Importantly, low H2 S concentrations induced the formation of a ROL barrier. In conclusion, the ROL barrier plays a key role in protecting against H2 S intrusion, and H2 S can act as an environmental signalling molecule for the induction of the barrier. This study demonstrates the multiple functions of the suberized/lignified outer part of the rice root beyond that of restricting ROL.

Post-embryonic function of GLOBULAR EMBRYO 4 (GLE4)/OsMPK6 in rice development.

In plants, mitogen activated protein kinases (MPKs) are involved in various signaling pathways that lead to biotic and abiotic responses as well as that regulate developmental processes. Among them, MPK6 and its closely related homologue, MPK3, act redundantly and are known to be involved in asymmetric cell divisions of meristemoid mother cells in stomata development and of zygotes in Arabidopsis. Loss-of-function mutants of GLE4/OsMPK6, which is an orthologue of MPK6 in rice, showed a defect in polarity establishment in early stage of embryogenesis. However, because of the embryo lethality of the mutations, the function of GLE4/OsMPK6 in post-embryonic development is not clarified. Here, we report the analysis of post embryonic function of GLE4/OsMPK6 in vegetative stage of rice using regenerated gle4/osmpk6 homozygous plants from tissue culture. The regenerated plants are dwarf and produce multiple shoots with small leaves. These shoots never develop into reproductive stage, instead, proliferate vegetative shoots repeatedly. Leaves of gle4/osmpk6 have small leaf blade at the tip and blade-sheath boundary become obscure. Stomata arrangement is also disturbed in gle4/osmpk6 leaf blade. The shape of shoot apical meristem of gle4/osmpk6 become disorganized. Thus, GLE4/OsMPK6 functions in shoot organization and stomata patterning in the post embryonic development in rice.

Time-series multispectral imaging in soybean for improving biomass and genomic prediction accuracy.

Multispectral (MS) imaging enables the measurement of characteristics important for increasing the prediction accuracy of genotypic and phenotypic values for yield-related traits. In this study, we evaluated the potential application of temporal MS imaging for the prediction of aboveground biomass (AGB) in soybean [Glycine max (L.) Merr.]. Field experiments with 198 accessions of soybean were conducted with four different irrigation levels. Five vegetation indices (VIs) were calculated using MS images from soybean canopies from early vegetative to early reproductive stage. To predict the genotypic values of AGB, VIs at the different growth stages were used as secondary traits in a multitrait genomic prediction. The prediction accuracy of the genotypic values of AGB from MS and genomic data largely outperformed that of the genomic data alone before the flowering stage (90% of accessions did not flower), suggesting that it would be possible to determine cross-combinations based on the predicted genotypic values of AGB. We compared the prediction accuracy of a model using the five VIs and a model using only one VI to predict the phenotypic values of AGB and found that the difference in prediction accuracy decreased over time at all irrigation levels except for the most severe drought. The difference in the most severe drought was not as small as that in the other treatments. Only the prediction accuracy of a model using the five VIs in the most severe droughts gradually increased over time. Therefore, the optimal timing for MS imaging may depend on the irrigation levels.

Tandemly duplicated genes encoding polygalacturonase inhibitors are associated with bruchid (Callosobruchus chinensis) resistance in moth bean (Vigna aconitifolia).

Bruchids are stored-grain insect pests responsible for serious seed loss in legume crops. A previous study using an F2 population (F2OA) derived from a cross between wild moth-bean (Vigna aconitifolia [Jacq.] Maréchal) accession TN67 (resistant) and cultivated moth-bean accession ICPMO056 (susceptible) revealed that resistance to the azuki bean weevil (Callosobruchus chinensis L.) in TN67 was regulated by a single gene located in the major quantitative trait locus-qVacBrc2.1. In this study, qVacBrc2.1 was finely mapped and candidate genes in this locus were identified using F2OA and another large F2 population (F2NB) derived from the cross mentioned previously. In contrast to the previous study, segregation analysis in the F2NB population revealed that resistance against this pest was controlled by two genes. Furthermore, the addition of novel markers to qVacBrc2.1 and reanalysis of the QTL in the F2OA population demonstrated that qVacBrc2.1 constituted two linked QTLs-qVacBrc2.1-A and qVacBrc2.1-B. The presence of qVacBrc2.1-B was verified using the population F2NB. Comparative genomics using three Vigna spp. strongly suggested the presence of two tandemly duplicated genes, VacPGIP1 and VacPGIP2, which encoded polygalacturonase inhibitors (polygalacturonase-inhibiting proteins) as the candidates for conferring resistance, but only VacPGIP1 could be successfully cloned and sequenced. The alignment of VacPGIP1 coding sequences of TN67 and ICPMO056 revealed eight single nucleotide polymorphisms, three of which altered the amino-acid sequence of the predicted domains of polygalacturonase inhibitors in ICPMO056. Overall, these findings indicate that VacPGIP1 and VacPGIP2 regulated C. chinensis resistance in TN67.

Modeling-based age-dependent analysis reveals the net patterns of ethylene-dependent and -independent aerenchyma formation in rice and maize roots.

Plants require oxygen for the functioning of roots, and thus the establishment of a long-distance diffusion path from above-water tissues to the submerged roots is essential to survive flooding. Rice (Oryza sativa) constitutively forms aerenchyma (gas spaces) under aerobic conditions, and induces its formation in response to low-oxygen conditions. Constitutive aerenchyma formation in rice roots is regulated by the phytohormone auxin, whereas ethylene stimulates inducible aerenchyma formation. However, the net patterns of the ethylene-dependent and -independent (auxin-dependent) aerenchyma formation remain unclear. In the present study, we used a modeling approach to determine age-dependent aerenchyma formation in the wild-type rice and reduced culm number 1 mutant, in which ethylene production is reduced, to reveal the net patterns of ethylene-dependent and -independent aerenchyma formation. Subsequent comparison of age-dependent aerenchyma formation between rice and maize roots suggested that more rapid induction of ethylene-dependent aerenchyma formation and more aerenchyma in rice roots are essential to achieve efficient oxygen diffusion under low-oxygen conditions. Moreover, rice roots showed rapid increase in the expression levels of ethylene biosynthesis and responsive genes, suggesting that the local ethylene production at an early time point after root-cell emergence contributes to the rapid induction of the ethylene-dependent aerenchyma formation in rice. DATA AVAILABILITY: All data included in this study are available upon request by contact with the corresponding author.

Genomic Prediction of Green Fraction Dynamics in Soybean Using Unmanned Aerial Vehicles Observations.

With the widespread use of high-throughput phenotyping systems, growth process data are expected to become more easily available. By applying genomic prediction to growth data, it will be possible to predict the growth of untested genotypes. Predicting the growth process will be useful for crop breeding, as variability in the growth process has a significant impact on the management of plant cultivation. However, the integration of growth modeling and genomic prediction has yet to be studied in depth. In this study, we implemented new prediction models to propose a novel growth prediction scheme. Phenotype data of 198 soybean germplasm genotypes were acquired for 3 years in experimental fields in Tottori, Japan. The longitudinal changes in the green fractions were measured using UAV remote sensing. Then, a dynamic model was fitted to the green fraction to extract the dynamic characteristics of the green fraction as five parameters. Using the estimated growth parameters, we developed models for genomic prediction of the growth process and tested whether the inclusion of the dynamic model contributed to better prediction of growth. Our proposed models consist of two steps: first, predicting the parameters of the dynamics model with genomic prediction, and then substituting the predicted values for the parameters of the dynamics model. By evaluating the heritability of the growth parameters, the dynamic model was able to effectively extract genetic diversity in the growth characteristics of the green fraction. In addition, the proposed prediction model showed higher prediction accuracy than conventional genomic prediction models, especially when the future growth of the test population is a prediction target given the observed values in the first half of growth as training data. This indicates that our model was able to successfully combine information from the early growth period with phenotypic data from the training population for prediction. This prediction method could be applied to selection at an early growth stage in crop breeding, and could reduce the cost and time of field trials.

WUSCHEL-related homeobox family genes in rice control lateral root primordium size.

The development of a plastic root system is essential for stable crop production under variable environments. Rice plants have two types of lateral roots (LRs): S-type (short and thin) and L-type (long, thick, and capable of further branching). LR types are determined at the primordium stage, with a larger primordium size in L-types than S-types. Despite the importance of LR types for rice adaptability to variable water conditions, molecular mechanisms underlying the primordium size control of LRs are unknown. Here, we show that two WUSCHEL-related homeobox (WOX) genes have opposing roles in controlling LR primordium (LRP) size in rice. Root tip excision on seminal roots induced L-type LR formation with wider primordia formed from an early developmental stage. QHB/OsWOX5 was isolated as a causative gene of a mutant that is defective in S-type LR formation but produces more L-type LRs than wild-type (WT) plants following root tip excision. A transcriptome analysis revealed that OsWOX10 is highly up-regulated in L-type LRPs. OsWOX10 overexpression in LRPs increased the LR diameter in an expression-dependent manner. Conversely, the mutation in OsWOX10 decreased the L-type LR diameter under mild drought conditions. The qhb mutants had higher OsWOX10 expression than WT after root tip excision. A yeast one-hybrid assay revealed that the transcriptional repressive activity of QHB was lost in qhb mutants. An electrophoresis mobility shift assay revealed that OsWOX10 is a potential target of QHB. These data suggest that QHB represses LR diameter increase, repressing OsWOX10 Our findings could help improve root system plasticity under variable environments.

Effects of irrigation on root growth and development of soybean: A 3-year sandy field experiment.

Increasing the water use efficiency of crops is an important agricultural goal closely related to the root system -the primary plant organ for water and nutrient acquisition. In an attempt to evaluate the response of root growth and development of soybean to water supply levels, 200 genotypes were grown in a sandy field for 3 years under irrigated and non-irrigated conditions, and 14 root traits together with shoot fresh weight and plant height were investigated. Three-way ANOVA revealed a significant effect of treatments and years on growth of plants, accounting for more than 80% of the total variability. The response of roots to irrigation was consistent over the years as most root traits were improved by irrigation. However, the actual values varied between years because the growth of plants was largely affected by the field microclimatic conditions (i.e., temperature, sunshine duration, and precipitation). Therefore, the best linear unbiased prediction values for each trait were calculated using the original data. Principal component analysis showed that most traits contributed to principal component (PC) 1, whereas average diameter, the ratio of thin and medium thickness root length to total root length contributed to PC2. Subsequently, we focused on selecting genotypes that exhibited significant improvements in root traits under irrigation than under non-irrigated conditions using the increment (I-index) and relative increment (RI-index) indices calculated for all traits. Finally, we screened for genotypes with high stability and root growth over the 3 years using the multi-trait selection index (MTSI).Six genotypes namely, GmJMC130, GmWMC178, GmJMC092, GmJMC068, GmWMC075, and GmJMC081 from the top 10% of genotypes scoring MTSI less than the selection threshold of 7.04 and 4.11 under irrigated and non-irrigated conditions, respectively, were selected. The selected genotypes have great potential for breeding cultivars with improved water usage abilities, meeting the goal of water-saving agriculture.

Mechanisms of lysigenous aerenchyma formation under abiotic stress.

Lysigenous aerenchyma is a gas space created by cortical cell death that enables efficient oxygen diffusion within plants and reduces the energy costs associated with root cells. Recent studies have demonstrated the benefits of aerenchyma formation under flooding, drought, and nutrient deficiency, although further studies will be necessary to understand the mechanisms involved.

Radial Oxygen Loss from Plant Roots-Methods.

In flooded soils, an efficient internal aeration system is essential for root growth and plant survival. Roots of many wetland species form barriers to restrict radial O2 loss (ROL) to the rhizosphere. The formation of such barriers greatly enhances longitudinal O2 diffusion from basal parts towards the root tip, and the barrier also impedes the entry of phytotoxic compounds produced in flooded soils into the root. Nevertheless, ROL from roots is an important source of O2 for rhizosphere oxygenation and the oxidation of toxic compounds. In this paper, we review the methodological aspects for the most widely used techniques for the qualitative visualization and quantitative determination of ROL from roots. Detailed methodological approaches, practical set-ups and examples of ROL from roots with or without barriers to ROL are included. This paper provides practical knowledge relevant to several disciplines, including plant-soil interactions, biogeochemistry and eco-physiological aspects of roots and soil biota.

Genetic regulation of root traits for soil flooding tolerance in genus Zea.

Flooding stress caused by excessive precipitation and poor drainage threatens upland crop production and food sustainability, so new upland crop cultivars are needed with greater tolerance to soil flooding (waterlogging). So far, however, there have been no reports of highly flooding-tolerant upland crop cultivars, including maize, because of the lack of flooding-tolerant germplasm and the presence of a large number of traits affecting flooding tolerance. To achieve the goal of breeding flooding-tolerant maize cultivars by overcoming these difficulties, we chose highly flooding-tolerant teosinte germplasm. These flooding-tolerance-related traits were separately assessed by establishing a method for the accurate evaluation of each one, followed by performing quantitative trait locus (QTL) analyses for each trait using maize × teosinte mapping populations, developing introgression lines (ILs) or near-isogenic lines (NILs) containing QTLs and pyramiding useful traits. We have identified QTLs for flooding-tolerance-related root traits, including the capacity to form aerenchyma, formation of radial oxygen loss barriers, tolerance to flooded reducing soil conditions, flooding-induced adventitious root formation and shallow root angle. In addition, we have developed several ILs and NILs with flooding-tolerance-related QTLs and are currently developing pyramided lines. These lines should be valuable for practical maize breeding programs focused on flooding tolerance.