Changes in Leaf Structural and Functional Characteristics when Changing Planting Density at Di ﬀ erent Growth Stages Alters Cotton Lint Yield under a New Planting Model

: Manipulation of planting density and choice of variety are e ﬀ ective management components in any cropping system that aims to enhance the balance between environmental resource availability and crop requirements. One-time fertilization at ﬁrst ﬂower with a medium plant stand under late sowing has not yet been attempted. To ﬁll this knowledge gap, changes in leaf structural (stomatal density, stomatal length, stomata width, stomatal pore perimeter, and leaf thickness), leaf gas exchange, and chlorophyll ﬂuorescence attributes of di ﬀ erent cotton varieties were made in order to change the planting densities to improve lint yield under a new planting model. A two-year ﬁeld evaluation was carried out on cotton varieties—V 1 (Zhongmian-16) and V 2 (J-4B)—to examine the e ﬀ ect of changing the planting density (D 1 , low, 3 × 10 4 ; D 2 , moderate, 6 × 10 4 ; and D 3 , dense, 9 × 10 4 ) on cotton lint yield, leaf structure, chlorophyll ﬂuorescence, and leaf gas exchange attribute responses. Across these varieties, J-4B had higher lint yield compared with Zhongmian-16 in both years. Plants at high density had depressed leaf structural traits, net photosynthetic rate, stomatal conductance, intercellular CO 2 uptake, quenching (qP), actual quantum yield of photosystem II ( Φ PSII), and maximum quantum yield of PSII ( F v / F m ) in both years. Crops at moderate density had improved leaf gas exchange traits, stomatal density, number of stomata, pore perimeter, length, and width, as well as increased qP, Φ PSII, and F v / F m compared with low- and high-density plants. Improvement in leaf structural and functional traits contributed to 15.9%–10.7% and 12.3%–10.5%


Introduction
Cotton (Gossypium hirsutum L.) is a natural white fiber and cash crop that is grown globally [1]. The cotton plant is characterized by indeterminate growth habits and shows morphological and physiological adaptation to a wide range of environmental and management practices, including planting density and cultivar. An expanding population necessitates global efforts to increase crop production, especially those fulfilling food and fiber needs. Currently, numerous management practices have been introduced for cotton production systems, but lint production per unit area has remained stagnant [2]. High input costs combined with multiple management and material inputs have threatened cotton productivity. [2]. An efficient agricultural production system characterized by moderate planting density with one-time fertilization under a short growing season can reduce inputs without yield loss [3,4].
Planting density and choice of cultivar are important agronomic practices that have the potential to optimize the canopy photosynthetic rate and crop productivity of any cropping system [5]. Changes in plant architecture and canopy dynamics in response to planting density can have impacts on disease incidence, water use, canopy temperature, and enzymatic activity of assimilate metabolism [6]. Manipulations of planting density in cotton have remarkable impacts on biomass partitioning, nutrient uptake, boll distribution, changes in the light spectrum, and crop production [3,5,7,8], which can influence yield and profits for producers [9]. Plants at high density can minimize evaporation and irrigation frequency, as well as increase the utilization of irrigation water [10,11]. In contrast, high-density planting can slow down leaf appearance and reduce open boll density [12], boll weight, and boll number [7]. It also delays leaf senescence [13] and decreases nitrogen use efficiency and nitrogen recovery efficiency [14]. A planting density of up to nine plants m −2 has been reported to sustain leaf photosynthetic rate and reproductive organ biomass formation by increasing plant potassium uptake at various developmental stages. However, a sowing density of >10 plants m −2 and subsequent shading can result in disease infestation, small boll size, fruit shedding, delayed maturity, and decreased individual plant development [4]. Dense planting can also delay crop maturity by promoting vegetative growth and can substantially depress net photosynthetic rate [4] due to decreased RuBP carboxylase activity and chlorophyll content [15]. High planting density can increase the auxin (IAA) content and enhance auxin polar transport by increasing the expression of the auxin biosynthesis gene (GhYUC5) and the auxin polar transport gene (GhPIN1). It can also inhibit vegetative branching by decreasing IAA, cytokinin, gibberellic acid, and brassinosteroid contents, followed by increased strigolactone content due to differential expression of hormone-associated genes in the tips of vegetative branches [15]. Optimal plant density can ensure healthy plant development by maintaining a core population of plants synchronizing boll number and fiber quality to achieve optimal yield [16].
Leaf morphological and physiological attributes are important players in photosynthetic regulation [4] and can provide a structural framework for gas exchange as well as optimize the photosynthetic function [17]. Cotton leaf surface characteristics, including cuticular thickening, wax layer, and trichomes, play critical roles in the variability of optical properties [18]. Generally, leaves developed under high sunlight can have thicker and smaller leaves with well-developed plastid tissues, greater stomatal density, and smaller granal stacks than shade leaves [19]. Plants under low density planting have a lower chlorophyll content and a higher electron transfer rate and ribulose-1,5-bi-phasphate carboxylase/oxygenase compared with high-density planting [20,21]. Leaves developed under lower density (sun leaves) are tolerant to strong light; conversely, shade leaves have weak photoprotection potential and are more sensitive to high light [17,21].
Studies regarding cotton growth and lint yield in response to diverse populations are common [22][23][24]. However, we are the first to report the effects of changing the planting density on cotton lint yield, leaf structure, chlorophyll fluorescence, and leaf gas exchange characteristics in subtropical regions. The objectives of this study were to investigate leaf structural and functional characteristics in response to different planting densities and varieties. It also explored optimal plant density and variety for improved lint yield in subtropical regions.

