Identification and Validation of Aerobic Adaptation QTLs in Upland Rice

The aerobic adaptation of upland rice is considered as the key genetic difference between upland rice and lowland rice. Genetic dissection of the aerobic adaptation is important as the basis for improving drought tolerance and terrestrial adaptation by using the upland rice. We raised BC1-BC3 introgression lines (ILs) in lowland rice Minghui 63 (MH63) background. The QTLs of yield and yield-related traits were detected based on ILs under the aerobic and lowland environments, and then the yield-related QTLs were identified in a backcrossed inbred population of BC4F5 under aerobic condition. We further verified phenotypes of QTL near-isogenic lines. Finally, three QTLs responsible for increasing yield in aerobic environment were detected by multiple locations and generations, which were designated as qAER1, qAER3, and qAER9 (QTL of aerobic adaptation). The qAER1 and qAER9 were fine-mapped. We found that qAER1 and qAER9 controlled plant height and heading date, respectively; while both of them increased yields simultaneously by suitable plant height and heading date without delay in the aerobic environment. The phenotypic differences between lowland rice and upland rice in the aerobic environment further supported the above results. We pyramided the two QTLs as corresponding molecular modules in the irrigated lowland rice MH63 background, and successfully developed a new upland rice variety named as Zhongkexilu 2. This study will lay the foundation for using aerobic adaptation QTLs in rice breeding programs and for further cloning the key genes involved in aerobic adaptation.


Introduction
Rice (Oryza sativa L.) as a staple food has fed almost half of the world population, but there are many challenges in rice production because of notably global climate change and water deficit. So far, drought stress has become the main limited factor for rice production [1]. To cope with the shortage of water resources, breeders aim to increase or maintain rice yields during drought periods or in arid areas by using natural genetic variation [2]. Rice including upland rice and lowland rice ecotypes has separately adapted to either aerobic or paddy cultivation. The upland rice represents a predominant ecotype grown under aerobic and rain-fed conditions in the mountainous areas of Southwest China, fine mapping. Finally, a new upland rice variety named as Zhongkexilu 2 was bred in a lowland rice background by molecular breeding methods that introduced clear QTLs as a molecular module.

Introgression of Aerobic Adaptation from the Upland Rice to Lowland Rice
To understand the effect of introgression from the upland to lowland rice, the elite improved upland rice cultivar B6144F-MR-6 was used as the donor and the high-yield lowland variety MH63 as the recurrent parent in an advanced backcross selection program from the 2005 winter season at the Sanya Breeding Station, Sanya, Hannan Province, China. We obtained not less than 150 BC 1 F 1 individuals, and 111 of which were backcrossed as BC 2 F 1 family lines. From each BC 2 F 1 family line, we randomly selected one plant to backcross to produce BC 3 F 1 family lines. We harvested BC 1 F 1 as BC 1 F 2 populations and harvested BC 2 F 1 and BC 3 F 1 as BC 2 F 2 and BC 3 F 2 , respectively, according to their family lines. Crosses and backcrosses were carried out at the Sanya breeding station with three generations per year. By the 2006 winter season, one BC 1 F 4 population, 111 BC 2 F 3 and 96 BC 3 F 2 family lines were obtained. These materials were screened in the dry season in the aerobic environment of the Sanya Breeding Station (SY-DS, November to April) and in the rainy season of the Menglian Breeding Station, Menglian, Yunnan Province, China (ML-RS, mid-May to September) starting in 2007-2008. In SY-DS, mobile sprinkler irrigation facilities were used to maintain a humid aerobic environment. There was sufficient rainfall in ML-RS so that there was no drought stress during the whole growth period; all materials were grown and managed according to the local protocol. Heading date, vegetative vigor, and phenotypic acceptability (PAcp) were assessed in each screening according to the Standard Evaluation System for Rice (International Rice Research Institute, 2002). Visual selection was made based mainly upon PAcp. For duplicate-based screening, 33 BC 1 F 8 , 23 BC 2 F 7 , and 6 BC 3 F 6 introgression lines (ILs) were selected from the BC 1 population and BC 2 , BC 3 family lines, respectively, based upon better PAcp compared with the recurrent parent in the aerobic condition.

Performance of the ILs under Lowland and Aerobic Conditions at Sanya and Menglian
To understand the genetic basis of the aerobic adaptation, the above mentioned ILs of different generations were used to investigate the differences in biomass, yield, harvest index, plant height, and heading dates between the aerobic and lowland condition (SY-DS and ML-RS, respectively) in 2009. All experiments were conducted under both aerobic and lowland condition with the same experimental design and analysis method, and each experiment was performed with three replications with the respective recurrent parent as control. For the aerobic treatment, we used direct sowing with 3-4 seeds per hole, and thinned to one seedling at the three-leaf stage. For the lowland treatment, sowing and transplanting single seedlings were done. Each plot was 2 × 0.75 m 2 , with three rows per plot, and planting density 25 × 20 cm between rows and plants. The five plants in the middle of the middle rows were analyzed for biomass yield, harvest index, heading date, and plant height. Analysis of variance by ANOVA was performed for each trait. Significant differences between ILs and MH63 were determined by multiple comparisons of LSD values compared with the recurrent parent.

QTLs Mapping of Traits Based upon the ILs
For QTL mapping, the ILs lines were genotyped using 256 SSR markers that were polymorphic between the two parents, which were uniformly distributed throughout the rice genome. Based on the performance of the ILs under lowland and aerobic conditions, the probabilities of markers linked to each trait were scored by a binomial test. The selection in 2006 to 2008 was made for aerobic performance compared with the lowland-type parent MH63, if the introgression rate of a marker was significantly higher than that of the theoretical prediction from the binomial distribution, the marker may be linked with the corresponding phenotype. In order to avoid the high probability of false QTL from multiple tests, a high significance level of 0.0001 was set [35,36].

Identification of QTLs in the BC 4 F 5 Backcross Inbred Lines
Based on the grain yield QTL mapping using the ILs, BC 3 F 8 introgression line IL-U315 (which carried four grain yield QTLs) was used to raise a BC 4 F 5 population by backcrossing with MH63 in 2010. More than 500 BC 4 F 4 backcross inbred lines (BILs) were obtained by the single-seed descent method in the 2013 winter season. The leaves of each BC 4 F 4 individuals were used for genotypic analysis and the corresponding BC 4 F 5 individual seeds were used for phenotypic evaluation. We randomly selected 238 BC 4 F 5 BILs for field phenotypic evaluation in the aerobic condition of ML-RS. Experiments were laid out in an alpha lattice design with three replications; each plot was 2 × 0.75 m 2 with three rows per plot and planting density 25 × 20 cm. The five plants in the middle of the middle rows were used for biomass, yield, heading date, and plant height measurements. The experiment was conducted under the aerobic condition with direct sowing, normal fertilizer and weed management. A total of 52 introgression SSR markers of IL-U135 were used to genotype the BC 4 F 4 BILs population. Composite interval mapping of QTLs was performed using IciMapping software version 3.2, combined with the phenotypes of BC 4 F 5 .

Planting and Phenotypic Verification of QTLs Near-Isogenic Lines (NILs)
Base on the BC 4 F 5 BILs population, NILs of relevant QTLs were selected with target region and less introgression. Phenotype verification of the NILs and recurrent parent MH63 were conducted under both aerobic and lowland conditions with the same evaluation method of ILs.

Fine-Mapping of Targeted QTLs
The strategy of substitution mapping was used for fine-mapping the targeted QTLs. The targeted single QTL NILs (BC 4 ) backcrossed with the recurrent parent MH63, and BC 5 F 2 populations were cultivated in the lowland field. Homozygous BC 5 F 2 individuals with different recombination lengths of the QTL region were selected according to genotype, and the phenotype was verified in the aerobic condition of the Menglian Breeding Station.

