- Research
- Open access
- Published:
Phenotypic variation and genome-wide association studies of main culm panicle node number, maximum node production rate, and degree-days to heading in rice
BMC Genomics volume 23, Article number: 390 (2022)
Abstract
Background
Grain yield is a complex trait that results from interaction between underlying phenotypic traits and climatic, edaphic, and biotic variables. In rice, main culm panicle node number (MCPNN; the node number on which the panicle is borne) and maximum node production rate (MNPR; the number of leaves that emerge per degree-day > 10°C) are primary phenotypic plant traits that have significant positive direct effects on yield-related traits. Degree-days to heading (DDTH), which has a significant positive effect on grain yield, is influenced by the interaction between MCPNN and MNPR. The objective of this research is to assess the phenotypic variation of MCPNN, MNPR, and DDTH in a panel of diverse rice accessions, determine regions in the rice genome associated with these traits using genome-wide association studies (GWAS), and identify putative candidate genes that control these traits.
Results
Considerable variation was observed for the three traits in a 220-genotype diverse rice population. MCPNN ranged from 8.1 to 20.9 nodes in 2018 and from 9.9 to 21.0 nodes in 2019. MNPR ranged from 0.0097 to 0.0214 nodes/degree day > 10°C in 2018 and from 0.0108 to 0.0193 nodes/degree-day > 10°C in 2019. DDTH ranged from 713 to 2,345 degree-days > 10°C in 2018 and from 778 to 2,404 degree-days > 10°C in 2019.
Thirteen significant (P < 2.91 x 10-7) trait-single nucleotide polymorphism (SNP) associations were identified using the multilocus mixed linear model for GWAS. Significant associations between MCPNN and three SNPs in chromosome 2 (S02_12032235, S02_11971745, and S02_12030176) were detected with both the 2018 and best linear unbiased prediction (BLUP) datasets. Nine SNPs in chromosome 6 (S06_1970442, S06_2310856, S06_2550351, S06_1968653, S06_2296852, S06_1968680, S06_1968681, S06_1970597, and S06_1970602) were significantly associated with MNPR in the 2019 dataset. One SNP in chromosome 11 (S11_29358169) was significantly associated with the DDTH in the BLUP dataset.
Conclusions
This study identifies SNP markers that are putatively associated with MCPNN, MNPR, and DDTH. Some of these SNPs were located within or near gene models, which identify possible candidate genes involved in these traits. Validation of the putative candidate genes through expression and gene editing analyses are necessary to confirm their roles in regulating MCPNN, MNPR, and DDTH. Identifying the underlying genetic basis for primary phenotypic traits MCPNN and MNPR could lead to the development of fast and efficient approaches for their estimation, such as marker-assisted selection and gene editing, which is essential in increasing breeding efficiency and enhancing grain yield in rice. On the other hand, DDTH is a resultant variable that is highly affected by nitrogen and water management, plant density, and several other factors.
Background
Grain yield is greatly influenced by interactions between underlying phenotypic traits and environmental variables. As a result, yield cannot be directly selected with a high degree of certainty. Selection must instead focus on traits that impact yield performance. In rice, tiller density and panicle mass-related traits, such as spikelet density, grain size, and grain number have been linked to with grain yield [1]. Genes that are associated with these traits have been identified, such as IPA1, MOC1, and FC1 for tiller number [2,3,4], APO1, DEP1, and Gn1a for grain number [5,6,7], GS3 for grain size [8], and WFP for panicle branching [9]. DNA markers for some of these genes have also been developed for use in trait selection [10, 11]. Studies have also shown that genes controlling heading date may influence grain yield in rice. QTLs for days to heading such as Ghd7, Ghd7.1, Ghd8, Hd1, qHd1, and RFT1 have pleiotropic effects on yield-related traits in rice [12,13,14,15,16,17,18,19]. Interactions between QTLs for heading date also affect rice yield-related traits [20, 21]. Unfortunately, most measured rice phenotypic traits, including those mentioned above, are resultant variables that can be greatly impacted by several climatic, edaphic, biotic, and management variables. The further a resultant trait is from determining yield, the greater its variability and less its predictive value. In contrast, the closer a trait is to a genotype’s underlying yield response, herein referred to as a primary phenotypic trait, the less it is impacted by other variables and the greater its potential use as a selection criterion.
Process-based simulation modeling were used by our team to identify four primary traits that when combined putatively produce an ultra-high-yielding rice ideotype for the Gulf Coast environment in Texas, United States; namely, increased node production rate, increased main culm panicle node, increased leaf mass, and increased spikelet density [22]. These findings were verified in a field experiment, where 10 out of 10 rice genotypes with the combination of all four traits yielded more than the check cultivar, while three out of four genotypes that did not have all these four traits yielded less [23]. Using correlation and path analyses from field experiments, Samonte et al. [24] determined that main culm panicle node number (MCPNN) and maximum node production rate (MNPR) have significant direct effects on degree-days to heading (DDTH), mass per panicle, panicle density, and stem mass at heading and harvest, which are traits that have direct effects on grain yield.
Main culm panicle node number has been identified as an important biomass component since each culm node produces a leaf [25]. The average number of leaves that appear on the main stem per unit of thermal time has been suggested as an important constituent trait that explains genotypic variation of early vigor in rice [26]. This statistic is the inverse of the maximum node production rate when calculation is restricted to early season node production data. Similarly, the phyllochron has been used as an index in rice root and shoot development studies [27] and is suggested as a significant factor in predicting heading in rice [28]. The maximum phyllochron value when calculated prior to when node production begins to decrease as plant biomass rapidly increases is equal to the maximum node production rate when expressed on a heat unit basis.
Genome-wide association studies (GWAS) usually analyze diverse populations or individuals from various geographical locations or origins, as these have the advantage of capturing historical recombination events that occur during development [29]. Rice is an ideal plant species for GWAS, as its self-pollination and its long history of artificial selection have allowed the fixing of favorable alleles with large effects, and the diverse environments where rice is planted has resulted in rice subpopulations with distinctive combinations of traits that are adapted to respective local environments [30].
To understand the genetic basis of primary phenotypic traits MCPNN and MNPR, and resultant variable DDTH, this study was conducted with the following objectives: 1). to assess the phenotypic diversity of these traits in a panel of diverse rice accessions; 2). to determine genomic regions with which these traits are associated; and 3). to identify candidate genes involved with these traits.
Results
Phenotypic Variation in MCPNN, MNPR, and DDTH
Wide variation was observed among the three traits (Table 1). MCPNN ranged from 8.1 to 20.9 nodes in 2018 and from 9.9 to 21.0 nodes in 2019. MNPR ranged from 0.0097 to 0.0214 nodes/degree-day > 10°C in 2018 and from 0.0108 to 0.0193 nodes/degree-day > 10°C in 2019. DDTH ranged from 713 to 2,345 degree-days > 10°C in 2018 and from 778 to 2,404 degree-days > 10°C in 2019. Presidio, one of the check cultivars, had an average of 14.1 nodes on the main culm, 0.0156 nodes-degree day-1, and heading at 1300 degree-days in 2018. In 2019, Presidio had an average of 14.6 nodes on the main culm, 0.0144 nodes-degree day-1, and heading at 1384 degree-days. Some of the rice accessions were consistently in the tails of the trait distributions. The accession C4-63 had the highest MCPNN in 2018 and 2019, while accessions Csornuj and Short Grain were the earliest and latest, respectively, in terms of DDTH in both years. Broad-sense heritability (H2) estimates were high for MCPNN (0.83 in 2018 and 0.90 in 2019) and DDTH (0.99 in 2018 and 0.98 in 2019), while MNPR had low to moderate H2 estimates (0.20 in 2018 and 0.60 in 2019).
The three traits followed a normal or slightly positively skewed distribution (Fig. 1). Based on the BLUPs estimated from the two-years' data, MCPNN and DDTH had a very strong, positive correlation (Pearson correlation coefficient (r) = 0.85), MCPNN and MNPR had a moderately positive correlation (r = 0.34), and MNPR and DDTH had a weak positive correlation (r = 0.097) (Fig. 1). While frequency distribution trends for MCPNN and DDTH were similar in both years, MNPR had a slightly negatively-skewed distribution in 2018 and a normal distribution in 2019. Some variations in correlation coefficients of the three traits were observed in 2018 and 2019 (Supplementary Figure 1). MCPNN and MNPR had moderately positive correlations in 2018 (r = 0.29) and 2019 (r = 0.32), MCPNN and DDTH had a strong positive correlation in 2018 (r = 0.78) and a moderately positive correlation in 2019 (r = 0.34). MNPR and DDTH had a weak positive correlation in 2018 (r = 0.039) and a weak negative correlation in 2019 (r = -0.028). Analysis of variance (ANOVA) showed significant variation due to year and genotype effects for MCPNN and DDTH. For MNPR, significant variation due to year but no significant effect for genotype (P = 0.1509) was observed (Table 2). A response surface analysis shows both MCPNN and MNPR have highly significant effects (P < 0.0001) on DDTH [DDTH = 1475.34 + (147.92*MCPNN) - (148512.8*MNPR)], explaining 77.54% of the variation (Figure 2).
Genome-wide Association Studies and Identification of Candidate Genes
GWAS was conducted in the 220 accessions with both genotype and phenotype data. Because there was a significant effect of year on all three traits (Table 2), GWAS were conducted separately for the 2018 and 2019 data, as well as the BLUP estimates based on the data from both years. Principal component analysis (PCA) was used to determine population structure. The top four principal components (PCs) explained 56.7% of the genetic variation, wherein the first PC (38.9% of the variation) divided most of the indicas and japonicas into separate groups, while the second PC (12.0% of the variation) further separated the japonica group into temperate japonicas and tropical japonicas. The admixed accessions can be found in between these three major subgroups (Fig. 3).
Genome-wide linkage disequilibrium (LD) decay was estimated to be at ~150,000 bp. The results are within the range estimated in previous findings, ranged from close to 100,000 bp to over 200,000 bp [31,32,33,34].
Thirteen significant trait-SNP associations were identified (Table 3 and Figure 4). Using the genome browser in RAP-DB (Nipponbare IRGSP Build 5), genes or gene models that include or are within 100 kilobase (kb) pairs to the SNPs declared to have significant associations with MCPNN, MNPR, and DDTH were identified (Supplementary Table 2). A total of 19 candidate genes were found to be in LD with the three SNPs significantly associated with MCPNN, 101 with the nine SNPs significantly associated with MNPR, and 21 with the SNP significantly-associated with DDTH.
Significant association (P < 2.91 x 10-7, after multiple testing) between MCPNN and three SNPs in chromosome 2 (S02_12032235, S02_11971745, and S02_12030176) were detected in the 2018 and BLUP datasets. The SNP S02_12030176 is within the Oryza sativa nicotianamine aminotransferase 1 (OsNAAT1) gene [35], in chromosome 2. The closest gene model to the SNP S02_11971745 is Os02g0305600, described as a spectrin repeat-containing protein. Nine SNPs in chromosome 6 (S06_1970442, S06_2310856, S06_2550351, S06_1968653, S06_2296852, S06_1968680, S06_1968681, S06_1970597, and S06_1970602) were significantly associated with MNPR in the 2019 dataset. The SNP S06_2296852 is within the SULTR-like phosphorus distribution transporter (SPDT) gene, also known as Oryza sativa sulphate transporter 3;4 (OsSULTR3;4) [36, 37]. Three other SNPs associated with MNPR, S06_1970442, S06_1970597, and S06_1970602, are close to the candidate gene Os06g0137100 in the RAP-DP genome browser. The counterpart of Os06g0137100 in the MSU-RGAP genome browser [38], LOC_Os06g04560, is described as a kinesin motor domain-containing protein. In addition, candidate genes related to plant development, cell division and elongation, as well as protein and carbohydrate metabolism, were located within 100,000 bp of the identified SNPs for MCPNN and MNPR. One SNP in chromosome 11 (S11_29358169) was significantly associated with DDTH in the BLUP dataset. The closest candidate genes are Os11g0672300, which is described as being similar to a kinase domain-containing protein, and Os11g0672400, a calponin homology domain containing protein in chromosome 11.
Discussion
Phenotypic Analyses
Phenotyping MCPNN and MNPR is relatively time-consuming, compared to the typical measuring of resultant variables such as plant height, days to heading, days to peak flowering, and other secondary traits measured in most rice breeding programs. However, to state measuring a new trait is “time-consuming” is a value judgement based on a paradigm that assumes shallow phenotyping of large numbers of genotypes is the most efficient way to eliminate lines because they are suspected to result in poor yield performance. It might in fact be more cost-effective to spend more time measuring traits that have a greater influence on yield performance and possibly even evaluate fewer genotypes and thereby decrease the average amount of effort required to release a higher yielding genotype. A major goal of the current research is to identify relationships between the measured traits, and to develop markers to more efficiently select these traits.
There is wide variation in MCPNN, MNPR, and DDTH across the rice accessions. The standard deviations for MCPNN and MNPR are comparable to those of Samonte et al. [24], while that of DDTH was higher. The ranges for MNPR and the inverse of development rate observed by Rebolledo et al. [39], were similar to the estimates we obtained.
Broad-sense heritability (H2) estimates for the three traits ranged from 0.20 for MNPR (2018) to 0.99 for DDTH (2019). Modifying the planting method in 2019 (i.e., using jiffy pots and not thinning the plants before observation as was conducted in 2018) may have reduced plant stress, which may have improved estimation of the genotype effects, and thereby improved the H2 estimates for MNPR in 2019. However, because plant density impacts plant to plant competition, not thinning the plants could have also contributed to the observed differences comparing the years. The timing of N fertilization may have also affected the heritability estimated for MNPR, where two splits were done in 2018 (at planting and six weeks after planting), and three splits were done in 2019 (at planting, two weeks after planting, and 11 weeks after planting). Khing et al. [40] observed that phyllochrons are shortened (i.e., faster node production rate) with higher N rates or wider spacing. Martinez-Eixarch et al. [41] also reported that plant density and the timing of N fertilization, in combination with water management, affect the rate of leaf emergence. The total N fertilizer applied in both years were similar, hence, MCPNN and DDTH, which are considered as end-point traits, had consistently high heritability estimates.
The correlations between MCPNN, MNPR, and DDTH in the BLUP estimates from this study confirm the findings of Samonte et al. [24]. The pairwise correlation coefficients among MCPNN and MNPR (r = 0.34 in this study, and r = 0.35 in [24]) as well as MCPNN and DDTH (r = 0.85 in this study, and r = 0.86 [24]) were very close. The correlation between MNPR and DDTH was positive but low in both studies.
Samonte et al.[24] reported MCPNN and MNPR have significant positive and negative direct effects, respectively, on DDTH based on path analysis. In this study, response surface regression shows both MCPNN and MNPR have highly significant effects on DDTH (R2 = 0.78), and follow the same trend as reported by Samonte et al. [24]. The high degree of explained variability suggests the model can accurately predict DDTH relatively early in a genotype’s development.
Genome-Wide Association Studies and Identification of Candidate Genes Related to MCPNN and MNPR
SNP-trait associations for MCPNN, MNPR, and DDTH were determined using the multi-locus mixed model approach, which included PCA for population structure and kinship matrix to minimize false associations and stepwise linear mixed-model regression with significantly-associated markers as cofactors. MCPNN was significantly associated with three SNPs in chromosome 2, which were detected in the BLUP dataset, which uses the data from both 2018 and 2019. In the individual experiments, only S02_12032235 is significantly associated with MCPNN. In the 2019 dataset, while no significantly associated SNPs were detected, these three SNPs were in the top 20, with the peak at S02_12032235. The OsNAAT1 gene harbors the SNP S02_12030176 and is in high linkage disequilibrium (r2 = 0.91) with S02_12032235. OsNAAT1 plays a role in iron (Fe) homeostasis in rice, where it is involved in Fe acquisition from the soil and transport within the plant [35, 42]. OsNAAT1 also flanks a QTL for grain yield in chromosome 2 in rice [43]. The gene model Os02g0305600 is described as a spectrin repeat-containing protein and is 983 kb from S02_11971745. Spectrin-like proteins have been detected in the plant nuclei and plasma membranes and are suggested to perform multiple roles, such as stability and ion transport [44,45,46].
MNPR is significantly associated with S06_2296852, which lies within the SPDT (OsSULTR3;4) gene. SPDT is found to be highly expressed in node 1, which connects the panicle to the flag leaf, and serves as a switch to distribute phosphorus to the grains [37]. Three other SNPs significantly associated with MNPR, S06_1970442, S06_1970597, and S06_1970602 are in high linkage disequilibrium (r2 = 0.73) with the gene model Os06g0137100. This gene model is similar to a predicted kinesin-like protein. Kinesins are microtubule motor proteins involved in cell division and growth. In rice, genes encoding kinesins confer traits related to plant height, grain length, and shape, and pollen partial sterility [47,48,49,50,51,52,53]. The gene model Os06g0137100 is also associated with grain number based on GWAS conducted by Huang et al. [33].
Some of the candidate genes within 100 kb of SNPs (Supplementary Table 2) are significantly associated with MCPNN and MNPR and are of interest as these have roles in plant development, cell division and elongation, as well as protein and carbohydrate metabolism. For instance, the gene models Os02g0305950 and Os06g0137400, which are close to the significantly-associated SNPs for MCPNN and MNPR, respectively, are putative small auxin-up RNA (SAUR) genes. SAUR genes encode auxin-responsive proteins that are primarily expressed in elongating tissues, and may play a role in regulating cell elongation [54].
S11_29358169, a SNP significantly associated with DDTH in the BLUP dataset, is in linkage disequilibrium (r2 = 0.96) with the gene model Os11g0672300, which is similar to a protein kinase domain-containing protein but needs to be characterized further to determine its specific function. In plants, the protein kinase superfamily is very broad and diverse, with roles in metabolic signaling, stress response, cell division regulation, and plant-specific functions, such as flowering [55].
Genes, QTLs, and chromosomal regions identified to be associated with MCPNN, MNPR, and DDTH in the literature were checked to determine if they overlap with the SNPs identified in this study. While none of the SNPs are co-located with reported or identified genes, some are close to previously reported QTL associated with heading date, such as Hd17, which is located at a 2.2 Mb region in chromosome 6 [56] and qHD11.3, located in a 28.7 Mb region in chromosome 11 [57]. Heuer et al. [58] reported that ZmMADS3 in maize is involved in node number. In this study, there is no significant SNP-trait association in the region where the closest homolog of ZmMADS3 in rice, OsMADS15 in chromosome 7, is located. For MNPR, two regions in Chromosome 7 are significantly associated with node development rate [39], while QTLs for phyllochron have been mapped in chromosomes 1, 2, 4, 9, 10, and 11 [59, 60]. However, these regions do not contain any SNPs that are significantly associated with MNPR in this study.
Conclusions
SNP markers significantly associated with MCPNN, MNPR, and DDTH were detected by association analyses using the multilocus linear mixed model. The identified SNPs are located within or in LD with gene models, which could be potential candidate genes for these traits. Further validation of the candidate genes through expression and gene editing analyses are needed to confirm if they are the causal agents for these traits. Identifying the genes regulating MCPNN and MNPR could lead to the development of fast and efficient approaches for their estimation, such as marker-assisted breeding or gene editing. Functional markers can be designed and could be used to select for early-generation breeding lines possessing higher MCPNN and faster MNPR. Gene editing can also be applied to alter the MCPNN and MNPR of a rice variety in order to achieve optimum yields. Understanding the underlying genetic basis for MCPNN, MNPR is an essential step in increasing breeding efficiency for enhancing grain yield in rice.
Methods
Plant Materials
A total of 220 rice accessions consisting of a diverse set of indica and japonica cultivars, landraces, inbred lines, and hybrids (Supplementary Table 1) were selected for their variation in degree-days to heading. These were planted at the Texas A&M AgriLife Research Center (Beaumont) in 2018 and 2019.
Field Experiment Set-up
In 2018, the entries were drill-seeded on April 19 in three-row plots that were 2.44 m long, with rows spaced 0.25 m apart, using a randomized complete block design with 4 replications. A one-meter segment in the middle row of each plot was marked with flags and thinned to 15 plants, for data collection. Urea was applied in two splits: 108.4 kg ha-1 N at planting and 128.9 kg ha-1 N at 6 weeks after planting.
In 2019, the entries were sown on April 16-17 in three-row plots that were 2.44 m long, with rows spaced 0.28 m apart. The slight change in row spacing for the two years is due to the different planters used in each year. The side rows of each plot were mechanically drill-seeded, and the middle rows were manually planted at the same time as the drill-seeding. The middle row of each plot was created by digging a 2.44-m long furrow (approximately 6 cm deep), and a strip of eighteen 5.7 cm square x 5.7 cm deep biodegradable Jiffy pots (Jiffy International; www.jiffygroup.com) each with one seed placed in the center of each row. The Jiffy pots were marked with flags, while the ends of the middle rows were seeded with the designated entry using a manual seeder. This was done to eliminate thinning the plants before data collection. The field was flash-flooded on April 17, after the drill-seeding and manual planting were completed for all accessions. Fifteen of the 18 plants were used for data collection, with the remaining three plants as backup. Urea was applied in 3 splits for a more efficient distribution of N: 59.4 kg N ha-1 at planting (17 April 2019), 128.9 kg N ha-1 at two weeks after planting (23 May 2019), and 47.1 kg N ha-1 at 11 weeks after planting (3 July 2019).
Phenotyping for MCPNN, MNPR, and DDTH
Data collection for MNPR and MCPNN are described in Appendix 1. The emergence of leaves on the main culms of 15 plants within a designated one-meter length was recorded every week, starting at the third leaf stage. MNPR was estimated through regression of the average leaf emergence data during the 3rd to 7th leaf stages with cumulative degree-days > 10℃ from planting, with degree-day accumulation estimated from ambient temperature recorded within 0.4 km of the research site. In 2018, MNPR was estimated using leaf count data starting at two weeks after thinning to minimize the effects of stress caused by thinning, while in 2019, thinning was not necessary because the plants were already spaced with the use of the Jiffy pots. Leaf counting was continued until heading, where MCPNN was estimated as the number of leaves on the main culm plus one for its panicle [24]. Days to 50% heading for each entry was estimated to occur when 50% of the panicles have exerted from tillers and used to estimate degree-days from sowing to 50% heading.
Phenotypic Data Analyses
The combined analysis of both years of data was conducted using SAS Version 9.4 (SAS 2016) for Windows. Fixed and random effects were estimated using PROC MIXED in SAS software, and these were used to estimate best linear unbiased predictions (BLUPs) in R version 3.6.1 [61]. Analyses of variance and broad-sense heritability for the measured variables for each year were estimated using the 'augmentedRCBD' package [62] in R, with Pearson's pairwise correlation coefficients also calculated using R. The relationships between MCPNN and MNPR and the response variable DDTH was estimated using the 'rsm' package [63] in R.
Marker Data
DNA was extracted from leaf tissues of 220 rice accessions (Supplementary Table 1), collected from the field at the Texas A&M AgriLife Research Center (Beaumont) in July 2018. Extraction was performed using standard protocol for leaf tissue with the Thermo Fisher Scientific KingFisher Flex (Thermo Fisher Scientific, Waltham, MA, USA). The DNA samples were sent to the Texas A&M AgriLife Genomics and Bioinformatics Service (TxGen) at College Station for genotyping-by-sequencing (GBS), with 1X coverage. The reference genome used was Oryza sativa ssp. japonica cultivar Nipponbare, International Rice Genome Sequencing Project (IRGSP) Build 5 [38]. The raw genotype data was filtered, selecting single nucleotide polymorphisms (SNPs) having less than 50% missing data and minimum allele frequency (MAF) >5%. After initial filtering, imputation was conducted using BEAGLE V4.0 [64] in 1,075,302 SNP markers. The genotype data was filtered a second time after imputation using TASSEL 5.2.61 [65], and 854,832 SNPs were used in the analyses after removing SNPs with less than 5% MAF and more than 5% missing data. Four genotypes were eliminated from the analysis as outliers due to the ratio of 2018 to 2019 MNPR rates being greater than the 75th quartile or less than the 25th quartile by more than 1.5 x (75th quartile - 25th quartile) [66].
Genome-wide Association Studies
GWAS for MCPNN, MNPR, and DDTH were conducted using the phenotype data for 2018 and 2019 separately, as well as the BLUPs estimated from data collected in both years. Factors that may cause false trait-SNP associations (i.e., population structure (Q) and genetic relatedness (K)) were controlled using principal component analysis (PCA) and kinship matrix, respectively. The most probable number of subpopulations was determined by plotting the number of principal components (PC) against the variance explained by the PCs, and the optimum number of PCs was selected when the decrease in variance has reached a plateau. In this study, the total number of PCs used to account for population structure was four. The VanRaden kinship algorithm [67] was used to construct a kinship matrix. Both PCA and kinship matrices were generated using the R package Genomic Association and Prediction Integrated Tool (GAPIT), Version 3 [68, 69]. Linkage disequilibrium decay was used to estimate the appropriate resolution for association mapping, with a sliding window of 50 markers calculated using TASSEL 5.2.61. Mean r2 was then computed every 10,000 base pairs (bp), with linkage disequilibrium decay determined as the distance in bp wherein the average r2 decreased to half its maximum value.
Association analyses were conducted using the multilocus mixed model (MLMM) [70] implemented in GAPIT Version 3. MLMM uses stepwise linear mixed-model regression that includes significantly associated markers as cofactors, with Q and K to account for false positives. Multiple testing was accounted for using the statistical program 'simpleM' [71, 72] implemented in R, which calculates the number of independent tests (Meff_G). The Meff_G was used to compute for the multiple testing threshold in a similar way as the Bonferroni correction method in which the significant threshold (α= 0.05) was divided by the Meff_G or P = α/ Meff_G. For this study, the multiple testing threshold to declare significant association was set to P = 2.91 x 10-7. Manhattan and quantile-quantile (Q-Q) plots were constructed using the R package 'qqman' [73]. Identification of genes that contain the significant SNPs was achieved using the Nipponbare IRGSP Build 5 genome browser (https://rapdb.dna.affrc.go.jp/viewer/gbrowse/build5/) in the Rice Annotation Project Database (RAP-DB) [74] and Michigan State University (MSU) Rice Genome Annotation Project (RGAP) genome browser (http://rice.plantbiology.msu.edu/cgi-bin/gbrowse/rice/) [38].
Availability of data and materials
All the data supporting the results of this article are provided within the article or in the additional files. The genotyping data has been deposited at Dryad (https://doi.org/10.5061/dryad.4qrfj6qbs).
References
Samonte SOP, Wilson LT, McClung AM. Path analyses of yield and yield-related traits of fifteen diverse rice genotypes. Crop Sci. 1998;38(5):1130–6.
Jiao Y, Wang Y, Xue D, Wang J, Yan M, Liu G, et al. Regulation of OsSPL14 by OsmiR156 defines ideal plant architecture in rice. Nat Genet. 2010;42(6):541–4.
Li X, Qian Q, Fu Z, Wang Y, Xiong G, Zeng D, et al. Control of tillering in rice. Nature [Internet]. 2003;422(6932):618–21. Available from: https://doi.org/10.1038/nature01518
Minakuchi K, Kameoka H, Yasuno N, Umehara M, Luo L, Kobayashi K, et al. FINE CULM1 (FC1) works downstream of strigolactones to inhibit the outgrowth of axillary buds in rice. Plant Cell Physiol. 2010;51(7):1127–35.
Ashikari M, Sakakibara H, Lin S, Yamamoto T, Takashi T, Nishimura A, et al. Plant science: Cytokinin oxidase regulates rice grain production. Science (1979). 2005;309(5735):741–5.
Huang X, Qian Q, Liu Z, Sun H, He S, Luo D, et al. Natural variation at the DEP1 locus enhances grain yield in rice. Nat Genet. 2009;41(4):494–7.
Ikeda K, Ito M, Nagasawa N, Kyozuka J, Nagato Y. Rice ABERRANT PANICLE ORGANIZATION 1, encoding an F-box protein, regulates meristem fate. Plant J. 2007;51(6):1030–40.
Fan C, Xing Y, Mao H, Lu T, Han B, Xu C, et al. GS3, a major QTL for grain length and weight and minor QTL for grain width and thickness in rice, encodes a putative transmembrane protein. Theor Appl Genet. 2006;112(6):1164–71.
Miura K, Ikeda M, Matsubara A, Song XJ, Ito M, Asano K, et al. OsSPL14 promotes panicle branching and higher grain productivity in rice. Nat Genet. 2010;42(6):545–9.
Kim SR, Ramos J, Ashikari M, Virk PS, Torres EA, Nissila E, Hechanova SL, Mauleon R, Jena KK. Development and validation of allele-specific SNP/indel markers for eight yield-enhancing genes using whole-genome sequencing strategy to increase yield potential of rice, Oryza sativa L. Rice. 2016;9(1):1–17.
Kim SR, Ramos JM, Hizon RJ, Ashikari M, Virk PS, Torres EA, Nissila E, Jena KK. Introgression of a functional epigenetic OsSPL14WFP allele into elite indica rice genomes greatly improved panicle traits and grain yield. Scientific reports. 2018;8(1):1–12.
Xue W, Xing Y, Weng X, Zhao Y, Tang W, Wang L, et al. Natural variation in Ghd7 is an important regulator of heading date and yield potential in rice. Nat Genet. 2008;40(6):761–7.
Wei X, Xu J, Guo H, Jiang L, Chen S, Yu C, et al. DTH8 suppresses flowering in rice, influencing plant height and yield potential simultaneously. Plant Physiol. 2010;153(4):1747–58.
Yan W, Liu H, Zhou X, Li Q, Zhang J, Lu L, et al. Natural variation in Ghd7.1 plays an important role in grain yield and adaptation in rice. Cell Res. 2013;23:969–71.
Yan WH, Wang P, Chen HX, Zhou HJ, Li QP, Wang CR, et al. A major QTL, Ghd8, plays pleiotropic roles in regulating grain productivity, plant height, and heading date in rice. Mol Plant. 2011;4(2):319–30.
Zhang ZH, Wang K, Guo L, Zhu YJ, Fan YY, Cheng SH, Zhuang JY. Pleiotropism of the photoperiod-insensitive allele of Hd1 on heading date, plant height and yield traits in rice. PloS one. 2012;7(12):e52538.
Chen JY, Zhang HW, Zhang HL, Ying JZ, Ma LY, Zhuang JY. Natural variation at qHd1 affects heading date acceleration at high temperatures with pleiotropism for yield traits in rice. BMC plant biology. 2018;18(1):1–11.
Chen JY, Guo L, Ma H, Chen YY, Zhang HW, Ying JZ, et al. Fine mapping of qHd1, a minor heading date QTL with pleiotropism for yield traits in rice (Oryza sativa L.). Theor Appl Genet. 2014;127(11):2515–24.
Zhu YJ, Fan YY, Wang K, Huang DR, Liu WZ, Ying JZ, Zhuang JY. Rice Flowering Locus T 1 plays an important role in heading date influencing yield traits in rice. Scientific reports. 2017;7(1):1–10.
Endo-Higashi N, Izawa T. Flowering time genes heading date 1 and early heading date 1 together control panicle development in rice. Plant Cell Physiol. 2011;52(6):1083–94.
Zhang ZH, Zhu YJ, Wang SL, Fan YY, Zhuang JY. Importance of the interaction between heading date genes Hd1 and Ghd7 for controlling yield traits in rice. International journal of molecular sciences. 2019;20(3):516.
Wilson L, Wu G, Samonte O, McClung A, Park W, Pinson S, Stansel JW. Identifying optimal phenotypic trait sets using physiologically‐based modeling. In: Proceedings of the 27th Rice Technical Working Group, Texas Agricultural Experiment Station, College Station, Texas. 1998;70.
Wilson LT, Wu G, Pinson S, McClung A. Rice Physiological Modeling. Texas Rice Special Section, Highlighting Research. Texas A&M University System Agricultural Research and Extension Center, Beaumont, Texas. 2001;13. Available online at https://beaumont.tamu.edu/eLibrary/Newsletter/2001_Highlights_in_Research.pdf.
Samonte SOPB, Wilson LT, Tabien RE. Maximum node production rate and main culm node number contributions to yield and yield-related traits in rice. Field Crop Res. 2006;96(2–3):313–9.
Eizenga GC, Edwards JD, Yeater KM, McCouch SR, McClung AM. Transgressive variation for yield components measured throughout the growth cycle of jefferson rice (Oryza sativa) × O. rufipogon introgression lines. Crop Science. Crop Science. 2016;56(5):2336–47.
Rebolledo MC, Dingkuhn M, Péré P, Mcnally KL, Luquet D. Developmental Dynamics and Early Growth Vigour in Rice. I. Relationship Between Development Rate (1/Phyllochron) and Growth. Journal of Agronomy and Crop Science. 2012;198(5):374–84.
Nemoto K, Morita S, Baba T. Shoot and Root Development in Rice Related to the Phyllochron. Crop Sci. 1995;35(1):24–9.
Itoh Y, Sano Y. Phyllochron dynamics under controlled environments in rice (Oryza sativa L.). Euphytica. 2006;150(1–2):87–95.
Lipka AE, Kandianis CB, Hudson ME, Yu J, Drnevich J, Bradbury PJ, et al. From association to prediction: Statistical methods for the dissection and selection of complex traits in plants. Curr Opin Plant Biol. 2015;24:110–8.
McCouch SR, Wright MH, Tung CW, Maron LG, McNally KL, Fitzgerald M, et al. Open access resources for genome-wide association mapping in rice. Nat Commun. 2016;4:7.
Mather KA, Caicedo AL, Polato NR, Olsen KM, McCouch S, Purugganan MD. The extent of linkage disequilibrium in rice (Oryza sativa L.). Genetics. 2007;177(4):2223–32.
McNally KL, Childs KL, Bohnert R, Davidson RM, Zhao K, Ulat VJ, et al. Genomewide SNP variation reveals relationships among landraces and modern varieties of rice. Proceedings of the National Academy of Sciences [Internet]. 2009 Jul 28;106(30):12273. Available from: http://www.pnas.org/content/106/30/12273.abstract
Huang X, Zhao Y, Wei X, Li C, Wang A, Zhao Q, et al. Genome-wide association study of flowering time and grain yield traits in a worldwide collection of rice germplasm. Nat Genet. 2012;44(1):32–9.
Courtois B, Audebert A, Dardou A, Roques S, Ghneim-Herrera T, Droc G, Frouin J, Rouan L, Gozé E, Kilian A, Ahmadi N. Genome-wide association mapping of root traits in a japonica rice panel. PloS one. 2013;8(11):e78037.
Inoue H, Takahashi M, Kobayashi T, Suzuki M, Nakanishi H, Mori S, et al. Identification and localisation of the rice nicotianamine aminotransferase gene OsNAAT1 expression suggests the site of phytosiderophore synthesis in rice. Plant Mol Biol. 2008;66(1–2):193–203.
Buchner P, Takahashi H, Hawkesford MJ. Plant sulphate transporters: co-ordination of uptake, intracellular and long-distance transport. Journal of experimental botany. 2004;55(404):1765–73.
Yamaji N, Takemoto Y, Miyaji T, Mitani-Ueno N, Yoshida KT, Ma JF. Reducing phosphorus accumulation in rice grains with an impaired transporter in the node. Nature. 2017;541(7635):92–5.
Kawahara Y, de la Bastide M, Hamilton JP, Kanamori H, Mccombie WR, Ouyang S, et al. Improvement of the oryza sativa nipponbare reference genome using next generation sequence and optical map data. Rice. 2013;6(1):3–10.
Rebolledo MC, Dingkuhn M, Courtois B, Gibon Y, Clément-Vidal A, Cruz DF, et al. Phenotypic and genetic dissection of component traits for early vigour in rice using plant growth modelling, sugar content analyses and association mapping. J Exp Bot. 2015;66(18):5555–66.
Khing EE, Kyaw Win K, Min D, Win S. Effects of Different Nitrogen Fertilizer Levels and Spacing on Phyllochron, Yield and Yield Components of Rice (Oryza sativa L.) under System of Rice Intensification (SRI). J Agric Res. 2019;6:70–9.
Martínez-Eixarch M, Zhu DF, Catalá-Forner M del M, Pla-Mayor E, Tomás-Navarro N. Water, Nitrogen and Plant Density Affect the Response of Leaf Appearance of Direct Seeded Rice to Thermal Time. Rice Science. 2013;20(1):52–60.
Nozoye T, Inoue H, Takahashi M, Ishimaru Y, Nakanishi H, Mori S, et al. The expression of iron homeostasis-related genes during rice germination. Plant Mol Biol. 2007;64(1–2):35–47.
Kotla A, Agarwal S, Yadavalli VR, Vinukonda VP, Chakravarthi Dhavala VN, Neelamraju S. Quantitative trait loci and candidate genes for yield and related traits in Madhukar x Swarna RIL population of rice. J Crop Sci Biotechnol. 2013;16(1):35–44.
Faraday CD, Spanswick RM. Evidence for a membrane skeleton in higher plants: A spectrin‐like polypeptide co‐isolates with rice root plasma membranes. FEBS letters. 1993;318(3):313–6.
de Ruijter NCA, Ketelaar T, Blumenthal SSD, Emons AMC, Schel JHN. Spectrin-like proteins in plant nuclei. Cell Biol Int. 2000;24(7):427–38.
Pérez-Munive C, de la Espina SMD. Nuclear spectrin-like proteins are structural actin-binding proteins in plants. Biol Cell. 2011;103(3):145–57.
Fang J, Yuan S, Li C, Jiang D, Zhao L, Peng L, et al. Reduction of ATPase activity in the rice kinesin protein Stemless Dwarf 1 inhibits cell division and organ development. Plant J. 2018;96(3):620–34.
Kitagawa K, Kurinami S, Oki K, Abe Y, Ando T, Kono I, et al. A novel kinesin 13 protein regulating rice seed length. Plant Cell Physiol. 2010;51(8):1315–29.
Li J, Jiang J, Qian Q, Xu Y, Zhang C, Xiao J, et al. Mutation of rice BC12/GDD1, which encodes a kinesin-like protein that binds to a GA biosynthesis gene promoter, leads to dwarfism with impaired cell elongation. Plant Cell. 2011;23(2):628–40.
Ran Q, Akhter D, Chengcong Y, Nath UK, Eshag J, Xiaoli J, et al. SRG1, Encoding a Kinesin-4 Protein, Is an Important Factor for Determining Grain Shape in Rice. Rice Sci. 2018;25(6):297–307.
Wu T, Shen Y, Zheng M, Yang C, Chen Y, Feng Z, et al. Gene SGL, encoding a kinesin-like protein with transactivation activity, is involved in grain length and plant height in rice. Plant Cell Rep. 2014;33(2):235–44.
Zhang M, Zhang B, Qian Q, Yu Y, Li R, Zhang J, et al. Brittle Culm 12, a dual-targeting kinesin-4 protein, controls cell-cycle progression and wall properties in rice. Plant J. 2010;63(2):312–28.
Zhou S, Wang Y, Li W, Zhao Z, Ren Y, Wang Y, et al. Pollen semi-sterility1 encodes a kinesin-1-like protein important for male meiosis, anther dehiscence, and fertility in rice. Plant Cell. 2011;23(1):111–29.
Gil P, Green PJ. Regulatory activity exerted by the SAUR-AC1 promoter region in transgenic plants, vol. 34. Plant Molecular Biology: Kluwer Academic Publishers; 1997.
Lehti-Shiu MD, Shiu SH. Diversity, classification and function of the plant protein kinase superfamily. Philosophical Transactions of the Royal Society B: Biological Sciences. 2012;367(1602):2619–39.
Matsubara K, Kono I, Hori K, Nonoue Y, Ono N, Shomura A, et al. Novel QTLs for photoperiodic flowering revealed by using reciprocal backcross inbred lines from crosses between japonica rice cultivars. Theor Appl Genet. 2008;117(6):935–45.
Yu SB, Li JX, Xu CG, Tan YF, Li XH, Zhang Q. Identification of quantitative trait loci and epistatic interactions for plant height and heading date in rice. Theor Appl Genet. 2002;104(4):619–25.
Heuer S, Hansen S, rg Bantin J, Brettschneider R, Kranz E, Lö rz H, et al. The Maize MADS Box Gene ZmMADS3 Affects Node Number and Spikelet Development and Is Co-Expressed with ZmMADS1 during Flower Development, in Egg Cells, and Early Embryogenesis 1 [Internet]. Vol. 127, Plant Physiology. 2001. Available from: www.plantphysiol.org
Miyamoto N, Goto Y, Matsui M, Ukai Y, Morita M, Nemoto K. Quantitative trait loci for phyllochron and tillering in rice. Theor Appl Genet. 2004;109(4):700–6.
Morita M, Tang DQ, Miyamoto N, Goto Y, Ukai Y, Nemoto K. Quantitative trait loci for rice phyllochron in lemont x IR36 cross. Plant Production Science. 2005;8(2):199–202.
R Core Team. R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2021.
Aravind J, Mukesh Sankar S, Wankhede DP, Kaur V. Augmented RCBD: Analysis of augmented randomised complete block designs. R package version 0.1. 2020;2.
Lenth R v. Response-Surface Methods in R, Using rsm [Internet]. Vol. 32, JSS Journal of Statistical Software. 2009. Available from: http://www.jstatsoft.org/
Browning SR, Browning BL. Rapid and accurate haplotype phasing and missing-data inference for whole-genome association studies by use of localized haplotype clustering. Am J Hum Genet. 2007;81(5):1084–97.
Bradbury PJ, Zhang Z, Kroon DE, Casstevens TM, Ramdoss Y, Buckler ES. TASSEL: Software for association mapping of complex traits in diverse samples. Bioinformatics. 2007;23(19):2633–5.
Tukey JW. Exploratory Data Analysis. Reading: Addison-Wesley; 1977.
VanRaden PM. Efficient methods to compute genomic predictions. J Dairy Sci. 2008;91(11):4414–23.
Lipka AE, Tian F, Wang Q, Peiffer J, Li M, Bradbury PJ, et al. GAPIT: Genome association and prediction integrated tool. Bioinformatics. 2012;28(18):2397–9.
Wang J, Zhang Z. GAPIT Version 3: Boosting Power and Accuracy for Genomic Association and Prediction. Genomics, Proteomics & Bioinformatics. 2021;19(4):629–40.
Segura V, Vilhjálmsson BJ, Platt A, Korte A, Seren Ü, Long Q, et al. An efficient multi-locus mixed-model approach for genome-wide association studies in structured populations. Nat Genet. 2012;44(7):825–30.
Gao X, Becker LC, Becker DM, Starmer JD, Province MA. Avoiding the high bonferroni penalty in genome-wide association studies. Genet Epidemiol. 2010;34(1):100–5.
Johnson RC, Nelson GW, Troyer JL, Lautenberger JA, Kessing BD, Winkler CA, O'Brien SJ. Accounting for multiple comparisons in a genome-wide association study (GWAS). BMC genomics. 2010;11(1):1–6.
Turner S. qqman: an R package for visualizing GWAS results using Q-Q and manhattan plots. Journal of Open Source Software. 2018;3(25):731.
Sakai H, Lee SS, Tanaka T, Numa H, Kim J, Kawahara Y, Wakimoto H, Yang CC, Iwamoto M, Abe T, Yamada Y. Rice Annotation Project Database (RAP-DB): an integrative and interactive database for rice genomics. Plant and Cell Physiology. 2013;54(2):e6.
Acknowledgements
Chersty Harper, Leon Holgate, Eric Christensen, and Leanna Martin are acknowledged for their technical assistance in conducting this study.
Funding
The authors appreciate the funding provided by the Texas Rice Research Foundation (124214-94410 and 114214-94410) and Texas A&M AgriLife Research (203463-94410).
Author information
Authors and Affiliations
Contributions
DLS and SOPBS conceptualized the study, designed the experiments, and wrote the original draft. DLS, SOPBS, YY, JBBA, and PAC carried out the phenotypic measurements. DLS conducted the data analyses. SOPBS, MJT, EMS, RET, LTW, and YY acquired the funding for this project. RET and ZY provided breeding lines that were used in the mapping population. KYM and MJT were instrumental in providing the genotypic marker data. JBBA, PAC, KYM, YY, LTW, MJT, and EMS edited the manuscript. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
The rice germplasm used in this study were obtained from the United States Department of Agriculture (USDA) National Plant Germplasm System (NPGS) and from the inbred and hybrid rice breeding programs of the Texas A&M AgriLife Research at Beaumont. The collection and use of rice germplasm, as well as the methods carried out in this study comply with relevant institutional, national, and international guidelines and legislation.
Consent for publication
Not applicable.
Competing Interests
The authors declare that they have no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional File 1.
Supplementary Figures
Additional File 2.
Supplementary Tables
Additional File 3.
Appendix 1. Phenotyping for Maximum Node Production Rate (MNPR) and Main Culm Panicle Node Number (MCPNN)
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Sanchez, D.L., Samonte, S.O.P., Alpuerto, J.B.B. et al. Phenotypic variation and genome-wide association studies of main culm panicle node number, maximum node production rate, and degree-days to heading in rice. BMC Genomics 23, 390 (2022). https://doi.org/10.1186/s12864-022-08629-y
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12864-022-08629-y