Next Article in Journal
The Seed Development Factors TT2 and MYB5 Regulate Heat Stress Response in Arabidopsis
Next Article in Special Issue
Low Temperature Antioxidant Activity QTL Associate with Genomic Regions Involved in Physiological Cold Stress Tolerance Responses in Rice (Oryza sativa L.)
Previous Article in Journal
Whole-Exome Sequencing, Proteome Landscape, and Immune Cell Migration Patterns in a Clinical Context of Menkes Disease
Previous Article in Special Issue
Diverse Roles of MAX1 Homologues in Rice
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Novel QTL Associated with Shoot Branching Identified in Doubled Haploid Rice (Oryza sativa L.) under Low Nitrogen Cultivation

Department of Southern Area Crop Science, National Institute of Crop Science, RDA, Miryang 50424, Korea
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work as first authors.
Submission received: 12 April 2021 / Revised: 5 May 2021 / Accepted: 12 May 2021 / Published: 14 May 2021
(This article belongs to the Special Issue Genetic Diversity in Rice Cultivars)

Abstract

:
Shoot branching is considered as an important trait for the architecture of plants and contributes to their growth and productivity. In cereal crops, such as rice, shoot branching is controlled by many factors, including phytohormones signaling networks, operating either in synergy or antagonizing each other. In rice, shoot branching indicates the ability to produce more tillers that are essential for achieving high productivity and yield potential. In the present study, we evaluated the growth and development, and yield components of a doubled haploid population derived from a cross between 93-11 (P1, indica) and Milyang352 (P2, japonica), grown under normal nitrogen and low nitrogen cultivation open field conditions. The results of the phenotypic evaluation indicated that parental lines 93-11 (P1, a high tillering indica cultivar) and Milyang352 (P2, a low tillering japonica cultivar) showed distinctive phenotypic responses, also reflected in their derived population. In addition, the linkage mapping and quantitative trait locus (QTL) analysis detected three QTLs associated with tiller number on chromosome 2 (qTNN2-1, 130 cM, logarithm of the odds (LOD) 4.14, PVE 14.5%; and qTNL2-1, 134 cM, LOD: 6.05, PVE: 20.5%) and chromosome 4 (qTN4-1, 134 cM, LOD 3.92, PVE 14.5%), with qTNL2-1 having the highest phenotypic variation explained, and the only QTL associated with tiller number under low nitrogen cultivation conditions, using Kompetitive Allele-Specific PCR (KASP) and Fluidigm markers. The additive effect (1.81) of qTNL2-1 indicates that the allele from 93-11 (P1) contributed to the observed phenotypic variation for tiller number under low nitrogen cultivation. The breakthrough is that the majority of the candidate genes harbored by the QTLs qTNL2-1 and qTNN4-1 (here associated with the control of shoot branching under low and normal nitrogen cultivation, respectively), were also proposed to be involved in plant stress signaling or response mechanisms, with regard to their annotations and previous reports. Therefore, put together, these results would suggest that a possible crosstalk exists between the control of plant growth and development and the stress response in rice.

1. Introduction

The global population is increasing at a relatively high growth rate per annum, and is projected to reach about 9.8 billion people by 2050 [1]. The current status of the food and nutrition security in the world question on how to provide sufficient, healthy and nutritious food, at all times to the growing population. To achieve this, several initiatives have been taken to combat malnutrition and food insecurity, and alleviate poverty and eradicate hunger, in line with the United Nations Sustainable Development Goals (SDGs), through research and development, among others.
Since the domestication of rice, this Poaceae is the second most important cereal crop after corn and remains today the only cereal crop solely cultivated for human consumption, in addition to being the staple food crop for more than half of the world’s population, despite the increasing preference for various foods and the diversification of diet. According to the Food and Agriculture Organization of the United Nations (FAO) statistics (http://www.fao.org/faostat/en/#data/QC/visualize, accessed on 5 March 2021), nearly 755.5 million metric tons of paddy rice were harvested in 2019, and around 162.1 million hectares were occupied by rice (upland, wetland, and irrigated) fields across the globe.
To achieve optimum plant growth and development, and realize their potential in terms of productivity and quality, plant crops require nitrogen (N), among other essential macronutrients [2]. The role of nitrogen in the metabolism of plants is important in promoting seed germination, cell division, plant growth and development, formation of amino acids and proteins, and plant productivity [3,4,5,6]. Rice cultivation requires abundant use of nitrogenous fertilizers [7], particularly in high yielding varieties [8]. However, excessive application of nitrogenous fertilizers has been associated with the emission of greenhouse gases (GHG), such as methane (CH4) and nitrous oxide (N2O) [9,10,11,12].
In the recent years, the trend in plant breeding revealed a growing interest in developing rice varieties with an optimized and balanced productivity and quality, while exhibiting a high nitrogen use efficiency (NUE), which may help reduce the use of large amounts of fertilizers, particularly the N-rich fertilizers or a combined mineral–organic fertilization regime [13,14,15,16], as part the efforts to reducing GHG emission from agricultural activities, and to promoting the use of ecofriendly agricultural practices. To achieve this, scientists have employed various strategies that include the use of advanced plant breeding techniques and the available and emerging technologies to develop rice varieties with an improved nitrogen use efficiency (NUE), which includes nitrogen uptake, transport, assimilation, and remobilization. It has been shown that an increase in nitrogen (N) fertilizers is associated with a large decrease in NUE [17]. Therefore, there is a pressing need to improve agricultural practices in order to contribute to the reduction of GHG emissions from irrigated, rainfed, and flood-prone rice ecosystems, to predict the future emissions and develop rice technologies that will reduce GHG emissions, and yet will meet the required increase in rice production [18].
In rice, shoot branching is considered as one of the major agronomic traits with a high potential to contribute to the yield and productivity of rice [19,20]. It is well established that phytohormones play fundamental roles in the control of shoot branching and bud outgrowth in plants, operating in a synergetic or antagonistic manner [21,22,23,24,25], in addition to their roles in seed dormancy and germination, plant growth and development, flowering and organogenesis, seed formation and maturation, fruit ripening, and signaling during abiotic or biotic stress occurrence [26]. Several studies have proposed a number of quantitative trait loci (QTLs) controlling tiller number in rice, and mapped to almost all rice chromosomes (10–11) using different sorts of mapping populations, including doubled haploid (DH) lines [27]. Many genes have been identified; however, only a few genes have been cloned and functionally characterized [28,29,30,31,32,33,34,35,36,37]. A mapping population derived from a cross between indica and japonica would offer the best cross combination to investigate QTLs associated with rice plant architecture, due to the morphological difference between indica and japonica subspecies.
The use of fixed homozygous populations, such as doubled haploid (DH) [38,39] and recombinant inbred lines (RILs) [40] have been favored when investigating QTLs for various agronomic traits in rice. The indica and japonica rice subspecies, are the most cultivated rice, and are said to be independently domesticated from the wild ancestor [41,42,43].
Therefore, from this perspective, this study performed a linkage mapping analysis on a doubled haploid (DH) developed through anther culture of a cross between 93-11 (P1) and Milyang352 (P2), typical indica and japonica rice cultivars, respectively, grown under normal nitrogen application and low nitrogen-level field conditions, to investigate putative QTLs controlling rice tiller number as well as agronomic traits contributing to the yield of rice using Kompetitive Allele-Specific PCR (KASP) and Fluidigm markers. Additionally, this study identified candidate genes associated with tiller number in rice under low nitrogen cultivation, with interesting annotated molecular functions, and suggested to have the potential to regulate shoot branching and panicle number in rice.

2. Materials and Methods

2.1. Mapping Population, Growth Conditions, and Phenotypic Measurements

A total of 117 doubled haploid (DH) rice lines developed through another culture of the cross between Oryza sativa L. ssp. indica cv. 93-11 (P1) and Oryza sativa L. ssp. japonica cv. Milyang352 (P2), respectively, were used as the mapping population. Prior to sowing, seeds were soaked in 0.7% nitric acid (HNO3) CAS: 7697-37-2, Lot No. 2016B3902; Junsei Chemical Co. Ltd., Tokyo, Japan) for about 24 h to break the dormancy, followed by incubation at 32 °C for 48 h to induce germination. Then, germinated seeds were sown in 50-well trays containing an enriched soil (Doobaena Plus, Nong Kyung Ltd., Yeongcheon-si, Korea) until transplanting time. Seedlings with uniform height were transplanted in the experimental field located at the National Institute of Crop Science (NICS) (altitude: 11 m, 35°29′31.4″ N and 128°44′30.0″ E), Rural Development Administration (RDA), Miryang, Republic of Korea. The experiments were initially conducted during the rice cropping season (June–November) of 2019 and repeated in 2020 to validate the phenotypes prior to performing QTL analysis. One set of the mapping population and parental lines were transplanted on a plot with a normal nitrogen cultivation regime (90 kg ha−1 N, 45 kg ha−1 P2O5, and 57 kg ha−1 K2O), and the other on a plot with a low nitrogen level (45 kg ha−1 N, 22.5 kg ha−1 P2O5, and 28.5 kg ha−1 K2O) maintained over 15 years. Each rice line had at total of 100 seedlings transplanted in four rows, with 25 plants per row and the spacing between and within the lines of 30 cm × 15 cm, respectively. Parental lines were planted after every 20 DH lines, starting from the initial rows. The weather parameters and conditions during the experimental period can be found in Figure S1.
The soil chemical analysis was performed following the method described by Park et al. [44] and the basic soil chemical properties are given in Table 1.
The phenotypic measurements of the studied traits, which included plant height (PH), tiller number (TN), culm length (CL), panicle length (PL), number of panicles per plant (PN), root length (RL), shoot dry weight (SDW), and root dry weight (RDW), were done on ten individual plants in triplicate randomly pooled from the inside rows, excluding the border rows to avoid the border effects on the traits studied and competition between lines. The plant height and culm length were measured using a 1-m long wooden stick ruler at maximum tillering stage. The panicle length was measured with a regular ruler soon after harvesting. The number of tillers per plant was counted at maximum tillering stage.
SDW was measured from an oven-dried (at about 45 °C overnight) separated shoot from the roots at the collar, cooled down, and weighed on a precision scale (CASKOREA Co. Ltd, MW2-6004, 0.02–60 g, Seongnam-si, Korea). Similarly, RDW was measured after washing off any loose soil particles, followed by blotting of samples with a soft paper towel to remove any free surface moisture. Samples were oven-dried, weighed, and the results expressed in grams.

2.2. Frequency Distribution, Quantile–Quantile Plots, Principal Component Analysis, and Correlation Analysis