Pyramiding Aerobic Adaptation Molecular Modules for Breeding
BC 4 BILs carried two aerobic adaptation QTLs were chosen for breeding by introduced the QTLs into the lowland rice MH63 background. The multi-site trials were carried out in Hainan and Yunnan Provinces of China from 2016 to 2018 to seek provincial approval of new varieties. In Hainan Province, the multi-site trails were carried out in three sites from 2016 to 2017. Since the test breeding lines was only Zhongkexilu 2 and the control variety B6144F-MR-6, the experimental area of each variety was 0.5 hectare. In Yunnan Province, the trials were carried out in five sites from 2017 to 2018; three breeding lines participated in the trial with the upland rice variety Yunlu 140 as control. The trials were performed for three repetitions with a randomized block, each block was 13.34 m 2 (0.02 mu). All trials were conducted under upland condition. Sowing at each site according to the beginning of the local rainy season, direct sowing were used with 3-4 seeds per hole, and planting density 25 × 20 cm. Water management depended on natural rainfall without irrigation, and normal fertilizer and weed management were manipulated. The yield and related characters were measured. Significant differences among (or between) breeding lines and control were determined.

Performance of the Introgression Lines under Lowland and Rain-Fed Aerobic Condition
We produced the introgression lines (ILs) by crossing and backcrossing the elite upland rice cultivar B6144F-MR-6 with lowland rice cultivar MH63. We tested the performance of three different generations of the ILs in lowland and aerobic environments. For the three experiments, analysis of variance (ANOVA) was performed for each trait (Table 1). In the lowland conditions, the same trend was found at both sites. There were no significant differences among the ILs in yield and biomass, but there were significant differences in harvest index and plant height, except for BC 3 harvest index at the Menglian site. Heading date showed significant differences at the Sanya site, but not at the Menglian site. In the aerobic conditions, there were significant differences in yield, harvest index and plant height at both sites and in heading date at Sanya, but no significant differences were found in biomass, except for BC 2 at Menglian. These data indicated that it was possible to introgress aerobic adaptation gene(s) from upland genotype to lowland genotype via backcross and selection, which improved the aerobic adaptation of lowland genotype. The traits with significant differences were analyzed by multiple comparisons with the recurrent parent MH63. Significant and non-significant phenotypes in each generation ILs were divided into phenotypic classes of bimorph and were used for further QTL analysis (Table 1).

QTL Mapping of Traits Based upon the ILs
Among the ILs that showed significant differences from their recurrent parent MH63, the probabilities of markers linked to particular traits were scored by a binomial test. QTLs were detected for each trait under both aerobic and lowland condition at Sanya and Menglian, respectively ( Table 2). In the aerobic condition, a plant height QTL (qPH1) was detected on chromosome 1; four heading date QTLs (qHD3, qHD4, qHD7, qHD12) were detected on chromosome 3, 4, 7, and 12, respectively. Four grains yield QTLs (qGY1, qGY3, qGY7, qGY12) and harvest index QTLs (qHI1, qHI3, qHI7, qHI12) were detected in the same region, respectively. In the lowland condition, no QTLs for biomass, yield or harvest index were detected since there were no significant differences between the ILs and the parent MH 63; only a few significant ILs were selected. For plant height, the QTLs detected under irrigation were the same as in the aerobic condition, which indicated that plant height showed consistent genetic influence in both lowland and aerobic condition. Similarly, heading date QTLs were almost the same under irrigation compared with the aerobic condition, except qHD4. Thus, QTLs for different traits varied among sites, condition, and generations, but the QTLs for grain yield in the aerobic condition were consistent across sites and generations, and overlapped with plant height and heading date QTLs in the aerobic condition. These results suggest that plant height and heading date were the important traits for aerobic adaptation underlying the increase of yield in the aerobic condition.

Identification of QTLs in the BC 4 F 5 Backcross Inbred Lines
The line IL-U315 (BC 3 F 8 ) carried four grain yield QTLs was used to produce the BC 4 F 5 BILs population (n = 238). Phenotypic variations of the BILs population and the parents were measured for grain yield, plant height, and heading date under aerobic conditions. The descriptive statistics of phenotypic variation were shown in Table 3. Significant phenotypic differences were observed between the parents in this cross; the heading date of IL-U135 was earlier than MH63, while the yield and plant height were higher than MH63. The frequency distributions of the traits in the BILs were plotted; and similar continuous bimodal or multimodal distributions were found for heading date, plant height, and yield ( Figure 1). The phenotypic correlations showed that grain yield was highly significantly and negatively correlated with heading date (correlation coefficient = −0.578 ** ), and highly significantly correlated with plant height (correlation coefficient = 0.45 ** ). Combined 52 introgressions and polymorphic SSR markers between IL-U135 and MH 63, we identified QTLs for grain yield, plant height and heading date in the BC 4 F 5 BILs population. In the RM212-RM543 interval of chromosome 1, QTLs were detected for plant height and grain yield that accounted for 24.14% of plant height phenotypic variations and 18.7% of yield variations. In the RM218-RM232 interval of chromosome 3, QTLs were detected for heading date and grain yield that accounted for 30% and 19.21% of the phenotypic variations, respectively. In the RM444-RM1328 interval of chromosome 9, one QTL that acted on the heading date was detected, accounting for 19.19% of the heading date phenotypic variations; this QTL also acted on the yield, and accounted for 16.71% of the phenotypic variance (Table 3). These results confirmed the two grain yield QTLs on chromosomes 1 and 3 that were detected in the aerobic condition based upon the ILs, and a new QTL was found on chromosome 9. We named them as Aerobic Adaptation QTLs (qAER1, qAER3, and qAER9).  (Table 3). These results confirmed the two grain yield QTLs on chromosomes 1 and 3 that were detected in the aerobic condition based upon the ILs, and a new QTL was found on chromosome 9. We named them as Aerobic Adaptation QTLs (qAER1, qAER3, and qAER9).

NILs Planting and Phenotype Evaluation
Based on the BC4F5 BILs population, near-isogenic lines (NILs) for Aerobic Adaptation QTLs were selected from the target region with introgression of smaller segments. Homozygous single-QTL NILs for qAER1, qAER9, double-QTL pyramided lines for qAER1 + qAER9, qAER3 + qAER9, and triple-QTL pyramided lines qAER1 + qAER3 + qAER9 were screened, but no NILs carried qAER3 alone or qAER1 + qAER3 double-QTL pyramided lines were detected. Phenotypic verification of NILs was conducted under aerobic and lowland conditions. For heading dates, there were no significant difference between all NILs and MH63 except for qAER1 + qAER9 in lowland condition, but all NILs were significantly earlier than MH63. For plant height, the NILs carried qAER1 was significantly higher than that of MH63, qAER9, and qAER3 + qAER9 in both environments. Similarly, Maximum tillers number per plant of the NILs carried qAER1 was significantly lower than MH63 in both conditions. Interestingly, MH63 maximum tillers number was more in aerobic condition than in lowland. The effective panicles number was consistent with maximum tillers number, but all NILs was significantly more than MH63. The pyramided lines qAER1 + qAER3 + qAER9 showed more spikelets per panicle in both conditions. There was no significant difference between near isogenic lines and MH63 for the biomass in lowland condition. The NILs carried qAER1 + qAER3 + qAER9, qAER3 + qAER9, and qAER1 were significantly higher than MH63 in aerobic condition. The grain setting percentage, harvest index, and yield showed consistent performance; there was no significant difference between NILs and MH63 in lowland, but were significantly higher than

