Niegel La Borde
John Rajewski
Ismail Dweikat- Department of Agronomy and Horticulture, University of Nebraska, Lincoln, NE, United States
Sorghum (Sorghum bicolor L.) a drought tolerant staple crop for half a billion people in Africa and Asia, an important source of animal feed throughout the world and a biofuel feedstock of growing importanceorghum’s originated from tropical regions rendering the crop to be cold sensitive. Low temperature stresses such as chilling and frost greatly affect the agronomic performance of sorghum and limit its geographical distribution, posing a major problem in temperate environments when sorghum is planted early. Understanding the genetic basis of wide adaptability and of sorghum would facilitate molecular breeding programs and studies of other C4 crops. The objective of this study is to conduct quantitative trait loci analysis using genotying by sequencing for early seed germination and seedling cold tolerance in two sorghum recombinant inbred lines populations. To accomplish that, we used two populations of recombinant inbred lines (RIL) developed from crosses between cold-tolerant (CT19, ICSV700) and cold-sensitive (TX430, M81E) parents. The derived RIL populations were evaluated for single nucleotide polymorphism (SNP) using genotype-by-sequencing (GBS) in the field and under controlled environments for their response to chilling stress. Linkage maps were constructed with 464 and 875 SNPs for the CT19 X TX430 (C1) and ICSV700 X M81 E (C2) populations respectively. Using quantitative trait loci (QTL) mapping, we identified QTL conferring tolerance to chilling tolerance at the seedling stage. A total of 16 and 39 total QTL were identified in the C1 and C2 populations, respectively. Two major QTL were identified in the C1 population, and three major QTL were mapped in the C2 population. Comparisons between the two populations and with previously identified QTL show a high degree of similarity in QTL locations. Given the amount of co-localization of QTL across traits and the direction of allelic effect supports that these regions have a pleiotropic effect. These QTL regions were also identified to be highly enriched for genes encoding chilling stress and hormonal response genes. This identified QTL can be useful in developing tools for molecular breeding of sorghums with improved low-temperature germinability.
1 Introduction
Sorghum (Sorghum bicolor L.) is the fifth most important cereal crop in the world and shows high production performance in a wide range of adverse conditions (Espitia-Hernández et al., 2020; Rashwan et al., 2021). A multi-purpose crop, sorghum has been traditionally used as grain and straw, and is also a promising crop for bioenergy production from starch, sugar, or cellulose on marginal lands with limited water and other (Teferra et al., 2019; Espitia-Hernández et al., 2020; Kimani et al., 2020; Palacios et al., 2021). It has been used as a staple for millions of people in several countries, mainly in Africa and Asia. Itis a tropical C4 grass possessing many advantageous qualities to be a designated energy crop. Sorghum has a high biomass yield, high drought tolerance, high sugar content in the stalk, a short life cycle, and high adaptability to various soils and climates. Unfortunately, to be optimally utilized in the temperate US Great Plains, sorghum must develop the ability to withstand early-season chilling temperatures. The Great Plains generally have a frost-free growing season from May through September. Sorghum must reach maturity during this brief period and be agronomically competitive with other biofuel crops.
Sorghum is susceptible to chilling injury as a tropical plant, and its development is negatively affected by temperatures of 20°C and below (Peacock, 1982; Anda and Pinter, 1994). Chilling temperature stress is typical when periods of wet and cool weather overlap with planting and is exacerbated on poorly draining heavier soils. Poor seed establishment, decreased emergence rate, decreased growth after emergence, and increased susceptibility to seedling pathogens (especially Pythium spp.) have been observed in low temperatures (Soujeole Miller 1984). Poor tolerance to cold in the early season reduces yields of both biomass and grain (Saballos, 2008). The photosynthetic machinery is susceptible to chilling under field conditions in the spring, where chilling temperatures occur along with high light intensities (Jompuk et al., 2005). Root development and architecture are crucial for chilling tolerance in sorghum (Bekele et al., 2014). Although it originates in the tropics of Africa, the remarkable scope of genetic diversity among the different subspecies has conferred an extraordinarily broad adaptability and a highly versatile range of end uses (Boyles et al., 2019). Genetic variation has been observed in early and late plantings in both fields and controlled environments (McBee and Miller, 1982; Brar and Stewart, 1994). Identifying and characterizing genetics responsible for this observed variation is critical to improving sorghum performance in temperate climates. Sowing sorghum earlier in the spring would allow for better utilization of winter moisture (Patane et al., 2006). Genetic markers associated with early season seedling emergence and germination under cold stress were studied under uncontrolled field and controlled indoor conditions (Cisse and Gebisa, 2003; Yu et al., 2004; Ejeta and Knoll, 2007; Burow et al., 2010; Bekele et al., 2014; Fielder et al., 2016; Marla et al., 2019; Schaffasz et al., 2019). These studies provided valuable insights into the chromosomal regions (QTL) underlying sorghum adaption/vigor in chilling conditions, yet some of the interpretations may be clouded by QTL-by-environment interactions (Burow et al., 2010), lack correlation to field-based screenings (Yu et al., 2004; Franks et al., 2006), and limited validation in multiple genetic materials. Therefore, the present study aimed to conduct a QTL analysis of physiological traits and identify potential candidate genes underlying chilling tolerance in sorghum seedlings under field and controlled environments. The identification of significant QTL-marker associations can be used to facilitate indirect selection for cold tolerance in sorghum breeding. To accomplish this, two segregating (for chilling tolerance) recombinant inbred (RIL) populations were evaluated for pre- and post-seedling development using genotyping by sequencing method.
Genotyping-by-sequencing (GBS) provides a low-cost high-density genotyping approach that avoids ascertainment bias (Elshire et al., 2011). Genotyping-by-sequencing has been used to increase genomic resources rapidly in many crops (Poland et al., 2012; Hu et al., 2015; Lasky et al., 2015; Felderhoff et al., 2016; Wallace and Mitchell, 2017). In sorghum, large-scale GBS data were generated to characterize genomic diversity patterns and map genomic loci underlying complex trait variation (Lasky et al., 2015; Yu et al., 2016; Bouchet et al., 2017; Hu et al., 2019).
The present study aimed to: 1) investigate the physiological responses of sorghum seedlings responses to chilling stress using two strongly segregating (for chilling tolerance) recombinant inbred (RIL) populations. 2) The identification of significant QTL-marker associations of physiological traits and identify potential candidate genes underlying chilling tolerance in sorghum seedlings under field and controlled environments.
2 Materials and methods
2.1 Environmental conditions
Typical sorghum planting dates in Nebraska are mid-to-late May when soil temperatures average about 20°C. Early season sorghum planting was during mid-April at Lincoln and Havelock. The cumulative precipitation in each environment in 2014 and 2012 was higher than in the environments in 2013 (Supplementary Table S1). The environment in 2013 also experienced snowfall twice during the evaluation period. In 2013, the experimental locations had fewer days with maximum temperatures above 15°C (29 days in 2012 and 23 days in 2014, while there were 18 days in 2013).
2.2 Germination tests
For this study, two populations were used, CT19 X TX430 (C1) and ICSV700 X M81 E (C2) populations, respectively The variation in cold and optimal temperature seed germination of the parents and RILs, were examined according to procedures set forth by the Association of Official Seed Analysts (AOAC, 2019). Three replications of a completely randomized design of 50 seeds were imbibed on polystyrene Petri dishes (10 cm) with filter paper (Whatman No. 1) moistened with distilled water (4 mL). The seed was incubated in the dark for 7 days at optimal temperature (28°C) and 30 days at chilling temperature (10°C). Total germination was scored daily visually by the protrusion of radicle (1 mm) from the seed at 7 days for optimal conditions and 7, 15, 22, and 30 days after sowing for cold treatment. The cumulative germination data was used as input into a curve fitting module (GERMINATOR) developed by Joosen et al. (2010) to extract germinability parameters. The module determined the maximum percentage of germination (Gmax), the onset of germination, based on 4% germination (T4) germination speed measured as time to 50% germination (T50) uniformity of germination, based on the time interval between 14% and 86% seed germination (U8416), and the area under the curve after 336 and 720 h of imbibition under optimal and chilling respectively (AUC).
2.3 Field experiment
Typical sorghum planting dates in Nebraska are mid-to-late May when soil temperatures average about 20°C. Early season sorghum planting was during mid-April at Lincoln and Havelock. The cumulative precipitation in each environment in 2014 and 2012 was higher than in the environments in 2013 (Supplementary Table S1). The environment in 2013 also experienced snowfall twice during the evaluation period. In 2013, the experimental locations had fewer days with maximum temperatures above 15°C (29 days in 2012 and 23 days in 2014, while there were 18 days in 2013.
To explore the variation in chilling/optimal field emergence, 50 seeds from the parents and RILs were sown in five × .75 m plots under rain-fed conditions at Havelock and Lincoln, Nebraska, from 2012 to 2013 and 2014. Both locations were planted in mid-April, about a month earlier than the average sorghum planting time in the Great Plains. The C2 population was sown only in the Lincoln, NE, location in 2013 and at both locations in 2014. All but the C2 in the 2013 experiments were laid in an alpha lattice incomplete block design with 16 unfinished blocks of fifteen plots (16 × 15αlattice) per replication with two replications per environment. Due to limited seed, the C2 population was sown in an augmented incomplete block design with five blocks of 35 plots, with the parents of both people serving as controls in each block. Emergence was assessed by counting, on alternate days, the total number of plants having emerged 30 days after sowing. Emergence data was also evaluated by GERMINATOR software (Joosen et al., 2010). Maximum emergence was coded Emax. All other parameter labels remained the same.
2.4 Soil-based assay
Seed emergence and seedling development in soil were observed in an alpha lattice incomplete block design with two replications. Each block was planted in 4 × 200-cell cone-trainer flats (Stewart and Sons, Inc. Corvallis, OR) with one seed per cell. Each flat was divided into quarters with 50 cells per RIL. Plants were grown under 16-h days at optimum and chilling temperatures. The plants were watered every other day, and humidity was set at 70%.
The seedlings sown at optimum condition were harvested 14 days after sowing. The seedlings grown in chilling conditions were harvested 30 days post-sowing. Total emergence was recorded every other day until harvest for optimal and cold experiments. The cumulative emergence data were input into the curve fitting module developed by Joosen et al. (2010). Seedling biomass development was determined by recording stem length (STML), root length (RTL), stem fresh weight (STMWW), stem dry weight (STMDW), root fresh weight (ROOTWW), and root dry weight (ROOTDW). At 14 days after planting, ten seedlings per RIL were taken out of the soil and rinsed with water to remove soil from the root. Shoot and root were separated from one another, and roots were wrapped in a paper towel for 2 hours to absorb excess moisture before determining lengths and fresh root weight. Roots and shoots were then dried for 7 days at 85°C, after which dry weights were determined. Separated root and shoot measurements were pooled over the ten seedlings to prevent measurement errors on individual sample weights. The greenness of the third leaf was recorded as the mean of three measuring points on ten random seedlings using a SPAD-502 plus chlorophyll meter (Konica Minolta Sensing Inc., Osaka, Japan).
2.5 DNA extraction
Extraction of high-quality DNA from leaf tissue was carried out utilizing a modified protocol proposed by (Xin and Chen, 2012; Jiao et al., 2018). Leaf tissue was harvested from 14- day-old seedlings in the greenhouse. Approximately 1 g of fresh tissue was harvested and placed on ice in a 2 mL microcentrifuge tube. The harvested samples were transferred to a −80°C freezer overnight. The next day, the samples were transferred to a freeze dryer and allowed to lyophilize for 2 days. A 4.5 mm steel zinc-plated pellet (Daisy Outdoor Products, Rogers, Arkansas) was added to each microcentrifuge tube containing the dried leaf samples. The tubes were sealed and transferred into Tissuelyser II grinding racks. Forty-eight tubes were ground, by TissueLyser II, at 28 strokes per minute for 1 minute. Tubes were then rotated and ground for an additional minute.
To each tube, 750 µL of extraction buffer (100 mM Tris (pH8), 20 mM EDTA, 2%CTAB, 1.2M NaCl, 0.1% B-mercapoethanol) is added. Tubes were then sealed and vortexed for 20 s and then incubated at 60°C, on a dry block, for an hour. Tubes were allowed to cool on benchtop for 5 minutes then 750 µL of chloroform:isoamylalchohol (24:1 v/v) was added to each tube. Tubes were vortexed for 20 s and then centrifuged at 3000 g for 15 min. A new set of 2 mL microcentrifuge tubes were labeled and filled with 500 µL of the aqueous layer from centrifuged tubes. To these tubes, 1 mL of dilution buffer (100 mM Tris (pH8), 20 mM EDTA 0,2%CTAB) was added. Tubes were vortexed for 20 s and incubated at 60°C for 30 min. Tubes were then centrifuged at 3000 g for 13 min at 4°C. Each tube’s supernatant was carefully poured off, and 1 mL of washing buffer (30% ethanol 70% TE (10 mM Tris, 1 mM EDTA)) was mixed gently by inversion and allowed to rest for 30 min. Tubes were then centrifuged at 3000 g for 15 min at 4°C. The supernatant of each tube was gently poured off, and 100 µL of high salt TE (10 mM Tris (pH 8), 2 mM EDTA, 1M NaCl) and 5 µL of RNAse was added to each tube. Tubes were gently mixed by inversion and incubated for 30 min at 60°C. The content of each microcentrifuge tube was transferred to 96-well PCR plates. To each well of the PCR plate, 5 µL of MagAttract suspension G solution (Qiagen, Venlo, Netherlands) was added along with 120 µL of absolute ethanol. A silicone mat was used to seal the PCR plate. The plates were mixed by inversion and allowed to rest on the bench for 5 minutes. The plates were then placed on 96-well magnetic plate, and ethanol was removed by inverting the plate. The magnetic beads in each well were washed three times by alternatively adding 200 µL of washing buffer to each well, mixing the plate by inversion, placing a plate on a magnetic plate, and discarding the washing buffer. After the third wash, plates were allowed to air dry on the lab bench for 10 min. To each well, 100 µL of TE buffer has been added. The plates were sealed, placed in an incubator, and incubated for 5 minutes at 60°C. Plates were gently mixed by inversion and then placed on a magnetic plate, and the TE solution was then transferred to a new 96-well plate.
2.6 DNA quantification and quality
Genomic DNA (2 µL/sample) was quantified using Thermo Scientific Nanodrop 8,000 spectrophotometer instrument (Fisher Scientific, Waltham, MA). The quality of DNA was then examined by digesting genomic DNA (2 µL per sample) with Hind III restriction endonuclease. DNA pre- and post-digestion were visualized on 1% agarose gel in 1x TBE stained with ethidium bromide, utilizing λ standard DNA dilution series as a control. Samples were stored at −20°C until shipped to Institute for Genomic Diversity at Cornell University (Ithaca, NY).
2.7 Statistical analysis
Germination and emergence performance can be interpreted by extracting relative parameters from a time curve as described by El-Kassaby et al. (2008). Where the four-parameter Hill function is utilized to fit the germination/emergence data:
Where y is the cumulative germination/emergence at x (time in hours), y0 is the y-intercept, greater than or equal to zero, a is the maximum cumulative germination/emergence percentage, b controls the shape and steepness of the curve and c is the time for 50% germination/emergence (T50). The initial a and c values are extracted from the cumulative germination/emergence count, and b is set to 20. Utilizing the Germinator curve fitting module (Joosen et al., 2010), five emergence/germination parameters were estimated: time to onset of germination/emergence [time to 4% germination, in hrs. (T4)], rate of germination/emergence (time to 50% germination/emergence, in hours), uniformity of germination/emergence (time between 16% and 84% germination/emergence, (U8416) maximum germination/emergence (Gmax,%) and area under the curve between time zero and time until the last sampling. Greater AUC values are indicative of an earlier and greater germination/emergence.
Descriptive statistics collected for traits from both the indoor and outdoor screenings were subjected to an analysis of variance (ANOVA) statistical analysis, to determine the significance of genotypic effect to cold tolerance. The means and significant differences of means, for each phenotype, were determined, for each individual environment, using the Agricolae package (Mendiburu and Simon. 2015), in the R statistical sotware (R Core Team, 2012). The alpha lattice incomplete block design model for each screening was:
Where µ is the overall population mean, gi is the genotypic effect, rj is the replication effect, bkj is the random block effect nested within the replication effect, and εij is the residual effect with a ∼ N (0, σ2b). Where the genotypes were treated as fixed effects and replication and blocks were considered random effects. The augmented partially balanced incomplete block design model for the 2013 C2 screening was:
Where µ is the overall population mean, ρij is the block effect, τj is the genotypic effect, and εij is the residual effect with a ∼ N (0,σ2b). The ANOVA results were used to estimate the heritability: h^2 = V_l/(V_l + V_le + V_e), respectively, the estimates of the variance components of lines (genotype), line × environment interaction, and environments (Basford et al., 2004). Pearson’s correlation coefficients among the traits were calculated on a least square mean basis.
2.8 Genotyping-by-sequencing (GBS)
Genotyping of the C1 and C2populations, along with their progenitors, was performed by the Institute for Genomic Diversity according to the genotyping-by-sequencing (GBS) protocol described by Elshire et al. (2011). Four 96-well microtiter plates (two per population, containing 35 µL of the previously extracted DNA were sent to the Institute for Genomic Diversity at Cornell University. Where briefly the DNA samples were digested using the APeKI restriction enzyme, 96x multiplexed libraries were assembled and sequenced via Illumina Genome Analyzer IIx (Illumina, San Diego, CA). To extract the single nucleotide polymorphisms (SNP) genotypes, the raw reads provided by the Institute were analyzed in the TASSEL software GBS pipeline (Lu et al., 2012), 86 bp raw reads were aligned to the BTX623 sorghum reference genome (Sorghum bicolor v1.4). Loci polymorphisms were detected by comparison of consensus sequences from all samples. As RILs were utilized in the construction of the libraries, loci with heterozygotes >10% of total RIL were discarded to reduce false positive results. Only loci with less than 20% missing data were used in the mapping.
2.9 Linkage map construction
Linkage maps for the two populations were created using the SNP data from GBS using R/ASMap (Taylor and Rowley, 1971). The Kosambi function (Kosambi, 1944) was utilized to convert recombination fractions into centiMorgans (cM). To detect segregation distortion chi-square (χ2) was calculated using R/ASMap. Unlinked and highly distorted markers were excluded from the analysis. The family-wise error rate, arising from multiple χ two tests, was controlled with the Bonferroni correction. Linkage maps for both C1 and C2 populations were visualized graphically with Mapchart 2.3 (Voorrips, 2002).
2.10 QTL analysis
QTL mapping was performed on all the germination/emergence parameters, along with SPAD, root and shoot length, and fresh and dry weight of both roots and shoots. QTL mapping was based on the average phenotypic values recorded in each environment and data averaged across all environments. QTL and their positions were determined using the inclusive composite interval mapping (ICIM) function of the IciMapping 4.0 program (Li et al., 2007). A minimum logarithm of odds (LOD) of 3.0 was set to declare the significance of detected QTL. QTL with a LOD score between two and three were considered putative QTL. The “walking speed” of 1 cM and a “window size” of 10 cM were set for the ICIM. A permutation threshold of 0.01 for each trait using 1,000 permutations (Doerge and Churchill, 1996) was set. 1-LOD support intervals were determined as described by Wang et al. (2011). The percentage of variation (R2) and the additive effects are determined at their peak LOD value for each trait. The QTL maps were graphically visualized using MapChart 2.3 software (Voorrips, 2002).
3 Results
3.1 Germination assay
The germination test confirmed that the viability of both RIL populations was relatively high; the C1 population had a lower germination rate (75%) than the C2 population (88%). Seed germination of sorghum was significantly affected by temperature; the germination decreased to 31% (C1) and 45% (C2) at 10°C (Tables 1, 2). As expected, the cold-tolerant genotypes were observed to germinate earlier (T4), at a faster rate, more uniform, and with a more significant percentage than the cold-susceptible genotypes. In the RIL population, observed germination (Gmax) was more critical, earlier (lower T4), at an increased rate (lower T50), and more uniform (lower U8416) under greenhouse conditions than in the cold room for both the C1 and the C2 populations (Table 3). The variance analysis indicated that RIL’s genotypic effect was significant in both population.
TABLE 1. Early field emergence data for C1 Population at Lincoln, NE and Havelock, NE during springs 2012–2014.
TABLE 2. Early field emergence data for C2 Population at Lincoln, NE and Havelock, NE during springs 2012–2014.
TABLE 3. Indoor germination data for C1 and C2 population.
3.2 Field screening of emergence
The mean performance of field emergence rates for both the C1 and C2 populations under chilling and optimal conditions are presented in Tables 1, 2. In both population, the cold-tolerant parent (CT19 and ISCV700) emerged earlier (lower T4) and had a faster rate of germination (lower T50) across all locations and years (Table 1). In the C1 population, seedling emergence was greater (↑ Emax), faster (↓ T50), and earlier (↓ T4) at the Lincoln location than at the Havelock location. Generally, seedling emergence was most significant at the Lincoln location due to warmer temperatures. Still, in 2014 there was an herbicidal carryover from a previous season that reduced germination rates in the Havelock, NE, area. Heritability for all emergence-related traits was high in both populations and environments. The cold-tolerant parent exhibited greater tolerance to the emergence in the presence of chilling stress than the cold-susceptible parent. The emergence parameters analyzed by the hill curve are presented in Tables 1, 2.
3.3 Soil-based indoor screening
The cold-tolerant parental lines exhibited greater tolerance, as described by the emergence parameters recorded in both cold and optimal environments, than the cold-susceptible parent. Both RIL populations showed considerable variation in observed traits. In both RIL populations, emergence, as expected, was greater (↑ Emax), faster (↓ T50), (↓ T4), and more uniform (↓ U8416) under greenhouse conditions compared to the cold room (Table 3). Likewise, 14-day seedlings were observed to be greener (SPAD), taller (STML), and with longer roots under ideal conditions versus chilling temperatures (Table 4). The narrow sense heritability was moderately high for all phenotypical traits. Under chilling stress, heritability was generally lower than under optimal conditions.
TABLE 4. Indoor physiological chilling response data for C1 and C2 population.
3.4 Correlations between controlled and field screening
Correlation analysis among germination trails and field emergence screenings revealed the relationships between indoor and outdoor parameters. As anticipated, germination in optimal and chilling conditions exhibited significant correlations in both populations (Figures 1, 2). Indoor germination under chilling conditions also significantly correlated with field emergence traits—especially pooled emergence traits. While correlations were significant, the correlation coefficient was very low for both populations. Field emergence traits also correlated significantly with other field triats in both populations.
FIGURE 1. Pearson’s Correlation between indoor germination and field emergence test for CT19 x TX430 (non-significant correlations have been left blank).
FIGURE 2. Correlation between indoor germination and field emergence test for ICSV700 x M81E (non-significant correlations have been left blank.
3.5 Linkage map
For both the C1 and C2 populations, linkage maps were constructed using GBS-SNP markers (Supplementary Figures S1, S2) with less than 30% missing data. For the C1 population, 464 polymorphic GBS-SNPs were mapped to 12 linkage groups. The total map length was 2080.1 cM and covered all ten chromosomes. The lengths of the individual linkage groups ranged from 57.3 cM to 295.8 cM, with an average marker density of 4.6 cM. Eight hundred seventy-five polymorphic SNP markers were mapped to 14 linkage groups (1,515.2 cM total length) in the C2 population. Markers covered all ten chromosomes; the linkage groups ranged from 15.03 cM to 242 cM with an average marker density of 1.8 cM per marker.
3.6 QTL analysis
QTL mapping was carried out to elucidate genomic regions underlying seedling response to chilling. The results of the QTL analysis for all germination/emergence-related traits, along with SPAD, root and stem length, and their fresh and dry weights in both C1 and C2 populations, are highlighted in Tables 5, 6, Figures 3A,B, Figures 4A,B. QTL were considered stable if they appeared in more than one location for a specific trait and consistent if they appeared in more than 1 year/season. QTL identified are highlighted below.
TABLE 5. Main Effect QTL positions for C1 population.
TABLE 6. Main effect QTL positions for C2 population.
FIGURE 3. (A, B) QTL mapping for population C1 using 464 SNP markers with traits data obtained from the field (during 2012.2013.2014) and the combined controlkd indoor experiments. The traits considered are seedling chlorophyll content (wad). phnt height(M). Shoot wet weight (shootww). Shoot dry inter! (shootdw). and root dry weight (rootdw). Emergence percentage lemp). time to 4% emergence (T4), time to 50% emergence (T50), mean germination time (mgt), and uniformity of emergence (U8416).
FIGURE 4. (A, B) QTL mapping for population C2 using 875 SNP markers with traits data obtained from the field (during 2012.2013.2014) and the combined controlled indoor experiment The traits considered and seedlling Chlorophyll content (spad). plant height (ht), shoot wet weight (shootww). Shoot dry weight (shootdw). and root dry weight (rootdw). Emergence percentage (emp). time to 4% emergence (T4), time to 50% emergence (T50), and uniformily of emerengence (U8416).
3.6.1 Field experiment
In population C1, ICIM detected ten additive QTL for seedling emergence and growth under chilling stress across 3 years and two environments (Figures 3A,B). In 2012, at the Lincoln site, two significant QTL were identified on SB05. One for time to 50% emergence (T50) and the other for time to onset of emergence (T4); the former explained 9.95 the latter 22.21% phenotypic variation. At the Havelock, NE site, one significant QTL for time to cessation of emergence (T4) was identified on SB_09, explaining 6.43 percent of the observed variation (Table 5). At the Lincoln location, in 2013, five QTL were detected. Two QTL describing 10.71% and 5.38% of the variation in onset to seedling emergence were observed on SBI-05 and SBI-04. Two QTL for uniformity of emergence (U8416) were identified on SB_01, one accounting for 12.4 and the other for 7.18% of the variability in uniformity. Finally, a QTL for time to 10% emergence was identified on SB_05, describing a 6.56% variation in time. At Havelock location, two QTL were observed on SB_09. One for maximum emergence percentage (Emax), describing 6.43% of the variation in emergence percentage, and the other for the area under the curve (AUC) was responsible for 6.92% variation in the area under a curve. On SBI-06, a QTL for root length was identified, describing 6.04% of the variation in length. In population C2, 19 additive QTL were identified for seedling emergence and growth (Table 6; Figures 4A,B). At Lincoln, NE, in 2013, two QTL for uniformity of emergence were identified. One was located on SBI-03, the other on SB_09, accounting for 9.63% and 7.96% of the variation. A QTL for time to 4% seedling emergence was observed on SB_07, which explained 9.42% of the observed differences in time to 4% germination. Two QTL describing the maximum emergence percentage was located on SB_08, accounting for 8.38% of the variation. Six additive QTL were detected for the Lincoln 2014 experiment. Three QTL for the velocity of emergence were identified, two of which resided on SB_09, accounting for 43.12% of the variability in observed velocity. The final velocity QTL was observed on SBI-07, contributing 6.46% to variability in emergence velocity. Two QTL for time to 50% germination were identified, one on SB_09 and the other on SBI-2.2, accounting for 11.53% and 6.26% variation in time, respectively. One QTL for maximum emergence percentage was identified on SB9, explaining 12.41% of the variability in emergence. At the Havelock site in 2014, two additive QTL for maximum emergence time were identified on SBI-01 and SBI-04, explaining 7.11% and 7.8% of phenotypic variation. A QTL for emergence percentage was identified on SBI-09, explaining a 6.67% variation. A velocity of emergence QTL, explaining 7.52% of the variation in the emergence rate, was identified on SBI-06. A QTL for 4% emergence was located on SBI-09, accounting for a 6.99% variation in time to 4% emergence. A QTL for uniformity was identified on SBI-03, accounting for 6.81% of the observable variation.
3.6.2 Cold room experiment
In the C1 population, three QTL were identified in the linkage group SB_2.2. One for time to 50% germination, another for seedling shoot dry weight and seedling root dry weight. Each QTL accounted for 19.1, 15.45, and 5.45% of phenotypic variation, respectively.
Eight additive QTL were observed from the cold room experiment in the C2population. Two QTL for seedling root lengths were identified on SB_04 and SB_06, explaining 28.01% and 8.65% of the respective differences in root lengths. Two QTL for seedling stem dry weight were observed on SB_03, and SB_06 explained 13.52% and 9.1% of the variation in dry weight. Two QTL for seedling wet weight were detected on SB_04, and SB_1.2 accounted for 11.41% and 9.07% of the observed variation. One QTL was identified on SB_02 for SPAD (leaf greenness), accounting for 5.2% of the observed variation. Finally, a QTL on SB_04 was detected, accounting for 5.4% of stem length variability.
Five QTL were observed for germination in the cold room. One for time to 4% germination was identified on SB_03, explaining 12.2 of the variability on germination time. Two QTL for uniformity of germination (U8416) were identified on SB_10 and SB_03, accounting for 10.6% and 9.2% of the variation in uniformity observed. On SB_02, a QTL for AUC, which accounted for 8.6% of phenotypic variation, was identified. A QTL for Gmax was identified on SB_07, describing a 5.1% of the variation in germination percentage.
3.6.3 Co-localization of QTL
All the recorded germination, emergence, and vigor phenotypes for both populations were identified to co-localize on unique regions across the ten linkage groups (Figures 3A,B; Figures 4A,B). These locations are similar to those identified in previous QTL mapping experiments.
4 Discussion
Early season chilling is a significant constraint in sorghum production in the temperate US Great Plains. To establish/expand sorghum as a viable crop, good seed germination, emergence, and high seedling vigor in the early-stage chilling condition is necessary. Cold tolerance in sorghum at the germination stage can be defined as the ability of the seeds to germinate at temperatures that are usually lower than the optimum ones, which is about 10°C. Selection of an appropriate screening temperature for cold tolerance is crucial. Too high a temperature will not allow differentiation in development among lines, and below 10°C temperature will stop all growth (Hope et al., 1992). Significant variation of cold tolerance was also reported in the literatures in commercial sorghum hybrids under controlled low temperature in the laboratory (Yu et al., 2004; Maulana et al., 2017; Marla et al., 2019, Rutayisire et al., 2021). Sorghum planted under chilling conditions typically results in reduced germination, emergence, and uneven plant stands within rows. Planting sorghum earlier in the spring will allow producers to take advantage of increased rainfall and cooler temperature, allowing for pollination and fertilization to avoid the hot and dry summer.
In this study, seed germination/emergence was evaluated in two different S. bicolor populations: genotypes from the respective RIL populations from a cross of tolerant and susceptible sorghum lines were used to elucidate genomic regions associated with seed germination, emergence, and seedling vigor under cold and optimum conditions. Using field and growth chambers for phenotyping allowed a greater understanding of chilling stress’s genetic complexity. QTL regions for seedling germination, emergence, and physiological traits were detected for chilling and non-chilling stress. The vast number of QTL detected further illuminates the complex and polygenic nature of seed germination, emergence, and seedling physiology under cold conditions.
4.1 Genetic variation in chilling tolerance and related traits
During early-season planting, sorghum seedlings may be exposed to chilling stress, which often decreases field establishment. From our study, seedling soil emergence onset, duration, and speed were hindered by chilling temperatures, indicating that environmental factors play a crucial role in seedling development. Also, the filter-based germination test proved an excellent indicator of field-based emergence under chilling temperatures. Under both chilling and optimal conditions, the cold-tolerant parents showed early onset, more uniformity, and faster germination/emergence. These parents were also observed with greater root and shoot length and biomass. These results are similar to those previously reported by Podder et al. (2020) and Salas Fernandez et al. (2014). These results suggest that germination rates may be used to screen germplasm for tolerance to early-season chilling. These chilling tolerant parents were also observed with more significant root and shoot length and biomass than non-tolerant parents. This was similar to what was observed by Moghimi et al. (2019), suggesting that seedling root and shoot vigor are indicators of chilling tolerance. The high variability in traits between the genotypes in both populations demonstrates the possibility of breeding to improve cold tolerance. The high heritability estimates allow for the possible selection of high cold tolerance. The C1 and C2 populations have varying degrees of cold tolerance, which can be exploited to improve sorghum to early season chilling temperatures.
4.2 SNP markers and linkage mapping
Most previously constructed linkage maps for intra-specific sorghum have utilized SSR markers. The identification of polymorphic SSR markers is both time and labor-intensive. As a result, previous studies (Bekele et al., 2014; Burow et al., 2010; Fiedler et al., 2016, 2014; Knoll et al., 2008). have been limited to a few hundred markers at best. In recent years, next-generation sequencing (NGS) has effectively constructed high marker-density genetic maps with respect to time, labor, and costs. A GBS approach was utilized to produce a high-density, low-cost genetic map (Gore et al., 2014), allowing the identification of SNPs and genotyping at the same time. We identified a total of 1,339 novel SNPs between the two populations. In comparison to previous GBS studies (Gore et al., 2014; Zheng et al., 2011), in other plant species, we identified fewer SNPs in the present experiment. This may be attributed to the intrinsic disadvantage of GBS-SNP: a large amount of missing data due to narrow sequencing depth.
Additionally, other SNPs may have needed to be recovered in attempts to capture high-quality data by removing SNPs with >30% missing data. Sequencing errors may have also expanded the genetic distance between markers in the two linkage maps. These extended regions may hinder the discovery of QTL in such regions (Kumawat et al., 2012). Compared to previous studies, our linkage maps fell into the ranges of those previously published (Ejeta and Knoll, 2007; Burow et al., 2010). In the C1 population, chromosomes Sb02 and Sb08 and Sb01, Sb05, Sb09 (Figures 3A,B), and Sb10 in the C2 (Figures 4A,B) populations were fragmented into two linkage groups. Hiremath et al. (2012) reported linkage group fragmentation resulting from high-density genetic mapping and its lack of impact on the fidelity of QTL mapping. From allele analysis, the need for splitting linkage groups was attributed to distorted allele frequencies in the Sb02 (C1), Sb01 (C2), and Sb10 (C2). Large blocks of homogeneity between parental alleles caused other fragments.
4.3 QTL associated with cold tolerance
The four parental lines’ early season germination and seedling growth vigor varied considerably when exposed to chilling temperatures. Unlike other DNA markers, SNPs allow for the development of gene haplotyping systems (Rafalski, 2002), mainly if several defined haplotypes exist in the region of interest. Due to the absence of allelic ambiguity, co-dominant SNPs are more informative than dominant markers (Kumar et al., 2012). Unfortunately, to our knowledge, no QTL studies utilizing SNPs for early-season chilling tolerance in sorghum exist. For QTL mapping studies, populations of 100–500 individuals are suggested. Thus, the 189 genotypes for C1 and C2 were acceptable for distinguishing QTL, conferring chilling tolerance and vigor traits.
Previous mapping experiments have identified several loci contributing to early-season chilling tolerance in sorghum. Knoll et al. (2008) identified one QTL for low-temperature germination on SB_03 and two for early/late season emergence, one on SB_01 and the other on SB_02. Burrow et al. (2010) uncovered four QTL: one on SB_02 for germination in both cold and optimal conditions, one for field emergence on SB_09, and two on SB_01 for late field emergence. Bekele et al. (2014) identified three QTL for cold germination: one on SB_01, one on SB_03, and one on SB_06 (Fiedler et al. (2012) reported three QTLs for emergence percentage on SB_01, two on SB_03, and one on both SB_06 and SB_07. Finally, in an association experiment, Upadhyaya et al. (2015) detected a single QTL for cold germination on SB_07. Moghimi et al. (2019) reported markers for germination on SB_02, SB_03, SB_08, and SB_09. Emergence related markers were also detected on SB_04 and SB_07. Upadhyaya et al. (2015) detected a single QTL for cold germination on SB_07. Moghimi et al. (2019) reported markers for germination on SB_02, SB_03, SB_08, and SB_09. Emergence related markers were also detected on SB_04 and SB_07.
The complex and polygenic nature of early-season chilling tolerance was confirmed by the number of QTL detected. Sixteen and 39 QTL were identified in the C1 and C2 populations (Figures 3A,B; Figures 4A,B). Of these, four QTL in the C1 and three QTL in the C2 population have major effects (PVE 14%). Despite the different environmental settings and many genotypes, several QTL were detected across multiple years and environments. QTL detected across multiple environments are considered stable and may prove beneficial for marker-assisted plant breeding. Regions with co-locating QTL possibly housed many tightly linked genes conferring chilling tolerance and seedling vigor-related traits with pleiotropic effects. SNP’s linked to the QTL identified in this study will be used as selection tools in breeding early season vigor and cold tolerance from the cold-tolerant parents into elite sorghum breeding lines. Validation of the effects of these molecular markers across other populations and environments will be an additional, but very important, step in developing a marker-assisted breeding effort. SNP’s which can be further validated will be introgressed into elite breeding lines, which can then be evaluated in hybrid combinations.
4.4 QTL for physiological traits related to cold tolerance
Several seedling vigor traits were evaluated under chilling and optimal conditions to identify genomic regions associated with seedling vigor. Six QTL regions in the C1 population and 18 QTL for the C2 population were identified, confirming the polygenic nature of seedling vigor traits. Of the QTL identified, one major QTL was identified in the C1 population, and three major QTL for the C2 population. In both populations, the majority of the seedling vigor QTL were centered on SB_01, SB_02, SB_03, SB_04, and SB_05, which coincided with previously mentioned studies.
On chromosome on Sb-01, two major stem-related QTL were identified in the C2 population, one for dry stem weight under optimal conditions where the M81e alleles conferred extracellular mass to seedlings’ stem. The other QTL identified stem length under optimal conditions. Investigating the intervals spanning these QTL, numerous putative genes were identified (Supplementary Table S2). Genes in the underlying stem dry weight QTL interval included a response to cytokinin and desiccation (Sb01g001660). Genes identified in the underlying stem length QTL interval included a response to growth hormones (Sb01g031060; Sb01g030930; Sb01g031050), a response to cold (Sb01g031870), etc.
SB_02 contained one major QTL for shoot dry weight, in optimal conditions, in the C1 population (Figures 3A,B) and one major QTL for root dry weight in the C2 population (Figures 4A,B). The QTL for shoot dry weight conferred a reduction in stem mass. This QTL co-localized with QTL for germination and emergence in chilling temperatures, as described earlier. The underlying intervals contained the same putative genes. Under optimal conditions, the QTL for root dry weight led to a reduction in root mass in seedlings possessing alleles. The interval spanned by the QTL contained putative genes for response to abiotic stress (Sb02g026490; Sb02g026500) and membrane functions (Sb02g026380; Sb02g026460).
Under optimal conditions, a major QTL for root length was discovered on SB_03 in the C2 population (Figures 4A,B). The allele led to a reduction in the length of seedling roots. Root length QTL was previously mapped to SB_03 (Bekele et al., 2014). However, these were different intervals we reported earlier. In our QTL interval, putative genes are highlighted in the Supplementary Table S2. Two major QTL for dry root weight and root length were detected on SB_04. The alleles at both QTL reduced the seedlings’ stem weight and root length. Moghimi et al. (2019) also reported a two QTLs for root biomass in the same region. One QTL for stem length under optimal conditions was identified on SB_05. The QTL for root dry weight has not previously been detected; upon an investigation of the underlying interval, some putative genes are highlighted in Supplementary Table S2. The QTL for root length co-localized with a QTL previously identified by Bekele et al. (2014).
5 Conclusion
The study was conducted based on the hypothesis that sorghum plants can tolerate low temperatures during germination and must contain genes that can be mapped that contribute to this tolerance compared to plants that will not germinate under those conditions. Early season chilling stress significantly limits sorghum productivity in the Great Plains. Chilling stress is a complex trait to study. While a great deal of the physiological and biological responses to chilling is well characterized in sorghum, the genetics underlying chilling tolerance are still being gleaned. The present study has added to understanding the genetics and physiology governing early season chilling tolerance through elucidation of QTL for germination, emergence, and seedling vigor and assigning functions to said QTL. The identified areas may serve as starting points for the enhancement of molecular breeding programs by allowing for an improvement in selection efficiency through a selection of QTL that allows for sorghum to improvegermination, emergence percentage, and respective rates while decreasing the time needed for both in the presence of chilling temperatures. As more NGS-based mapping and association studies are conducted, we anticipate the fine mapping of major QTL hubs and the discovery of new areas conferring chilling tolerance to sorghum. The results provide important new insights for adaptive crop breeding in the face of climate change and the expansion of sorghum production to different regions. This will facilitate sorghum from being a “crop of the future” to transforming into a real-life important agricultural alternative. Further studies are needed to confirm these QTLs by using near isogenic lines (NILs) and backcrossing methods. With the identification of increasing numbers of favorable alleles at QTLs for cold tolerance by QTL analysis, the pyramiding approach and MAS strategy could become a promising approach for improving cold tolerance.
Data availability statement
The authors acknowledge that the data presented in this study must be deposited and made publicly available in an acceptable repository, prior to publication. Frontiers cannot accept a manuscript that does not adhere to our open data policies.
Author contributions
ID conceived and designed the experiments. NL and JR performed the experiments. NL analyzed the data. NL and ID wrote the article. All authors contributed to the article and approved the submitted version.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2023.1129460/full#supplementary-material
References
Anda, A., and Pinter, L. (1994). Sorghum germination and development as influenced by soil temperature and water content. Agron. J. 86, 621–624. doi:10.2134/agronj1994.00021962008600040008x
AOAC (2019). Association of Official analytical chemists. Official methods of analysis of AOAC international. 21st ed. Washington: AOAC International. MD (2019).
Bekele, W., Fiedler, K., Shiringani, A., Schnaubelt, D., Windpassinger, S., Uptmoor, R., et al. (2014). Unravelling the genetic complexity of sorghum seedling development under low-temperature conditions. Cell & Environ. 37, 707–723. doi:10.1111/pce.12189
Bouchet, S., Olatoye, M. O., Marla, S. R., Perumal, R., Tesso, T., Yu, J., et al. (2017). Increased power to dissect adaptive traits in global sorghum diversity using a nested association mapping population. Genetics 206 (2), 573–585. doi:10.1534/genetics.116.198499
Boyles, R. E., Brenton, Z. W., and Kresovich, S. (2019). Genetic and genomic resources of sorghum to connect genotype with phenotype in contrasting environments. Plant J. 97 (1), 19–39. doi:10.1111/tpj.14113
Brar, G. S., and Stewart, B. A. (1994). Germination under controlled temperature and field emergence of 13 sorghum cultivars. Crop Sci. 34, 1336–1340. doi:10.2135/cropsci1994.0011183x003400050036x
Burow, G., Burke, J. J., Xin, Z., and Franks, C. (2010). Genetic dissection of early-season cold tolerance in sorghum (Sorghum bicolor (L.) Moench). Mol. Breed. 28, 391–402. doi:10.1007/s11032-010-9491-4
Cisse, N. D., and Gebisa, E. (2003). Genetic variation and relationships among seedling vigor traits in sorghum. Crop Sci. 43, 824–825. doi:10.2135/cropsci2003.0824
Collins, N. C., Tardieu, F., and Tuberosa, R. (2008).Quantitative trait loci and crop performance under abiotic stress: Where do we stand? Plant Physiol. 147, 469–486. doi:10.1104/pp.108.118117
Doerge, R. W., and Churchill, G. A. (1996). Permutation tests for multiple loci affecting a quantitative character. Genetics 142, 285–294. doi:10.1093/genetics/142.1.285
Ejeta, G., and Knoll, J. E. (2007). Marker-assisted selection in sorghum in Genomics-assisted crop improvement. Editors R. K. Varshney, and R. Tuberosa (Dordrecht: Springer). doi:10.1007/978-1-4020-6297-1
El-Kassaby, Y. A., Moss, I., Kolotelo, D., and Stoehr, M. (2008). Seed germination: Mathematical representation and parameters extraction. For. Sci. 54 (2), 220–227.
Elshire, R. J., Jeffrey, C., Glaubitz, J., Sun, Q., Poland, J., Kawamoto, K., et al. (2011). A robust, simple genotyping-by-sequencing (GBS) approach for high diversity species. PloSOne 6 (5), e19379. doi:10.1371/journal.pone.0019379
Espitia-Hernandez, P., Gonz ́alez, M. L. C., Ascacio-Vald ́es, J. A., D ́avila-Medina, D., Flores-Naveda, A., and Silva, T. (2020). Sorghum (Sorghum bicolor L.) as a potential source of bioactive substances and their biological properties. Crit. Rev. Food Sci. Nutr., 62 1–12. doi:10.1080/10408398.2020.1852389
Felderhoff, T. J., McIntyre, L. M., Saballos, A., and Vermerris, W. (2016). Using genotyping by sequencing to map two novel anthracnose resistance loci in Sorghum bicolor. G3 (Bethesda) 6, 1935–1946. doi:10.1534/g3.116.030510
Fiedler, K., Bekele, W. A., Matschegewski, C., Snowdon, R., Wieckhorst, S., Zacharias, A., et al. (2016). Cold tolerance during juvenile development in sorghum: A comparative analysis by genomewide association and linkage mapping. Plant Breed. 135 (5), 598–606. doi:10.1111/pbr.12394
Franks, C. D., Burow, G., and Burke, J. (2006). A comparison of U.S. And Chinese sorghum germplasm for early season cold tolerance. Crop Sci. 46, 1371–1376. doi:10.2135/cropsci2005.08-0279
Gore, M. A., Fang, D. D., Poland, J. A., Zhang, J., Percy, R. G., Cantrell, R. G., et al. (2014). Linkage map construction and quantitative trait locus analysis of agronomic and fiber quality traits in cotton. Plant Genome 7, 1–62. doi:10.3835/plantgenome2013.07.0023
Hiremath, P. J., Kumar, A., Penmetsa, R. V., Farmer, A., Schlueter, J. A., Chamarthi, S. K., et al. (2012). Large-scale development of cost-effective SNP marker assays for diversity assessment and genetic mapping in chickpea and comparative mapping in legumes. Plant Biotechnol. J. 10, 716–732. doi:10.1111/j.1467-7652.2012.00710.x
Hope, H. J., Maamari, R., Séguin, S., Hamilton, R. I., Dwyer, L. M., and White, R. P. (1992). Low temperature emergence potential of short season corn hybrids grown under controlled environment and plot conditions. Can. J. Plant Sci. 72 (1), 83–91. doi:10.4141/cjps92-009
Hu, Z., Mbacké, B., Perumal, R., Guèye, M. C., Sy, O., Bouchet, S., et al. (2015). Population genomics of pearl millet (Pennisetum glaucum (L.) R. Br.): Comparative analysis of global accessions and Senegalese landraces. BMC Genomics 16, 1048. doi:10.1186/s12864-015-2255-0
Hu, Z., Olatoye, M. O., Marla, S., and Morris, G. P. (2019). An integrated genotyping-by-sequencing polymorphism map for over 10,000 sorghum genotypes. Plant Genome 12 (1), 180044. doi:10.3835/plantgenome2018.06.0044
Hu, Z., Olatoye, M. O., Marla, S., and Morris, G. P. (2019). An integrated genotyping-by-sequencing polymorphism map for over 10,000 sorghum genotypes. Plant Genome 12, 180044.doi:10.3835/plantgenome2018.06.0044
Jiao, Y., Burow, G., Gladman, N., Acosta-Martinez, V., Chen, J., Burke, J., et al. (2018). Efficient identification of causal mutations through sequencing of bulked F2 from two allelic bloomless mutants of sorghum bicolor. Front. Plant Sci. 8, 2267. doi:10.3389/fpls.2017.02267
Jompuk, C., Fracheboud, Y., Stamp, P., and Leipner, J. (2005). Mapping of quantitative trait loci associated with chilling tolerance in maize (Zea mays L.) seedlings grown under field conditions. J. Exp. Bot. 56 (414), 1153–1163. doi:10.1093/jxb/eri108
Joosen, R. V., Kodde, J., Willems, L. A., Ligterink, W., van der Plas, L. H., and Hilhorst, H. W. (2010). Germinator: A software package for high-throughput scoring and curve fitting of arabidopsis seed germination. Plant J. 62, 148–159. doi:10.1111/j.1365-313X.2009.04116.x
Kimani, W., Zhang, L.-M., Wu, X.-Y., Hao, H.-Q., and Jing, H.-C. (2020). Genome-wide association study reveals that different pathways contribute to grain quality variation in sorghum (Sorghum bicolor). BMC Genomics 21, 112. doi:10.1186/s12864-020-6538-8
Kumar, S., Banks, T. W., and Cloutier, (2012). SNP discovery through next-generation sequencing and its applications. Int. J. Plant Genomics 2012, 831460. doi:10.1155/2012/831460
Kumawat, G., Raje, R. S., Bhutani, S., Pal, J. K., Mithra, A. S. V. C. R., Gaikwad, K., et al. (2012). Molecular mapping of QTLs for plant type and earliness traits in pigeonpea (Cajanus cajanL. Millsp.). BMC Genet. 13, 84. doi:10.1186/1471-2156-13-84
Lander, E. S., and Botstein, D. (1989). Mapping mendelian factors underlying quantitative traits using RFLP linkage maps. Genetics 121, 185–199. doi:10.1093/genetics/121.1.185
Lasky, J. R., Upadhyaya, H. D., Ramu, P., Deshpande, S., Hash, C. T., Hyma, K., et al. (2015). Genome–environment associations in sorghum landraces predict adaptive traits. Sci. Adv. 1 (6), e1400218. doi:10.1126/sciadv.1400218
Li, H., Ye, G., and Wang, J. (2007). A modified algorithm for the improvement of composite interval mapping. A Modified Algorithm Improv. Compos. Interval Mapping.” Genet. 175, 361–374. doi:10.1534/genetics.106.066811
Lu, F., Glaubitz, J., Harriman, J., Casstevens, T., and Elshire, R. (2012). TASSEL 3.0 universal network enabled analysis kit (UNEAK) pipeline documentation. White Pap. 2012, 1–12.
Marla, S. R., Burow, G., Chopra, R., Hayes, C., Olatoye, M. O., Felderhoff, T., et al. (2019). Genetic architecture of chilling tolerance in sorghum dissected with a nested association mapping population. G3: Genes, Genomes. Genetics 9, 4045–4057. doi:10.1534/g3.119.400353
Maulana, F., Weerasooriya, D., and Tesso, T. (2017). Sorghum landrace collections from cooler Re.gions of the world exhibit magnificent genetic differentiation and early season cold tolerance. Front. Plant Sci. 8, 756. doi:10.3389/fpls.2017.00756
McBee, G. G., and Miller, F. M. (1982). Carbohydrates in sorghum culms as influenced by cultivars, spacing, and maturity over a diurnal Period1. Crop Sci. 22, 381–385. doi:10.2135/cropsci1982.0011183x002200020042x
Mendiburu, F., and Simon, R. (2015). Agricolae -ten years of an open source statistical tool for experiments in breeding, agriculture and biology. Peer J. Pre Prints 3, e1404v1. doi:10.7287/peerj.preprints.1404v1
Moghimi, N., Desai, J. S., Bheemanahalli, R., Impa, S. M., Vennapusa, A. R., Sebela, D., et al. (2019). New candidate loci and marker genes on chromosome 7 for improved chilling tolerance in sorghum. J. Exp. Bot. 70 (12), 3357–3371. doi:10.1093/jxb/erz143
Palacios, C. E., Nagai, A., Torres, P., Rodrigues, J. A., and Salatino, A. (2021). Contents of tannins of cultivars of sorghum cultivated in Brazil, as determined by four quantification methods. Food Chem. 327, 127970. doi:10.1016/j.foodchem.2020.127970
Patane, C., Cavallaro, V., Avola, G., and D'Agosta, G. (2006). Seed respiration of sorghum [Sorghum bicolor (L.) Moench] during germination as affected by temperature and osmoconditioning. Seed Sci. Res. 16 (4), 251–260. doi:10.1017/SSR2006259
Peacock, J. M. (1982). “Response and tolerance of sorghum to temperature stress,” in Sorghum in the Eighties: Proceedings of the International Symposium on Sorghum, Patancheru. A.P. India, 2–7 November 1981.
Podder, S., Samarappuli, D., Anderson, J., and Berti, M. T. (2020). Phenotyping a diverse collection of forage sorghum genotypes for chilling tolerance. Agronomy, 10(8), 1074. doi:10.3390/agronomy10081074
Poland, J. A., Brown, P. J., Sorrells, M. E., and Jannink, J.-L. (2012). Development of high-density genetic maps for barley and wheat using a novel two-enzyme genotyping-by-sequencing approach. PLoS One 7 (2), e32253. doi:10.1371/journal.pone.0032253
Rafalski, A. (2002). Applications of single nucleotide polymorphisms in crop genetics. Curr. Opin. Plant Biol. 5, 94–100. doi:10.1016/s1369-5266(02)00240-6
Rashwan, A. K., Yones, H. A., Karim, N., Taha, E. M., and Chen, W. (2021). Potential processing technologies for developing sorghum-based food products: An update and comprehensive review. Trends Food Sci. Technol., 110, 168. doi:10.1016/j.tifs.2021.01.087
Saballos, A. (2008). Development and utilization of sorghum as a bioenergy crop in Genetic improvement of bioenergy crops. Editor W. Vermerris (New York, NY: Springer). doi:10.1007/978-0-387-70805-8
Salas Fernandez, M. G., Schoenbaum, G. R., and Goggi, A. S. (2014). Novel germplasm and screening methods for early cold tolerance in sorghum. Crop Sci. 54 (6), 2631–2638. doi:10.2135/cropsci2014.01.0025
Schaffasz, A., Windpassinger, S., Snowdon, R., and Wittkop, B. (2019). Reproductive cold stress tolerance in sorghum F1 hybrids is a heterotic trait. Agronomy 9, 508. doi:10.3390/agronomy9090508
Silvia, R., Casazza, A. P., Enrico, C. M., Havaux, M., Jennings, R. C., and Soave, C. (2006). Suppression of both ELIP1 and ELIP2 in arabidopsis does not affect tolerance to photoinhibition and photooxidative stress. Plant Physiol. 141, 1264–1273. doi:10.1104/pp.106.083055
Taylor, A. O., and Rowley, J. A. (1971). Plants under climatic stress: I. Low temperature, high light effects on photosynthesis. Plant Physiol. 47, 713–718. doi:10.1104/pp.47.5.713
Teferra, T. F., Amoako, D. B., Rooney, W. L., and Awika, J. M. (2019). Qualitative assessment of ‘highly digestible’ protein mutation in hard endosperm sorghum and its functional properties. Food Chem., 271, 561. doi:10.1016/j.foodchem.2018.08.014
Voorrips, R. E. (2002). MapChart: Software for the graphical presentation of linkage maps and QTLs. J. Hered. 93, 77–78. doi:10.1093/jhered/93.1.77
Wallace, J. G., and Mitchell, S. E. (2017). Genotyping-by-sequencing. Curr. Protoc. Plant Biol. 2, 64–77. doi:10.1002/cppb.20042
Wang, Z., Wang, J., Bao, Y., Wu, Y., and Zhang, H. (2011). Quantitative trait loci controlling rice seed germination under salt stress. Springer Neth. 178, 297–307. doi:10.1007/s10681-010-0287-8
Xin, Z., and Chen, J. (2012). A high throughput DNA extraction method with high yield and quality. Plant Methods 8, 26. doi:10.1186/1746-4811-8-26
Yu, J., Tuinstra, M. R., Claassen, M. M., Gordon, W. B., and Witt, M. D. (2004). Analysis of cold tolerance in sorghum under controlled environment conditions. Field Crops Res. 85, 21–30. doi:10.1016/s0378-4290(03)00125-4
Keywords: SNP’s, sorghum breeding, QTL, GBS, chilling tolerance
Citation: La Borde N, Rajewski J and Dweikat I (2023) Novel QTL for chilling tolerance at germination and early seedling stages in sorghum. Front. Genet. 14:1129460. doi: 10.3389/fgene.2023.1129460
Received: 22 December 2022; Accepted: 23 February 2023;
Published: 15 March 2023.
Edited by:
Edel Pérez-López, Laval University, CanadaReviewed by:
Javaid Akhter Bhat, Nanjing Agricultural University, ChinaKazutoshi Okuno, University of Tsukuba, Japan
Copyright © 2023 La Borde, Rajewski and Dweikat. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ismail Dweikat, idweikat2@unl.edu