The frequency distributions of each phenotypic trait under normal and low nitrogen cultivation and the correlation analysis were performed using GraphPad Prism software (version 7.00, © 1992–2016 GraphPad, San Diego, CA, USA). The normality of the distribution was assessed using the Shapiro–Wilk W-statistic for the test of normality in IciMapping (version 4.1.0.0, 2016, Chinese Academy of Agriculture Sciences, Beijing, China). In addition, the Quantile–Quantile (Q–Q) plots and the pairwise kinship matrix were generated using the GAPIT function: my_GAPIT <-GAPIT(Y = myY, G = myG, Model.selection = TRUE, SNP.MAF = 0.05) in RStudio (version 1.2.5042, © 2009–2020 RStudio, Inc., Boston, MA, USA). The cluster analysis was performed using the R packages ggplot2, tidyverse, cluster, and fviz_cluster R function. The principal component analysis (PCA) was performed using the R package ggplot2 and functions myPr, plot, and biplot [45,46].

2.3. Genotyping and Molecular Markers Analysis

For the molecular markers analysis, Kompetitive Allele-Specific PCR (KASP) markers and Fluidigm markers were employed. KASP markers were amplified from genomic DNA and the allelic discrimination was performed using the Nexar system (LGC, Biosearch Technologies, Hoddesdon, UK), while the Fluidigm markers for SNP genotypes were determined using the BioMarkTM DH system (Fluidigm, San Francisco, CA, USA) and 96.96 Dynamic Array IFC (96.96 IFC) chip as described earlier [47].

2.4. Construction of Linkage Maps and QTL Analysis