Plant Material and Experimental Site
Seeds of two cotton cultivars-V 1 (Zhongmian-16) and V 2 (J-4B)-were procured from the Cotton Research Institute, Chinese Academy of Agricultural Sciences and were grown under field conditions for two years. A replicated two-year (2017 and 2018 growing seasons) field experiment was conducted at Guangxi University, Nanning, China. The soil properties of the experimental field were sandy loam and yellowish, having a pH of 6.5; organic matter of 23.37 mg kg −1 ; and available nitrogen, phosphorus, and potassium content of 53.24, 77.58, and 6.30 mg kg −1 , respectively. The experimental design layout was a balanced split plot with three replications.

Crop Management and Experimental Design
Before sowing, the experimental field was ploughed, laser leveled approximately three weeks prior, and covered with plastic film to conserve moisture and suppress weed germination. The experiments were designed in a split plot arrangement with three replications of each of the six treatment combinations. Two cotton varieties (V 1 , Zhongmian-16; V 2 , J-4B), were randomly allocated to the main plots and three plant population levels (D 1 , low, 3 × 10 4 ; D 2 , medium, 6 × 10 4 ; and D 3 , dense, 9 × 10 4 ha −1 ) were randomized in subplots. By increasing the precision of comparisons, split plot arrangements were adopted. Seeds were sown on 5 June in double rows on each raised plot (3.0 m wide and 11 m long), with a total plot size of 33.0 m 2 . Each subplot was 11 m long and 1.5 m wide, consisting of four rows with narrow (10 cm) and wide (66 cm) row spaces for a total of eight rows on each main plot. Plant-to-plant spacing was controlled according to the corresponding population level. Crops were irrigated one day after sowing to ensure uniform germination. Cotton seedlings were hand-thinned at the third leaf stage to the target population level for each plot. A basal application of phosphorus (P 2 O 5 ) at 66 kg ha −1 , nitrogen (N) at 170 N kg ha −1 , and potash (K 2 O) at 190 kg ha −1 was applied using superphosphate (12% P 2 O 5 ), urea (46% N), and potassium chloride (59% K 2 O), respectively, during the pinhead stage. A plant growth regulator (i.e., mepiquat chloride) at the rate of 0.057 active ingredient ha −1 was sprayed to control vegetative growth. All the necessary field management practices were performed according to crop requirements during the whole crop cycle.

Data Collection
Data on leaf structure, chlorophyll fluorescence, leaf gas exchange attributes, cotton yield, and yield contributors were recorded for each treatment in three replications. The details of each measurement are given below.

Yield and Yield Components
To assess cotton yield, fully opened bolls were hand-picked at three times in each treatment. The harvested seed cotton was sun-dried to ≤11% moisture content [16]. The seed cotton was ginned to calculate seed cotton and lint yield. During the second picking, 100 mature bolls were manually picked to calculate single boll mass and lint percentage. Seed cotton yield of 100 bolls was divided by the number of bolls to assess individual boll weight. Lint % was determined using the lint yield of 100 bolls and divided by seed cotton mass.

Cotton Leaf Structure Attributes
Ten plants in each plot were randomly tagged to measure leaf structure and plant growth characteristics at the boll setting stage. Leaf thickness was determined on 10 fully expanded leaves from the upper part of three plants (functional leaves, i.e., upper fourth leaf). A hand-held micrometer (Mitutoyo Digital Micrometer Model 293-185, Kawasaki, Japan) with a digital display and a clutch that ensured uniform pressure [25] was used for leaf thickness assessment. A 5 × 8 mm leaf section was removed for each treatment. Samples were then added into 10 mL tubes containing 50%, 5%, and 5% alcohol solution, formaldehyde, and glacial acetic acid, respectively. Scanning electron microscopy was performed at Guanxi Medical University using a SUPRA 55VP (Carl Zeiss AG, Oberkochen, Germany). Image software was used to assess cotton leaf stomatal length, width, density, and pore perimeter according to the method reported in [26].  [27]. The effective quantum yield of PSII photochemistry of light-adapted leaves was determined by Φ PSII (F m ' − F)/F m ' [28]. Coefficient of photochemical quenching (qP) was assessed using the formula qP = [29]. Minimal fluorescence of light-adapted leaves (F o ') was calculated according to the equation [28]. Nonphotochemical quenching (NPQ) was recorded according to [6] as NPQ = (F m − F m ')/F m ', where F m represents the value of the predawn observations. The electron transport rate (ETR) was assessed using a leaf absorptance of 0.85 and half of the absorbed light was partitioned to each photosystem: ETR = PSII × PPFD × 0.85 × 0.5 [30].

Leaf Gas Exchange Attributes
At squaring, flowering, peak bloom, and boll setting stages, fully expanded leaves from the upper part of three plants (functional leaves, i.e., upper fourth leaf) were chosen to assess net photosynthetic rate (Pn), stomatal conductance (g s ), intercellular CO 2 concentration (C i ), and transpiration rate (E). Net rate of photosynthesis was measured from the six functional leaves of three plants in each treatment using a portable infrared gas exchange analyzer (Li-6400, Li-Cor, Lincoln, NE, USA). These observations were made on a clear day between 10:00 a.m. and 12:00 p.m. Beijing time in each experimental unit of four replications. Leaves in each plot followed the following adjustments: PAR, 1800 µmol m −2 s −2 ; air flow, 389.42 mmol −1 m −2 s −1 ; water vapor pressure into leaf chamber, 3.13 mbar; leaf temperature, 30 • C; ambient temperature, 33.69 • C; and ambient carbon dioxide concentration, 330-350 mol mol −1 .

Statistical Analysis
All the data were processed using Microsoft Excel 2016. Figures were plotted using Sigma Plot 14.00 software. Analysis of variance was implemented using SAS software (version 8.1, SAS Institute, Cary, NC, USA). The initial combined data showed no interactions with years. Therefore, the data were pooled and presented across the two years. Means of planting density were separated using the least significant difference (LSD) test at the 5% probability level. Both planting densities and cultivars were taken as main factors and fixed effects with cropping season as the repetitive measured factor with a fixed effect. Similarly, the interaction was taken as fixed effects and treatment × replication interaction, which was taken as a random effect. Differences among treatments imply statistical difference (p = 0.05).

Yield and Yield Components
The analysis of variance (Table 1) showed that effects of year, planting density, variety, and their interaction on cotton yield and yield contributors. The year effect was statistically significant, but the differences were not large. Planting density and variety did not affect lint percentage and boll weight. There were 14.5% and 7.1% more bolls m −2 with a 19% and 11.5% higher lint yield in moderate-density crops compared with low-and high-density crops, respectively. Under high-density conditions, a reduction of 9.6% and 2.3% was noted in boll weight and lint percentage, respectively, compared with low-and moderate-density crops. Across the varieties, J-4B produced 6% and 7.8% greater bolls m −2 and lint yield, respectively, compared with the Zhongmian-16 variety ( Figure 1A-D). Interaction was significant for density × variety across two years. Cotton plant individual boll weight, boll density m −2 , and lint yield were highest under moderate-to-high planting density for J-4B, while under low-density conditions, Zhongmian-16 had a higher boll weight.

Yield and Yield Components
The analysis of variance (Table 1) showed that effects of year, planting density, variety, and their interaction on cotton yield and yield contributors. The year effect was statistically significant, but the differences were not large. Planting density and variety did not affect lint percentage and boll weight. There were 14.5% and 7.1% more bolls m −2 with a 19% and 11.5% higher lint yield in moderate-density crops compared with low-and high-density crops, respectively. Under high-density conditions, a reduction of 9.6% and 2.3% was noted in boll weight and lint percentage, respectively, compared with low-and moderate-density crops. Across the varieties, J-4B produced 6% and 7.8% greater bolls m −2 and lint yield, respectively, compared with the Zhongmian-16 variety ( Figure 1A-D). Interaction was significant for density × variety across two years. Cotton plant individual boll weight, boll density m −2 , and lint yield were highest under moderate-to-high planting density for J-4B, while under low-density conditions, Zhongmian-16 had a higher boll weight.

Chlorophyll Fluorescence Traits
During both years, planting densities, varieties, and their interaction had significant impacts on chlorophyll fluorescence traits in different growth stages (Tables 3-5). Except the squaring stage, ΦPSII at first bloom, peak bloom, and boll setting stages were increased by moderate-density compared with lowand high-density crops, while the F v /F m yield was greater at all growth stages (Table 3). Across the varieties, J-4B had higher ΦPSII and F v /F m at peak bloom and boll setting stages than Zhongmian-16, respectively. The interaction between density × variety remained significant for ΦPSII and F v /F m at different growth stages. The J-4B variety with moderate crops had greater ΦPSII and F v /F m across the years.
Significant variation between planting densities, varieties, and years was found for photochemical quenching (qP) and nonphotochemical quenching (NPQ) of cotton at all growth stages (Table 4). Across densities, medium competitive plants yielded higher qP and NPQ rates. The variety J-4B resulted in higher qP at squaring and boll setting stages, while Zhongmian-16 had higher NPQ rates at first and full bloom stages. The interaction showed that J-4B had a higher qP under moderate density at different growth stages. J-4B had higher values for NPQ at low density compared with Zhongmian-16, followed by moderate density for the same variety at the peak bloom stage.
Significant differences existed between years, densities, and varieties for the ETR at four growth stages (Table 5). Interaction between density × variety revealed substantial variation between varieties to planting density at all growth stages. Increased planting density substantially reduced ETR at all growth stages in both years ( Table 5). The low-density plants improved ETR at squaring, first, peak bloom, and boll setting stages, followed by moderate density, while there was a lower ETR in high-density crops. A higher ETR was noted for the variety Zhongmian-16 at squaring, first, and peak bloom stages compared with J-4B; however, J-4B had a higher ETR at the boll setting stage than Zhongmian-16. ETR values were substantially reduced under high density for both varieties. Lower planting density had higher ETR values for Zhongmian-16 or J-4B during both years.

Leaf Gas Exchange Attributes
Cotton leaf gas attributes were significantly influenced by plant density, variety, and growing year (Tables 6 and 7). Under moderate-density conditions, net photosynthetic rate (Pn) was increased at all growth stages except squaring, while stomatal conductance (g s ) was higher at the first bloom and boll setting stages. Plants under high density had significantly lower Pn and g s compared with low and moderate density ( Table 6). J-4B had higher Pn and g s compared with Zhongmian-16 under moderate density. Interaction between density × variety was significant only at full bloom and boll setting for Pn and at the peak bloom stage for g s . J-4B under low-to-moderate planting density had a higher Pn at squaring and first bloom stages, while it was higher in Zhongmian-16 at the peak bloom and boll set stages. A higher g s under moderate planting density was noted in J-4B at the peak stage than Zhongmian-16 at low or high density.  Increasing planting density significantly reduced C i in cotton leaves for both varieties. Plants with moderate density had higher C i uptake at first bloom and peak bloom stages compared with lowand high-density crops, respectively (Table 7). Plants under low density resulted in a higher rate of E during first bloom, peak bloom, and boll set stages compared with moderate-and high-density crops, respectively (Table 7). Across the varieties, Zhongmian-16 yielded higher for both C i uptake and E rates compared with J-4B. Interaction between density × variety remained significant at all growth stages for C i . The transpiration rate was decreased in both varieties when the planting density increased.

Discussion
The current study has provided new data on the common perception that high planting density significantly decreases leaf structural characteristics, such as stomatal density, length, width, pore perimeter, and leaf thickness, as well as functional traits (leaf gas exchange and chlorophyll fluorescence traits), which leads to lint yield loss. However, we found that improved leaf functional and structural traits for J-4B under moderate density had a higher lint yield. Under high-density treatment, reductions in lint yield for Zhongmian-16 were associated with repression in leaf structural and functional attributes, which in turn caused depression in leaf photosynthetic capacity due to nutrient competition. The difference between varieties from changing planting density might be associated with canopy architecture and genetic variation. Therefore, these changes in varieties might have significant impacts on leaf structural and functional attributes and, ultimately, on yield formation.
High planting density responses to cotton lint yield, growth, biomass production, nutrient uptake, and fiber quality have been extensively investigated [3,4,13,22]. The mechanisms of interplant competitiveness under low-to-high planting density on leaf structure, chlorophyll fluorescence, and leaf gas exchange attributes for optimal cotton lint yield have not yet been reported. Across densities, the moderate population had a higher boll number m −2 with improved lint yield for J-4B compared with Zhongmian-16 across two years. High-density plants substantially reduced yield and yield components in both years, probably due to competition for nutrients. The phenomenon of increased lint yield under moderate density can be associated with improved leaf structural and chlorophyll fluorescence traits and higher leaf photosynthetic capacity, which resulted in higher boll density m -2 compared with other densities.
Moderate density favors dry matter partitioning to the reproductive structures rather than vegetative organs [31] and less fruit shedding compared with denser plants. The reductions in lint yield under high density can be attributed to decreased leaf structural and physiological traits, which were observed in this study. The differences that existed between varieties for yield when changing planting density might be attributable to canopy architecture. Differences in plant canopy architectural traits among varieties have an impact on growth characteristics and lint yield. These data further confirmed that an appropriate selection of variety and optimal density can contribute to successful cotton production. Reducing population density may also have other implications, such as decreased frequency and insecticide inputs per season without any yield loss to increase profit. Moreover, high plant density can substantially depress leaf structural and physiological attributes, which in turn cause a severe yield penalty.
Plants respond to ambient and management interventions via architectural and structural changes. Plant growth and leaf morphological attributes, including stomatal density, size, number of pores, width, length, and leaf thickness features, are pivotal windows regulating leaf photosynthetic capacity [10,25] and offer a structural framework for CO 2 exchange and optimization of photosynthetic activities, which in turn can improve crop yield [17]. In this study, high planting density substantially decreased leaf thickness, stomatal density, leaf length, width, and number of stomatal pores. Limitations in these attributes disrupted the photosynthetic capacity of plants by restricting entry of CO 2 to the mesophyll through the stomata of leaves, which is extremely responsive to light environments. Thus, the exchange of CO 2 by means of stomata might be restricted [32]. Higher stomatal density, thicker leaves, and rapid metabolite transfer between the mesophyll and bundle sheet cells can favor higher leaf photosynthetic capacity [33]. Increasing planting density has been proposed to decrease the stomatal density of wheat leaves [34]. A greater stomatal size can facilitate CO 2 distribution into the leaf due to its conductance being proportional to the square of the effective radius of the stomatal pore, resulting in increased stomatal conductance [35]. However, the responses of leaf structural attributes vary under different abiotic stresses in different plant species or varieties [36]. These data suggest that plants under high-density conditions have significantly decreased leaf morphological characteristics, which might be particularly responsible for depressing leaf photosynthetic capacity.
Chlorophyll fluorescence is a nondestructive evaluation of PSII activity. In plant physiology, this technique is commonly used and has become a classical method for crop improvement, screening of beneficial traits, and linking genomic knowledge to phenological response. Due to the sensitivity of PSII to undesirable ambient conditions, this is a useful method for understanding photosynthetic mechanisms and a good indicator of how plants respond to ambient change [37,38].
ΦPSII is a measure of light energy capture efficiency, which reflects the actual primary sunlight energy conversion efficiency of the PSII reaction center [15]. In this study, ΦPSII substantially declined under high-density conditions. Probably, a lower ΦPSII value under high-density conditions did not efficiently convert photon energy to chemical energy; however, this phenomenon needs further exploration. Under shading conditions, a low ΦPSII may be responsible for depressing Pn due to the adjustment in photochemical reaction centers [39], which was observed in our study. The efficient use of limited light energy and the degree of the PSII reaction center openness can increase, resulting in improved energy conversion efficiency. This is associated with the increase of F v /F m , ΦPSII, and qP at early shading [40]. The maximal photochemical efficiency of PSII (F v /F m ) determines the potential quantum efficiency of PSII [41]. In this study, F v /F m had higher values under low rather than high planting density, which is consistent with [32], and reductions in F v /F m values might be due to the lower values of F m and increased values of F o . The ETR is an important chlorophyll fluorescence attribute affected by the external light environment. The rate of ETR declined from low to high density in this study, which corresponds with [40], and shading can significantly decrease ETR values by affecting PSII photochemical reaction centers and consequently diminish the primary stable quinine acceptor of PSII, leading to a decrease in the activity of photosynthetic electron transport efficiency via PSII [27,33]. NPQ can have critical roles in the nonradiative dissipation of surplus light energy [42]. A low-light environment can cause a reduction in NPQ, possibly associated with reduced light energy [32,43]. In this study, a severe decline in NPQ values was noted under high-density compared with low-density crops. This can be explained as the decreased NPQ being associated with the decreased efficiency of photochemical reactions through the reduced fraction of incident light in photochemical energy utilization, which resulted in lower thermal dissipation in PSII [44]. The rate of photochemical quenching (qP) under dense crops showed a substantial reduction compared with low and moderate densities. Probably, a low-light environment can cause reductions in the amount of pigment and the efficiency of photochemical energy conversion, resulting in the depressed quantum yield of PSII and decreased qP. The qP reflects the efficiency of light quantum harvested by PSII to chemical energy and represents the openness degree of the PSII reaction center, and a greater qP results in greater activity of electron transfer in PSII.
Leaf gas exchange traits can play a central role in biomass formation and the prime determination of cotton lint yield [45]. High planting density results in rapid canopy closure and an increase in radiation interception, which reduces weed competition [46], but this impedes leaf gas exchange traits, leading to yield loss [47]. In the current study, cotton leaf gas exchange parameters were substantially depressed under close planting at different growth stages. Accordingly, high-density conditions resulted in reductions in leaf stomatal density, length, width, pores, and leaf thickness, probably due to mutual shading, which may be responsible for depressing stomatal conductance (g s ) and CO 2 uptake through the stomata, which in turn suppressed the photosynthetic capacity. Plants under high-density conditions can significantly decrease g s and Ci, which can negatively influence the photosynthetic system [38]. The CO 2 concentration plays a central role in net photosynthetic rate (Pn), but this varies across species and ambient conditions [48,49]. The g s might respond to alterations in Pn and thus prevents C i near saturation. The primary function of stomata is to avoid desiccation and enable the passage of CO 2 . Stomata induce a substantial disruption in the CO 2 assimilation rate, which reduces more in C 4 than C 3 plants. The stomatal limitation of Pn is the role of stomatal resistance to contribute to "resistance" to CO 2 uptake and stomatal limitation in spite of a decline in C i [50]. The higher transpiration (E) rates in low-density conditions may have been due to low mutual shading, which allowed rapid stomata opening. Our data showed that high plant density substantially decreased leaf thickness, stomatal density, width, length, and stomatal pores and resulted in lower C i and g s , which in turn depressed leaf photosynthetic capacity.

Conclusions
In the present study, planting densities and varieties significantly influenced lint yield by affecting leaf stomatal density, thickness, width, length, pore perimeter, leaf gas exchange, and chlorophyll fluorescence characteristics. The J-4B variety in the moderate-density condition produced a higher lint yield due to improved leaf structure, leaf gas exchange, and chlorophyll fluorescence attributes compared with low or high planting densities. Plants at high density substantially depressed leaf stomatal density, thickness, width, length, and pore perimeter, probably due to more competition for nutrients compared with low and moderate planting densities in both varieties. The offset in these attributes further disrupted ΦPSII, F v /F m , ETR, and NPQ, which in turn reduced leaf photosynthetic capacity and consequently, lint yield loss. Conclusively, J-4B and Zhongmian-16 grown under mediumand lower-density conditions may be a promising option based on improved leaf structural and functional traits in subtropical regions. Our data will substantially contribute to cotton breeding programs in subtropical environments in the future.