NILs Planting and Phenotype Evaluation
Based on the BC 4 F 5 BILs population, near-isogenic lines (NILs) for Aerobic Adaptation QTLs were selected from the target region with introgression of smaller segments. Homozygous single-QTL NILs for qAER1, qAER9, double-QTL pyramided lines for qAER1 + qAER9, qAER3 + qAER9, and triple-QTL pyramided lines qAER1 + qAER3 + qAER9 were screened, but no NILs carried qAER3 alone or qAER1 + qAER3 double-QTL pyramided lines were detected. Phenotypic verification of NILs was conducted under aerobic and lowland conditions. For heading dates, there were no significant difference between all NILs and MH63 except for qAER1 + qAER9 in lowland condition, but all NILs were significantly earlier than MH63. For plant height, the NILs carried qAER1 was significantly higher than that of MH63, qAER9, and qAER3 + qAER9 in both environments. Similarly, Maximum tillers number per plant of the NILs carried qAER1 was significantly lower than MH63 in both conditions. Interestingly, MH63 maximum tillers number was more in aerobic condition than in lowland. The effective panicles number was consistent with maximum tillers number, but all NILs was significantly more than MH63. The pyramided lines qAER1 + qAER3 + qAER9 showed more spikelets per panicle in both conditions. There was no significant difference between near isogenic lines and MH63 for the biomass in lowland condition. The NILs carried qAER1 + qAER3 + qAER9, qAER3 + qAER9, and qAER1 were significantly higher than MH63 in aerobic condition. The grain setting percentage, harvest index, and yield showed consistent performance; there was no significant difference between NILs and MH63 in lowland, but were significantly higher than MH63 in aerobic condition. We compared QTLs NILs in both lowland and aerobic environments and further verified the phenotype of QTLs. The results indicated that these QTLs increased yield by affecting the plant height (appropriate plant type), not delaying the heading dates, and increased the productive tiller percentage, grain setting percentage, and harvest index in aerobic environment (Figure 2).

Fine-Mapping of the Targeted QTLs
We carried out fine mapping for two main effects-QTLs (qAER1, qAER9). A NIL for qAER1 carried a 6 cM integration segment (RM6333-RM12045) was backcrossed with parent MH63 to produce the BC 5 F 2 population. Based on 3210 BC 5 F 2 individuals, nine polymorphic SSR markers were used to investigate individual genotypes. Eight homozygous lines (BC 5 F 3 ) with different length overlapping segments were screened for phenotypic investigation under aerobic condition (Figure 3a). The phenotypes of the overlapping segment lines showed that lines L5, L6, and L7 had significantly higher plant height and grain yield, while the other lines showed the same phenotype as the parent MH63 (Figure 3b,c). By substitution mapping, qAER1 was fine-mapped into a 134-kb region between RM11974 and RM11982 (Figure 3a).
Life 2020, 10, x FOR PEER REVIEW 9 of 17 MH63 in aerobic condition. We compared QTLs NILs in both lowland and aerobic environments and further verified the phenotype of QTLs. The results indicated that these QTLs increased yield by affecting the plant height (appropriate plant type), not delaying the heading dates, and increased the productive tiller percentage, grain setting percentage, and harvest index in aerobic environment ( Figure 2).

Figure 2.
Phenotypic comparison of the NILs with MH63 under aerobic and lowland condition respectively. The asterisk * and ** indicate significant difference from MH63 at p < 0.05 and p < 0.01, respectively.