A total of 240 markers, including KASP markers earlier reported [48] to be specific for detecting polymorphism between japonica subspecies and Fluidigm markers [49] specific for detecting polymorphism between indica and japonica subspecies, and the phenotypic data of the mapping population consisting of 117 doubled haploid (DH) rice lines and their parental lines (93-11, P1 and Milyang352, P2) were used to perform QTL analysis in order to identify putative QTLs associated with plant growth, tiller number, panicle length, panicle number, root length, and shoot and root dry weights in rice, under low nitrogen cultivation. The genotype and phenotype raw data were initially formatted in RStudio using the fread, pheno.raw, geno.raw, and cbind functions. The QTL analysis was performed and the linkage maps constructed with IciMapping software v.4.1.0.0, for a biparental population using position mapping and Kosambi mapping functions [50]. The permutation (1000 times) parameters, which explains the probability for detecting statistically significant QTLs associated with the target traits, were selected.
The rice genome annotation database publicly available online browser (http://rice.plantbiology.msu.edu/, accessed on 22 February 2021) was used to identify candidate genes pooled from the detected QTL qTNL2-1 (684.4 kbp, Chr2:26706424..27390847) region flanked by KASP markers KJ02_057 (left, closest marker) and id2012773 (right marker) on chromosome 2, associated to the control of the number of tillers under low nitrogen cultivation.

3. Results

3.1. Differential Growth Patterns between Doubled Haploid Lines under Normal and Low Nitrogen Conditions

We evaluated the phenotypic response of the mapping population and their parental lines under normal and low nitrogen cultivation conditions. The results revealed a normal distribution for plant height (Figure 1A,B, Table S3, under normal and low nitrogen cultivation), number of panicles (Figure 1I, Table S3, under normal nitrogen cultivation), and root length (Figure 1K,L, Table S3, under normal and low nitrogen cultivation). However, for culm length (Figure 1F, under low nitrogen cultivation) a negative skewness was observed. In contrast, a positive skewness was recorded for tiller number (Figure 1C,D, Table S3), culm length (Figure 1E, Table S3, under normal nitrogen cultivation), panicle length (Figure 1G,H, Table S3, under normal and low nitrogen cultivation), panicle number (Figure 1J, Table S3, under low nitrogen cultivation), shoot dry weight (Figure 1M,N, Table S3, under both normal and low nitrogen cultivation), and root dry weight (Figure 1O,P, Table S3, under normal and low nitrogen cultivation). Similarly, the Quantile–Quantile (Q–Q) plots (observed quantiles plotted against theoretical quantile or z-scores) reflected the same distribution patterns of the traits observed earlier (Figure S2A–P).
Our data indicate that about 8.5% and 91.5% of the DH lines had relatively short plants (P1-like phenotype) and tall plants (P2-like phenotype), respectively, under normal nitrogen conditions. Whereas, under low nitrogen cultivation, nearly 39.3% and 60.7% had P1-like (relatively short plants) and P2-like (relatively tall plants) plant height, respectively (Figure 2A,I,J). In addition, about 68.4% of the 117 DH lines grown under normal nitrogen conditions recorded an increase in number of tillers per plant similar to P1 (>13 tillers per plant) against 31.6% with P2-like tiller number (≤13 tillers per plant), while about 69.2% and 30.8% of the DH lines cultivated under low nitrogen level had P1-like and P2-like tiller number, respectively (Figure 2B). In the same way, about 83.8% of the mapping population exhibited an increase in culm length (P1-like, >58 cm) against 14.5% that showed a P2-like culm length (≤58 cm) under normal nitrogen cultivation. However, when the DH lines were grown under low nitrogen conditions, nearly 82.1% and 17.9% showed P1-like (>59 cm) and P2-like (≤59 cm) culm length, respectively (Figure 2C). Additionally, of the 117 DH lines cultivated under normal nitrogen level, about 37.6% and 62.4% had long panicles (P1-like) and short panicles (P2-like), respectively (Figure 2D), while under low nitrogen cultivation, 39.3% and 60.7% had long panicles (P1-like) and short panicles (P2-like), respectively. Furthermore, from the set of the mapping population grown under normal range of nitrogen level, we observed that about 93.2% exhibited an increase in number of panicles per plant (P1-like), and 6.8% had P2-like number of panicles. Similarly, when grown under low nitrogen cultivation conditions, 98.3% and 1.7% of the DH lines recorded an increase in number of panicles per plant, P1-like and P2-like, respectively (Figure 2E).
We were also interested to see the changes in the root growth and the biomass dry weight of the mapping population under low nitrogen cultivation compared to the normal range nitrogen supply. Therefore, we measured the root length, shoot dry weight, and root dry weight. The results in Figure 2F indicate that only 3.4% of the DH lines had P1-like root growth pattern (relatively long roots), while nearly 96.6% of the same population grown under normal nitrogen cultivation exhibited a P2-like root growth pattern (relatively short roots). When grown under low nitrogen level, we observed that only 7.7% of the DH lines showed a P1-like root growth pattern against 92.3% with a P2-like root growth pattern (Figure 2F). Regarding the shoot dry weight of the mapping population cultivated under normal nitrogen conditions, our data indicate that about 96.6% and 3.4% had P1-like (high) and P2-like (low) shoot dry weight (SDW), while about 87.2% and 12.8% recorded a high SDW (P1-like) and low SDW (P2-like) under low nitrogen cultivation conditions (Figure 2G). As for root dry weight (RDW), about 81.2% and 14.5% recorded a high (P1-like) and low (P2-like) RDW under normal nitrogen cultivation, respectively; whereas, about 39.3% and 60.7% had a P1-like and P2-like RDW under low nitrogen cultivation (Figure 2H).
In the same perspective, the root–shoot ratio (R:S) results indicate that about 46.2% and 38.5% of the DH lines recorded a high R:S (P1-like, R:S > 0.42 and >0.35 based on length (Figure 2I) and dry weight (Figure 2J), respectively), while 61.5% and 53.8% (P2-like, R:S ≤ 0.42 and ≤0.35 based on length and dry weight, respectively) showed a low R:S (P2-like). Meanwhile, under low nitrogen conditions, about 19.7% and 62.4% (R:S based on length and dry weight, respectively) had a P1-like R:S (>0.35 and >0.37, based on length and dry weight, respectively) against 37.6% and 59% (R:S ≤ 0.37 and ≤0.35, based on length and dry weight, respectively).
In addition, the phenotypic response of the parental lines (Figure 2A–L) indicates that both 93-11 (P1) and Milyang352 (P2) experienced a reduction in almost all the traits when grown under low nitrogen cultivation, compared to the normal nitrogen supply conditions, except a slight increase in culm length recorded in 93-11 (1.1% increase) and Milyang352 (0.4% increase), and a significant increase in number of panicles per plant observed in 93-11 (P1) against 16.6% reduction recorded in Milyang352 (P2) under low nitrogen cultivation.

3.2. Relatedness, Principal Component Analysis, and Correlation between Traits

In Figure 3A, our data show the KASP marker-based kinship matrix, also known as the co-ancestry or half relatedness, revealing the distribution of coefficients of co-ancestry of the mapping population. Additionally, the DH lines were grouped into two distinct clusters with regard to their number of tillers pattern (P1-type or P2-type) under low nitrogen cultivation (Figure 3B). The principal component analysis (PCA) proposes a correlation between traits under normal and low nitrogen cultivation conditions. Principal component 1 (PC1), 2 (PC2), 3 (PC3), and 4 (PC4) explained 27.3%, 19.5%, 15.2%, and 14.1% of the proportion of variance of the observed phenotypes, respectively, of the DH lines grown under normal nitrogen conditions. Similarly, under low nitrogen cultivation, PC1, PC2, PC3, and PC4 explained 27.3%, 19.1%, 15.6%, and 14.4% of the proportion of variance of the observed phenotypes.
Initially, we were interested in investigating a possible relationship between tiller number and other quantitative traits of interest recorded in the mapping population. For this reason, we performed a linear regression analysis. The results shown in Figure 3C,D suggest a strong positive correlation between tiller number and panicle number (R2 = 0.603*** and R2 = 0.560***, under normal and low nitrogen cultivation, respectively). In contrast, plant height showed a negative correlation with tiller number (R2 = −0.465*** and R2 = −0.279***, under normal and low nitrogen cultivation, respectively). However, culm length (R2 = 0.085ns and R2 = 0.198ns), panicle length (R2 = −0.199* and R2 = −0.222ns), root length (R2 = 0.046ns and R2 = 0.003ns), and root dry weight (R2 = 0.173ns and R2 = 0.288**) showed a nonexistent or very low correlation with tiller number (Tables S1 and S2). In the same perspective, root length and root dry weight had a positive but weak correlation under normal nitrogen cultivation (R2 = 0.209*, Table S1), while showing a stronger positive correlation under low nitrogen cultivation (R2 = 0.440***, Table S2). Likewise, shoot dry weight and root dry weight are predicted to have a strong positive correlation. Furthermore, it was interesting to see that shoot dry weight (SDW) of the mapping population showed a strong positive correlation with tiller number (R2 = 0.583***, Figure 3F, Table S2) under low nitrogen cultivation conditions, while under normal nitrogen supply, a weak positive correlation was observed (R2 = 0.317***, Figure 3E, Table S1). It is then thought that a unit increase in the number of tillers would cause an increase in the shoot dry weight of rice by the same measure. Tables S1 and S2 provide details information about the Person r correlation between traits. Data in these two tables also revealed that the correlation level between tiller number and other traits differ depending on each nitrogen condition.

3.3. Quantitative Trait Loci (QTLs) Associated with Plant Growth and Shoot Branching under Normal and Low Nitrogen Cultivation Were Mapped to Different Chromosomes

We used 240 KASP makers and Fluidigm markers, and the phenotype data of 117 DH lines and parental lines, cultivated under normal nitrogen and low nitrogen levels, to perform linkage mapping and QTL analysis. The results revealed that, in total, 28 QTLs (all traits considered) associated with plant growth and shoot branching were mapped to all rice chromosomes, except chromosomes 6 and 11 (Table 2, Figure 4A–J). In essence, five QTLs (qPHN1-1, LOD: 4.7; qPHN3-1, LOD: 4.88; qPHN5-1, LOD: 5.69; qPHN10-1, LOD: 5.3; and qPHN12-1, LOD: 5.13) were associated with plant height of rice plants grown under normal nitrogen cultivation (Table 2, Figure 4A,C,E,I,J). The additive effect recorded by qPHN1 (2.08) and qPHN10 (2.14), and those of the QTLs qPHN3 (−2.1), qPHN5 (−2.2), qPHN12 (−2), indicate that the alleles from both 93-11 and Milyang352 contributed to the observed plant height phenotypic variation under normal nitrogen cultivation.
In the case of tiller number, two QTLs (qTNN2-1, LOD: 4.14; and qTNN4-1, LOD: 3.92) were detected in DH lines cultivated under normal nitrogen conditions, accounted for about 14.5% of the phenotypic variance explained (PVE) (Table 2, Figure 4B,D). However, solely one QTL (qTNL2-1, LOD:6.05) mapped to chromosome 2, flanked by KJ02_057 and id2012773, was associated with tiller number in DH lines grown under low nitrogen cultivation conditions, and accounted for about 20.5% of the PVE (Table 2, Figure 4B). The recorded additive effect (1.81) indicates that the allele from 93-11 (P1, indica) contributed to the observed phenotypic variation for tiller number under low nitrogen cultivation.
In addition, seven QTLs associated with culm length, of which number five were detected in the mapping population cultivated under normal nitrogen conditions (qCLN2-1, LOD: 4.42; qCLN2-2, LOD: 13.45; qCLN3-1, LOD: 11.45; qCLN9-1, LOD:4.56; and qCLN10-1, LOD: 3.18) (Table 2, Figure 4B,C,H,I). The QTL qCLN2-1 (for culm length) and the qTNL2-2 (for tiller number) coincided within the same region and are flanked by the same markers (KJ02_057 and id20112773), which would suggest that qTNL2-2 could be associated with both tiller number and culm length. Furthermore, two QTLs (qCLL1-1, LOD: 3.95; and qCLL9-1, LOD: 4.63) were detected in DH lines grown under low nitrogen level (Table 2, Figure 4A,H). The recorded additive effects indicate that both 93-11 and Milyang352 contributed to the observed culm length phenotypic variation.
In the same way, a set of five QTLs (qPLN2-1, LOD: 4.81; qPLN3-1, LOD: 3.92; qPLN3-2, LOD: 8.84; qPLN7-1, LOD: 5.94; and qPLN8-1, LOD: 8.06) associated with panicle length were only detected in the mapping population cultivated under normal nitrogen conditions (Table 2, Figure 4B,C,F,G). The additive effects recorded for both normal and low nitrogen cultivation conditions show that 93-11 and Milyang352, P1 and P2, respectively, contributed to the observed panicle length phenotypic variation. For panicle number trait, a total of four QTLs were detected, of which number three (qPNN2-1, LOD: 3.49; qPNN4-1, LOD: 9.63; and qPNN8-1, LOD: 4.34) were associated with normal nitrogen cultivation regime (Table 2, Figure 4B,D,G), and one QTL (qPNL4-1, LOD: 2.95) with low nitrogen cultivation (Table 2, Figure 4D). The additive effects indicate that, for both normal and low nitrogen cultivation conditions, the alleles from the indica parental line (93-11, P1) contributed to the observed phenotypic variation for panicle number (Table 3).
A successful establishment of plants is expected when they develop a robust rooting system for efficient nutrients and water acquisition. Our results indicate that a single QTL (qRLL4-1, LOD:3.41) associated with root growth under low nitrogen cultivation was detected, and mapped to chromosome 4, and flanked by KJ04_053 and KJ04_057 (closest marker) KASP markers. The additive effect (−1.1) suggests that the allele from Milyang352 (P2, japonica) contributed to the observed phenotypic variance explained (PVE, 11.98%) (Table 2, Figure 4D). Biomass dry weight is affected by many factors, including the acquisition of mineral nutrients and crop productivity [51]. It was also proposed that there is an interdependent relationship between the shoot and roots [52,53]. The QTL analysis results showed that one QTL associated with shoot dry weight (qSDWN5-1, LOD: 3.74) (Table 2, Figure 4E), and two QTLs (qRDW1-1, LOD: 4.06; and qRDW4-1, LOD: 3.22) (Table 2, Figure 4A,D).
Although many QTLs associated with different traits analyzed were identified, we opted to unveil candidate genes found within qTNL2-1 (684.4 kbp, Chr2:26706424..27390847) essentially due to the fact that this major QTL was the only one associated with the control of shoot branching under low nitrogen cultivation (Table 3), which may help reduce excessive nitrogen-rich fertilizers application in rice cultivation, and contribute to reducing methane (NH4) emissions from paddy fields. An additional list of candidate genes with interesting annotated functions, pooled from the qTNN4-1 region (1.193 Mbp) and associated with tiller number under normal nitrogen cultivation conditions, is attached as Table S4.
Shoot branching has been shown to be regulated at different levels of the plant cellular metabolism. It is well established that phytohormones signaling pathways are key regulators of shoot branching in plants, operating in synergetic or antagonistic manner. Among the identified candidate genes (Table 3), we found two genes encoding GRAS (GIBBERELILIN-INSENSITIVE (GAI), REPRESSOR of GA1-3 (RGA), and SCARECROW (SCR)) transcription factor (TF) (Os02g44360 and Os2g44370). Proteins in the GRAS domain family are known as key players in gibberellin (GA) signaling, which regulates various aspects of plant growth and development, by inhibiting the proliferation and expansion of cell-mediated plant growth [54]. Another gene found in the same region is OsGPX3, a glutathione peroxidase putative expressed earlier reported to be involved shoot and root development in rice and stress response [55], as well as abscisic acid (ABA) response in rice [56]. Similarly, we could also find many other genes with interesting predicted functions, such as OsDhn1 (Os02g44870), earlier suggested to be involved in abiotic stress response in rice [57,58], OsRALF8 (rapid alkalinization factor, Os02g44940) gene, one of which FERONIA (a receptor-like kinase (RLK) member of CrRLK1L, known as a major plant cell growth modulator in distinct tissues) is known to be a receptor, and previously reported to play a critical role in the regulation of plant defense and development (inhibitory or suppression of growth) [59].
Reports indicated that carbohydrates regulate different aspects of plant growth via the modulation of cell division and expansion. For instance, a gene encoding DUF581 domain was reported to be a generic SnRK1 (sucrose nonfermenting-1 (SNF1)-related protein kinase) interaction module and co-expressed with SnRK1 during plant cell signaling [60]. We also identified a DUF encoding gene OsFBDUF13 (F-box-like domain and domain of unknown function (DUF295) containing protein) within the qTNL2-1 region. In the same way, a gene identified as being involved in carbohydrates metabolic process, among other biological processes, was found (OsHUGT1, Os02g44510). The other set of genes encoding the Aquaporin (OsPIP1 (plasma membrane intrinsic protein 1) or OsRWC1 (water channel protein), Os02g44630), the cytochrome P450 (Os02g44654) [61,62], zinc finger RING-type domain (Os02g44700), OsSID1 or OsIDD4 (Os02g45054) [63], and MYB (Myeloblastosis) transcription factor encoding gene (Os02g4580) was found within the detected QTL, among other genes (Table 3).

4. Discussion

4.1. Tiller Number Positively Correlated with Shoot Dry Weight, While Shoot and Root Dry Weights Showed Positive Relationship

In rice, tillers represent the branches that develop from the leaf axils at each unelongated node of the main shoot or from other tillers during vegetative growth. It is said that tiller number determines the number panicles per plant [39,64,65], which happens to be an important yield component of rice. The number of tillers is influenced by environmental factors [66], agricultural practices, including nitrogen fertilization [67,68,69,70,71,72], or genetic factors [65,73,74,75,76,77]. It is also expected that an increase or decrease in tiller number per plant could cause significant changes in the shoot biomass. Here, the results of the linear regression analysis showed that shoot dry weight positively correlated with tiller number (R2 = 0.583*** and R2 = 0.317***, under low and normal nitrogen cultivation, respectively) and plant height (R2 = 0.374*** and R2 = 0.350***, under low and normal nitrogen cultivation, respectively). Therefore, these findings suggest that a unit increase in tiller number would result in the increase in shoot dry weight by the same measure. In the same way, we recorded that panicle number showed a strong positive correlation with tiller number (R2 = 0.560*** and R2 = 0.603***, under low and normal nitrogen cultivation, respectively). It could then be said that rice plants with a high tiller number would eventually develop more panicles, which would result in an increased productivity.
Roots and shoot are said to be positionally and developmentally related in various ways [78,79]. Therefore, the recorded positive correlation between shoot dry weight and root dry weight (R2 = 0.491*** and R2 = 0.309***, under low and normal nitrogen cultivation, respectively), coupled with the respective phenotypes of the mapping population concerning SDW and RDW (Figure 2G,H: 96.6% and 87.2% of the DH lines showed a high SDW and RDW, respectively, similar to 93-11 (P1, indica and high tillering cultivar)), would suggest that a robust rooting system may contribute to achieving an optimum shoot growth and development in rice. In addition, it is thought that the observed difference of the correlation level between tiller number and other traits under normal and low nitrogen conditions would also give insights into the recorded N-effects on the phenotypic traits analyzed, and on the genetic loci difference controlling tiller number.

4.2. Quantiative Trait Loci (QTLs) Associated with Shoot Branching under Low Nitrogen Cultivation Conditions Are Identified

Several agronomic traits of rice have been suggested to be controlled by multiple quantitative trait loci (QTLs) [27,80,81]. A number of studies have proposed a large number of putative QTLs controlling tiller number and panicle number in rice. Most of the reported QTLs related to tiller number or panicle number are mapped across the 10–11 chromosomes of rice, using various methods, marker-types on different sorts of rice populations, at various developmental stages [27,39,82]. In this study, of the 28 QTLs identified for all the growth related and productivity related traits, only three were associated with tiller number. Of this number, two (qTNN2-1 and qTNN4-1) were detected in the set of DH lines grown under normal nitrogen conditions, and one (qTNL2-1) in that grown under low nitrogen conditions, were mapped on chromosomes 2 and 4. Meanwhile, four QTLs related to panicle number were detected (Table 3), with qPNL4-1 being the only QTL associated with panicle number identified in DH lines grown under low nitrogen conditions, and mapped on chromosome 4. A study conducted by Liu et al. [82] identified 12 QTLs associated with tiller number mapped on chromosomes 2 (two QTLs: 12.5 cM and 23.1 cM), 3 (three QTLs: 6.0 cM, 17.2 cM, and 20.4 cM) and 6–9 (three QTLs on chromosome 6: 8.5 cM, 8.8 cM, and 28 cM; one on chromosome 7: 24.5 cM; one on chromosome 8: 8.3 cM; one on chromosome 9: 17.0 cM, and one on chromosome 12: 2.6 cM) using rice microsatellites (RM) markers on single segment substitution lines (SSSLs). In a similar perspective, Yan et al. [39] reported 15 QTLs associated with tiller number at different growth stages, and mapped on chromosomes 1–4, 6–8, and 12. In a converse approach, Liu et al. [27] identified 27 QTLs associated with tiller number at different development stages in rice, and mapped on chromosomes 1–12. However, Zhu and his colleagues reported a QTL associated with panicle number in rice, mapped on chromosome 1 (qPN1) [83]. None of the previously reported QTLs associated with tiller number or panicle number overlapped with the identified QTLs by this study. However, some of them are mapped on similar chromosomes but different positions.

4.3. QTLs qTNL2-1 and qTN4-1 Harbor Candidate Genes Proposed to Be Involved in the Control of Shoot Branching in Rice under Low and Regular Nitrogen Cultivation

In the perspective of unveiling the identity of candidate genes located within the genetic region covered by the QTL qTNL2-1 associated with the control of shoot branching under low nitrogen cultivation, we used the available maker positions to explore this region. Interestingly, a substantial number of genes were proposed to have transcription factor binding activity. This is the case for two GRAS transcription factor (TF) encoding genes (Os02g44360 and Os2g44370), which have a sequence similarity of about 86.6%. Proteins in the GRAS domain family are known as key players in gibberellin (GA) signaling [84,85], which regulates various aspects of plant growth and development, by inhibiting the proliferation and expansion of cell-mediated plant growth [54]. Similarly, OsHMG1 (Os02g44930), OsSID1 (Os02g45054), also known as OsIDD4 (a Zinc finger C2H2-type domain containing protein), and MYB TF encoding gene (Os02g45080), are proposed to code for transcription factors.
Another set of interesting genes, such as OsGPX3, was earlier reported to be involved in the development of the roots and shoot in rice, in the stress response mechanisms [55], and abscisic acid (ABA) response in rice [56]; OsDhn1, earlier suggested to be involved in abiotic stress response in rice [57,58]; OsRALF8 [86], previously reported to play a critical role in the regulation of plant defense and development (inhibition or suppression of growth) [59], could play important roles in the control of shoot branching, while maintaining a balanced plant growth and development under low nitrogen cultivation in rice. Nietzsche et al. [60] supported that carbohydrates regulate different aspects of plant growth via the modulation of cell division and expansion. Here, a DUF encoding gene OsFBDUF1 or OsMAIF23 found within the qTNL2-1 region, and OsHUGT1 known as UDP-glucose glycoprotein glucosyltransferase 1 precursor, predicted as being involved in carbohydrates metabolic process would contribute to the observed tiller number phenotype under low nitrogen cultivation.
In the same way, Aquaporin (OsPIP1 or OsRWC1), and OsDhn1 or OsLip9 [57,58], associated with cytochrome P450 (earlier suggested to be involved in auxin signaling by Chaban et al. [87] or encoded by OsMAX1, a strigolactone biosynthesis genes according to Marek et al. [61,88] or encoded by OsSLR1/OsGAI (SLENDER RICE 1 encodes a DELLA protein belonging to GRAS [89]), which modulate gibberellin response in rice [62]), are believed to actively contribute to modulating water movements in and outside the cell to facilitate transport of nutrients, including nitrogen, as well as plant growth and development under low nitrogen conditions. From another perspective, Zou et al. [90,91] identified and characterized the OsHTD1 (HIGH-TILLERING DWARF1, Os04g46470) gene encoding a carotenoid cleavage dioxygenase (CCD7 or D7 ortholog of Arabidopsis MAX3, involved in the strigolactone biosynthesis). The OsHTD1 was reported to be required for negative regulation of the outgrowth of axillary buds. Similarly, Wang et al. [29] demonstrated that a partial loss-of-function allele of OsHTD1/CCD7 gene increased tiller number and improved rice grain yield. Here, our finding revealed that the OsHTD1 gene (Chr4:27567824..27570926) is located within the QTL qPNL4-1 mapped on chromosome 4, identified by this study as being associated with the control of panicle number under low nitrogen cultivation. The qPNL4-1 is flanked by ad04009559 (left marker, 26443283 kbp) and ah04001252 (right marker, 28788134 kbp).
The breakthrough is that many of the candidate genes found either in qTNL2-1 or qTNN4-1 are also proposed to be involved in plant stress signaling or response mechanisms. A recent study investigating the transcriptional regulation of hormonal biosynthesis or signaling pathways genes under drought stress conditions, proposed that shoot branching related genes could be involved in the adaptive response mechanisms of plants towards abiotic stress tolerance in Arabidopsis [92]. We could speculate that under normal growth conditions, plants utilize their resources for their growth and development. However, under a reduced level of nitrogen, which could also be perceived as a stress, plants would activate abiotic stress responsive genes and/or growth related genes, in order to either escape or avoid the stress or alleviate the adverse effects that would be caused by the stress, while tending to maintain a balanced growth. Owing to the above, a possible crosstalk between plant growth and development and defense mechanisms could be suggested.

5. Conclusions

Shoot branching is an important agronomic trait contributing to the yield and productivity of cereal food crops, such as rice. In plants, shoot branching is controlled by a complex hormonal signaling network, operating either in synergetic or antagonistic manner. It is well established that most of the genes involved in the regulation of shoot branching, through promotion or inhibition, are also involved in the adaptive response mechanisms of plants towards stress tolerance or resistance (abiotic or biotic stress). In this study, we identified QTLs proposed to be associated with tiller number and panicle number under low nitrogen cultivation, mapped to chromosomes 2 (qTNL2-1) and chromosome 4 (qPNL4-1). Another QTL associated with tiller number was detected in the doubled haploid population grown under normal nitrogen cultivation conditions (qTNN4-1) and mapped on chromosome 4. These QTLs harbor candidate genes with interesting annotated functions, and proposed for playing important roles in the adaptive response mechanisms towards stress (abiotic and biotic) in rice, in addition to being associated with shoot branching. Candidate genes encoding GRAS protein, cytochrome P450, OsRALF8, Aquaporin and OsDhn1, among others, would play important roles in controlling shoot branching under low nitrogen cultivation with regard to their functions and their similarity with previously reported genes. Therefore, this study suggests a possible crosstalk between the regulation of plant growth and development, and the stress response mechanisms. Functional analysis would help identify major genes and elucidate their function.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/genes12050745/s1, Figure S1: Changes in weather parameters during the rice cropping season of 2020, Figure S2: Quantile–Quantile plots, Table S1: Person r correlation between traits under normal nitrogen cultivation, Table S2: Person r correlation between traits under low nitrogen cultivation, Table S3: Shapiro Wilk W-statistic for the test of normality of the distribution, Table S4: List of candidate genes located within the qTNN4-1.

Author Contributions

Conceptualization, J.-H.L., J.-Y.L. and D.-J.S.; methodology, J.-H.L., J.-W.K. and N.-R.K.; validation, J.-M.K., J.-H.L. and D.-J.S.; formal analysis, Y.-H.K., S.-M.L., N.-R.K. and J.-K.C.; investigation, Y.-H.K., J.-H.L., J.-W.K. and S.-M.L.; resources, J.-M.K., J.-H.L., J.-H.C. and J.-Y.L.; data curation, N.-R.K.; writing—original draft preparation, Y.-H.K. and N.-R.K.; writing—review and editing, J.-H.L. and D.-J.S.; visualization and supervision, J.-H.L. and J.-M.K.; project administration, J.-H.L., J.-H.C. and J.-M.K.; funding acquisition, J.-M.K. and J.-H.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Cooperative Research Program for Agriculture Science and Technology Development (Project No. PJ01506901).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

This study was supported by 2021 the RDA Fellowship Program of NICS, Rural Development Administration, Korea.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. DESA UN. World Population Projected to Reach 9.8 Billion in 2050, and 11.2 Billion in 2100. 2019, Volume 2. Available online: https://www.un.org/development/desa/en/news/population/world-population-prospects-2017.html (accessed on 8 March 2021).
  2. Olson, R.; Kurtz, L. Crop nitrogen requirements, utilization, and fertilization. Nitrogen Agric. Soils 1982, 22, 567–604. [Google Scholar] [CrossRef]
  3. Leghari, S.J.; Wahocho, N.A.; Laghari, G.M.; Hafeez Laghari, A.; Mustafa Bhabhan, G.; Hussain Talpur, K.; Bhutto, T.A.; Wahocho, S.A.; Lashari, A.A. Role of nitrogen for plant growth and development: A review. Adv. Environ. Biol. 2016, 10, 209–219. [Google Scholar]
  4. Causin, H.F. The central role of amino acids on nitrogen utilization and plant growth. J. Plant Physiol. 1996, 149, 358–362. [Google Scholar] [CrossRef]
  5. Novoa, R.; Loomis, R.S. Nitrogen and plant production. Plant Soil 1981, 58, 177–204. [Google Scholar] [CrossRef]
  6. Kirkby, E. Plant growth in relation to nitrogen supply. Ecol. Bull. 1981, 33, 249–267. [Google Scholar]
  7. Van Keulen, H. Nitrogen Requirements of Rice with Special Reference to Java; CRIA: Bogor, Indonesia, 1977; Volume 30. [Google Scholar]
  8. Fageria, N.; Baligar, V.C. Lowland rice response to nitrogen fertilization. Commun. Soil Sci. Plant Anal. 2001, 32, 1405–1429. [Google Scholar] [CrossRef]
  9. Mohanty, S.R.; Bodelier, P.L.; Floris, V.; Conrad, R.J.A. Differential effects of nitrogenous fertilizers on methane-consuming microbes in rice field and forest soils. Appl. Environ. Microbiol. 2006, 72, 1346–1354. [Google Scholar] [CrossRef] [Green Version]
  10. Dubey, S.K. Microbial ecology of methane emission in rice agroecosystem: A review. Appl. Ecol. Environ. Res. 2005, 3, 1–27. [Google Scholar] [CrossRef]
  11. Minamikawa, K.; Sakai, N.; Yagi, K.J.M. Methane emission from paddy fields and its mitigation options on a field scale. Microbes Environ. 2006, 21, 135–147. [Google Scholar] [CrossRef] [Green Version]
  12. Schimel, J.J.N. Rice, microbes and methane. Nature 2000, 403, 375–377. [Google Scholar] [CrossRef]
  13. Namai, S.; Toriyama, K.; Fukuta, Y. Genetic variations in dry matter production and physiological nitrogen use efficiency in rice (Oryza sativa L.) varieties. Breed. Sci. 2009, 59, 269–276. [Google Scholar] [CrossRef] [Green Version]
  14. Fei, W.; Peng, S.-B. Yield potential and nitrogen use efficiency of China’s super rice. J. Integr. Agric. 2017, 16, 1000–1008. [Google Scholar] [CrossRef] [Green Version]
  15. Raun, W.R.; Johnson, G.V. Improving nitrogen use efficiency for cereal production. Agron. J. 1999, 91, 357–363. [Google Scholar] [CrossRef] [Green Version]
  16. Chen, G.; Chen, Y.; Zhao, G.; Cheng, W.; Guo, S.; Zhang, H.; Shi, W. Do high nitrogen use efficiency rice cultivars reduce nitrogen losses from paddy fields? Agric. Ecosyst. Environ. 2015, 209, 26–33. [Google Scholar] [CrossRef] [Green Version]
  17. Pires, M.V.; da Cunha, D.A.; de Matos Carlos, S.; Costa, M.H. Nitrogen-use efficiency, nitrous oxide emissions, and cereal production in Brazil: Current trends and forecasts. PLoS ONE 2015, 10, e0135234. [Google Scholar] [CrossRef] [Green Version]
  18. Neue, H.-U.; Sass, R.L. Trace gas emissions from rice fields. In Global Atmospheric-Biospheric Chemistry; Springer: Berlin/Heidelberg, Germany, 1994; pp. 119–147. [Google Scholar] [CrossRef]
  19. Yasuno, N.; Takamure, I.; Kidou, S.; Tokuji, Y.; Ureshi, A.; Funabiki, A.; Ashikaga, K.; Yamanouchi, U.; Yano, M.; Kato, K. Rice shoot branching requires an ATP-binding cassette subfamily G protein. New Phytol. 2009, 182, 91–101. [Google Scholar] [CrossRef]
  20. Ariyaratne, M.; Takamure, I.; Kato, K. Shoot branching control by reduced culm number 4 in rice (Oryza sativa L.). New Phytol. 2009, 176, 744–748. [Google Scholar] [CrossRef]
  21. Leyser, O. The control of shoot branching: An example of plant information processing. Plant Cell Environ. 2009, 32, 694–703. [Google Scholar] [CrossRef]
  22. Yaish, M.; Guevara, D.; El-Kereamy, A.; Rothstein, S. Axillary shoot branching in plants. In Plant Developmental Biology-Biotechnological Perspectives; Springer: Berlin/Heidelberg, Germany, 2010; pp. 37–52. [Google Scholar]
  23. Umehara, M.; Hanada, A.; Yoshida, S.; Akiyama, K.; Arite, T.; Takeda-Kamiya, N.; Magome, H.; Kamiya, Y.; Shirasu, K.; Yoneyama, K.J.N. Inhibition of shoot branching by new terpenoid plant hormones. Nature 2008, 455, 195–200. [Google Scholar] [CrossRef]
  24. Gomez-Roldan, V.; Fermas, S.; Brewer, P.B.; Puech-Pagès, V.; Dun, E.A.; Pillot, J.-P.; Letisse, F.; Matusova, R.; Danoun, S.; Portais, J.-C.J.N. Strigolactone inhibition of shoot branching. Nature 2008, 455, 189–194. [Google Scholar] [CrossRef]
  25. Barbier, F.F.; Dun, E.A.; Kerr, S.C.; Chabikwa, T.G.; Beveridge, C.A. An update on the signals controlling shoot branching. Trends Plant Sci. 2019, 24, 220–236. [Google Scholar] [CrossRef]
  26. Peleg, Z.; Blumwald, E. Hormone balance and abiotic stress tolerance in crop plants. Curr. Opin. Plant Biol. 2011, 14, 290–295. [Google Scholar] [CrossRef]
  27. Liu, G.; Li, M.; Wen, J.; Du, Y.; Zhang, Y.-M. Functional mapping of quantitative trait loci associated with rice tillering. Mol. Genet. Genom. 2010, 284, 263–271. [Google Scholar] [CrossRef]
  28. Zhao, J.; Wang, T.; Wang, M.; Liu, Y.; Yuan, S.; Gao, Y.; Yin, L.; Sun, W.; Peng, L.; Zhang, W.J.P.; et al. DWARF3 participates in an SCF complex and associates with DWARF14 to suppress rice shoot branching. Plant Cell Physiol. 2014, 55, 1096–1109. [Google Scholar] [CrossRef] [Green Version]
  29. Wang, Y.; Shang, L.; Yu, H.; Zeng, L.; Hu, J.; Ni, S.; Rao, Y.; Li, S.; Chu, J.; Meng, X.J.; et al. A strigolactone biosynthesis gene contributed to the green revolution in rice. Mol. Plant 2020, 13, 923–932. [Google Scholar] [CrossRef]
  30. Gao, Z.; Qian, Q.; Liu, X.; Yan, M.; Feng, Q.; Dong, G.; Liu, J.; Han, B. Dwarf 88, a novel putative esterase gene affecting architecture of rice plant. Plant Mol. Biol. 2009, 71, 265–276. [Google Scholar] [CrossRef]
  31. Liu, W.; Wu, C.; Fu, Y.; Hu, G.; Si, H.; Zhu, L.; Luan, W.; He, Z.; Sun, Z.J.P. Identification and characterization of HTD2: A novel gene negatively regulating tiller bud outgrowth in rice. Planta 2009, 230, 649–658. [Google Scholar] [CrossRef] [PubMed]
  32. Kulkarni, K.P.; Vishwakarma, C.; Sahoo, S.P.; Lima, J.M.; Nath, M.; Dokku, P.; Gacche, R.N.; Mohapatra, T.; Robin, S.; Sarla, N.; et al. A substitution mutation in OsCCD7 cosegregates with dwarf and increased tillering phenotype in rice. J. Genet. 2014, 93, 389–401. [Google Scholar] [CrossRef] [PubMed]
  33. Lo, S.-F.; Yang, S.-Y.; Chen, K.-T.; Hsing, Y.-I.; Zeevaart, J.A.; Chen, L.-J.; Yu, S.-M. A novel class of gibberellin 2-oxidases control semidwarfism, tillering, and root development in rice. Plant Cell 2008, 20, 2603–2618. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Zhang, Q.Q.; Wang, J.G.; Wang, L.Y.; Wang, J.F.; Wang, Q.; Yu, P.; Bai, M.Y.; Fan, M. Gibberellin repression of axillary bud formation in Arabidopsis by modulation of DELLA-SPL9 complex activity. J. Integr. Plant Biol. 2020, 62, 421–432. [Google Scholar] [CrossRef] [Green Version]
  35. Zou, X.; Wang, Q.; Chen, P.; Yin, C.; Lin, Y. Strigolactones regulate shoot elongation by mediating gibberellin metabolism and signaling in rice (Oryza sativa L.). J. Plant Physiol. 2019, 237, 72–79. [Google Scholar] [CrossRef]
  36. Wang, M.; Le Moigne, M.-A.; Bertheloot, J.; Crespel, L.; Perez-Garcia, M.-D.; Ogé, L.; Demotes-Mainard, S.; Hamama, L.; Davière, J.-M.; Sakr, S. BRANCHED1: A key hub of shoot branching. Front. Plant Sci. 2019, 10, 76. [Google Scholar] [CrossRef] [Green Version]
  37. Komatsu, K.; Maekawa, M.; Ujiie, S.; Satake, Y.; Furutani, I.; Okamoto, H.; Shimamoto, K.; Kyozuka, J. LAX and SPA: Major regulators of shoot branching in rice. Proc. Natl. Acad. Sci. USA 2003, 100, 11765–11770. [Google Scholar] [CrossRef] [Green Version]
  38. Yoshida, S.; Ikegami, M.; Kuze, J.; Sawada, K.; Hashimoto, Z.; Ishii, T.; Nakamura, C.; Kamijima, O. QTL analysis for plant and grain characters of sake-brewing rice using a doubled haploid population. Breed. Sci. 2002, 52, 309–317. [Google Scholar] [CrossRef] [Green Version]
  39. Yan, J.; Zhu, J.; He, C.; Benmoussa, M.; Wu, P. Quantitative trait loci analysis for the developmental behavior of tiller number in rice (Oryza sativa L.). Theor. Appl. Genet. 1998, 97, 267–274. [Google Scholar] [CrossRef]
  40. Onishi, K.; Horiuchi, Y.; Ishigoh-Oka, N.; Takagi, K.; Ichikawa, N.; Maruoka, M.; Sano, Y. A QTL cluster for plant architecture and its ecological significance in Asian wild rice. Breed. Sci. 2007, 57, 7–16. [Google Scholar] [CrossRef] [Green Version]
  41. Khush, G.S. Origin, dispersal, cultivation and variation of rice. Plant Mol. Biol. 1997, 35, 25–34. [Google Scholar] [CrossRef]
  42. Zhu, Q.; Ge, S. Phylogenetic relationships among A-genome species of the genus Oryza revealed by intron sequences of four nuclear genes. New Phytol. 2005, 167, 249–265. [Google Scholar] [CrossRef]
  43. Kato, S. On the affinity of rice varieties as shown by fertility of hybrid plants. Bull. Sci. Fac. Agric. Kyushu Univ. 1928, 3, 132–147. [Google Scholar]
  44. Park, H.J.; Kim, S.U.; Jung, K.Y.; Lee, S.; Choi, Y.D.; Owens, V.N.; Kumar, S.; Yun, S.W.; Hong, C.O. Cadmium phytoavailability from 1976 through 2016: Changes in soil amended with phosphate fertilizer and compost. Sci. Total Environ. 2021, 762, 143132. [Google Scholar] [CrossRef]
  45. Thorpe, R.S. Multiple group principal component analysis and population differentiation. J. Zool. 1988, 216, 37–40. [Google Scholar] [CrossRef]
  46. R Development Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2013. [Google Scholar]
  47. Kang, J.-W.; Kabange, N.R.; Phyo, Z.; Park, S.-Y.; Lee, S.-M.; Lee, J.-Y.; Shin, D.; Cho, J.H.; Park, D.-S.; Ko, J.-M. Combined Linkage Mapping and Genome-Wide Association Study Identified QTLs Associated with Grain Shape and Weight in Rice (Oryza sativa L.). Agronomy 2020, 10, 1532. [Google Scholar] [CrossRef]
  48. Cheon, K.-S.; Baek, J.; Cho, Y.-I.; Jeong, Y.-M.; Lee, Y.-Y.; Oh, J.; Won, Y.J.; Kang, D.-Y.; Oh, H.; Kim, S.L.; et al. Single nucleotide polymorphism (SNP) discovery and kompetitive allele-specific PCR (KASP) marker development with Korean japonica rice varieties. Plant Breed. Biotechnol. 2018, 6, 391–403. [Google Scholar] [CrossRef]
  49. Seo, J.; Lee, G.; Jin, Z.; Kim, B.; Chin, J.H.; Koh, H.-J. Development and application of indica–japonica SNP assays using the Fluidigm platform for rice genetic analysis and molecular breeding. Mol. Breed. 2020, 40, 1–16. [Google Scholar] [CrossRef] [Green Version]
  50. Meng, L.; Li, H.; Zhang, L.; Wang, J. QTL IciMapping: Integrated software for genetic linkage map construction and quantitative trait locus mapping in biparental populations. Crop. J. 2015, 3, 269–283. [Google Scholar] [CrossRef] [Green Version]
  51. Nada, R.M.; Abogadallah, G.M. Restricting the above ground sink corrects the root/shoot ratio and substantially boosts the yield potential per panicle in field-grown rice (Oryza sativa L.). Physiol. Plant. 2016, 156, 371–386. [Google Scholar] [CrossRef]
  52. Andrews, M.; Sprent, J.; Raven, J.; Eady, P.E. Relationships between shoot to root ratio, growth and leaf soluble protein concentration of Pisum sativum, Phaseolus vulgaris and Triticum aestivum under different nutrient deficiencies. Plant Cell Environ. 1999, 22, 949–958. [Google Scholar] [CrossRef]
  53. Xu, W.; Cui, K.; Xu, A.; Nie, L.; Huang, J.; Peng, S. Drought stress condition increases root to shoot ratio via alteration of carbohydrate partitioning and enzymatic activity in rice seedlings. Acta Physiol. Plant. 2015, 37, 9. [Google Scholar] [CrossRef]
  54. Pysh, L.D.; Wysocka-Diller, J.W.; Camilleri, C.; Bouchez, D.; Benfey, P.N. The GRAS gene family in Arabidopsis: Sequence characterization and basic expression analysis of the SCARECROW-LIKE genes. Plant J. 1999, 18, 111–119. [Google Scholar] [CrossRef]
  55. Passaia, G.; Fonini, L.S.; Caverzan, A.; Jardim-Messeder, D.; Christoff, A.P.; Gaeta, M.L.; de Araujo Mariath, J.E.; Margis, R.; Margis-Pinheiro, M. The mitochondrial glutathione peroxidase GPX3 is essential for H2O2 homeostasis and root and shoot development in rice. Plant Sci. 2013, 208, 93–101. [Google Scholar] [CrossRef]
  56. Paiva, A.L.S.; Passaia, G.; Jardim-Messeder, D.; Nogueira, F.C.; Domont, G.B.; Margis-Pinheiro, M. The mitochondrial isoform glutathione peroxidase 3 (OsGPX3) is involved in ABA responses in rice plants. J. Proteom. 2021, 232, 104029. [Google Scholar] [CrossRef] [PubMed]
  57. Lee, S.-C.; Kim, S.-H.; Kim, S.-R. Drought inducible OsDhn1 promoter is activated by OsDREB1A and OsDREB1D. J. Plant Biol. 2013, 56, 115–121. [Google Scholar] [CrossRef]
  58. Lee, S.-C.; Lee, M.-Y.; Kim, S.-J.; Jun, S.-H.; An, G.; Kim, S.-R. Characterization of an abiotic stress-inducible dehydrin gene, OsDhn1, in rice (Oryza sativa L.). Mol. Cells 2005, 19, 212–218. [Google Scholar]
  59. Sharma, A.; Hussain, A.; Mun, B.-G.; Imran, Q.M.; Falak, N.; Lee, S.-U.; Kim, J.Y.; Hong, J.K.; Loake, G.J.; Ali, A.; et al. Comprehensive analysis of plant rapid alkalization factor (RALF) genes. Plant Physiol. Biochem. 2016, 106, 82–90. [Google Scholar] [CrossRef]
  60. Nietzsche, M.; Schießl, I.; Börnke, F. The complex becomes more complex: Protein-protein interactions of SnRK1 with DUF581 family proteins provide a framework for cell-and stimulus type-specific SnRK1 signaling in plants. Front. Plant Sci. 2014, 5, 54. [Google Scholar] [CrossRef] [Green Version]
  61. Zhang, Y.; Van Dijk, A.D.; Scaffidi, A.; Flematti, G.R.; Hofmann, M.; Charnikhova, T.; Verstappen, F.; Hepworth, J.; Van Der Krol, S.; Leyser, O.; et al. Rice cytochrome P450 MAX1 homologs catalyze distinct steps in strigolactone biosynthesis. Nat. Chem. Biol. 2014, 10, 1028–1033. [Google Scholar] [CrossRef]
  62. Luo, A.; Qian, Q.; Yin, H.; Liu, X.; Yin, C.; Lan, Y.; Tang, J.; Tang, Z.; Cao, S.; Wang, X.J.P.; et al. EUI1, encoding a putative cytochrome P450 monooxygenase, regulates internode elongation by modulating gibberellin responses in rice. Plant Cell Physiol. 2006, 47, 181–191. [Google Scholar] [CrossRef]
  63. Deng, L.; Li, L.; Zhang, S.; Shen, J.; Li, S.; Hu, S.; Peng, Q.; Xiao, J.; Wu, C. Suppressor of rid1 (SID1) shares common targets with RID1 on florigen genes to initiate floral transition in rice. PLoS Genet. 2017, 13, e1006642. [Google Scholar] [CrossRef]
  64. Xing, Y.; Zhang, Q. Genetic and molecular bases of rice yield. Ann. Rev. Plant Biol. 2010, 61, 421–442. [Google Scholar] [CrossRef]
  65. Zhen, X.; Li, X.; Yu, J.; Xu, F. OsATG8c-mediated increased autophagy regulates the yield and nitrogen use efficiency in rice. Int. J. Mol. Sci. 2019, 20, 4956. [Google Scholar] [CrossRef] [Green Version]
  66. Krishnan, P.; Ramakrishnan, B.; Reddy, K.R.; Reddy, V. High-temperature effects on rice growth, yield, and grain quality. Adv. Agron. 2011, 111, 87–206. [Google Scholar] [CrossRef]
  67. Wang, Y.; Lu, J.; Ren, T.; Hussain, S.; Guo, C.; Wang, S.; Cong, R.; Li, X. Effects of nitrogen and tiller type on grain yield and physiological responses in rice. AoB Plants 2017, 9. [Google Scholar] [CrossRef]
  68. Zhang, J.; Tong, T.; Potcho, P.M.; Huang, S.; Ma, L.; Tang, X.J.A. Nitrogen Effects on Yield, Quality and Physiological Characteristics of Giant Rice. Agronomy 2020, 10, 1816. [Google Scholar] [CrossRef]
  69. Tanaka, A.; Garcia, C. Studies of the relationship between tillering and nitrogen uptake of the rice plant: 2. Relation between tillering and nitrogen metabolism of the plant. Soil Sci. Plant Nutr. 1965, 11, 31–37. [Google Scholar] [CrossRef] [Green Version]
  70. Tian, G.; Gao, L.; Kong, Y.; Hu, X.; Xie, K.; Zhang, R.; Ling, N.; Shen, Q.; Guo, S. Improving rice population productivity by reducing nitrogen rate and increasing plant density. PLoS ONE 2017, 12, e0182310. [Google Scholar] [CrossRef] [Green Version]
  71. Huang, M.; Shan, S.; Cao, J.; Fang, S.; Tian, A.; Liu, Y.; Cao, F.; Yin, X.; Zou, Y. Primary-tiller panicle number is critical to achieving high grain yields in machine-transplanted hybrid rice. Sci. Rep. 2020, 10, 1–7. [Google Scholar] [CrossRef] [Green Version]
  72. Badshah, M.; Naimei, T.; Zou, Y.; Ibrahim, M.; Wang, K. Yield and tillering response of super hybrid rice Liangyoupeijiu to tillage and establishment methods. Crop. J. 2014, 2, 79–86. [Google Scholar] [CrossRef] [Green Version]
  73. Wang, J.; Lu, K.; Nie, H.; Zeng, Q.; Wu, B.; Qian, J.; Fang, Z.J.R. Rice nitrate transporter OsNPF7.2 positively regulates tiller number and grain yield. Rice 2018, 11, 1–13. [Google Scholar] [CrossRef] [Green Version]
  74. Liu, Y.; Wang, H.; Jiang, Z.; Wang, W.; Xu, R.; Wang, Q.; Zhang, Z.; Li, A.; Liang, Y.; Ou, S.J.N. Genomic basis of geographical adaptation to soil nitrogen in rice. Nature 2021, 1–6. [Google Scholar] [CrossRef]
  75. Liao, Z.; Yu, H.; Duan, J.; Yuan, K.; Yu, C.; Meng, X.; Kou, L.; Chen, M.; Jing, Y.; Liu, G.; et al. SLR1 inhibits MOC1 degradation to coordinate tiller number and plant height in rice. Nat. Commun. 2019, 10, 1–9. [Google Scholar] [CrossRef] [Green Version]
  76. Lacchini, E.; Kiegle, E.; Castellani, M.; Adam, H.; Jouannic, S.; Gregis, V.; Kater, M.M. CRISPR-mediated accelerated domestication of African rice landraces. PLoS ONE 2020, 15, e0229782. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  77. Wen, X.; Sun, L.; Chen, Y.; Xue, P.; Yang, Q.; Wang, B.; Yu, N.; Cao, Y.; Zhang, Y.; Gong, K.; et al. Rice dwarf and low tillering 10 (OsDLT10) regulates tiller number by monitoring auxin homeostasis. Plant Sci. 2020, 297, 110502. [Google Scholar] [CrossRef] [PubMed]
  78. Groff, P.A.; Kaplan, D.R. The relation of root systems to shoot systems in vascular plants. Bot. Rev. 1988, 54, 387–422. [Google Scholar] [CrossRef]
  79. Qi, Y.; Wei, W.; Chen, C.; Chen, L. Plant root-shoot biomass allocation over diverse biomes: A global synthesis. Glob. Ecol. Conserv. 2019, 18, e00606. [Google Scholar] [CrossRef]
  80. Yano, M.; Sasaki, T. Genetic and molecular dissection of quantitative traits in rice. In Oryza: From Molecule to Plant; Springer: Berlin/Heidelberg, Germany, 1997; pp. 145–153. [Google Scholar] [CrossRef]
  81. Yamamoto, T.; Yonemaru, J.; Yano, M. Towards the understanding of complex traits in rice: Substantially or superficially? DNA Res. 2009, 16, 141–154. [Google Scholar] [CrossRef] [Green Version]
  82. Liu, G.; Zhu, H.; Zhang, G.; Li, L.; Ye, G. Dynamic analysis of QTLs on tiller number in rice (Oryza sativa L.) with single segment substitution lines. Theor. Appl. Genet. 2012, 125, 143–153. [Google Scholar] [CrossRef]
  83. Zhu, J.; Zhou, Y.; Liu, Y.; Wang, Z.; Tang, Z.; Yi, C.; Tang, S.; Gu, M.; Liang, G. Fine mapping of a major QTL controlling panicle number in rice. Mol. Breed. 2011, 27, 171–180. [Google Scholar] [CrossRef]
  84. Bolle, C. The role of GRAS proteins in plant signal transduction and development. Planta 2004, 218, 683–692. [Google Scholar] [CrossRef]
  85. Xu, K.; Chen, S.; Li, T.; Ma, X.; Liang, X.; Ding, X.; Liu, H.; Luo, L. OsGRAS23, a rice GRAS transcription factor gene, is involved in drought stress response through regulating expression of stress-responsive genes. BMC Plant Biol. 2015, 15, 141. [Google Scholar] [CrossRef] [Green Version]
  86. Zhang, X.; Yang, Z.; Wu, D.; Yu, F. RALF-FERONIA signaling: Linking plant immune response with cell growth. Plant Commun. 2020, 100084. [Google Scholar] [CrossRef]
  87. Chaban, C.; Waller, F.; Furuya, M.; Nick, P. Auxin responsiveness of a novel cytochrome P450 in rice coleoptiles. Plant Physiol. 2003, 133, 2000–2009. [Google Scholar] [CrossRef] [Green Version]
  88. Marzec, M.; Situmorang, A.; Brewer, P.B.; Brąszewska, A.J.G. Diverse roles of MAX1 homologues in rice. Genes 2020, 11, 1348. [Google Scholar] [CrossRef]
  89. Itoh, H.; Shimada, A.; Ueguchi-Tanaka, M.; Kamiya, N.; Hasegawa, Y.; Ashikari, M.; Matsuoka, M. Overexpression of a GRAS protein lacking the DELLA domain confers altered gibberellin responses in rice. Plant J. 2005, 44, 669–679. [Google Scholar] [CrossRef]
  90. Zou, J.; Zhang, S.; Zhang, W.; Li, G.; Chen, Z.; Zhai, W.; Zhao, X.; Pan, X.; Xie, Q.; Zhu, L. The rice HIGH-TILLERING DWARF1 encoding an ortholog of Arabidopsis MAX3 is required for negative regulation of the outgrowth of axillary buds. Plant J. 2006, 48, 687–698. [Google Scholar] [CrossRef]
  91. Zou, J.; Chen, Z.; Zhang, S.; Zhang, W.; Jiang, G.; Zhao, X.; Zhai, W.; Pan, X.; Zhu, L. Characterizations and fine mapping of a mutant gene for high tillering and dwarf in rice (Oryza sativa L.). Planta 2005, 222, 604–612. [Google Scholar] [CrossRef]
  92. Rolly, N.K.; Mun, B.-G.; Yun, B.-W. Insights into the Transcriptional Regulation of Branching Hormonal Signaling Pathways Genes under Drought Stress in Arabidopsis. Genes 2021, 12, 298. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Frequency distribution of traits. Frequency distribution of plant height of a DH population grown under a normal (A) and low (B) nitrogen (N) cultivation regimes. Frequency distribution of tiller number under normal (C) and low (D) nitrogen levels. Frequency distribution of culm length (E,F) and panicle length (G,H) under normal and low N cultivation. Frequency distribution of panicle number (I,J) and root length (K,L) under normal and low N cultivation. Frequency distribution of shoot (M,N) and root (O,P) dry weights under normal and low N levels.
Figure 1. Frequency distribution of traits. Frequency distribution of plant height of a DH population grown under a normal (A) and low (B) nitrogen (N) cultivation regimes. Frequency distribution of tiller number under normal (C) and low (D) nitrogen levels. Frequency distribution of culm length (E,F) and panicle length (G,H) under normal and low N cultivation. Frequency distribution of panicle number (I,J) and root length (K,L) under normal and low N cultivation. Frequency distribution of shoot (M,N) and root (O,P) dry weights under normal and low N levels.
Genes 12 00745 g001
Figure 2. Genotype–phenotype correlation of the doubled haploid population and parental lines under normal and low nitrogen cultivation conditions. (A) The box plot displays the effects of a low nitrogen (N) cultivation regime on the growth of 93-11 (P1, indica), Milyang352 (P2, japonica), and that of a doubled haploid (DH) rice population (n = 117) compared to the normal N application, (B) box plot showing the number of tillers per plant, (C) culm length, (D) panicle length, (E) number of panicles per plant, (F) root length, (G) shoot dry weight, (H) root dry weight, (I) root–shoot (R:S) ratio calculated from the root length and plant height, and (J) based on the root dry weight and shoot dry weight of the same population, and (K,L) phenotypes of 93-11 (P1) and Milyang352 (P2) under the same conditions.
Figure 2. Genotype–phenotype correlation of the doubled haploid population and parental lines under normal and low nitrogen cultivation conditions. (A) The box plot displays the effects of a low nitrogen (N) cultivation regime on the growth of 93-11 (P1, indica), Milyang352 (P2, japonica), and that of a doubled haploid (DH) rice population (n = 117) compared to the normal N application, (B) box plot showing the number of tillers per plant, (C) culm length, (D) panicle length, (E) number of panicles per plant, (F) root length, (G) shoot dry weight, (H) root dry weight, (I) root–shoot (R:S) ratio calculated from the root length and plant height, and (J) based on the root dry weight and shoot dry weight of the same population, and (K,L) phenotypes of 93-11 (P1) and Milyang352 (P2) under the same conditions.
Genes 12 00745 g002
Figure 3. Population genome background and principal component analysis results. (A) The density map of pairwise kinship matrix values generated based on 240 KASP markers on a DH population, (B) cluster plot showing the distinctive phenotypic response between P1 (93-11) and P2 (Milyang352) based on the number of tillers under low N cultivation, (C,D) 2D principal component analysis (PCA) indicating the relationship between traits, and (E,F) linear regression analysis proposing the existence of a positive correlation between number of tillers and shoot dry weight in rice grown under a normal and low N cultivation regime, respectively.
Figure 3. Population genome background and principal component analysis results. (A) The density map of pairwise kinship matrix values generated based on 240 KASP markers on a DH population, (B) cluster plot showing the distinctive phenotypic response between P1 (93-11) and P2 (Milyang352) based on the number of tillers under low N cultivation, (C,D) 2D principal component analysis (PCA) indicating the relationship between traits, and (E,F) linear regression analysis proposing the existence of a positive correlation between number of tillers and shoot dry weight in rice grown under a normal and low N cultivation regime, respectively.
Genes 12 00745 g003
Figure 4. Linkage maps and identified QTLs associated with various agronomic traits in rice under normal and low nitrogen cultivation. Five (5) QTLs (qPHN-1, 43 cM; qPHN3-1, 166 cM; qPHN5-1, 23 cM; qPHN10-1, 1 cM; and qPHN12-1, 130 cM) associated with plant height of rice were mapped on five (5) different chromosomes, for normal nitrogen cultivation. The QTLs associated with tiller number were mapped on chromosome 2 (qTNN-1, 130 cM for normal nitrogen cultivation, and qTNL2-1, 134 cM for low nitrogen cultivation) and chromosome 4 (qTNN4-1, 173 cM). For both normal and low nitrogen cultivation conditions, five QTLs (qCLN2-1, 132 cM; qCLN2-2, 187 cM; qCLN3-1, 128 cM; qCLN9-1, 9 cM; and qCLN10-1, 17 cM) and two (2) QTLs (qCLL1-1, 108 cM; and qCLL9-1, 72 cM) were mapped on four and two chromosomes, respectively. Five QTLs (qPLN2-1, 74 cM; qPLN3-1, 123 cM; qPLN3-2, 201 cM; qPLN7-1, 89 cM; and qPLN8-1, 101 cM) associated with panicle length under normal nitrogen cultivation were mapped on four chromosomes. Three QTLs (qPNN2-1, 53 cM; qPNN4-1, 169 cM; qPNN8-1, 84 cM) and one QTL (qPNL4-1, 154 cM) associated with panicle number were mapped on three chromosomes and one chromosome, for normal nitrogen and low nitrogen cultivation, respectively. A unique QTL (qRLL4-1, 59 cM) associated with root length was mapped on chromosome 4. Finally, one QTL (qSDWN-51, 46 cM) and two QTLs (qRDW1-1, 27 cM; and qRDW4-1, 21 cM) associated with shoot dry weight (SDW) and root dry weight, respectively, under normal nitrogen cultivation, were mapped to chromosomes 5 (for SDW) and 1 and 4 (for RDW). (AJ) Chromosomes 1, 2, 3, 4, 5, 7, 8, 9, 10, and 12. LOD score, logarithm of the odds indicating the significance of genetic linkage between markers and the detected QTLs above a threshold.
Figure 4. Linkage maps and identified QTLs associated with various agronomic traits in rice under normal and low nitrogen cultivation. Five (5) QTLs (qPHN-1, 43 cM; qPHN3-1, 166 cM; qPHN5-1, 23 cM; qPHN10-1, 1 cM; and qPHN12-1, 130 cM) associated with plant height of rice were mapped on five (5) different chromosomes, for normal nitrogen cultivation. The QTLs associated with tiller number were mapped on chromosome 2 (qTNN-1, 130 cM for normal nitrogen cultivation, and qTNL2-1, 134 cM for low nitrogen cultivation) and chromosome 4 (qTNN4-1, 173 cM). For both normal and low nitrogen cultivation conditions, five QTLs (qCLN2-1, 132 cM; qCLN2-2, 187 cM; qCLN3-1, 128 cM; qCLN9-1, 9 cM; and qCLN10-1, 17 cM) and two (2) QTLs (qCLL1-1, 108 cM; and qCLL9-1, 72 cM) were mapped on four and two chromosomes, respectively. Five QTLs (qPLN2-1, 74 cM; qPLN3-1, 123 cM; qPLN3-2, 201 cM; qPLN7-1, 89 cM; and qPLN8-1, 101 cM) associated with panicle length under normal nitrogen cultivation were mapped on four chromosomes. Three QTLs (qPNN2-1, 53 cM; qPNN4-1, 169 cM; qPNN8-1, 84 cM) and one QTL (qPNL4-1, 154 cM) associated with panicle number were mapped on three chromosomes and one chromosome, for normal nitrogen and low nitrogen cultivation, respectively. A unique QTL (qRLL4-1, 59 cM) associated with root length was mapped on chromosome 4. Finally, one QTL (qSDWN-51, 46 cM) and two QTLs (qRDW1-1, 27 cM; and qRDW4-1, 21 cM) associated with shoot dry weight (SDW) and root dry weight, respectively, under normal nitrogen cultivation, were mapped to chromosomes 5 (for SDW) and 1 and 4 (for RDW). (AJ) Chromosomes 1, 2, 3, 4, 5, 7, 8, 9, 10, and 12. LOD score, logarithm of the odds indicating the significance of genetic linkage between markers and the detected QTLs above a threshold.
Genes 12 00745 g004
Table 1. Soil mineral elements composition and hydrogen potential.
Table 1. Soil mineral elements composition and hydrogen potential.
FieldpHTN (%)ECOM (%)P2O5 (mg kg−1)K+
(cmolc kg−1)
Ca2+ (cmolc kg−1)Mg2+ (cmolc kg−1)Na+ (cmolc kg−1)
Low N5.50.270.752.7444.770.844.761.010.18
Normal N5.71.381.383.7996.870.271.177.821.27
pH: potential of hydrogen; TN: total nitrogen; EC: cationic exchange; OM: organic matter; P2O5: phosphate; K+: potassium; Ca2+: calcium; Mg2+: magnesium; Na+: sodium, cmolc kg−1: centimoles of charges per kilogram.
Table 2. Detected QTLs associated with growth and development, and yield components of rice under normal and low cultivation.
Table 2. Detected QTLs associated with growth and development, and yield components of rice under normal and low cultivation.
Traits (a)QTL (b)Chr.(c)Position (cM) (d)Left Marker (e)Right Marker (f)LOD (g)PVE (%) (h)Add (i)
PHNqPHN1-1143Os01_id1028304Os01_ad010208244.7038.7452.08
qPHN3-13166Os03_Hd6-1ATOs03_id30154534.8838.798−2.1
qPHN5-1523Os05_KJ05_019Os05_KJ05_0175.6949.915−2.2
qPHN10-1101Os10_id10002069Os10_KJ10_0035.3059.1472.14
qPHN12-112130Os12_cmb1224_0Os12_cmb1226_05.1348.836−2
TNLqTNL2-12134Os02_KJ02_057Os02_id20127736.05220.451.81
TNNqTNN2-12130Os02_KJ02_053Os02_KJ02_0574.14114.511.25
qTNN4-14173Os04_cmb0432_2Os04_cmb0434_13.92214.491.33
CLLqCLL1-11108Os01_id1015984Os01_ah010018433.94710.573.3
qCLL9-1972Os09_KJ09_075Os09_TAC1_CT4.63211.59−2.8
CLNqCLN2-12132Os02_KJ02_057Os02_id20127734.426.509−1.9
qCLN2-22187Os02_id2014575Os02_cmb0236_613.4524.813.77
qCLN3-13128Os03_id3010700Os03_ad0301390511.4519.513.37
qCLN9-199Os09_cmb0909_6Os09_KJ09_0274.5566.805−2.2
qCLN10-11017Os10_KJ10_003Os10_KJ10_0193.1785.4011.89
PLNqPLN2-1274Os02_id2002293Os02_KJ02_0094.8117.985−0.7
qPLN3-13123Os03_id3010700Os03_ad030139053.9176.7680.61
qPLN3-23201Os03_ah03002520Os03_cmb0336_58.84416.060.93
qPLN7-1789Os07_KJ07_021Os07_id70011555.93811.140.78
qPLN8-18101Os08_GW8-AGOs08_id80077648.05914.35−0.9
PNLqPNL4-14154Os04_ad04009559Os04_ah040012522.95411.991.18
PNNqPNN2-1253Os02_id2000007Os02_ad020005123.4877.2690.63
qPNN4-14169Os04_id4009823Os04_cmb0432_29.63323.760.94
qPNN8-1884Os08_ae08007378Os08_id80067514.349.3380.55
RLLqRLL4-1459Os04_KJ04_053Os04_KJ04_0573.24111.99−1.1
SDWNqSDWN5-1546Os05_KJ05_011Os05_id50024973.74812.88−0.6
RDWNqRDWN1-1127Os01_KJ01_125Os01_KJ01_1274.06911.74−0.4
qRDWN4-1421Os04_KJ04_029Os04_KJ04_0473.2188.654−0.3
(a) Evaluated rice traits for which QTLs were detected: plant height under normal nitrogen cultivation (PHN), tiller number under low nitrogen availability (TNL), tiller number under normal nitrogen availability (TNN), culm length under low nitrogen level (CLL), culm length under normal nitrogen level (CNN), panicle length under normal nitrogen cultivation (PLN), panicle number under low nitrogen level (PNL), panicle number under normal nitrogen cultivation (PNN), root length under low nitrogen level (RLL), shoot dry weight under normal nitrogen cultivation (SDWN), and root dry weight under normal nitrogen cultivation; (b) detected QTLs names; (c) chromosome number; (d) absolute position of the QTLs from top of the linkage map in centimorgan (cM); (e) left flanking markers; (f) right flanking markers; (g) logarithm of the odds scores; (h) phenotypic variance explained (PVE) by the QTLs, expressed in percentage; (i) additive effect: the negative value indicates that the allele from Milyang352 (P2) increased the corresponding trait value, while the positive value shows that the allele from 93-11 (P1) increased the trait value.
Table 3. Candidate genes pooled from the qTNL2-1 associated with tiller number under low nitrogen availability.
Table 3. Candidate genes pooled from the qTNL2-1 associated with tiller number under low nitrogen availability.
No.Gene NameLocusAnnotationMolecular FunctionCellular Component
1OsPBF19Os02g44220Peroxisomal biogenesis factor 19, putative, expressedProtein bindingCytosol; peroxisome
2OsSc11/Os02g44360Scarecrow transcription factor family protein, putative, expressed; GRAS: GIBBERELLIN-INSENSITIVE (GAI), REPRESSOR of GA1-3 (RGA), and SCARECROW (SCR)Sequence-specific DNA binding transcription factor activityPlastid, nucleus
3OsGRASOs02g44370Scarecrow, putative, expressed; GRAS transcription factor domain containing protein NSENSITIVE (GAI), REPRESSOR of GA1-3 (RGA), and SCARECROW (SCR)Sequence-specific DNA binding transcription factor activityPlastid, nucleus
4OsGpx3Os02g44500Glutathione peroxidase, putative, expressed; OsGpx3|GPX3|OsGPx03Glutathione peroxidase activity; peroxidase activity; oxidoreductase activityMitochondrion; plasma membrane; plastid; cytosol; chloroplast;
5OsHUGT1Os02g44510UDP-glucose glycoprotein glucosyltransferase 1 precursor, putative, expressedUDP-glucose:glycoprotein glucosyltransferase activityEndoplasmic reticulum
6OsSub19Os02g44520Putative Subtilisin homologue, expressed; Peptidase S8 and S53, subtilisin, kexin, sedolisin domain containing proteinHydrolase activity; serine-type endopeptidase activityVacuole; cytoplasm; ribosome; cytosol; plastid; vacuolar membrane; chloroplast
7OsMalicOs02g44550NADP-dependent malic enzyme, putative, expressedCatalytic activity; malic enzyme activity; oxidoreductase activity, acting on the CH-OH group of donors Cytosol
8 Os02g44570Mitochondrial carrier protein, putative, expressed; Mitochondrial substrate carrier family proteinbindingMitochondrial inner membrane
9OsPIP1;1|RWC1Os02g44630Aquaporin. Plasma membrane intrinsic proteins subfamily Transporter activityPlasma membrane; mitochondrion; chloroplast envelop
10MAP3K.10—STEOs02g44642STE_MEKK_ste11_MAP3K.10—STE kinases include homologs to sterile 7, sterile 11 and sterile 20 from yeast, expressed; Similar to MAP3KaTransferring phosphorus-containing groups; ATP binding; protein kinase activity
11CYP450Os02g44654Cytochrome P450, putative, expressed; Similar to Cytochrome P450Heme binding; oxidoreductase activity
12 Os02g44700Zinc finger, C3HC4 type domain containing protein, expressed; Zinc finger, RING-type domain containing proteinBinding; catalytic activity; zinc ion binding; ubiquitin–protein transferase activity; metal ion binding; protein binding
13OsDhn1/OsLip9Os02g44870Dehydrin, putative, expressed; Dehydrin family protein BindingCytosol; plasma membrane; plastid
14OsHMG1Os02g44930HMG1/2, putative, expressed; Similar to HMGc1 proteinSequence-specific DNA binding transcription factor activity; structural molecule activityIntracellular
15OsRALFL8Os02g44940Rapid Alkalinization Factor RALF family protein precursor, expressed; Rapid Alkalinization Factor family proteinSignal transducer activityExtracellular region
15 Os02g44980Transmembrane amino acid transporter protein, putative, expressed; Similar to amino acid transport proteinTransporter activityMembrane; integral component of membrane
16OsFBDUF1/OsMAIF13Os02g44990F-box and DUF domain containing protein, expressed; Cyclin-like F-box domain containing protein.Protein binding
17bHLHOs02g45010Ethylene-responsive protein related, putative, expressed; Basic helix-loop-helix dimerization region bHLH domain containing proteinSequence-specific DNA binding transcription factor activity; protein dimerization activityNucleus
18SID1|OsIDD4Os02g45054ZOS2-15—C2H2-type domain containing protein, expressed; Zinc fingerNucleic acid and protein binding; sequence-specific DNA binding transcription factor activity; metal ion bindingNucleus
19MYBOs02g45080MYB family transcription factor, putative, expressed, containing a myb_SHAQKYF class DNA-binding domain at the C-terminal end of the motifDNA binding; sequence-specific DNA binding transcription factor activity
20 Os02g45130Protein kinase, putative, expressedprotein serine/threonine kinase activity; ATP bindingNucleus
21OsALMT7Os02g45160Aluminum-activated malate transporter, putative, expressed; Similar to ALMT3Malate transport
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kwon, Y.-H.; Kabange, N.-R.; Lee, J.-Y.; Lee, S.-M.; Cha, J.-K.; Shin, D.-J.; Cho, J.-H.; Kang, J.-W.; Ko, J.-M.; Lee, J.-H. Novel QTL Associated with Shoot Branching Identified in Doubled Haploid Rice (Oryza sativa L.) under Low Nitrogen Cultivation. Genes 2021, 12, 745. https://0-doi-org.brum.beds.ac.uk/10.3390/genes12050745

AMA Style

Kwon Y-H, Kabange N-R, Lee J-Y, Lee S-M, Cha J-K, Shin D-J, Cho J-H, Kang J-W, Ko J-M, Lee J-H. Novel QTL Associated with Shoot Branching Identified in Doubled Haploid Rice (Oryza sativa L.) under Low Nitrogen Cultivation. Genes. 2021; 12(5):745. https://0-doi-org.brum.beds.ac.uk/10.3390/genes12050745

Chicago/Turabian Style

Kwon, Young-Ho, Nkulu-Rolly Kabange, Ji-Yun Lee, So-Myeong Lee, Jin-Kyung Cha, Dong-Jin Shin, Jun-Hyeon Cho, Ju-Won Kang, Jong-Min Ko, and Jong-Hee Lee. 2021. "Novel QTL Associated with Shoot Branching Identified in Doubled Haploid Rice (Oryza sativa L.) under Low Nitrogen Cultivation" Genes 12, no. 5: 745. https://0-doi-org.brum.beds.ac.uk/10.3390/genes12050745

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop