1. Introduction
Upland cotton (Gossypium hirsutum L.) is one of the most important cash crops in the United States and worldwide. Although Upland cotton is relatively heat- and drought-resilient compared to some other major crops, enhancements are needed to maintain production, further improve sustainability, and confront global climate change. Most of the cotton-producing regions frequently face shortages of irrigation water due to variable rainfall and access to irrigation or irrigation water, commonly coupled with higher evapotranspiration demands [1,2]. The global surface temperature has increased by approximately 0.6–0.85 °C since 1880 and is projected to continue, increasing 0.2 °C per decade [3]. In addition to decreasing and uneven precipitation, warmer days with a decreased diurnal temperature incidence are projected to occur more often in the future [4,5]. A combination of too high temperatures and drought events could reduce major food crops’ yields up to 60% worldwide ([6,7]. Furthermore, a comparison of reported average crop yields with historical record yields in the USA-grown crops is one third to one seventh of record, due to unfavorable environmental conditions [8].
Due to the indeterminate growth habits of cotton, the stored soil moisture and local microclimatic conditions strongly influence new branching nodes, flowers, bolls, and the partitioning of resources between source and sink [1,2]. In the southern United States, reproductive and boll formation stages of cotton coincide with high temperatures and low precipitation. A rise in temperature above the optimum (>30 °C) in July-August combined with high humidity and low rainfall results in significant impairment to the plant growth of cotton, changes leaf biochemical and biophysical properties and reduces boll formation and boll size due to higher evapotranspiration demand under field conditions [2,9,10,11]. While small variations in temperature may not cause significant damage to seed weight/boll weight, such a narrow spike in temperatures during anther dehiscence causes a substantial reduction in fertility/seed number ([12,13,14]. Short periods of heat events (>30 °C) are reported to significantly limit cotton productivity by disrupting a series of reproductive events [12,15,16]. Frequent heat spikes during flowering impair reproduction by inhibiting male [17] and female [18] gametophyte development, pollen germination [19], and pollen tube length growth [20]. Similar to other crops, a decline in boll set and the abortion of squares was observed upon exposure to daytime temperatures exceeding 30 °C for about 13 h both in Upland and Pima cotton [2,12,21,22,23].
Genetic variability in physiological, reproductive traits, lint yield, and fiber quality in cotton is evident in cotton, indicating the potential to increase stress tolerance [24,25]. To date, selections have been focused on yield and fiber quality traits to develop progenitors of G. hirsutum from its wild relatives [26] for non-stress conditions. Wild resources have been utilized in the breeding program to improve the stress resilience of Upland cotton lines by stacking progenitor alleles. G. barbadense (superior fiber quality) and G. tomentosum (heat-resistant species of the genus, [27] have been extensively used to explore beneficial alleles for useful traits in the breeding programs [26,28,29,30,31,32]. For example, G. barbadense (Pima cotton) is being used as a donor for fiber quality [26,33]. However, even with the continued genetic gain (2015–2019), lint yields of Upland cotton average 600 kg ha−1 more than Pima cotton [34]. This difference in cotton yield could be due to differences in stress tolerance [12,13,35].
High temperature and drought trigger morph-physiological, biochemical, and molecular alterations leading to weaker plant growth and productivity, and ultimately yield in many crops, including cotton [35,36]. Among numerous traits, plants’ ability to maintain robust growth and development under stress signifies the greater carboxylation efficiency, a cooler canopy, and greater enzyme activities at mitochondria and chloroplasts (high chlorophyll and membrane stability). The impacts of heat and drought stress on morphological (leaf area and biomass), physiological (gas exchange and pigments), biophysical (membrane fluidity or stability and canopy temperature), and reproductive (pollen germination ability) traits have been widely studied across crops [37,38,39]. For instance, gas exchange (photosynthesis and stomatal conductance), canopy temperature, and cellular or membrane stability traits are integrated attributes that reflect the morphological, anatomical, physiological, and biochemical features at the root and leaf [36,37,38,39]. Genetic variability in these traits has been effectively utilized to understand mechanisms, discriminate between tolerant and sensitive genotypes, and for crop improvement under different stress conditions [37,38,39].
The effects of these abiotic stresses are often inter-related. Previous analyses of CS lines have associated many novel traits with substituted chromosomes and chromosome segments from the wild species and unadapted germplasm [29,30,31,32]. However, the existing collection of CS lines have not been tested for physiological and reproductive traits associated with enhanced agronomic performance under field conditions. The long-term goal of screening cotton breeding lines and new germplasm for abiotic stress tolerance, such as heat and drought, is to identify a new genetic variation that minimizes yield losses in common cotton-producing environments. Genetic gains would complement farm management efforts that partially alleviate stress effects and stabilize cotton yields across increasingly variable environmental conditions. Given that genomes of CS lines are mostly of Upland origin and true-breeding, CS lines identified as stress-resilient can be directly incorporated into breeding programs aiming to improve sustainability and guard against the climate crisis, which has significantly projected more frequent episodes of higher temperatures and drought intensities. This research’s immediate objective was to identify cotton substitution lines with improved traits against the two most important challenges (high temperature and drought) in sustainable cotton production.
We utilized CS lines developed through hypoaneuploidy-based backcrossing, using “TM-1” as a recurrent parent, to quantify selected physiological and reproductive characteristics associated with line yield. The present study investigated CS lines’ performance in carrying a wild genomic composition of G. barbadense and G. tomentosum in TM-1 background for leaf gas exchange, pigments, biophysical and reproductive traits under field conditions. We classified the CS lines based on the degree of tolerance to heat and drought using vegetative, physiological, and reproductive parameters. Finally, we associate the assessed phenotype with exotic substitution chromosome/s to identify potentially favorable CS lines for future cotton breeding. 2. Results and Discussion
The present study characterized CS lines under field conditions to identify potential CS lines with differences in leaf gas exchange, pigments, biophysical and reproductive traits compared with TM-1. Each CS-B or CS-T line was bred to replace a pair of chromosomes or chromosome arms of G. hirsutum inbred TM-1, with the respective pair from G. barbadense doubled-haploid 3–79 or G. tomentosum, respectively [28,31,40]. These backcrossed CS lines are nearly isogenic to the recurrent parent TM-1 for 25 chromosome pairs and similar to each other for 24 chromosome pairs [26,28]. Therefore, this uniform genetic background provides us with opportunities for comparing these lines and an accentuated possibility of detecting a genetic effect(s) of the substituted chromosomes and chromosome arms on specific traits inherited from the alien species.
Planting was done to coincide with the reproductive stage of the CS lines to occur specifically between late-July and mid-August. The range in air temperature (minimum and maximum) and rainfall during the cropping season is visualized in Figure 1. We observed the maximum air temperature combined with the lowest rainfall at the experiment site during this period. The average maximum air temperature varied from 34.4 °C (SD ± 1.9) in July to 35.3 °C (SD ± 2.8) (Figure 1), which is 4.4 to 5.3 °C above the optimum temperature of 30 °C for successful growth, development and reproduction in cotton [2,12,21,22,23]. On the other hand, on nearly all days in July (30 days) and August (29 days), the maximum daily temperature was above 30 °C. The same hotter window coincided with the year’s lowest rainfall; on average, about 58.5% lower rainfall was witnessed compared to the early vegetative stage (157.4 mm) (see Figure 1). The corresponding total rainfall was 68.7 mm in July and 63.7 mm in August. In the southern United States, a rise in temperature above the optimum (>30 °C) in July-August combined with low rainfall is known to induce damage to leaf gas exchange, biochemical and biophysical properties, and reduces boll formation and boll size due to higher evapotranspiration demand under field conditions [2,9,10,11]. In this study, CS lines differed significantly for many measuresd physiological and reproductive characteristics, and various response indexes (Table 1).
The relative performance of CS lines compared with TM-1, the recurrent parent of the CS lines, has been given in Table S1. The Dunnett’s test was performed using trait means to test the difference between each CS line and TM-1 and the overall result from the stress response index analysis for different traits used to identify CS lines that exhibit increased tolerance to high temperature, drought, or a combination of both abiotic stresses (Table 2).
2.1. Leaf Gas Exchange and Pigment Parameters
The eleven CS-B and CS-T lines varied significantly for various leaf gas exchange and pigment parameters (Table 1). Five CS lines had photosynthesis (Pn) significantly greater than TM-1, 13.7% on average (Table 1). These included three CS-T lines, namely CS-T07, CS-T18, and CS-T15sh, and two CS-B lines, CS-B04 and CS-B15sh (Table 2). Significant variation also occurred for stomatal conductance (Gs, mmol m−2 s−1). Two CS-T lines, CS-T07 and CS-T15sh, and one CS-B line, CS-B15sh, exhibited significantly greater Gs than TM-1, whereas one was considerably lower, i.e., CS-T04 (Table 2). Somewhat in parallel, lines CS-T15sh and CS-T04 also exhibited the maximum and minimum transpiration rates (T, mmol H2O m−2 s−1); the T for CS-T15sh (24%) was significantly higher than TM-1 (Table 1 and Table S1), and the T for CS-T04 was significantly lower than the average of all CS lines (11.48). In this study, Pn had significant correlations with Gs (r = 0.88, p < 0.001) and T (r = 0.86, p < 0.001), indicating the strong influence of stomatal conductance to CO2 on final leaf photosynthesis [42,43]. Thus, two CS lines (CS-T15sh and CS-T07) with increased Gs and T could be used as potential donors for improving heat tolerance (an indication of evaporative cooling) or cotton yield potential in a warmer environment [44,45].
The water use efficiency of several crops is linked with Pn, T, or both [46]. Values of instantaneous water-use efficiency (iWUE) varied by a factor of 1.2 among the CS lines, with CS-T04 and CS-T15sh being the highest and lowest, respectively (Table 1). The iWUE of CS-T04 was significantly higher (16.6%) than the iWUE of TM-1, while CS-T15sh was nearly significantly lower than TM-1 (Table 1 and Table S1). The iWUE was negatively correlated with Pn, Gs, and T, with intensities increasing in that order (Figure 2). The strongest negative correlation for iWUE, however, was with Ci/Ca (r = −0.79, p < 0.001) (Figure 2). The Ci/Ca ratio was fairly strongly correlated with T (r = 0.57) (Figure 2). These results seem to suggest that the final carbon assimilation in CS lines was most strongly influenced by Gs. A higher iWUE in CS-T04 indicates a lesser Gs. Though CS-T07 and CS-T18 had a combination of high Pn and Gs, the iWUEs of these are not significantly different compared to TM-1 (Table 2). That variation among the CS lines for gas exchange traits was statistically significant, which suggests that genetic differences among the CS lines affected these traits.
Genetic difference in total chlorophyll content (TCHL) at the reproductive stage was significant at p < 0.05, and no CS line had greater TCHL than either TM-1 or the CS line average (45.6) for this trait (Table 1). Nonetheless, the fact that CS-T07 maintained greater TCHL (13%) and Pn (19%) at flowering compared to CS-B07 could reflect a differential genetic ability of CS-T07 to maintain chlorophyll content and higher photosynthetic rate (Table 1). Similarly, soybean with higher chlorophyll content showed a higher Pn at different growth stages [47]. During the reproductive phase in the current work, carotenoids (Caro, µg cm−2) showed a narrow range from about 9.3 to 11.0 and no significant variance due to differences among the lines (Table 1). However, there was a positive (r = 0.87) correlation between Caro and THCL (Figure 2). Because a plant with relatively high Pn value is associated with delayed senescence (or the ‘stay trait’), the CS-T07 line with greater leaf pigment values can reflect more efficient Pn as well as light-harvesting capacity, similar to studies in other species [48,49]. The combination of a relatively high Pn and TCHL (greener canopy at flowering) observed for CS-T07 and CS-T18 (Table 1) could indicate that these lines have a delayed senescence stay-green trait for improving agriculturally elite types of cotton.
2.2. Biophysical and Temperature-Tolerant Indices
Measurements of canopy temperature depression (CTD, °C) reported the difference between ambient air and canopy temperatures. They indicated the ability of a line to lower canopy temperature through transpiration-based cooling. Variation of CTD among CS lines was significant and ranged from −2.52 to 5.05, with three CS lines being significantly cooler than TM-1, namely CS-T15sh, CS-B04, and CS-T08sh (Table 1 and Table 2). CS-T04 had the hottest canopy, but differences from TM-1 were nonsignificant. The significance of internal leaf cooling in association with Gs driven by T and yield while under stress has been highlighted in G. barbadense Pima cotton [44]. Values for CTD were correlated negatively with transpiration (T) (r = −0.65, p < 0.05), and positively with iWUE (r = 0.62, p < 0.05) (Figure 2). Studies in Arizona on the impact of high temperature on cotton yields indicated that tolerant cultivars transpired more water, which resulted in cooler canopies that led to flowers and young bolls being less exposed to high temperatures [45,50]. In those studies, the high-temperature tolerant cultivars produced a greater yield than temperature-sensitive cultivars. Breeding cultivars tolerant to high- and low temperatures would be an alternative approach to minimize future cotton production’s negative impact. Low-temperature tolerance would enable farmers to plant early. Thus, flowering could occur before the mid-summer high temperatures that often limit the fruit set [51].
The resilience of cell membranes to elevated temperatures and greater chlorophyll stability values are correlated with heat and drought tolerance in cotton [41]. A higher cell membrane and chlorophyll stability value signifies a CS line plant’s ability to withstand stress due to greater cellular level tolerance, which would lead to increased photosynthesis, more dry matter production, and improved yield performance [52]. The variance of cell membrane temperature stability (CMTS) due to line effects was highly significant. While the average of CS lines (39.5%) was similar to that of TM-1 (39.8%), the values of the three CS lines, CS-B04, CS-T01, and CS-T18, were significantly higher than TM-1, and the importance of one line, CS-T08sh, was substantially lower than TM-1 (Table 1 and Table 2). Plants with high CMTS were used as a proxy to select high yielding cotton under high-temperature conditions [53]. While CMTS and CTD were slightly negatively correlated overall (Figure 2), CS-T04 exhibited the highest CMTS and the highest CTD (Table 1), indicating the existence of complex genetic interactions. The CSI is related to cotton plant tolerance to environmental stresses [41] and differed significantly among CS lines (Table 2).
Interestingly, none of the tested CS lines achieved greater CSI than TM-1, and the average CSI of CS lines (86.4%) was well below the observed TM-1 value (94%). While CSI values need not be consistent across studies, the present study’s values were within the range reported for 38 cotton cultivars [41], whereas the TM-1 value was spuriously high. Among all CS lines, CS-B04 had relatively high CMTS and CSI, the second lowest CTD, and among the highest Pn and T (Table 1).
Specific leaf area (SLA, cm2 g−1) is a critical plant functional trait, and variation in SLA is a useful indicator for leaf mass allocation under stress conditions [54,55]. The SLA values varied significantly and ranged from 117 cm2 g−1 in CS-B01 to 138 cm2 g−1 in CS-T01, with an average of 126 cm2 g−1 (Table 1). As compared to TM-1, CS-T01 recorded 7% greater SLA, whereas CB-B01 recorded 9.3% lower SLA (Table 1 and Table S1). SLA at the reproductive stage was most closely correlated with CMTS (r = 0.65) (Figure 2) and indicated that CS lines with high SLA might also maintain high CMTS under non-stressed conditions. In previous studies, low SLA has been used as an indicator of drought-tolerant lines in field crops [56,57,58] and in cotton [45] because thicker leaves (low SLA with high chlorophyll density) usually maintain higher photosynthetic capacity and WUE [59]. Thus, CS-T07, which had low SLA, high TCHL, and the highest Pn capacity, would seem to be prospectively useful as a source for drought tolerance attributes (Table S1).
2.3. Pollen Viability and Germination Parameters
Compromised pollen viability and function have been implicated as yield-influencing casualties of high temperatures in cotton [19] and other crops [60]. In this study, line-to-line differences significantly affected all of the pollen performance traits (Table 1). Pollen viability (PV), based on TTC staining, ranged from 33% in CS-T04 to 52% in CS-B07, with an average of 45%, roughly 80% of the observed TM-1 percentage (55%). Likewise, pollen germination (PG30) of CS lines at 30 °C was approximately lower by 20% than TM-1 (52%). However, CS-T15sh recorded high pollen germination (49%) at 30 °C, similar to TM-1 (Table 2). For all lines, the frequencies of PG30 were just a few percent lower than the PV rates, and patterns of variation among lines for PV and PG30 were nearly identical, resulting in a near-perfect correlation (Table 1 and Figure 2). Almost all of the pollens judged to be viable were capable of germination at 30 °C (non-stress condition). Five CS-lines had negative/lower pollen viability and germination percentages than the adapted Upland TM-1 cotton line (Table 2). The pollen viability and germination rates observed in the present study are similar to or slightly higher than observed previously in cultivars with similar methods [19,41] but lower than those observed using different methods [61]. The observed lower PV by 25% and PG by 27% (at 30 °C) compared to TM-1 in the present study might indicate stronger reproductive incompatibility, which could be due to interspecific and conspecific competitions between TM-1 and wild progenitors. These observations promote follow up studies to unravel the mystery of reproductive incompatibility between wild accessions and adapted cotton lines.
It is interesting to note that chromosomes 1, 4, and 18 from CS-B and CS-T lines in our previous studies generally showed a negative additive effect on lint yield in cultivars [29,30]. However, this report showed that chromosomes 1 and 4 from CS-B lines showed more significant additive effects in fiber quality traits than fiber strength in homologs. A previous study conducted in controlled-temperature sunlit plant growth chambers at the flowering stage showed that when Upland cotton plants were placed at 29 °C (long-term July temperature at Starkville, Mississippi, the USA plus 2 °C), all the flowers that were set after 5 days of stress gradually declined [51]. On the contrary, all the bolls abscised within five days after flowering in the sunlit plant growth chambers kept at the long-term July mean temperature plus 7 °C. About 200 to 250 bolls m−2 were retained in all temperatures; however, the number of bolls retained was strongly influenced by temperature. These observations suggest a need for follow up investigation to unravel the mystery of reproductive incompatibility at different temperatures between wild accessions and adopted Upland cotton lines and its effects on yield.
Production of sterile pollen under high temperatures has been associated with the early degeneration of the anther tapetal layer [62]. High day temperatures that exceeded 33 °C during heat-sensitive micro and macro sporogenesis resulted in reduced pollen viability and pollen number in many field crops [60]. Along the same line, decreased pollen viability could significantly affect pollen germination on the stigma [63,64], depending on the growth temperature [65]. Therefore, pollen germination inhibition under heat stress [64,66,67,68,69] negatively impacts the seed/fruit set in crops [69,70,71]. Upland cotton lines have been bred for heat tolerance by selecting progenies developed from surviving pollen grains pre-exposed to 35 °C for 15 min [72]. This result suggests that higher pollen viability germination could be used as a proxy to screen cotton cultivars for high-temperature tolerance. Burke et al. [61] recorded sharp declines in the percentages of cotton pollen germination at 37 (71%) to 40 °C (23%) and little germination at 43 °C. In this study, pollen grains showed a significant (p ≤ 0.001) decrease (23.1%) in the pollen germination percentage at 38 °C compared to 30 °C (Table 1). The PG38 for the CS lines varied by only about 11% in absolute terms, from nearly 27 or 28% for CS-T04 and CS-T01 to nearly 37 or 38% for CS-B04 and CS-T08sh. Thus, no line had substantially different PG38 than TM-1 (33.1%) (Table 1). While the overall line-to-line variance PG38 was significant statistically, most differences between individual lines and TM-1 were small relative to the LSD value of 5.9% and statistically significant for CS-T04 only (p < 0.05). The largest decrease in the pollen germination response (PGR) to elevated temperature was recorded in TM-1 (37%). However, this merely followed the overall trend that the variation in PGR values among lines largely reflected the variation in PG30 (r = −0.84), rather than the variation in PG38 (r = −0.32) (Figure 2). Additional data will be required to determine whether the desirable variations in PG38 and PGR, e.g., CS-T04 and CS-T01, point to potentially useful genes for heat tolerance.
2.4. Promising CS Lines for Future Studies
Cotton, like some other plant species, produces dynamic responses at the morpho-physiological and cellular levels to tolerate unfavorable temperature and drought conditions. Cotton is cultivated under a wide range of very diverse environments and extremely variable temperatures throughout the world. A recent report showed that increased high temperature and drought stress episodes could reduce crop yields by 50% [73]. Using an index-based classification system for drought-, heat- and combined-resistance, similar to Singh et al. [41], each CS line was categorized as stress tolerance towards heat, drought, and combined heat and drought (Table 3). The heat stress response index (HSRI) varied significantly among the 11 CS lines and ranged from 10.9 (heat-sensitive, CS-T08sh) to 12.4 (heat tolerant, CS-B15sh). Based on the drought stress response index (DSRI), CS-T08sh recorded the lowest value of 9.6, while three CS lines, CS-T07, CS-B15sh, and CS-T18, exhibited a greater value by 14.8% (p < 0.05) compared with the drought-sensitive CS-T08sh (Table 3). Based on the cumulative heat and drought stress response index (CHDSRI) value and an increment of SD, four CS lines (CS-T01, CS-T07, CS-B15sh, and CS-T18) were identified as heat- and drought-sensitive (Table 3). Among the tested lines, CS-B07 and CS-T08sh exhibited low drought, low heat, and combined low heat and drought stress tolerance (Table 3). Two lines from G. tomentosum, i.e., CS-T07 and CS-T18, showed the highest indices for drought and heat tolerance individually and under a combination of drought and heat tolerance. As CHDSRI had strong positive associations with HSRI (R2 = 0.97) and DSRI (R2 = 0.98), CS-T07, CS-T18, and CS-B15sh were identified as the most drought- and heat-tolerant CS lines (Table 3). These CS lines had significantly greater Pn (14%), Gs (29%), T (13%), and PGR (23.6%), as compared to TM-1 under field conditions (Table S1). Future studies should aim towards validating a relationship between total plant variation (cumulative indices) and heat/drought using a diversity panel. The identified heat- and drought-tolerant lines represent a promising means to develop novel climate-ready cotton genotypes after validating under gradient temperature and drought conditions.
2.5. Comparative Analysis of CS Lines Physiological and Pollen Phenotype and Modern Advanced Upland Cotton Cultivars
Our comparative analysis among these 11 CS lines and previous observations across 38 cultivars [41] shows that CS lines have increased leaf Pn, Gs, and T by 17%, 23%, and 12% compared to modern advanced Upland cotton cultivars, respectively (Table 1). CS lines have altered the magnitude of these traits. On the other hand, the larger TCHL (25%) and Caro (21%) differences were also observed between the CS lines and modern advanced cotton cultivars at the reproductive stage (Table 1). It was shown that the selection of Pima cotton with higher yield potential and heat resistance was linked with smaller leaf areas [45]. Thus, it could be expected that the decreased leaf characteristic (low SLA) in CS lines could account for genetic differences in gas exchange parameters, yield potential, and heat resistance. The CS lines (CS-T07, CS-B15sh, and CS-T15sh) maintained a considerably greater Gs than the mean for 38 cotton cultivars (0.76 vs. 0.51 mol m−2 s−1; see [41]), whereas TM-1 exhibited intermediate behavior (Table 1). The difference in CTD between the advanced CS lines and cotton cultivars were observed under field conditions. Averaged across CS lines and cultivars, the CS lines exhibited ~1.3 °C cooler canopy than the advanced Upland cotton cultivars (Table 1). The data show that CMTS in the CS lines is more significant (66%) than the mean of 38 cotton cultivars under field conditions. These findings indicate that the increases in CMTS and lower CTD could be causally related to the increases in Gs in the CS lines. The combination of such trait relationships in cotton showed higher yield potential under hotter climates [45]. Lower CTD (negative values) and higher Gs appear to have provided CS lines with a heat resilient mechanism during pollination, which favored maintaining higher pollen germination under heat stress conditions (38 °C, Table 1). The PGR was comparatively high in CS lines by 19 and 12%, compared to TM-1 and other cotton cultivars, respectively (Table 1). Phenotypic comparison between CS lines and domesticated Upland cotton cultivars made by Singh et al. [41] showed an increased phenotypic magnitude (Table 1).
The available data suggest that these CS lines will be useful for diversifying and enhancing heat and drought tolerance traits in elite Upland cotton. The findings indicate that additional studies are needed to define the physiological potential of CS lines under various combinations of abiotic stress conditions involving high and low temperatures, drought, and nutrient deficit conditions, and various salt conditions. CS line responses to various CO2 levels should also be characterized, given that the CO2 levels are projected to keep increasing and can affect many of these traits [74,75,76].
2.6. Compositional Analysis
Each CS line used in this analysis was bred previously by a pedigreed modified backcross-inbreeding scheme with the maintenance of hypoaneuploid-enforced hemizygosity during the backcross generations and the recovery of a disomic homozygote after self-pollination of the BC5F1 hemizygote. The procedure takes advantage of low-to-moderate transmission of monosomy and segmental monosomy (e.g., segmental deletion heterozygote) via an isogenic hypoaneuploid TM-1-like Upland cotton seed parent and the euploid-only transmission of the donor chromosome via pollen from the backcross hemizygotes. While hemizygous, a chromosome or segment is denied an opportunity for recombination, and in this case, assures that pollen-transmitted donor chromosomes are genetically intact. Meanwhile, recombination and segregation affect the rest of the genome for each generation, i.e., the heterozygous non-target chromosomes and segments. Expectations are that the end-products will be substituted for the targeted chromosome or chromosome segment, plus any donor segments inadvertently retained during backcrossing and rendered homozygous upon selfing and euploid recovery. The presence of donor segments at any location, of course, presents opportunities for the discovery of donor germplasm effects, subsequently refined localization and use in breeding programs. As part of a more extensive single nucleotide polymorphism (SNP) study for all CS lines that will be reported in greater detail separately, we (I) assessed each CS line used in this study for the substitution status of the respective target chromosome or segment with both SSR and SNP markers [77,78], and (II) conducted a genome-wide survey for donor germplasm based on SNP content of the CottonSNP63K array, according to the methods described in previous studies [28,31,77].
Molecular marker analysis confirmed targeted substitutions for 9 of the 11 CS lines, while two of them lacked the intended substitution but contained non-target donor segments (Table 4). As expected, the genomes of all of the CS lines were substituted for donor germplasm at non-target sites varying in size and number (Table 4). These were characterized into two classes, equal to and larger versus smaller than 5 Mbp, and identified as to their chromosomal association (Table 4). Additional details about each donor segment in CS lines will be communicated as part of an upcoming comprehensive analysis of all CS lines. To put relative sizes into perspective, the average Gossypium AD-genome (n = 26) chromosome is roughly 100 Mbp. The overall genome was flowing cytometrically, estimated to be about 2.4 Gbp [79], about 3500 cM recombinationally [77], and the most recent public sequence assembly is 2.3 Gbp [80]. Thus, each chromosome would, on average, be about 4% of the genome, whereas each 5 Mbp segment would about 0.2% of the genome and about 5% of one chromosome (Table 4).
The presence of moderate numbers of non-targeted alien segments in CS lines is more advantageous than disadvantageous. Still, the detection of their presence underscores the importance of being guarded about the location(s) of genes causing experimentally detected effects. The non-targeted segmental substitutions increase the coverage of the donor genome and thus increase the chances of seeing donor genetic effects of interest. Moreover, the availability of relatively inexpensive high-density genome-wide SNP genotyping as well as very affordable non-destructive simplex genotyping of seed and seedlings makes it relatively simple to track individual segments of hybrids, i.e., after chromosome assortment or homology-based recombination. 2.7. Outlook
High temperatures and drought often occur coincidently, and the stresses they cause are often interrelated. Individually or combined, they cause a wide range of phenotype changes that adversely affect cotton development and productivity. How plants deal with heat and drought stresses can be categorized as escape, avoidance, tolerance, and recovery strategies [81]. Drought or heat tolerance provides plants with the ability to endure severe dehydration or heat effects through specific physiological and biochemical activities, such as osmotic adjustment via osmoprotectants [82]. Drought or heat escape is plants’ ability to regulate their growth period or lifecycle, such as the cotton variety with a short life cycle, to get out of the seasonal drought stress or heat period [83]. Drought recovery of plants is a genotype’s capacity to resume growth and yield after exposure to severe drought or heat stresses. It is critical to understand the genetic mechanisms associated with abiotic stresses by detecting the phenotype of the complex cotton morpho-physiological process related to such abiotic stress. We used the indices for heat, drought, and heat and drought combined to identify the CS lines tolerant against such stresses by using a method that was reported in detail in our previous research [35,41,84,85]. It is interesting to note that CS-T07 and CS-T18 showed high tolerance against heat, drought, and heat and drought based on the indices (Table 3). While it remains most likely that the G. tomentosum chromosome 7 and 18 substitutions in CS-T07 and CS-T18, respectively, carry the genes responsible for the observed drought and heat tolerance, genes in the other non-targeted donor segments in these lines could be accountable.
It is possible that G. tomentosum will be a valuable source of multiple abiotic resilience traits, in that it grows in dry, rocky, or clay coastal areas in the main islands of Hawai at elevations below 390 feet. [86]. Cotton originated in a wild perennial habitat and currently grows in environments that suffer from periodic drought and heat stress [87]. Nonetheless, G. tomentosum was reported to be the most heat-resistant genus [27]. Modern cotton elite lines resulting in the subsequent intensive selection are from different latitudes, geographies, and climates to produce high lint yield and improved fiber quality for hand and mechanical harvesting and textile processing. This domestication of elite Upland cotton cultivars has unintentionally narrowed genetic diversity among the cotton lines for several essential traits such as heat tolerance and water use efficiency [88,89]. The primary challenges to the genetic improvement of Upland cotton include (i) low genetic diversity among the growing cultivars, (ii) inadequate information about genomic regions or genes controlling traits such as drought and heat tolerance, and (iii) suitable genetic material for economically important characteristics to improve Upland cotton germplasm. Drought and heat are two of the most critical factors limiting productivity in many cotton-growing areas because cotton is susceptible to heat and water stress [90]. Previous research showed the negative impacts of these stresses on fiber quality and yield [91,92,93]. Identification of heat- and drought-tolerant cotton lines requires correct phenotyping of the associated complex morpho-physiological traits. The innovative method used in this research provided an opportunity to discover several heat- and drought-tolerant cotton CS lines that will be useful in sustainable cotton production against the severe threat of climate crisis from global warming. We advocate future research to estimate genotype-specific anatomical or cellular level tolerance at temporal and spatial scale to better relate total plant variation indices or health with heat and drought stress tolerance.
3. Materials and Methods
The overall methods of detecting the stress-resilient cotton lines were based on the stress response indices method used in previous studies [35,41,84,85,94].
3.1. Planting Materials
Seeds of eleven CS lines bred for chromosome or chromosome segment substitutions from G. barbadense and G. tomentosum in the genetic background of G. hirsutum (TM-1) were used in the experiment. Each line, e.g., CS-B01, was bred previously by crossing a developed TM-1 hypoaneuploid, e.g., TM-1 monosomic for chromosome-1, as seed parent with the donor species, selecting cytogenetically for the corresponding hypoaneuploid interspecific BC0F1 hybrid; additional introgression was by analogous backcrossing and recovery of hypoaneuploid BCnF1 hybrid, i.e., maintaining the targeted chromosome in a hemizygous state. After the fifth backcross (BC5F1) or higher, the hypoaneuploid BCnF1 was inbred and progeny cytogenetically selected for euploidy to establish a homozygous disomic substitution line, e.g., BC5S1, for seed increase. Five CS lines were bred from crosses with the G. barbadense donor doubled haploid line 3–79, namely CS-B01, CS-B15sh, CS-B04, CS-B07, and CS-B18, whereas six CS lines were bred from a cross with G. tomentosum donor, namely CS-T04, CS-T15sh, CS-T07, CS-T18, CS-T01, and CS-T08sh (Table 1). These CS lines are, as expected, nearly isogenic with each other and TM-1 and differ only for the targeted chromosome or chromosome segment in each line, Upland cotton (G. hirsutum). TM-1 is an inbred line bred from the commercial cultivar Deltapine 14 (Delta and Pine Land Co.; Scott, MS, USA) and maintained over 40 generations by self-pollination, as described by Kohel et al. [95]. Comparative analysis of the true-breeding homozygous CS lines with each other and/or TM-1 provides opportunities to associate genetic and trait variations with the lines and substitutions therein. Seeds were sown in an experimental field at the Rodney Foil Plant Science Research Center, Mississippi State University, Mississippi State (lat. 33°28′ N, long. 88°47′ W) during 2010 following a randomized complete block design (RCBD) with four replications. Plots were fertilized with 22.7 kg N pre-plant and 34.0 kg N side-dress 35 days after planting (DAP). The soil type at the Mississippi State University location is a Leeper silt clay loam (fine, smectitic, nonacid, thermic Vertic Epiaquept). Plants were grown under standard production practices for commercial production.
3.2. Pollen Viability under Field Conditions and Germination Performance under High-Temperature Stress
Twenty to thirty fresh flowers at anthesis were randomly collected from 10–15 plants for each CS line between 08:00 and 10:00 h from the first fruiting position during flowering, 55 to 60 days after planting. Pollen grains were collected by gently tapping the flowers over Petri dishes that were used to test pollen viability and germination. A pollen viability assay was performed using 2% (aq) 2,3,5-triphenyl tetrazolium chloride (TTC) stain in deionized water, as described by Aslam et al. [96]. The TTC solution stains the viable pollen with a reddish-purple color due to insoluble red formazan formation. A drop of TTC solution was added to the pollen dispersed on a microscope slide. Thirty minutes after staining, the total number of pollen grains and the number of stained pollen grains were counted in two microscopic fields of 2.4 mm2 containing more than 100 pollen grains from each microscopic field in all replications.
Taylor’s [97] pollen growth medium, consisting of 2 g agar, 30 g C12H22O11, 5.3 mg KNO3, 51.7 mg MnSO4, 10.3 mg H3BO3, 10.3 mg MgSO4.7H2O made up to 100 mL with deionized water, was used to study pollen germination [19]. The prepared germination medium was poured into three replicate Petri dishes for each CS line and cooled for 15 min for agar solidification. The collected pollen grains were distributed uniformly on the solidified germination medium using a tiny, clean, bristle paintbrush. The Petri dishes were then covered and incubated in an incubator (Precision Instruments, New York, NY, USA) at 30 and 38 °C for 24 h. Each Petri dish per TM-1, CS-line, and temperature treatment combination was considered as a replicate. Total pollen grains and the number of pollen germinated pollen grains were counted using a Nikon SMZ 800 microscope (Nikon Alphaphot YS microscope; Nikon Instrument, Kanagawa, Japan). A pollen grain was considered germinated when the pollen tube length was at least equal to or greater than the grain diameter [94]. Germination percentage was determined as the ratio of the number of germinated pollen grains total number of pollens per field of view. Further, pollen germination response (PGR) was calculated to evaluate the reproductive responses of TM-1 and each CS line using the following equation.
Pollen germination response (PGR) = (PG38 °C/PG30 °C)
where PG38 °C and PG30 °C are the pollen germination percentages of each line at temperatures of 38 and 30 °C, respectively.
3.3. Gas Exchange Parameters Leaf photosynthesis (Pn), transpiration (T), and stomatal conductance (Gs) were measured between 10.00 and 14.00 h on a cloud-free day, using the third or fourth leaves from the top using a portable photosynthesis system (Model LI-6200, LI-COR, Lincoln, NE, USA) 55 to 60 DAP. The measurements were taken at 1500 µmoL photon m−1 s−1 photosynthetically active radiation (PAR), the chamber temperature was set at 30 °C, and carbon dioxide concentration was maintained at 400 ppm. The internal CO2 concentration (Ci) normalized for external CO2 concentration (Ca), and the ratio of leaf internal to ambient carbon dioxide concentration (Ci/Ca) were also measured. Instantaneous water-use efficiency (iWUE) was calculated as the ratio between Pn/T, which indicates the amount of net CO2 assimilation per unit of water loss. 3.4. Leaf Pigment Content and Chlorophyll Stability Index (CSI)
Two sets of leaf samples were collected from five fully expanded leaves during the flowering stage. Five leaf discs, each with an area of 2.0 cm2, from each sample, were collected randomly and placed in vials containing four milliliters of dimethyl sulphoxide (DMSO) for chlorophyll (Chl) extraction. The sample vials were incubated at room temperature in the dark for 24 h, the extract decanted into a cuvette, and the absorbance was measured using a Bio-Rad ultraviolet/VIS spectrophotometer (Bio-Rad Laboratories, Hercules, CA, USA) at 470, 648, and 663 nm to calculate concentrations of Chl a, Chl b, and carotenoid content [98]. The total leaf Chl (TCHL) was estimated by summing Chl a and Chl b values.
The chlorophyll stability index (CSI) was assessed according to the procedure described by Murphy (1962). Another set of leaf discs, each with an area of 2.0 cm2, was collected similarly from each cotton line and incubated at 56 °C in a temperature-controlled water bath for 1 h. The set of tubes was brought to 25 °C, and the Chl content was measured from the heat-treated samples as described previously. The CSI was estimated as the percentage of TCHL in heated leaf (56 °C) relative to that in fresh leaf using the following equation.
Chlorophyll stability index (CSI, %) = (TCHL in heat-treated)/(TCHL in control) × 100
3.5. Leaf Cell Membrane Thermostability (CMTS)
The leaf CMTS in TM-1 and CS-lines was assessed following the procedure described by Singh et al. [41] in each of the three replicates. Two sets (control and treatment) of 10 leaf discs, each disc with a 1.3 cm2 area, were cut from five fully expanded third or fourth leaves from the top of the stem axis from each line. The leaf discs were thoroughly rinsed three times with deionized water to remove electrolytes, both adhering to the leaf surface and leaching from the leaf discs’ cut surfaces. All the test tubes with leaf discs were filled with 10 mL of deionized water and capped with aluminum foil to prevent water evaporation. After incubation for 20 min at 55 °C, the sets of test tubes were brought to 25 °C, and initial measurement of the conductance of the control (CEC1) and the treatment (TEC1) was measured by an electrical conductivity meter (Corning Checkmate II; Corning Inc., Corning, NY, USA) at room temperature. Tubes were then autoclaved at 0.1 MPa for 12 min to kill tissues completely, releasing all the electrolytes. Tubes were then cooled to 25 °C, and final conductance was measured (CEC2 and TEC2). Finally, the ratio of the initial to the last electrolyte leakage caused by the elevated temperature was estimated, giving a measure of the extent of damage to cellular membranes using the following equation.
CMTS(%)=[1 − (TEC1/TEC2)]/Leaf CMT (%) [1 − (CEC1/CEC2)] × 100
3.6. Canopy Temperature Depression
Canopy temperature depression (CTD) was measured during the reproductive stage between 12:00 and 13:00 h (cloudless, bright days) using a handheld infrared thermometer (Model OS533E-OMEGASCOPE; OMEGA Engineering, Inc., Stamford, CT, USA). Canopy temperature depression was estimated using the equation below
Canopy temperature depression (CTD,°C)=Tc−Ta
where Ta and Tc refer to the ambient air and canopy temperatures of the target cotton line, respectively.
3.7. Specific Leaf Area Leaf areas of three individual leaves (used for gas exchange) were measured using a leaf area meter (LI-3100: Li-COR, Lincoln, NB, USA). The dry weight and measured leaf area of all lines and cultivar TM1 were used to calculate the specific leaf area (SLA), expressed as cm2 g−1 d. Wt. 3.8. SSR Analysis
The overall method of SSR marker genotyping was as described by Saha et al. [78]. Two to three young leaves from each CS line and TM-1 were used for DNA extraction using DNAeasy cotton genomic DNA extraction protocol (Qiagen kit), and extracted genomic DNAs were checked with Nanodrop and diluted to appropriate concentrations before using for PCR. The PCR-amplifications were performed in a 10 µL reaction mix containing 1 µL 109 PCR buffer, 0.2 µL dNTPs (10 mM each), 1.2 µL 25 mM MgCl2, 0.3 µL 5 lM labeled primers (FAM, HEX, VIC, PET), 0.1 µL AmpliTaq Gold DNA polymerase (Applied Biosystems, Foster City, CA, USA), and 5 ng µL−1 genomic DNA. With labeled primers, polymorphisms at amplified SSR loci among chromosome substitution lines were detected in denaturing capillary electrophoresis according to manufacturer’s protocol on the Genetic Analyzer 3130xl (Applied Biosystems, Foster City, CA, USA). The size of the amplified products was detected using Gene Mapper 4.0 (Applied Biosystems, Foster City, CA, USA).
3.9. SNP Analysis
CS lines were SNP genotyped by (i) extracting DNA using a commercially available kit (Qiagen DNeasy Plant, Macherey-Nagel NucleoSpin Plant II, or OPS SYNERGY Plant DNA), (ii) testing the DNA quality and concentration, (iii) adjusting concentration, and (iv) analyzing with the DNA with the CottonSNP63K Array at the Texas Institute for Genome Sciences and Society (TIGSS) Lab. Genotypes were called with the published cluster file; then, genotypic data were purged of SNPs yielding less-than high uniformity and had more than 10 percent missing data [77]. Linkage groups were identified using a previously published G. hirsutum × G. barbadense F2 linkage map [77], and SNPs were ordered similarly. Processed genotype data were then scanned to find chromosomes or segments (multiple contiguous marker bins) homozygous or heterozygous for the donor or non-TM1 genotypes. The results were used first to assess if the targeted chromosome or chromosome segment was indeed substituted, and secondarily, to scan the genome for the presence of any inadvertent chromosome or segment substitutions which involved non-targeted chromosomes. The latter donor segment was categorized and tabulated according to physical size: more versus less than 5 Mbp, based on the most recent Upland cotton (TM-1) genome sequence assembly [80].
3.10. Stress Response Indices and Statistical Analysis
The indices for heat, drought, and heat and drought combined were developed following the approach that was reported in previous studies [35,41,84,85]. In brief, the relative response index of the CS line was calculated as the value for a CS line divided by the phenotypic value of the Upland recurrent parent, cultivar TM-1. Then, the cumulative indices, namely cumulative heat stress response indices (HSRI), cumulative drought response indices (DRSI), and cumulative heat and drought response indices (CHDSRI) calculated for each line, were used for the classification of high and low stress-tolerant lines. The DSRI was calculated as the summation of the individual relative index of five gas exchange (Pn + Gs + T + iWUE + Ci/Ca), two pigments (TCHL + Caro), two cell and chlorophyll stability indices (CMTS + CSI), and two biophysical (CTD + SLA) traits.
DSRI = (Pn + Gs + T + iWUE + Ci/Ca + TCHL + Caro + CMTS + CSI + CTD + SLA)
Likewise, the HSRI was calculated as the sum of individual relative gas exchange (Pn + Gs + T + iWUE+ Ci/Ca), pigments (TCHL + Caro), cell and chlorophyll stability indices (CMTS + CSI), biophysical (CTD), and reproductive response (PV + PG38/PG30) traits.
HSRI = (Pn + Gs + T + iWUE + Ci/Ca + TCHL + Caro + CMTS + CSI + CTD + PV + PGR)
The CHDSRI of each CS line was calculated as the sum of all 13- traits relative indices used for the phenotyping of heat tolerance and drought tolerance.
CHDSRI = (Pn + Gs + T + iWUE + Ci/Ca + TCHL + Caro + CMTS + CSI + CTD + SLA + PV + PGR)
The 11 CS-lines were classified as low, moderate, and tolerant based on minimum HSRI or DSRI or CHDSI values and an increment of one standard deviation, respectively, as Singh et al. [41] described.
Data were analyzed statistically using ANOVA in RStudio 3.6.1 (https://rstudio.com/) to estimate the significance of variation among lines. Pearson’s correlation coefficient was carried out between all traits using the library (“corrplot”). Dunnett’s test compared the difference between the CS line and TM-1 using the library (DescTools) in RStudio, and a difference detected at p < 0.05 by ANOVA was considered significant.
4. Conclusions In summary, replicated analyses of multiple isogenic CS lines at flowering enabled detection of differential responses to various physiological, biophysical, and reproductive traits relevant to performance stability, sustainability, and global warming. The CS lines that imparted significant desirable effects on traits will be useful in breeding programs for targeted introgression of abiotic stress resistance genes of drought and high-temperature conditions in Upland cotton improvement programs, for which additional genetic diversity is needed. For a few traits, the phenotypic extreme(s) among these 11 CS lines exceeded previous observations across 38 Upland cotton cultivars, which suggests that these CS lines could be compelling for diversifying and enhancing these traits in elite Upland cotton. In addition, all or most of the CS lines imparted statistically significant beneficial effects on traits of interest; these lines are likely sources of desirable genes. Upon introgression into elite breeding germplasm lines by marker-assisted breeding, these genes will expectedly extend the phenotypic range of elite cultivars beyond what is currently known. However, there is an urgent need to investigate how cotton cultivars and CS lines will perform under future projected elevated CO2 conditions combined with other abiotic stress conditions.
Enlarge this image.Figure 1. Daily weather parameters (temperature and rainfall) during crop growing season. Data were obtained from the PRISM weather data (http://www.prism.oregonstate.edu) for the experiment site.
Enlarge this image.Figure 2. Cotton chromosome substitution (CS) line trait correlations of the recorded parameters. Dark color values represent strong correlations, and light color values represent weaker correlations. Values close to zero indicate no correlation between the two variables; whereas, values close to ± indicate strong correlation (positive-dark blue and negative-dark red) between two traits. Leaf gas exchange (Pn-photosynthesis, Gs-stomatal conductance, T-transpiration, Ci/Ca-leaf internal to ambient carbon dioxide concentration ratio and iWUE-instantaneous water use efficiency), cell and chlorophyll stability indices (CMTS-cell membrane stability index, and CSI chlorophyll stability index), biophysical (CTD-canopy temperature depression, and SLA-specific leaf area), and reproductive (PV-pollen viability, PG30-pollen germination at 30 °C and PGR-pollen germination response). Correlation values ≥ ±0.55 are significant at p < 0.05.
| CS Line | Pn | Gs | T | iWUE | Ci/Ca | TCHL | Caro | CTD | CMTS | CSI | SLA | PV | PG30 | PG38 | PGR |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| CS-B01 | 33.3 | 0.51 | 10.4 | 3.20 | 0.71 | 43.5 | 9.8 | 3.40 | 38.5 | 92.0 | 117.0 | 45.0 | 42.2 | 32.0 | 75.9 |
| CS-B04 | 38.1 | 0.70 | 12.5 | 3.10 | 0.73 | 44.3 | 9.9 | −1.70 | 48.1 | 88.0 | 131.0 | 50.0 | 47.2 | 37.5 | 79.1 |
| CS-B07 | 30.8 | 0.50 | 10.1 | 3.10 | 0.70 | 42.9 | 9.3 | 2.90 | 35.4 | 80.0 | 124.0 | 52.0 | 49.3 | 32.7 | 66.3 |
| CS-B15sh | 37.2 | 0.83 | 12.3 | 3.00 | 0.71 | 47.5 | 10.3 | 3.50 | 38.9 | 82.0 | 125.0 | 50.0 | 47.1 | 33.9 | 71.9 |
| CS-B18 | 34.0 | 0.56 | 10.8 | 3.20 | 0.71 | 49.8 | 11.0 | 1.40 | 35.7 | 87.0 | 127.0 | 50.0 | 46.5 | 37.2 | 79.8 |
| CS-T01 | 33.4 | 0.63 | 11.7 | 2.90 | 0.74 | 43.3 | 9.7 | 0.90 | 45.5 | 83.0 | 138.0 | 36.0 | 33.3 | 28.0 | 84.0 |
| CS-T04 | 32.8 | 0.47 | 9.5 | 3.50 | 0.68 | 45.9 | 10.1 | 5.00 | 39.5 | 85.0 | 124.0 | 33.0 | 29.6 | 27.0 | 91.3 |
| CS-T07 | 38.1 | 0.86 | 12.4 | 3.10 | 0.71 | 49.4 | 10.4 | 3.70 | 35.3 | 86.0 | 121.0 | 43.0 | 40.2 | 33.2 | 74.6 |
| CS-T08sh | 35.4 | 0.64 | 11.5 | 3.10 | 0.72 | 42.7 | 9.9 | −0.70 | 29.0 | 94.0 | 123.0 | 49.0 | 45.8 | 37.8 | 82.4 |
| CS-T15sh | 37.4 | 0.84 | 13.3 | 2.80 | 0.74 | 44.6 | 10.2 | −2.50 | 42.0 | 90.0 | 123.0 | 50.0 | 46.9 | 33.1 | 70.7 |
| CS-T18 | 37.8 | 0.76 | 11.8 | 3.20 | 0.67 | 47.4 | 10.2 | 1.90 | 46.8 | 83.0 | 132.0 | 38.0 | 34.6 | 30.4 | 88.4 |
| TM-1 | 33.1 | 0.63 | 10.7 | 3.10 | 0.70 | 45.2 | 9.9 | 4.30 | 39.8 | 94.0 | 129.0 | 55.0 | 52.4 | 33.1 | 63.3 |
| LSD (p = 0.05) | 2.9 | 0.1 | 1.7 | 0.28 | ns | 4.7 | ns | 3.3 | 3.8 | 4.4 | 6.9 | 6.5 | 6.52 | 5.79 | 9.7 |
| Significance | *** | *** | ** | ** | * | *** | *** | *** | *** | *** | *** | ** | *** | ||
| CS lines mean + | 35.3 | 0.66 | 11.5 | 3.11 | 0.71 | 45.6 | 10.1 | 1.62 | 39.5 | 86.4 | 125.9 | 45.1 | 42.1 | 33.0 | 78.6 |
| Advanced Upland cotton cultivars # | |||||||||||||||
| Minimum | 17 | 0.17 | 4.5 | 1.60 | 0.43 | 30.4 | 6.7 | 1.2 | 15.6 | 72 | 155.9 | 43.8 | 25.6 | 15.4 | 42.8 |
| Maximum | 38.6 | 0.85 | 14.4 | 4.40 | 0.84 | 38.1 | 10.2 | 4.3 | 40.4 | 92.9 | 199.6 | 66.3 | 49.9 | 41.8 | 94.4 |
| Mean | 29.3 | 0.51 | 10.1 | 3.20 | 0.67 | 34.0 | 8.00 | 2.9 | 26.1 | 82.1 | 184 | 60.2 | 40.1 | 27.7 | 68.9 |
+ mean of 11 CS lines. # Phenotypic data of 38 diverse Upland cotton cultivars were obtained from an earlier field study [41] in order to compare results to those for the 11 CS lines. * p < 0.05, ** p < 0.01, *** p < 0.001, ns-nonsignificant.
| CS Line | Pn | Gs | T | iWUE | CTD | CMTS | CSI | SLA | PV | PG30 | PGR |
|---|---|---|---|---|---|---|---|---|---|---|---|
| CS-B01 | −11.9 | −10.2 | −10.2 | +12.6 | |||||||
| CS-B04 | +4.9 | −6.0 | +8.3 | +15.8 | |||||||
| CS-B07 | −13.6 | ||||||||||
| CS-B15sh | +4.1 | +0.20 | −11.3 | +8.6 | |||||||
| CS-B18 | −7.2 | +16.5 | |||||||||
| CS-T01 | +5.7 | −11.1 | −19.0 | −19.0 | +20.6 | ||||||
| CS-T04 | +0.42 | −9.0 | −22.8 | −22.8 | +27.9 | ||||||
| CS-T07 | +4.9 | +0.23 | −7.6 | −12.1 | −12.1 | +11.3 | |||||
| CS-T08sh | −5.0 | −10.9 | +19.1 | ||||||||
| CS-T15sh | +4.2 | +0.21 | +2.60 | −6.8 | |||||||
| CS-T18 | +4.6 | +6.9 | −10.6 | −17.8 | −17.8 | +25.1 |
Dunnett’s test was performed on trait means to examine differences between each of the CS lines and TM-01. +, significant positive performance; −, significant negative performance compared with TM-1.
| Heat Tolerant | Drought Tolerant | Combined Heat and Drought Tolerant | ||||||
|---|---|---|---|---|---|---|---|---|
| Low | Moderate | High | Low | Moderate | High | Low | Moderate | High |
| CS-T01 (11.2) | TM-1 (12) | CS-B15sh (12.4) | CS-B04 (10.2) | CS-T04 (10.8) | CS-T07 (11.4) | CS-T01 (12.3) | CS-T04 (12.8) | CS-T07 (13.4) |
| CS-T15sh (11.1) | CS-T04 (11.9) | CS-T07 (12.4) | CS-T15sh (10.1) | CS-B01 (10.4) | CS-B15sh (11.3) | CS-T15sh (12.1) | CS-B18 (12.5) | CS-B15sh (13.4) |
| CS-B07 (11) | CS-B01 (11.5) | CS-T18 (12.1) | CS-B07 (9.9) | CS-B18 (10.3) | CS-T18 (11.1) | CS-B07 (11.9) | CS-B01 (12.4) | CS-T18 (13.1) |
| CS-T08sh (10.9) | CS-B18 (11.5) | CS-T08sh (9.6) | CS-T01 (10.3) | TM-1 (11) | CS-T08sh (11.8) | CS-B04 (12.4) | TM-1 (13) | |
| CS-B04 (11.4) | ||||||||
The relative performance index of chromosome substitution (CS) lines was calculated as the value for a CS line divided by the phenotypic value of TM-1. The drought stress response index (DSRI) is the sum of the individual relative index of five gas exchange (Pn + Gs + T + iWUE + Ci/Ca), two pigments (TCHL + Caro), two cell and chlorophyll stability indices (CMTS + CSI), and two biophysical (CTD + SLA) traits. The heat stress response index (HSRI) is the sum of individual gas exchange (Pn + Gs + T + iWUE + Ci/Ca), pigments (TCHL + Caro), cell and chlorophyll stability indices (CMTS + CSI), biophysical (CTD), and pollen vitality response (PV + PG38/PG30) traits. The combined heat and drought stress response index (CHDSRI) is the sum of all 13 individual response indexes (see Materials and methods). The different CS lines were classified as low, moderate, and tolerant based on minimum HSRI or DSRI or CHDSI values, and an increment of one standard deviation, respectively, described by Singh et al. [41].
| CS Line | SSR Markers * | Non Target (>5 Mb) ** | Non Target (<5 Mb) *** | Target **** |
|---|---|---|---|---|
| CS-B01 | BNL2921-158, BNL3580-207, BNL3848-197, BNL-3888-197, CIR009-233, NAU2437-235 | c06, c25 | c09, c25 | Substituted |
| CS-B04 | BNL2821 = 192, BNL3089-124, BNL3988-120, CIR222-276, CIR249-190 | - | c18; c20 | Substituted |
| CS-B07 | - | c11, c16, | c04, c16, c20, c21, c22 | TM-1 |
| CS-B15sh | - | c07 | c02, c10; c13, c20 | TM-1 |
| CS-B18 | BNL0193-111, BNL3479-251, CIR216-141, Gh501-201, TMB2762-205, | - | c12, c16, c17, c19, c24, c25, c26 | Substituted |
| CS-T01 | BNL1693-244, BNL2921-162, BNL3580 = 202, BNL3848-223, BNL3888-182, BNL2437-247 | c25 | c03, c07, c10, c20, c25 | Substituted |
| CS-T04 | - | - | c18 | Substituted |
| CS-T07 | BNL1395-163, BNL1597-228, CIR141-170, CIR169-135, Gh548-113, | c02, c21 | c05, c08 | Substituted |
| CS-T08sh | - | c12, c20 | c05, c06, c09, c13, c21, c25 | Substituted |
| CS-T15sh | - | c09, c12 | c08, c14, c17 | Substituted |
| CS-T18 | BNL0193-113, BNL1721-171, BNL243-135, BNL3479-244, BNL569-142, CIR216-141, CIR216-141 | c01 | c26 | Substituted |
* The marker name- marker size in bp [78]; ** single nucleotide polymorphism (SNP) markers associated with more than 5 Mb size [77]; *** SNP markers associated with less than 5 Mb size [77]; and **** confirmation on the cytological and molecular evidence concerning the substituted chromosome or chromosome segment.
Supplementary Materials
The following are available online at https://www.mdpi.com/2223-7747/9/12/1747/s1, Table S1: The relative performance of cotton chromosome substitution (CS) lines compared with parent Upland cultivar, TM-1.
Author Contributions
K.R.R.; conceptualization, project administration, methodology, supervision, funding acquisition, review, co-writing and editing, and data curation, R.B.; formal analysis, review, co-writing-original draft, editing, and data visualization. B.G., S.B.L. and K.S.; data collection and editing. S.S., J.J.R. and J.N.J.; Field experimental design, editing, and review. D.A.R. and R.N.V.; CS line development. A.M.H.-K.; SNP genotyping; L.M.D.S.; SNP genome alignment and analysis and summarization; D.M.S.; CS line development, data analysis, SNP genome analysis, co-writing, editing, and review. All authors have read and agreed to the published version of the manuscript.
Funding
We thank the Cotton Incorporated and the National Institute of Food and Agriculture, NIFA 2016-34263-25763, and MIS 043050 for partial funding to carry out this research.
Acknowledgments
We thank David Brand for technical assistance and graduate students at the Environmental Plant Physiology Laboratory for their help during data collection. Kelli Kochan and Drew Hillhouse at the TIGSS laboratory helped run the CottonSNP63K arrays.
Conflicts of Interest
The authors declare no conflict of interest.
Abbreviations
| Car | Carotenoids |
| Ci/Ca | Ratio of leaf internal to ambient carbon dioxide concentration |
| CMTS | Cell membrane thermostability |
| CSI | Chlorophyll stability index |
| CTD | Canopy temperature depression |
| DSRI | Drought stress response index |
| Gs | Stomatal conductance |
| HSRI | Heat stress response index |
| iWUE | Instantaneous leaf water use efficiency |
| PG | Pollen germination |
| PGR | Pollen germination response |
| Pn | Leaf net photosynthesis |
| PV | Pollen viability |
| SLA | Specific leaf area |
| SNP | Single nucleotide polymorphism |
| SSR | Simple sequence repeat |
| TCHL | Total leaf chlorophyll |
| T | Transpiration |
1. ICAC. Global warming and cotton production-Part 2. ICAC Rec. 2009, 27, 9-13.
2. Zahid, K.R.; Ali, F.; Shah, F.; Younas, M.; Shah, T.; Shahwar, D.; Hassan, W.; Ahmad, Z.; Qi, C.; Lu, Y. Response and tolerance mechanism of cotton Gossypium hirsutum L. to elevated temperature stress: A review. Front. Plant Sci. 2016, 7, 937.
3. IPCC. Climate Change 2014: Synthesis Report. Contribution of Working Groups I, II and III to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change; IPCC: Geneva, Switzerland, 2014; p. 151.
4. Lobell, D.B.; Schlenker, W.; Costa-Roberts, J. Climate trends and global crop production since 1980. Science 2011, 333, 616-620.
5. Suzuki, N.; Rivero, R.M.; Shulaev, V.; Blumwald, E.; Mittler, R. Abiotic and biotic stress combinations. New Phytol. 2014, 203, 32-43.
6. Mittler, R. Abiotic stress, the field environment and stress combination. Trends Plant Sci. 2006, 11, 15-19.
7. FAO. World Food and Agriculture; Food and Agriculture Organization: Rome, Italy, 2015.
8. Kukal, M.S.; Irmak, S. Climate-driven crop yield and yield variability and climate change impacts on the US Great Plains agricultural production. Sci. Rep. 2018, 8, 1-18.
9. Ehlig, C.; LeMert, R. Effects of Fruit Load, Temperature, and relative humidity on boll retention of cotton 1. Crop Sci. 1973, 13, 168-171.
10. Khan, A.; Pan, X.; Najeeb, U.; Tan, D.K.Y.; Fahad, S.; Zahoor, R.; Luo, H. Coping with drought: Stress and adaptive mechanisms, and management through cultural and molecular alternatives in cotton as vital constituents for plant stress resilience and fitness. Biol. Res. 2018, 51, 47.
11. Sabagh, A.E.; Hossain, A.; Islam, M.S.; Barutcular, C.; Ratnasekera, D.; Gormus, O.; Amanet, K.; Mubeen, M.; Nasim, W.; Fahad, S. Drought and heat stress in cotton (Gossypium hirsutum L.): Consequences and their possible mitigation strategies. In Agronomic Crops; Springer: Berlin/Hedelberg, Germany, 2020; pp. 613-634.
12. Reddy, K.R.; Hodges, H.F.; Reddy, V.R. Temperature effects on cotton fruit retention. Agron. J. 1992, 84, 26-30.
13. Pettigrew, W. The effect of higher temperatures on cotton lint yield production and fiber quality. Crop Sci. 2008, 48, 278-285.
14. Bheemanahalli, R.; Sathishraj, R.; Tack, J.; Nalley, L.L.; Muthurajan, R.; Jagadish, K.S. Temperature thresholds for spikelet sterility and associated warming impacts for sub-tropical rice. Agric. For. Meteorol. 2016, 221, 122-130.
15. Reddy, K.R.; Reddy, V.R.; Hodges, H.F. Temperature effects on early season cotton growth and development. Agron. J. 1992, 84, 229-237.
16. Soliz, L.M.A.; Oosterhuis, D.M.; Coker, D.L.; Brown, R.S. Physiological response of cotton to high night temperature. Am. J. Sci. Biotechnol. 2008, 2, 63-68.
17. Jain, M.; Prasad, P.V.; Boote, K.J.; Hartwell, A.L.; Chourey, P.S. Effects of season-long high temperature growth conditions on sugar-to-starch metabolism in developing microspores of grain sorghum (Sorghum bicolor L. Moench). Planta 2007, 227, 67-79.
18. Saini, H.; Sedgley, M.; Aspinall, D. Effect of heat stress during floral development on pollen tube growth and ovary anatomy in wheat (Triticum aestivum L.). Funct. Plant Biol. 1983, 10, 137-144.
19. Kakani, V.G.; Reddy, K.R.; Koti, S.; Wallace, T.; Prasad, P.V.; Reddy, V.R.; Zhao, D. Differences in in vitro pollen germination and pollen tube growth of cotton cultivars in response to high temperature. Ann. Bot. 2005, 96, 59-67.
20. Hedhly, A.; Hormaza, J.; Herrero, M. Effect of temperature on pollen tube kinetics and dynamics in sweet cherry, Prunus avium (Rosaceae). Am. J. Bot. 2004, 91, 558-564.
21. Reddy, V.R.; Reddy, K.R.; Baker, D.N. Temperature effect on growth and development of cotton during the fruiting period. Agron. J. 1991, 83, 211-217.
22. Reddy, K.R.; Hodges, H.F.; McKinion, J.M. A temperature model for cotton phenology. Biotronics 1993, 22, 47-59.
23. Reddy, K.; Hodges, H.; McKinion, J. Carbon dioxide and temperature effects on pima cotton growth. Agric. Ecosyst. Environ. 1995, 54, 17-29.
24. Ulloa, M.; De Santiago, L.M.; Hulse-Kemp, A.M.; Stelly, D.M.; Burke, J.J. Enhancing Upland cotton for drought resilience, productivity, and fiber quality: Comparative evaluation and genetic dissection. Mol. Genet. Genom. 2020, 295, 155-176.
25. Witt, T.W.; Ulloa, M.; Schwartz, R.C.; Ritchie, G.L. Response to deficit irrigation of morphological, yield and fiber quality traits of upland (Gossypium hirsutum L.) and Pima (G. barbadense L.) cotton in the Texas High Plains. Field Crop. Res. 2020, 249, 107759.
26. Saha, S.; Jenkins, J.N.; Wu, J.; McCarty, J.C.; Gutiérrez, O.A.; Percy, R.G.; Cantrell, R.G.; Stelly, D.M. Effects of chromosome-specific introgression in upland cotton on fiber and agronomic traits. Genetics 2006, 172, 1927-1938.
27. Percival, A.; Wendel, J.; Stewart, J. Taxonomy and Germplasm Resources; Agronomy Monographs; John Wiley and Sons, Inc.: New York, NY, USA, 1999.
28. Stelly, D.; Saha, S.; Raska, D.; Jenkins, J.; McCarty, J., Jr.; Gutierrez, O. Registration of 17 upland (Gossypium hirsutum) cotton germplasm lines disomic for different G. barbadense chromosome or arm substitutions. Crop Sci. 2005, 45, 263-265.
29. Jenkins, J.N.; McCarty, J.C.; Campbell, B.T.; Hayes, R.; Wu, J.; Saha, S.; Stelly, D. Genotypic comparisons of chromosomes 01, 04, and 18 from three tetraploid species of Gossypium in topcrosses with five elite cultivars of G. hirsutum L. Euphytica 2017, 213, 107.
30. Jenkins, J.N.; McCarty, J.C., Jr.; Campbell, B.T.; Hayes, R.; Wu, J.; Saha, S.; Stelly, D. Genetic effects of chromosomes 1, 4, and 18 from three tetraploid Gossypium species in topcrosses with five elite cultivars. Crop Sci. 2017, 57, 1338-1346.
31. Saha, S.; Wu, J.; Jenkins, J.N.; McCarty, J.C.; Campbell, B.T.; Hayes, R.W.; Stelly, D.M. Tri-species shuffling of chromosomes to study the effects on fiber traits using chromosome substitution lines. Crop Sci. 2017, 57, 1211-1226.
32. Awasthi, A.; Reddy, K.R.; Saha, S.; Jenkins, J.N.; Stelly, D.M. Morph-physiological responses of cotton interspecific chromosome substitution lines to low temperature and drought stresses. Euphytica 2018, 214, 218.
33. Zhai, H.; Gong, W.; Tan, Y.; Liu, A.; Song, W.; Li, J.; Deng, Z.; Kong, L.; Gong, J.; Shang, H. Identification of chromosome segment substitution lines of Gossypium barbadense introgressed in G. hirsutum and quantitative trait locus mapping for fiber quality and yield traits. PLoS ONE 2016, 11, e0159101.
34. USDA-ERS. USDA ERS-Background [WWW Document]. Available online: https://www.ers.usda.gov/data-products/cotton-wool-and-textile-data/cotton-and-wool-yearbook/ (accessed on 24 August 2020).
35. Singh, R.P.; Prasad, P.V.; Sunita, K.; Giri, S.; Reddy, K.R. Influence of high temperature and breeding for heat tolerance in cotton: A review. Adv. Agron. 2007, 93, 313-385.
36. Bita, C.; Gerats, T. Plant tolerance to high temperature in a changing environment: Scientific fundamentals and production of heat stress-tolerant crops. Front. Plant Sci. 2013, 4, 273.
37. Tuberosa, R. Phenotyping for drought tolerance of crops in the genomics era. Front. Physiol. 2012, 3, 347.
38. Tricker, P.J.; ElHabti, A.; Schmidt, J.; Fleury, D. The physiological and genetic basis of combined drought and heat tolerance in wheat. J. Exp. Bot. 2018, 69, 3195-3210.
39. Abdelraheem, A.; Esmaeili, N.; O'Connell, M.; Zhang, J. Progress and perspective on drought and salt stress tolerance in cotton. Ind. Crop. Prod. 2019, 130, 118-129.
40. Saha, S.; Stelly, D.M.; Raska, D.A.; Wu, J.; Jenkins, J.N.; McCarty, J.C.; Makamov, A.; Gotmare, V.; Abdurakhmonov, I.Y.; Campbell, B. Chromosome substitution lines: Concept, development and utilization in the genetic improvement of Upland cotton. In Plant Breeding; IntechOpen: London, UK, 2012.
41. Singh, K.; Wijewardana, C.; Gajanayake, B.; Lokhande, S.; Wallace, T.; Jones, D.; Reddy, K.R. Genotypic variability among cotton cultivars for heat and drought tolerance using reproductive and physiological traits. Euphytica 2018, 214, 57.
42. Kusumi, K.; Hirotsuka, S.; Kumamaru, T.; Iba, K. Increased leaf photosynthesis caused by elevated stomatal conductance in a rice mutant deficient in SLAC1, a guard cell anion channel protein. J. Exp. Bot. 2012, 63, 5635-5644.
43. Kondamudi, R.; Swamy, K.; Rao, Y.V.; Kiran, T.V.; Suman, K.; Rao, D.S.; Rao, P.R.; Subrahmanyam, D.; Sarla, N.; Kumari, B.R. Gas exchange, carbon balance and stomatal traits in wild and cultivated rice (Oryza sativa L.) genotypes. Acta Physiol. Plant. 2016, 38, 160.
44. Radin, J.W.; Lu, Z.; Percy, R.G.; Zeiger, E. Genetic variability for stomatal conductance in Pima cotton and its relation to improvements of heat adaptation. Proc. Natl. Acad. Sci. USA 1994, 91, 7217-7221.
45. Lu, Z.; Radin, J.W.; Turcotte, E.L.; Percy, R.; Zeiger, E. High yields in advanced lines of Pima cotton are associated with higher stomatal conductance, reduced leaf area and lower leaf temperature. Physiol. Plant. 1994, 92, 266-272.
46. Udayakumar, M.; Sheshshayee, M.; Nataraj, K. Why has breeding for water use efficiency not been successful? An analysis and alternate approach to exploit this trait for crop improvement. Curr. Sci. 1998, 74, 994-1000.
47. Ma, B.; Morrison, M.; Voldeng, H. Leaf greenness and photosynthetic rates in soybean. Crop Sci. 1995, 35, 1411-1414.
48. Elfeky, S.S.; Osman, M.E.; Hamada, S.M.; Hasan, A.M. Effect of salinity and drought on growth criteria and biochemical analysis of Catharanthus roseus shoot. Int. J. Bot. 2007, 3, 202-207.
49. Wang, L.; Deng, F.; Ren, W.-J. Shading tolerance in rice is related to better light harvesting and use efficiency and grain filling rate during grain filling period. Field Crop. Res. 2015, 180, 54-62.
50. Kittock, D.; Turcotte, E.; Hofmann, W. Estimation of heat tolerance improvement in recent American pima cotton cultivars. J. Agron. Crop Sci. 1988, 161, 305-309.
51. Reddy, K.R.; Hodges, H.F.; McKinion, J.M. A comparison of scenarios for the effect of global climate change on cotton growth and yield. Funct. Plant Biol. 1997, 24, 707-713.
52. Mohan, M.; Narayanan, S.L.; Ibrahim, S. Chlorophyll stability index (CSI): Its impact on salt tolerance in rice. Int. Rice Res. Notes 2000, 25, 38-39.
53. ur Rahman, H.; Malik, S.A.; Saleem, M. Heat tolerance of upland cotton during the fruiting stage evaluated using cellular membrane thermostability. Field Crop. Res. 2004, 85, 149-158.
54. Rebetzke, G.; Botwright, T.; Moore, C.; Richards, R.; Condon, A. Genotypic variation in specific leaf area for genetic improvement of early vigour in wheat. Field Crop. Res. 2004, 88, 179-189.
55. Kiran, T.; Rao, Y.; Subrahmanyam, D.; Rani, N.; Bhadana, V.; Rao, P.; Voleti, S. Variation in leaf photosynthetic characteristics in wild rice species. Photosynthetica 2013, 51, 350-358.
56. Nageswara Rao, R.; Talwar, H.; Wright, G. Rapid assessment of specific leaf area and leaf nitrogen in peanut (Arachis hypogaea L.) using a chlorophyll meter. J. Agron. Crop Sci. 2001, 186, 175-182.
57. Sheshshayee, M.; Bindumadhava, H.; Rachaputi, N.; Prasad, T.; Udayakumar, M.; Wright, G.; Nigam, S. Leaf chlorophyll concentration relates to transpiration efficiency in peanut. Ann. Appl. Biol. 2006, 148, 7-15.
58. Songsri, P.; Jogloy, S.; Holbrook, C.; Kesmala, T.; Vorasoot, N.; Akkasaeng, C.; Patanothai, A. Association of root, specific leaf area and SPAD chlorophyll meter reading to water use efficiency of peanut under different available soil water. Agric. Water Manag. 2009, 96, 790-798.
59. Richardson, A.D.; Duigan, S.P.; Berlyn, G.P. An evaluation of noninvasive methods to estimate foliar chlorophyll content. New Phytol. 2002, 153, 185-194.
60. Prasad, P.V.; Bheemanahalli, R.; Jagadish, S.K. Field crops and the fear of heat stress-opportunities, challenges and future directions. Field Crop. Res. 2017, 200, 114-121.
61. Burke, J.J.; Velten, J.; Oliver, M.J. In vitro analysis of cotton pollen germination. Agron. J. 2004, 96, 359-368.
62. Porch, T.; Jahn, M. Effects of high-temperature stress on microsporogenesis in heat-sensitive and heat-tolerant genotypes of Phaseolus vulgaris. Plant Cell Environ. 2001, 24, 723-731.
63. Prasad, P.V.; Craufurd, P.; Summerfield, R. Fruit number in relation to pollen production and viability in groundnut exposed to short episodes of heat stress. Ann. Bot. 1999, 84, 381-386.
64. Cross, R.; McKay, S.; McHughen, A.; Bonham-Smith, P. Heat-stress effects on reproduction and seed set in Linum usitatissimum L. (flax). Plant Cell Environ. 2003, 26, 1013-1020.
65. Chiluwal, A.; Bheemanahalli, R.; Kanaganahalli, V.; Boyle, D.; Perumal, R.; Pokharel, M.; Oumarou, H.; Jagadish, S.K. Deterioration of ovary plays a key role in heat stress-induced spikelet sterility in sorghum. Plant Cell Environ. 2020, 43, 448-462.
66. Peet, M.; Sato, S.; Gardner, R. Comparing heat stress effects on male-fertile and male-sterile tomatoes. Plant Cell Environ. 1998, 21, 225-231.
67. Prasad, P.V.; Boote, K.J.; Allen, L.H., Jr.; Thomas, J.M. Effects of elevated temperature and carbon dioxide on seed-set and yield of kidney bean (Phaseolus vulgaris L.). Glob. Chang. Biol. 2002, 8, 710-721.
68. Vara Prasad, P.; Boote, K.J.; Hartwell Allen, L., Jr.; Thomas, J.M. Super-optimal temperatures are detrimental to peanut (Arachis hypogaea L.) reproductive processes and yield at both ambient and elevated carbon dioxide. Glob. Chang. Biol. 2003, 9, 1775-1787.
69. Young, L.W.; Wilen, R.W.; Bonham-Smith, P.C. High temperature stress of Brassica napus during flowering reduces micro-and megagametophyte fertility, induces fruit abortion, and disrupts seed production. J. Exp. Bot. 2004, 55, 485-495.
70. Salem, M.A.; Kakani, V.G.; Koti, S.; Reddy, K.R. Pollen-based screening of soybean genotypes for high temperatures. Crop Sci. 2007, 47, 219-231.
71. Bheemanahalli, R.; Sunoj, V.J.; Saripalli, G.; Prasad, P.V.; Balyan, H.; Gupta, P.; Grant, N.; Gill, K.S.; Jagadish, S.K. Quantifying the impact of heat stress on pollen germination, seed set, and grain filling in spring wheat. Crop Sci. 2019, 59, 684-696.
72. Rodriguez-Garay, B.; Barrow, J.R. Pollen selection for heat tolerance in cotton. Crop Sci. 1988, 28, 857-859.
73. Lamaoui, M.; Jemo, M.; Datla, R.; Bekkaoui, F. Heat and drought stresses in crops and approaches for their mitigation. Front. Chem. 2018, 6, 26.
74. Reddy, A.R.; Reddy, K.R.; Padjung, R.; Hodges, H.F. Nitrogen nutrition and photosynthesis in leaves of Pima cotton. J. Plant Nutr. 1996, 19, 755-770.
75. Reddy, K.R.; Hodges, H.F.; Kimball, B.A. Crop ecosystem responses to climatic change: Cotton. In Climate Change and Global Crop Productivity; Reddy, K.R., Hodges, H.F., Eds.; CABI Publishing: Wallingford, UK, 2000; pp. 161-187.
76. Zhao, D.; Reddy, K.R.; Kakani, V.G.; Mohammed, A.R.; Read, J.J.; Gao, W. Leaf and canopy photosynthetic characteristics of cotton (Gossypium hirsutum) under elevated CO2 concentration and UV-B radiation. J. Plant Physiol. 2004, 161, 581-590.
77. Hulse-Kemp, A.M.; Lemm, J.; Plieske, J.; Ashrafi, H.; Buyyarapu, R.; Fang, D.D.; Frelichowski, J.; Giband, M.; Hague, S.; Hinze, L.L. Development of a 63K SNP array for cotton and high-density mapping of intraspecific and interspecific populations of Gossypium spp. G3 Genes Genomes Genet. 2015, 5, 1187-1209.
78. Saha, S.; Stelly, D.M.; Makamov, A.K.; Ayubov, M.S.; Raska, D.; Gutiérrez, O.A.; Manchali, S.; Jenkins, J.N.; Deng, D.; Abdurakhmonov, I.Y. Molecular confirmation of Gossypium hirsutum chromosome substitution lines. Euphytica 2015, 205, 459-473.
79. Hendrix, B.; Stewart, J.M. Estimation of the nuclear DNA content of Gossypium species. Ann. Bot. 2005, 95, 789-797.
80. Chen, Z.J.; Sreedasyam, A.; Ando, A.; Song, Q.; De Santiago, L.M.; Hulse-Kemp, A.M.; Ding, M.; Ye, W.; Kirkbride, R.C.; Jenkins, J. Genomic diversifications of five Gossypium allopolyploid species and their impact on cotton improvement. Nat. Genet. 2020, 52, 525-533.
81. Fang, Y.; Xiong, L. General mechanisms of drought response and their application in drought resistance improvement in plants. Cell. Mol. Life Sci. 2015, 72, 673-689.
82. Luo, L. Breeding for water-saving and drought-resistance rice (WDR) in China. J. Exp. Bot. 2010, 61, 3509-3517.
83. Manavalan, L.P.; Guttikonda, S.K.; Phan Tran, L.-S.; Nguyen, H.T. Physiological and molecular approaches to improve drought resistance in soybean. Plant Cell Physiol. 2009, 50, 1260-1276.
84. Gajanayake, B.; Trader, B.W.; Reddy, K.R.; Harkess, R.L. Screening ornamental pepper cultivars for temperature tolerance using pollen and physiological parameters. HortScience 2011, 46, 878-884.
85. Wijewardana, C.; Henry, W.B.; Reddy, K.R. Evaluation of drought tolerant maize germplasm to induced drought stress. J. Miss. Acad. Sci. 2017, 62, 316-329.
86. Wagner, W.; Herbst, D.; Sohmer, S. Manual of the Flowering Plants of HAWAII; Two Volumes; Bernice Pauahi Bishop Museum Special Publication; University of Hawaii Press: Honolulu, HI, USA, 1990; Volume 83.
87. Lee, J.A.; Fang, D.D. Cotton as a world crop: Origin, history, and current status. Cotton 2015, 57, 1-23.
88. Morrow, M.; Krieg, D. Cotton management strategies for a short growing season environment: Water-nitrogen considerations. Agron. J. 1990, 82, 52-56.
89. Bowman, D.T.; Gutiérrez, O.A. Sources of fiber strength in the US upland cotton crop from 1980 to 2000. J. Cotton Sci. 2003, 7, 64-69.
90. Loka, D.A.; Oosterhuis, D.M. Water Stress and Reproductive Development in Cotton. In Flowering and Fruiting in Cotton; Oosterhuis, D.M., Cothren, J.T., Eds.; Cotton Foundation: Memphis, TN, USA, 2012; pp. 51-58.
91. Grimes, D.W.; Yamada, H.; Dickens, W. Functions for cotton (Gossypium hirsutum L.) production from irrigation and nitrogen fertilization variables: I. Yield and evapotranspiration. Agron. J. 1969, 61, 769-773.
92. Johnson, R.M.; Downer, R.G.; Bradow, J.M.; Bauer, P.J.; Sadler, E.J. Variability in cotton fiber yield, fiber quality, and soil properties in a southeastern coastal plain. Agron. J. 2002, 94, 1305-1316.
93. Lokhande, S.; Reddy, K.R. Reproductive and fiber quality responses of upland cotton to moisture deficiency. Agron. J. 2014, 106, 1060-1069.
94. Kakani, V.; Prasad, P.; Craufurd, P.; Wheeler, T. Response of in vitro pollen germination and pollen tube growth of groundnut (Arachis hypogaea L.) genotypes to temperature. Plant Cell Environ. 2002, 25, 1651-1661.
95. Kohel, R.J.; Yu, J.; Park, Y.-H.; Lazo, G.R. Molecular mapping and characterization of traits controlling fiber quality in cotton. Euphytica 2001, 121, 163-172.
96. Aslam, M.; Brown, M.; Kohel, R. Evaluation of seven tetrazolium salts as vital pollen stains in cotton Gossypium hirsutum L. 1. Crop Sci. 1964, 4, 508-510.
97. Taylor, R.M. Germination of cotton (Gossypium hirsutum L.) Pollen on an artificial medium 1. Crop Sci. 1972, 12, 243-244.
98. Chapple, C.; Vogt, T.; Ellis, B.E.; Somerville, C.R. An Arabidopsis mutant defective in the general phenylpropanoid pathway. Plant Cell 1992, 4, 1413-1424.
Kambham Raja Reddy
1,*,
Raju Bheemanahalli
1,
Sukumar Saha
2,*,
Kulvir Singh
1,
Suresh B. Lokhande
1,
Bandara Gajanayake
1,
John J. Read
2,
Johnie N. Jenkins
2,
Dwaine A. Raska
3,
Luis M. De Santiago
3,
Amanda M. Hulse-Kemp
3,4,
Robert N. Vaughn
3 and
David M. Stelly
3
1Department of Plant and Soil Sciences, Mississippi State University, Mississippi State, MS 39762, USA
2USDA-ARS, Genetics and Sustainable Agriculture Research Unit, Mississippi State, MS 39762, USA
3Department of Soil and Crop Sciences, Texas A&M AgriLife Research, College Station, TX 77843, USA
4USDA-ARS, Genomics and Bioinformatics Research Unit, Raleigh, NC 27695, USA
*Authors to whom correspondence should be addressed.
© 2020. This work is licensed under http://creativecommons.org/licenses/by/3.0/ (the “License”). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.