Fine-Mapping of the Targeted QTLs
We carried out fine mapping for two main effects-QTLs (qAER1, qAER9). A NIL for qAER1 carried a 6 cM integration segment (RM6333-RM12045) was backcrossed with parent MH63 to produce the BC5F2 population. Based on 3210 BC5F2 individuals, nine polymorphic SSR markers were used to investigate individual genotypes. Eight homozygous lines (BC5F3) with different length overlapping segments were screened for phenotypic investigation under aerobic condition ( Figure  3a). The phenotypes of the overlapping segment lines showed that lines L5, L6, and L7 had significantly higher plant height and grain yield, while the other lines showed the same phenotype as the parent MH63 (Figure 3b,c). By substitution mapping, qAER1 was fine-mapped into a 134-kb region between RM11974 and RM11982 (Figure 3a). Using the same strategy, the qAER9 NIL was backcrossed with MH63 to produce a BC 5 F 2 population. Based on 940 BC 5 F 2 individuals, eight polymorphic SSR markers were used to investigate individual genotypes in the qAER9 interval (RM444-RM1328, 13.1 cM), and five BC 5 F 3 homozygous lines with overlapping segments covered the entire qAER9 interval were screened (Figure 4a). Phenotypic analysis of all the overlapping lines grown in the aerobic condition showed that the heading date of overlapping lines R3, R4, and R5 were earlier than MH63; the effective panicles, grain setting percentage and grain yield were higher than MH63, while the phenotypes of other overlapping lines were the same as that of MH63 (Figure 4b-e). The three overlapping lines (R3, R4, R5) all carried the introgression marker RM105. We narrowed down and fine-mapped the interval to flanking markers RM5526-RM23966 with a length of 3.0 cM. Further investigations of the phenotype of qAER9 showed that it resulted in no delay of heading date, but increased the effective panicles, setting percentage and grain yield.

Pyramiding the Aerobic Adaptation Molecular Modules for Breeding
To verify the usefulness of qAER1 and qAER9 in upland rice breeding, BC 4 QTL pyramided line U272, which carried two aerobic adaptation QTLs (qAER1 + qAER9), was chosen for breeding. The multi-site trials were carried out from 2016 to 2018 in Hainan and Yunnan Province of China under the provincial approval system for new upland varieties. The results showed that the yield in Hainan and Yunnan was significantly higher than that of the control (Figure 5a,b). We named the new variety as Zhongkexilu 2, which has been officially approved in Hainan (Approval Number: Qiongshendao 2018033), and the approval in Yunnan is pending. Using the same strategy, the qAER9 NIL was backcrossed with MH63 to produce a BC5F2 population. Based on 940 BC5F2 individuals, eight polymorphic SSR markers were used to investigate individual genotypes in the qAER9 interval (RM444-RM1328, 13.1 cM), and five BC5F3 homozygous lines with overlapping segments covered the entire qAER9 interval were screened (Figure 4a). Phenotypic analysis of all the overlapping lines grown in the aerobic condition showed that the heading date of overlapping lines R3, R4, and R5 were earlier than MH63; the effective panicles, grain setting percentage and grain yield were higher than MH63, while the phenotypes of other overlapping lines were the same as that of MH63 (Figure 4b-e). The three overlapping lines (R3, R4, R5) all carried the introgression marker RM105. We narrowed down and fine-mapped the interval to flanking markers RM5526-RM23966 with a length of 3.0 cM. Further investigations of the phenotype of qAER9 showed that it resulted in no delay of heading date, but increased the effective panicles, setting percentage and grain yield.

Pyramiding the Aerobic Adaptation Molecular Modules for Breeding
To verify the usefulness of qAER1 and qAER9 in upland rice breeding, BC4 QTL pyramided line U272, which carried two aerobic adaptation QTLs (qAER1 + qAER9), was chosen for breeding. The multi-site trials were carried out from 2016 to 2018 in Hainan and Yunnan Province of China under the provincial approval system for new upland varieties. The results showed that the yield in Hainan and Yunnan was significantly higher than that of the control (Figure 5a,b). We named the new variety as Zhongkexilu 2, which has been officially approved in Hainan (Approval Number: Qiongshendao 2018033), and the approval in Yunnan is pending. Figure 5. Yield comparison of multi-site trials of the upland rice. (a) Hainan Province, Students' t-test was used to assess the differences of yield and determine significant, * P < 0.05, ** P < 0.01; (b) Yunnan Province, significant differences among breeding lines and control were determined by

