QTL Mapping for Yield Traits in Upland Cotton (Gossypium hirsutum L.)

Genetic Diversity and Subspecific Races of Upland Cotton (Gossypium hirsutum L.)

Background/Objectives: The classification and phylogenetic relationships of Gossypium hirsutum L. landraces, despite their proximity to southern Mexico, remain unresolved. This study aimed to clarify these relationships using SSR markers and hybridization methods, focusing on subspecies and race differentiation within G. hirsutum L. Methods: Seventy polymorphic SSR markers (out of 177 tested) were used to analyze 141 alleles and calculate genetic distances among accessions. Phylogenetic relationships were determined using MEGA software (version 11.0.13) and visualized in a phylogenetic tree. ANOVA in NCSS 12 was used for statistical analysis. Over 1000 inter-race crosses were conducted to assess boll-setting rates. Results: Distinct phylogenetic patterns were identified between G. hirsutum subspecies and races, correlating with boll-setting rates. Latifolium, richmondii, and morilli showed no significant increase in boll-setting rates in reciprocal crosses. Cultivars Omad and Bakht, as paternal parents, yielded higher boll-setting rates. Religiosum and yucatanense displayed high boll- and seed-setting rates as maternal parents but low rates as paternal parents. Additionally, phylogenetic analysis revealed a close relationship between cultivars 'Omad' and 'Bakht' with G. hirsutum race richmondii, indicating their close evolutionary relationship. Conclusions: Reciprocal differentiation characteristics of G. hirsutum subspecies and races, particularly religiosum and yucatanense, should be considered during hybridization for genetic and breeding studies. Understanding the phylogenetic relationships among G. hirsutum taxa is crucial for exploring the genetic diversity of this economically important species.

Genomic and Cytogenetic Analysis of Synthetic Polyploids between Diploid and Tetraploid Cotton (Gossypium) Species

Cotton (Gossypium spp.) is the most important natural fiber source in the world. The genetic potential of cotton can be successfully and efficiently exploited by identifying and solving the complex fundamental problems of systematics, evolution, and phylogeny, based on interspecific hybridization of cotton. This study describes the results of interspecific hybridization of G. herbaceum L. (A1-genome) and G. mustelinum Miers ex Watt (AD4-genome) species, obtaining fertile hybrids through synthetic polyploidization of otherwise sterile triploid forms with colchicine (C22H25NO6) treatment. The fertile F1C hybrids were produced from five different cross combinations: (1) G. herbaceum subsp. frutescens × G. mustelinum; (2) G. herbaceum subsp. pseudoarboreum × G. mustelinum; (3) G. herbaceum subsp. pseudoarboreum f. harga × G. mustelinum; (4) G. herbaceum subsp. africanum × G. mustelinum; (5) G. herbaceum subsp. euherbaceum (variety A-833) × G. mustelinum. Cytogenetic analysis discovered normal conjugation of bivalent chromosomes in addition to univalent, open, and closed ring-shaped quadrivalent chromosomes at the stage of metaphase I in the F1C and F2C hybrids. The setting of hybrid bolls obtained as a result of these crosses ranged from 13.8-92.2%, the fertility of seeds in hybrid bolls from 9.7-16.3%, and the pollen viability rates from 36.6-63.8%. Two transgressive plants with long fiber of 35.1-37.0 mm and one plant with extra-long fiber of 39.1-41.0 mm were identified in the F2C progeny of G. herbaceum subsp. frutescens × G. mustelinum cross. Phylogenetic analysis with 72 SSR markers that detect genomic changes showed that tetraploid hybrids derived from the G. herbaceum × G. mustelinum were closer to the species G. mustelinum. The G. herbaceum subsp. frutescens was closer to the cultivated form, and its subsp. africanum was closer to the wild form. New knowledge of the interspecific hybridization and synthetic polyploidization was developed for understanding the genetic mechanisms of the evolution of tetraploid cotton during speciation. The synthetic polyploids of cotton obtained in this study would provide beneficial genes for developing new cotton varieties of the G. hirsutum species, with high-quality cotton fiber and strong tolerance to biotic or abiotic stress. In particular, the introduction of these polyploids to conventional and molecular breeding can serve as a bridge of transferring valuable genes related to high-quality fiber and stress tolerance from different cotton species to the new cultivars.

Inheritance, QTLs, and Candidate Genes of Lint Percentage in Upland Cotton

Cotton (Gossypium spp.) is an important natural fiber plant. Lint percentage (LP) is one of the most important determinants of cotton yield and is a typical quantitative trait with high variation and heritability. Many cotton LP genetic linkages and association maps have been reported. This work summarizes the inheritance, quantitative trait loci (QTLs), and candidate genes of LP to facilitate LP genetic study and molecular breeding. More than 1439 QTLs controlling LP have been reported. Excluding replicate QTLs, 417 unique QTLs have been identified on 26 chromosomes, including 243 QTLs identified at LOD >3. More than 60 are stable, major effective QTLs that can be used in marker-assisted selection (MAS). More than 90 candidate genes for LP have been reported. These genes encode MYB, HOX, NET, and other proteins, and most are preferentially expressed during fiber initiation and elongation. A putative molecular regulatory model of LP was constructed and provides the foundation for the genetic study and molecular breeding of LP.

GBS-SNP and SSR based genetic mapping and QTL analysis for drought tolerance in upland cotton

A recombinant inbred line mapping population of intra-species upland cotton was generated from a cross between the drought-tolerant female parent (AS2) and the susceptible male parent (MCU13). A linkage map was constructed deploying 1,116 GBS-based SNPs and public domain-based 782 SSRs spanning a total genetic distance of 28,083.03 cM with an average chromosomal span length of 1,080.12 cM with inter-marker distance of 10.19 cM.A total of 19 quantitative trait loci (QTLs) were identified in nine chromosomes for field drought tolerance traits. Chromosomes 3 and 8 harbored important drought tolerant QTLs for chlorophyll stability index trait while for relative water content trait, three QTLs on chromosome 8 and one QTL each on chromosome 4, 12 were identified. One QTL on each chromosome 8, 5, and 7, and two QTLs on chromosome 15 linking to proline content were identified. For the nitrate reductase activity trait, two QTLs were identified on chromosome 3 and one on each chromosome 8, 13, and 26. To complement our QTL study, a meta-analysis was conducted along with the public domain database and resulted in a consensus map for chromosome 8. Under field drought stress, chromosome 8 harbored a drought tolerance QTL hotspot with two in-house QTLs for chlorophyll stability index (qCSI01, qCSI02) and three public domain QTLs (qLP.FDT_1, qLP.FDT_2, qCC.ST_3). Identified QTL hotspot on chromosome 8 could play a crucial role in exploring abiotic stress-associated genes/alleles for drought trait improvement.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01041-y.

Genome-wide dissection of hybridization for fiber quality- and yield-related traits in upland cotton

An evaluation of combining ability can facilitate the selection of suitable parents and superior F₁ hybrids for hybrid cotton breeding, although the molecular genetic basis of combining ability has not been fully characterized. In the present study, 282 female parents were crossed with four male parents in accordance with the North Carolina II mating scheme to generate 1128 hybrids. The parental lines were genotyped based on restriction site-associated DNA sequencing and 306 814 filtered single nucleotide polymorphisms were used for genome-wide association analysis involving the phenotypes, general combining ability (GCA) values, and specific combining ability values of eight fiber quality- and yield-related traits. The main results were: (i) all parents could be clustered into five subgroups based on population structure analyses and the GCA performance of the female parents had significant differences between subgroups; (ii) 20 accessions with a top 5% GCA value for more than one trait were identified as elite parents for hybrid cotton breeding; (iii) 120 significant single nucleotide polymorphisms, clustered into 66 quantitative trait loci, such as the previously reported Gh_A07G1769 and GhHOX3 genes, were found to be significantly associated with GCA; and (iv) identified quantitative trait loci for GCA had a cumulative effect on GCA of the accessions. Overall, our results suggest that pyramiding the favorable loci for GCA may improve the efficiency of hybrid cotton breeding.

Detection of favorable alleles for yield and yield components by association mapping in upland cotton

Association mapping based on linkage disequilibrium provides a promising tool for dissecting the genetic basis underlying complex traits. To reveal the genetic variations of yield and yield components traits in upland cotton, 403 upland cotton accessions were collected and analyzed by 560 genome-wide simple sequence repeats (SSRs). A diverse panel consisting of 403 upland cotton accessions was grown in six different environments, and the yield and yield component traits were measured, and 560 SSR markers covering the whole genome were mapped. Association studies were performed to uncover the genotypic and phenotypic variations using a mixed linear model. Favorable alleles and typical accessions for yield traits were identified. A total of 201 markers were polymorphic, revealing 394 alleles. The average gene diversity and polymorphism information content were 0.556 and 0.483, respectively. Based on a population structure analysis, 403 accessions were divided into two subgroups. A mixed linear model analysis of the association mapping detected 43 marker loci according to the best linear unbiased prediction and in at least three of the six environments(- lgP > 1.30, P < 0.05). Among the 43 associated markers, five were associated with more than two traits simultaneously and nine were coincident with those identified previously. Based on phenotypic effects, favorable alleles and typical accessions that contained the elite allele loci related to yield traits were identified and are widely used in practical breeding. This study detected favorable quantitative trait loci's alleles and typical accessions for yield traits, these are excellent genetic resources for future high-yield breeding by marker-assisted selection in upland cotton in China.

Inclusive Composite Interval Mapping of QTL by Environment Interactions in Biparental Populations

Identification of environment-specific QTL and stable QTL having consistent genetic effects across a wide range of environments is of great importance in plant breeding. Inclusive Composite Interval Mapping (ICIM) has been proposed for additive, dominant and epistatic QTL mapping in biparental populations for single environment. In this study, ICIM was extended to QTL by environment interaction (QEI) mapping for multi-environmental trials, where the QTL average effect and QEI effects could be properly estimated. Stepwise regression was firstly applied in each environment to identify the most significant marker variables which were then used to adjust the phenotypic values. One-dimensional scanning was then conducted on the adjusted phenotypic values across the environments in order to detect QTL with either average effect or QEI effects, or both average effect and QEI effects. In this way, the genetic background could be well controlled while the conventional interval mapping was applied. An empirical method to determine the threshold of logarithm of odds was developed, and the efficiency of the ICIM QEI mapping was demonstrated in simulated populations under different genetic models. One actual recombinant inbred line population was used to compare mapping results between QEI mapping and single-environment analysis.

Cotton QTLdb: a cotton QTL database for QTL analysis, visualization, and comparison between Gossypium hirsutum and G. hirsutum × G. barbadense populations

KEY MESSAGE

A specialized database currently containing more than 2200 QTL is established, which allows graphic presentation, visualization and submission of QTL. In cotton quantitative trait loci (QTL), studies are focused on intraspecific Gossypium hirsutum and interspecific G. hirsutum × G. barbadense populations. These two populations are commercially important for the textile industry and are evaluated for fiber quality, yield, seed quality, resistance, physiological, and morphological trait QTL. With meta-analysis data based on the vast amount of QTL studies in cotton it will be beneficial to organize the data into a functional database for the cotton community. Here we provide a tool for cotton researchers to visualize previously identified QTL and submit their own QTL to the Cotton QTLdb database. The database provides the user with the option of selecting various QTL trait types from either the G. hirsutum or G. hirsutum × G. barbadense populations. Based on the user's QTL trait selection, graphical representations of chromosomes of the population selected are displayed in publication ready images. The database also provides users with trait information on QTL, LOD scores, and explained phenotypic variances for all QTL selected. The CottonQTLdb database provides cotton geneticist and breeders with statistical data on cotton QTL previously identified and provides a visualization tool to view QTL positions on chromosomes. Currently the database (Release 1) contains 2274 QTLs, and succeeding QTL studies will be updated regularly by the curators and members of the cotton community that contribute their data to keep the database current. The database is accessible from http://www.cottonqtldb.org.

Identification of associated SSR markers for yield component and fiber quality traits based on frame map and Upland cotton collections

Detecting QTLs (quantitative trait loci) that enhance cotton yield and fiber quality traits and accelerate breeding has been the focus of many cotton breeders. In the present study, 359 SSR (simple sequence repeat) markers were used for the association mapping of 241 Upland cotton collections. A total of 333 markers, representing 733 polymorphic loci, were detected. The average linkage disequilibrium (LD) decay distances were 8.58 cM (r2 > 0.1) and 5.76 cM (r2 > 0.2). 241 collections were arranged into two subgroups using STRUCTURE software. Mixed linear modeling (MLM) methods (with population structure (Q) and relative kinship matrix (K)) were applied to analyze four phenotypic datasets obtained from four environments (two different locations and two years). Forty-six markers associated with the number of bolls per plant (NB), boll weight (BW), lint percentage (LP), fiber length (FL), fiber strength (FS) and fiber micornaire value (FM) were repeatedly detected in at least two environments. Of 46 associated markers, 32 were identified as new association markers, and 14 had been previously reported in the literature. Nine association markers were near QTLs (at a distance of less than 1-2 LD decay on the reference map) that had been previously described. These results provide new useful markers for marker-assisted selection in breeding programs and new insights for understanding the genetic basis of Upland cotton yields and fiber quality traits at the whole-genome level.

A comparative meta-analysis of QTL between intraspecific Gossypium hirsutum and interspecific G. hirsutum × G. barbadense populations

KEY MESSAGE

Based on 1075 and 1059 QTL from intraspecific Upland and interspecific Upland × Pima populations, respectively, the identification of QTL clusters and hotspots provides a useful resource for cotton breeding. Mapping of quantitative trait loci (QTL) is a pre-requisite of marker-assisted selection for crop yield and quality. Recent meta-analysis of QTL in tetraploid cotton (Gossypium spp.) has identified regions of the genome with high concentrations of QTL for various traits called clusters and specific trait QTL called hotspots or meta-QTL (mQTL). However, the meta-analysis included all population types of Gossypium mixing both intraspecific G. hirsutum and interspecific G. hirsutum × G. barbadense populations. This study used 1,075 QTL from 58 publications on intraspecific G. hirsutum and 1,059 QTL from 30 publications on G. hirsutum × G. barbadense populations to perform a comprehensive comparative analysis of QTL clusters and hotspots between the two populations for yield, fiber and seed quality, and biotic and abiotic stress tolerance. QTL hotspots were further analyzed for mQTL within the hotspots using Biomercator V3 software. The ratio of QTL between the two population types was proportional yet differences in hotspot type and placement were observed between the two population types. However, on some chromosomes QTL clusters and hotspots were similar between the two populations. This shows that there are some universal QTL regions in the cultivated tetraploid cotton which remain consistent and some regions which differ between population types. This study for the first time elucidates the similarities and differences in QTL clusters and hotspots between intraspecific and interspecific populations, providing an important resource to cotton breeding programs in marker-assisted selection .

Favorable QTL alleles for yield and its components identified by association mapping in Chinese Upland cotton cultivars