Pyramiding the Aerobic Adaptation Molecular Modules for Breeding
To verify the usefulness of qAER1 and qAER9 in upland rice breeding, BC4 QTL pyramided line U272, which carried two aerobic adaptation QTLs (qAER1 + qAER9), was chosen for breeding. The multi-site trials were carried out from 2016 to 2018 in Hainan and Yunnan Province of China under the provincial approval system for new upland varieties. The results showed that the yield in Hainan and Yunnan was significantly higher than that of the control (Figure 5a,b). We named the new variety as Zhongkexilu 2, which has been officially approved in Hainan (Approval Number: Qiongshendao 2018033), and the approval in Yunnan is pending. Figure 5. Yield comparison of multi-site trials of the upland rice. (a) Hainan Province, Students' t-test was used to assess the differences of yield and determine significant, * P < 0.05, ** P < 0.01; (b) Yunnan Province, significant differences among breeding lines and control were determined by

Aerobic Adaptation of Upland Rice
During the origin and domestication of rice, human selection, and adaptation to diverse environments have resulted in different ecotypes. Rice is grown under diverse growing conditions such as lowland, rain-fed lowland, rain-fed upland, and flood-prone ecosystems [37]. Thus, the lowland rice and upland rice domesticated as the two major ecotypes of the Asian cultivated rice. In the process of domestication from aquatic (paludal) wild rice or paddy rice to terrestrial upland rice, the growth and development process should be consistent with the soil water supply to avoid drought, and ensure that the water-sensitive reproductive stage is consistent with the rainy season peak through photoperiod control. History and present application areas of upland rice in the Southwest China, Southeast Asia, and South Asia are mainly distributed in the moist upland area. The growth period coincides with the rainy season of the monsoon climate. Drought tolerance may not be the first driving force of upland rice domestication, but may be a second force. Therefore, we suggest that the aerobic adaptation of upland rice is the key genetic difference between the upland rice and lowland rice. Other studies have put forward similar views [11].

Phenotype Differences between the Lowland Rice and Upland Rice in the Aerobic Environment
How to define the phenotype variations adaptation to the aerobic environment between the upland rice and lowland rice is the key to dissecting the aerobic adaptation of upland rice. Those two ecotypes of rice have high genetic variations in morphological and physiological traits [11,14]. However, most studies focus on the difference of drought resistance (tolerance) between lowland and upland rice, and try to design the selection criteria of drought resistance by understanding these drought-related traits [38][39][40], or use marker-assisted selection (MAS) to improve breeding efficiency [41]. In fact, whether using drought-related traits as indicators or using MAS to improve rice drought resistance in breeding practice is not clearly clarified. As the relationship between drought-related traits and yield is not clear, their effects on rice yield under the aerobic environment still limit [42]. Moreover, the expression of major yield QTLs found under drought stress was not consistent in different genetic backgrounds and environments [43]. The aerobic rice breeding progress is still jogged [34]. Comparison the phenotypic differences of a series of typical upland rice and lowland rice verities in aerobic and no-aerobic environments to avoid the complex drought stress, we found that dwarfing plant height and delayed growth period of paddy rice varieties were intuitive phenotypes in the aerobic environment, while the phenotypic differences of typical upland rice varieties in both environments were invisible. Similar results have also been highlighted by some studies [6,8]. In our research, we confirmed two aerobic adaptation QTLs by multiple locations and generations in the aerobic environments without water stress, respectively, acting on plant height and heading date and increasing yield simultaneously in the aerobic environment. These results were consistent with the observed phenotypes.