Linkage disequilibrium based association mapping is a powerful tool for dissecting the genetic basis underlying complex traits. In this study, an association mapping panel consisting of 356 representative Upland cotton cultivars was constructed, evaluated in three environments and genotyped using 381 SSRs to detect molecular markers associated with lint yield and its components. The results showed that abundant phenotypic and moderate genetic diversities existed within this germplasm panel. The population could be divided into two subpopulations, and weak relatedness was detected between pair-wise accessions. LD decayed to the background (r² = 0.1182, P≤0.01), r² = 0.1 and r² = 0.2 level within 12–13 cM, 17–18 cM and 3–4 cM, respectively, providing the potential for association mapping of agronomically important traits in Chinese Upland cotton. A total of 55 marker-trait associations were detected between 26 SSRs and seven lint yield traits, based on a mixed linear model (MLM) and Bonferroni correction (P≤0.05/145, −log₁₀P≥3.46). Of which 41 could be detected in more than one environment and 17 markers were simultaneously associated with two or more traits. Many associations were consistent with QTLs identified by linkage mapping in previous reports. Phenotypic values of alleles of each loci in 41 stably detected associations were compared, and 23 favorable alleles were identified. Population frequency of each favorable allele in historically released cultivar groups was also evaluated. The QTLs detected in this study will be helpful in further understanding the genetic basis of lint yield and its components, and the favorable alleles may facilitate future high-yield breeding by genomic selection in Upland cotton.

QTL analysis for early-maturing traits in cotton using two upland cotton (Gossypium hirsutum L.) crosses

Making use of the markers linked closely to QTL for early-maturing traits for MAS (Marker-assisted selection) is an effective method for the simultaneous improvement of early maturity and other properties in cotton. In this study, two F2 populations and their F2:3 families were generated from the two upland cotton (Gossypium hirsutum L.) crosses, Baimian2 × TM-1 and Baimian2 × CIR12. QTL for early-maturing traits were analyzed using F2:3 families. A total of 54 QTL (31 suggestive and 23 significant) were detected. Fourteen significant QTL had the LOD scores not only > 3 but also exceeding permutation threshold. At least four common QTL, qBP-17 for bud period (BP), qGP-17a/qGP-17b (qGP-17) for growth period (GP), qYPBF-17a/qYPBF-17b (qYPBF-17) for yield percentage before frost (YPBF) and qHFFBN-17 for height of first fruiting branch node (HFFBN), were found in both populations. These common QTL should be reliable and could be used for MAS to facilitate early maturity. The common QTL, qBP-17, had a LOD score not only > 3 but also exceeding permutation threshold, explaining 12.6% of the phenotypic variation. This QTL should be considered preferentially in MAS. Early-maturing traits of cotton are primarily controlled by dominant and over-dominant effects.

Variations and transmission of QTL alleles for yield and fiber qualities in upland cotton cultivars developed in China

Cotton is the world’s leading cash crop, and genetic improvement of fiber yield and quality is the primary objective of cotton breeding program. In this study, we used various approaches to identify QTLs related to fiber yield and quality. Firstly, we constructed a four-way cross (4WC) mapping population with four base core cultivars, Stoneville 2B, Foster 6, Deltapine 15 and Zhongmiansuo No.7 (CRI 7), as parents in Chinese cotton breeding history and identified 83 QTLs for 11 agronomic and fiber quality traits. Secondly, association mapping of agronomical and fiber quality traits was based on 121 simple sequence repeat (SSR) markers using a general linear model (GLM). For this, 81 Gossypium hirsutum L. accessions including the four core parents and their derived cultivars were grown in seven diverse environments. Using these approaches, we successfully identified 180 QTLs significantly associated with agronomic and fiber quality traits. Among them were 66 QTLs that were identified via linkage disequilibrium (LD) and 4WC family-based linkage (FBL) mapping and by previously published family-based linkage (FBL) mapping in modern Chinese cotton cultivars. Twenty eight and 44 consistent QTLs were identified by 4WC and LD mapping, and by FBL and LD mapping methods, respectively. Furthermore, transmission and variation of QTL-alleles mapped by LD association in the three breeding periods revealed that some could be detected in almost all Chinese cotton cultivars, suggesting their stable transmission and some identified only in the four base cultivars and not in the modern cultivars, suggesting they were missed in conventional breeding. These results will be useful to conduct genomics-assisted breeding effectively using these existing and novel QTL alleles to improve yield and fiber qualities in cotton.