Does the Green Revolution Gene SD1 Play an Important Role in the Aerobic Adaptation of Upland Rice?
Interestingly, the qAER1 region contains the so-called Green Revolution dwarfing gene SD1. Previous study has comprehensively analyzed co-located QTLs for yield and other traits under drought stress [44]. The region of SD1 was considered to be the key to controlling drought resistance or aerobic adaptation of upland rice [45][46][47][48]. The Green Revolution involved the widespread adoption of semi-dwarf rice varieties carried sd1, and produced more than doubled the global rice yields [49,50]. The high yield potential of modern semi-dwarf varieties was mainly attributed to their improved harvest index, lodging resistance and response to high inputs (mainly nitrogen and water) [51,52]. But most modern semi-dwarf (sd1) rice cultivars performed poorly in the rain-fed ecosystems of Asia and Africa [6,53]. Recent studies have shown that SD1 haplotypes, unique to deep-water rice, promoted internode elongation and the plants adaptation to flooding conditions when submerged. Evolutionary analysis showed that the deep-water rice-specific haplotype was derived from standing variation in wild rice and selected for deep-water rice cultivation [54]. Study of japonica rice domestication also proved that the SD1 locus has been a target of human-mediated selection since prehistoric times [55]. The Green Revolution gene sd1 related to the plant height of rice; meanwhile, the yield and other traits were also affected [6]. Previous studies have shown that there were abundant genomic variations in SD1 [56,57]. During the long-term domestication of different eco-types of rice, selection of different haplotypes of SD1 led to different ecological adaptability. For example, "Green Revolution" haplotypes adapted to a high-input paddy field environment, and deep-water rice haplotypes adapted to the rapid elongation of internodes in a flooded environment. Building on the fine mapping of qAER1, we will further analyze the functions of SD1 haplotypes and other candidate genes in the qAER1 region of upland rice in the aerobic condition.

The Yield of Plants Largely Depends on the Successful Transition from Vegetative Growth to Reproductive Growth
Flowering is not only an indispensable part of reproduction but also a key stage of development that is vulnerable to environmental stress. However, it is becoming increasingly clear that changing flowering time is an evolutionary strategy for maximizing plant reproduction under different pressures [58]. The aerobic environment of upland rice has more complex environmental stress than that of paddy rice, notably a more erratic water supply. The typical upland rice varieties adapt to the aerobic environment by following mechanisms: Dehydration avoidance characterized by significantly higher growth rate and biomass, efficient partitioning characterized by improving harvest index, and drought escape by accelerating heading [40]. In contrast, most paddy rice varieties show poor adaptation to the aerobic environment, with only some spikelets flowering and fruiting, which results in a large number of ineffective tillers and reduces seed setting rate, in turn leading to reduced yield [8]. We identified qAER9 on chromosome 9, which acted on orderly heading in the aerobic environment, and improved effective panicle and seed setting rate, ultimately increased the yield. Recent studies have shown that the OsWOX13 gene enhanced rice drought resistance and induced early flowering, and might be a regulator of other rice drought resistance genes [59]. Further study still need to be conducted for elucidation the molecular function and mechanism of qAER9.
Our study detected and validated two aerobic adaptation QTLs in upland rice, and carried out fine mapping. Aerobic adaptation is a complex agronomic trait controlled by multiple genes. The gene regulatory networks for complex traits often present the characteristics of modularization, in which gene expression is regulated by other major genes and their interactions. Combining into a functional unit, such regulatory genes are responsible for the development of related functions and the formation of target traits as a whole [60,61]. Taking the work forward into breeding practice, we pyramided two QTLs as corresponding molecular modules and introduced them into a lowland rice background. In BC 4 , we observed that the pyramided lines had similar yield potential to the lowland rice parents in the upland environment. By submitting to the Hainan and Yunnan provincial upland rice certification procedures, the stable pyramided line became a variety named as Zhongkexilu 2 that has been approved through Hainan Province. Trials in Yunnan Province also showed excellent field performance. To verify the effects of the QTLs in different genetic backgrounds, we have also screened rice varieties mainly used in Myanmar, Bangladesh, and Nepal as the recipients. We have used molecular marker-assisted selection (MAS) to confirm the introduction of the molecular modules into different genetic backgrounds.