Search

Abstract Global water scarcity and extreme weather intensify drought stress, significantly reducing cotton yield and quality worldwide. Drought treatments are conducted using a population of chromosome segment substitution lines generated from E22 (G. hirsutum) and 3–79 (G. barbadense) as parental lines either show superior yields or fiber quality under both control and drought conditions. Fourteen datasets, covering 4 yields and 4 quality traits, are compiled and assessed for drought resistance using the drought resistance coefficient (DRC) and membership function value of drought resistance (MFVD). Genome‐wide association studies, linkage analysis, and bulked segregant analysis are combined to analyze the DR‐related QTL. A total of 121 significant QTL are identified by DRC and MFVD of the 8 traits. CRISPR/Cas9 and virus‐induced gene silencing techniques verified DRR1 and DRT1 as pivotal genes in regulating drought resistant of cotton, with hap3‐79 exhibiting greater drought resistance than hapE22 concerning DRR1 and DRT1. Moreover, 14 markers with superior yield and fiber quality are selected for drought treatment. This study offers valuable insights into yield and fiber quality variations between G. hirsutum and G. barbadense amid drought, providing crucial theoretical and technological backing for developing cotton varieties resilient to drought, with high yield and superior fiber quality.

Introduction: Sea-island cotton (Gossypium barbadense, Gb) is one of the major sources of high-grade natural fiber. Besides the common annual Gb cotton, perennial Gb cotton is also cultivated, but studies on perennial Gb cotton are rare. Objectives: We aimed to make a systematic analysis of perennial sea-island cotton and lay a foundation for its utilization in breeding, and try to identify the representative structural variations (SVs) in sea-island cotton, and to reveal the population differentiation and adaptive improvement of sea-island cotton. Methods: Through genome assembly of one perennial Gb cotton accession (named Gb_M210936) and comparative genome analysis, variations during Gb cotton domestication were identified by comparing Gb_M210936 with annual Gb accession 3–79 and with wild allotetraploid cotton G. darwinii. Six perennial Gb accessions combining with the resequenced 1,129 cotton accessions were used to conduct population and genetic analysis. Large haplotype blocks (haploblocks), generated from interspecific introgressions and intraspecific inversions, were identified and were used to analyze their effects on population differentiation and agronomic traits of sea-island cotton. Results: One reference genome of perennial sea-island cotton was assembled. Representative SVs in sea-island cotton were identified, and 31 SVs were found to be associated with agronomic traits. Perennial Gb cotton had a closer kinship with the wild-to-landrace continuum Gb cotton from south America where Gb cotton is originally domesticated. Haploblocks were associated with agronomic traits improvement of sea-island cotton, promoted sea-island cotton differentiation into three subgroups, were suffered from breeding selection, and may drive Gb cotton to be adapted to central Asian. Conclusion: Our study made up the lack of perennial Gb cotton genome, and clarified that exotic introgressions improved the traits of sea-island cotton, promoted the population differentiation, and drove sea-island cotton adaptive to central Asia, which will provide new insights for the genetic breeding improvement of sea-island cottons.

Abstract Background Cotton fiber yield is a complex trait, which can be influenced by multiple agronomic traits. Unravelling the genetic basis of cotton fiber yield-related traits contributes to genetic improvement of cotton. Results In this study, 503 upland cotton varieties covering the four breeding stages (BS1–BS4, 1911–2011) in China were used for association mapping and domestication analysis. One hundred and forty SSR markers significantly associated with ten fiber yield-related traits were identified, among which, 29 markers showed an increasing trend contribution to cotton yield-related traits from BS1 to BS4, and 26 markers showed decreased trend effect. Four favorable alleles of 9 major loci (R 2 ≥ 3) were strongly selected during the breeding stages, and the candidate genes of the four strongly selected alleles were predicated according to the gene function annotation and tissue expression data. Conclusions The study not only uncovers the genetic basis of 10 cotton yield-related traits but also provides genetic evidence for cotton improvement during the cotton breeding process in China.

Hybrid utilization has proficiently increased crop production worldwide. The cytoplasmic male sterility (CMS) system has emerged as an efficient tool for commercial hybrid cotton seed production. The restorer line with dominant Rf2 gene can restore the fertility of the CMS-D8 sterile line. However, the molecular mechanism of fertility restoration remains unclear in CMS-D8 cotton that limits wider utilization of three-line hybrid breeding. In our study, the Pacific Biosciences (PacBio) Iso-Seq technology was applied to understand fertility restoration mechanism of CMS-D8 cotton. In total, 228,106 full-length non-chimeric transcriptome sequences were obtained from anthers of developing flowering buds. The analysis results identified 3,174 novel isoforms, 2,597 novel gene loci, 652 long non-coding RNAs predicted from novel isoforms, 7,234 alternative splicing events, 114 fusion transcripts, and 1,667 genes with alternative polyadenylation. Specially, two novel genes associated with restoration function, Ghir_D05.742.1 and m64033_190821_201011/21103726/ccs were identified and showed significant higher levels of expression in restorer line than sterile and maintainer lines. Our comparative full-length transcriptome analysis provides new insights into the molecular function of Rf2 fertility restorer gene. The results of this study offer a platform for fertility restoration candidate gene discovery in CMS-D8 cotton.

Abstract Background Gossypium hirsutum L. (cotton) is one of the most economically important crops in the world due to its significant source of fiber, feed, foodstuff, oil and biofuel products. However, the utilization of cottonseed was limited due to the presence of small and darkly pigmented glands that contain large amounts of gossypol, which is toxic to human beings and non-ruminant animals. To date, some progress has been made in the pigment gland formation, but the underlying molecular mechanism of its formation was still unclear. Results In this study, we identified an AP2/ERF transcription factor named GhERF105 (GH_A12G2166), which was involved in the regulation of gland pigmentation by the comparative transcriptome analysis of the leaf of glanded and glandless plants. It encoded an ERF protein containing a converved AP2 domain which was localized in the nucleus with transcriptional activity, and showed the high expression in glanded cotton accessions that contained much gossypol. Virus-induced gene silencing (VIGS) against GhERF105 caused the dramatic reduction in the number of glands and significantly lowered levels of gossypol in cotton leaves. GhERF105 showed the patterns of spatiotemporal and inducible expression in the glanded plants. Conclusions These results suggest that GhERF105 contributes to the pigment gland formation and gossypol biosynthesis in partial organs of glanded plant. It also provides a potential molecular basis to generate ‘glandless-seed’ and ‘glanded-plant’ cotton cultivar.

As one of the most widely cultivated cotton species in China, upland cotton has moderate fiber quality and wide applicability, but its genetic basis is relatively narrow. To expand genetic diversity and improve fiber quality, in this study an introgression population (BC5S5) containing 107 lines was constructed by using G. hirsutum acc. 4105 as the recurrent parent and G. tomentosum as the donor parent. Using the specific-locus amplified fragment sequencing (SLAF-seq) strategy, 3157 high-throughput single nucleotide polymorphism (SNP) markers were obtained. Linkage analysis showed that a total of ninety-one QTLs related to fiber quality traits were detected in three environments, and the phenotypic variance explained (PVE) rates were 4.53–20.92%. Forty-six QTL (50.55%) synergistic genes were derived from G. tomentosum. Among them, qFS-A02-1 and qSCI-A02-1 were stably detected with a PVE of 9.8–16.71% and 14.78–20.92%, respectively. Within the candidate interval, Ghir_A02G012730, Ghir_A02G012790 and Ghir_A02G012830 were found to be possibly involved in cellulose and cell wall biosynthesis, with a relatively high expression during fiber development, 20 DPA and 25 DPA, which suggested that these three genes may be involved in the regulation of fiber strength traits, but their functions need further validation to determine the regulatory mechanism. Our research lays the foundation of fiber quality related to basic genetic research and breeding in cotton.

Abstract Background Meta-analysis of quantitative trait locus (QTL) is a computational technique to identify consensus QTL and refine QTL positions on the consensus map from multiple mapping studies. The combination of meta-QTL intervals, significant SNPs and transcriptome analysis has been widely used to identify candidate genes in various plants. Results In our study, 884 QTLs associated with cotton fiber quality traits from 12 studies were used for meta-QTL analysis based on reference genome TM-1, as a result, 74 meta-QTLs were identified, including 19 meta-QTLs for fiber length; 18 meta-QTLs for fiber strength; 11 meta-QTLs for fiber uniformity; 11 meta-QTLs for fiber elongation; and 15 meta-QTLs for micronaire. Combined with 8 589 significant single nucleotide polymorphisms associated with fiber quality traits collected from 15 studies, 297 candidate genes were identified in the meta-QTL intervals, 20 of which showed high expression levels specifically in the developing fibers. According to the function annotations, some of the 20 key candidate genes are associated with the fiber development. Conclusions This study provides not only stable QTLs used for marker-assisted selection, but also candidate genes to uncover the molecular mechanisms for cotton fiber development.

Global warming is limiting availability of water resources in arid and semi-arid regions, and so understanding water use efficiency (WUE) is increasingly important for agricultural production in those areas. As China is the largest cotton producing area, the problem of balancing WUE and efficient cotton production is a major issue. In this study, we used a natural population of 517 Upland cotton accessions to conduct a water-controlled trial in south and north of Xinjiang over two years. A total of 18 traits including agronomic traits, fiber yield indices and fiber quality indices, were investigated for broad-sense heritability and coefficient of variation. Appropriate water limitation was found to promote the establishment of favorable agronomic traits in cotton, associated with an increased cotton yield of 8.46% in Xinjiang, at the expense of a certain degree of fiber quality, such as decreased fiber length and an over-higher micronaire value. We detected 33 QTL related to response to water limitation using a drought resistance coefficient (DRC), and 6 QTL were found using a comprehensive indicator of CIDT (comprehensive index of drought tolerance) at the genetic level by integrating resequencing data. Two novel QTL-hotspots including six differentially expressed genes (DEGs) were further identified related to the drought response of cotton. These findings not only suggested a new approach to irrigation of cotton fields in Xinjiang, but also provided abundant genetic evidence for genetic breeders to study drought improvement of crops.

Structural variation is a major type of genetic variation that can potentially induce powerful genetic effects. In this study, we examined the Inv(A07)p1.09p2.23 genetic inversion in brown fibre cotton at the individual and population genetics levels. A dark-brown fibre mutant that resulted from a distant hybridization between Gossypium barbadense and G. hirsutum, and a natural population including 30 dark-brown, 70 light-brown and 21 white fibre cotton accessions were collected to perform a functional study of this micro-inversion. The results showed that Inv(A07)p1.09p2.23 can be detected by high-throughput resequencing method, and induce micro-deletion, gene disruption (Ghir_A07G000980) and abnormal gene expression in the breakpoint regions. Inv(A07)p1.09p2.23 existed in only dark-brown fibre cotton, had undergone negative selection in elite brown fibre cultivars, and was significantly associated with fibre colour and nine fibre traits. In the Inv(A07)p1.09p2.23 region, nucleotide diversity was lower, recombination was absent, and linkage disequilibrium was higher. Overall, this inversion event in dark-brown fibre cotton produced significant genetic effects, and this study will guide us to better understand the genetic effects of inversion events in dark-brown fibre cotton.

Abstract Background The low genetic diversity of Upland cotton limits the potential for genetic improvement. Making full use of the genetic resources of Sea-island cotton will facilitate genetic improvement of widely cultivated Upland cotton varieties. The chromosome segments substitution lines (CSSLs) provide an ideal strategy for mapping quantitative trait loci (QTL) in interspecific hybridization. Results In this study, a CSSL population was developed by PCR-based markers assisted selection (MAS), derived from the crossing and backcrossing of Gossypium hirsutum (Gh) and G. barbadense (Gb), firstly. Then, by whole genome re-sequencing, 11,653,661 high-quality single nucleotide polymorphisms (SNPs) were identified which ultimately constructed 1211 recombination chromosome introgression segments from Gb. The sequencing-based physical map provided more accurate introgressions than the PCR-based markers. By exploiting CSSLs with mutant morphological traits, the genes responding for leaf shape and fuzz-less mutation in the Gb were identified. Based on a high-resolution recombination bin map to uncover genetic loci determining the phenotypic variance between Gh and Gb, 64 QTLs were identified for 14 agronomic traits with an interval length of 158 kb to 27 Mb. Surprisingly, multiple alleles of Gb showed extremely high value in enhancing cottonseed oil content (SOC). Conclusions This study provides guidance for studying interspecific inheritance, especially breeding researchers, for future studies using the traditional PCR-based molecular markers and high-throughput re-sequencing technology in the study of CSSLs. Available resources include candidate position for controlling cotton quality and quantitative traits, and excellent breeding materials. Collectively, our results provide insights into the genetic effects of Gb alleles on the Gh, and provide guidance for the utilization of Gb alleles in interspecific breeding.

With increasing demand for high-quality cotton, it is desirable to identify genes involved in fiber development for molecular improvement of cotton. In this study, 780 differentially expressed genes (DEGs) were identified in developing fibers at 10 days post-anthesis (DPA) in Gossypium hirsutum acc. DH962 and G. hirsutum cv. Jimian 5 using RNA-seq. Of 15 stable QTL for fiber quality identified in the same two parents in previous studies, 4, 3, 6, 1, and 1 QTL were associated with fiber length (FL), fiber strength (FS), micronaire (MIC), fiber elongation (FE) and fiber length uniformity ratio (FU), respectively. Integration of DEGs and QTL allowed the identification of 31 genes in 9 QTL regions, of which 25 were highly expressed in fibers based on the transcriptome datasets and 9 were preferentially expressed in different stages of fiber development. Gh_A01G0453 (GhDTX19), Gh_D07G1581 and Gh_D04G0942 were expressed specifically in 5 and 10 DPA fibers, with Gh_D04G0942 showing low expression in other tissues except pistil. Gh_D07G1799 (GhGAUT9), Gh_D11G0326 (GhVPS29), Gh_D11G0333 (GhTCP14), and Gh_D11G0334 (GhNRP2) were preferentially expressed in 5 or 10 DPA fibers; Gh_A01G0397 (GhABCG10) and Gh_D07G1744 were expressed specifically in 20 and 25 DPA fibers. These results suggest candidate genes for molecular improvement of cotton fiber quality. Keywords: Upland cotton, Fiber quality, RNA-seq, QTL, Candidate genes

Plant architecture traits influence crop yield. An understanding of the genetic basis of cotton plant architecture traits is beneficial for identifying favorable alleles and functional genes and breeding elite cultivars. We collected 121 cotton accessions including 100 brown-fiber and 21 white-fiber accessions, genotyped them by whole-genome resequencing, and phenotyped them in multiple environments. This genome-wide association study (GWAS) identified 11 quantitative trait loci (QTL) for two plant architecture traits: plant height and fruit spur branch number. Negative-effect alleles were enriched in the elite cultivars. Based on these QTL, gene annotation information, and published QTL, candidate genes and natural genetic variations in four QTL were identified. Ghir_D02G017510 and Ghir_D02G017600 were identified as candidate genes for qD02-FSBN-1, and a premature start codon gain variation was found in Ghir_D02G017510. Ghir_A12G026570, the candidate gene of qA12-FSBN-2, belongs to the pectin lyase-like superfamily, and a significantly associated SNP, A12_105366045 (T/C), in this gene represents an amino acid change. The QTL, candidate genes, and associated natural variations in this study are expected to lay a foundation for studying functional genes and developing breeding programs for desirable architecture in brown-fiber cotton. Keywords: Brown-fiber cotton, Plant architecture traits, GWAS, QTL, Genetic variation

Environmental stresses severely impair cotton production worldwide. To identify the genetic basis of, and molecular markers associated with, environmental stresses such as salt, cold and Verticillium wilt, association mapping of salt-, cold-, and disease-tolerance traits was performed in a population of 503 upland cotton accessions using 179 polymorphic SSR markers and 11,975 array-derived SNP markers. Salt and cold tolerance was evaluated via the relative germination rate (RGR) of the seeds under seven and four environments, respectively. The disease index of Verticillium wilt was investigated for two years in the field in Xinjiang. These three traits showed large variation across environments. A genome-wide association study revealed that 31, 19, and 15 SSR markers were associated with RGR-Salt (the relative germination rates of seeds under salt stress), RGR-Cold (the relative germination rates of seeds under cold stress), and DIV (the disease index of Verticillium wilt), respectively. Six SNPs in seven environments and two SNPs in BLUP (best linear unbiased prediction) results were associated with RGR-Salt, and the phenotypic variance explained ranged from 3.96 to 5.00%. Two SNPs (i02237Gh, i02243Gh) on D01 were concluded to be stable genetic loci associated with RGR-Salt. A total of 223 genes were found in a candidate gene interval (D01, 37771–1942912). Of these four genes, GhPIP3A, GhSAG29, GhTZF4, and GhTZF4a, showed expression changes in sensitive and tolerant genotype accessions under salt stress, and were assigned as candidate genes associated with RGR-Salt. This study revealed the genetic basis of stress resistance in upland cotton and will facilitate stress-resistance breeding in cotton. Keywords: Upland cotton, Stress, Relative germination rate, SSR, SNP, Candidate gene

Abstract Background Ongoing molecular processes in a cell could target microsatellites, a kind of repetitive DNA, owing to length variations and motif imperfection. Mutational mechanisms underlying such kind of genetic variations have been extensively investigated in diverse organisms. However, obscure impact of ploidization, an evolutionary process of genome content duplication prevails mostly in plants, on non-coding DNA is poorly understood. Results Genome sequences of diversely originated plant species were examined for genome-wide motif imperfection pattern, and various analytical tools were employed to canvass characteristic relationships among repeat density, imperfection and length of microsatellites. Moreover, comparative genomics approach aided in exploration of microsatellites conservation footprints in Gossypium evolution. Based on our results, motif imperfection in repeat length was found intricately related to genomic abundance of imperfect microsatellites among 13 genomes. Microsatellite decay estimation depicted slower decay of long motif repeats which led to predominant abundance of 5-nt repeat motif in Gossypium species. Short motif repeats exhibited rapid decay through the evolution of Gossypium lineage ensuing drastic decrease of 2-nt repeats, of which, “AT” motif type dilapidated in cultivated tetraploids of cotton. Conclusion The outcome could be a directive to explore comparative evolutionary footprints of simple non-coding genetic elements i.e., repeat elements, through the evolution of genus-specific characteristics in cotton genomes.

Recombination is crucial for genetic evolution, which not only provides new allele combinations but also influences the biological evolution and efficacy of natural selection. However, recombination variation is not well understood outside of the complex species' genomes, and it is particularly unclear in Gossypium. Cotton is the most important natural fibre crop and the second largest oil-seed crop. Here, we found that the genetic and physical maps distances did not have a simple linear relationship. Recombination rates were unevenly distributed throughout the cotton genome, which showed marked changes along the chromosome lengths and recombination was completely suppressed in the centromeric regions. Recombination rates significantly varied between A-subgenome (At) (range = 1.60 to 3.26 centimorgan/megabase [cM/Mb]) and D-subgenome (Dt) (range = 2.17 to 4.97 cM/Mb), which explained why the genetic maps of At and Dt are similar but the physical map of Dt is only half that of At. The translocation regions between A02 and A03 and between A04 and A05, and the inversion regions on A10, D10, A07 and D07 indicated relatively high recombination rates in the distal regions of the chromosomes. Recombination rates were positively correlated with the densities of genes, markers and the distance from the centromere, and negatively correlated with transposable elements (TEs). The gene ontology (GO) categories showed that genes in high recombination regions may tend to response to environmental stimuli, and genes in low recombination regions are related to mitosis and meiosis, which suggested that they may provide the primary driving force in adaptive evolution and assure the stability of basic cell cycle in a rapidly changing environment. Global knowledge of recombination rates will facilitate genetics and breeding in cotton.

Upland cotton (Gossypium hirsutum L.), an important source of natural fiber, can tolerate relatively high salinity and drought stresses. In the present study, a plasma membrane Na+/H+ antiporter gene, GhSOS1, was cloned from a salt-tolerant genotype of G. hirsutum, Zhong 9807. The expression level of GhSOS1 in cotton roots was significantly upregulated in the presence of high concentrations of NaCl (200 mM), while its transcript abundance was increased when exposed to low temperature and drought stresses. Localization analysis using onion epidermal cells showed that the GhSOS1 protein was localized to the plasma membrane. The overexpression of GhSOS1 in Arabidopsis enhanced tolerance to salt stress, as indicated by a lower MDA content and decreased Na+/K+ ratio in transgenic plants. Moreover, the transcript levels of stress-related genes were significantly higher in GhSOS1 overexpression lines than in wild-type plants under salt treatment. Hence, GhSOS1 may be a potential target gene for enhancing salt tolerance in transgenic plants.

Cotton is a pioneer of saline land crop, while salt stress still causes its growth inhibition and fiber production decrease. Phenotype identification showed better salt tolerance of a wild diploid cotton species Gossypium klotzschianum. To elucidate the salt-tolerant mechanisms in G. klotzschianum, we firstly detected the changes in hormones, H2O2 and glutathione (GSSH and GSH), then investigated the gene expression pattern of roots and leaves treated with 300 mM NaCl for 0, 3, 12, 48 h, and each time control by RNA-seq on the Illumina-Solexa platform. Physiological determination proved that the significant increase in hormone ABA at 48 h, while that in H2O2 was at 12 h, likewise, the GSH content decrease at 48 h and the GSSH content increase at 48 h, under salt stress. In total, 37,278 unigenes were identified from the transcriptome data, 8,312 and 6,732 differentially expressed genes (DEGs) were discovered to be involved in salt stress tolerance in roots and leaves, respectively. Gene function annotation and expression analysis elucidated hormone biosynthesis and signal transduction, reactive oxygen species (ROS), and salt overly sensitive (SOS) signal transduction related genes revealed the important roles of them in signal transmission, oxidation balance and ion homeostasis in response to salinity stress. This is a report which focuses on primary response to highly salty stress (upto 300 mM NaCl) in cotton using a wild diploid Gossypium species, broadening our understanding of the salt tolerance mechanism in cotton and laying a solid foundation of salt resistant for the genetic improvement of upland cotton with the resistance to salt stress.

Two immortalized backcross populations (DHBCF1s and JMBCF1s) were developed using a recombinant inbred line (RIL) population crossed with the two parents DH962 and Jimian5 (as the males), respectively. The fiber quality and yield component traits of the two backcross populations were phenotyped at four environments (two locations, two years). One hundred seventy-eight quantitative trait loci (QTL) were detected including 76 for fiber qualities and 102 for yield components, explaining 4.08-17.79% of the phenotypic variation (PV). Among the 178 QTL, 22 stable QTL were detected in more than one environment or population. A stable QTL, qFL-c10-1, was detected in the previous F2 population, a RIL population in 3 environments and the current two BCF1 populations in this study, explaining 5.79-37.09% of the PV. Additionally, 117 and 110 main-effect QTL (M-QTL) and 47 and 191 digenic epistatic QTL (E-QTL) were detected in the DHBCF1s and JMBCF1s populations, respectively. The effect of digenic epistasis played a more important role on lint percentage, fiber length and fiber strength. These results obtained in the present study provided more resources to obtain stable QTL, confirming the authenticity and reliability of the QTL for molecular marker-assisted selection breeding and QTL cloning.

The tetraploid species Gossypium hirsutum is cultivated widely throughout the world with high yield and moderate fiber quality, but its genetic basis is narrow. A set of 107 introgression lines (ILs) was developed with an interspecific cross using G. hirsutum acc. 4105 as the recurrent parent and G. tomentosum as the donor parent. A specific locus amplified fragment sequencing (SLAF-seq) strategy was used to obtain high-throughput single nucleotide polymorphism (SNP) markers. In total, 3157 high-quality SNP markers were obtained and further used for identification of quantitative trait loci (QTLs) for fiber quality and yield traits evaluated in multiple environments. In total, 74 QTLs were detected that were associated with five fiber quality traits (30 QTLs) and eight yield traits (44 QTLs), with 2.02–30.15% of the phenotypic variance explained (PVE), and 69 markers were found to be associated with these thirteen traits. Eleven chromosomes in the A sub-genome (At) harbored 47 QTLs, and nine chromosomes in the D sub-genome (Dt) harbored 27 QTLs. More than half (44 QTLs = 59.45%) showed positive additive effects for fiber and yield traits. Five QTL clusters were identified, with three in the At, comprised of thirteen QTLs, and two in the Dt comprised of seven QTLs. The ILs developed in this study and the identified QTLs will facilitate further molecular breeding for improvement of Upland cotton in terms of higher yield with enhanced fiber quality.

A population of 178 recombinant inbred lines (RILs) was developed using a single seed descendant from a cross between G. hirsutum. acc DH962 and G. hirsutum. cv Jimian5, was used to construct a genetic map and to map QTL for fiber and yield traits. A total of 644 polymorphic loci were used to construct a final genetic map, containing 616 loci and spanning 2016.44 cM, with an average of 3.27 cM between adjacent markers. Statistical analysis revealed that segregation distortion in the intraspecific population was more serious than that in the interspecific population. The RIL population and the two parents were phenotyped under 8 environments (two locations, six years), revealing a total of 134 QTL, including 64 for fiber qualities and 70 for yield components, independently detected in seven environments, explaining 4.40-15.28% of phenotypic variation (PV). Among the 134 QTL, 9 common QTL were detected in more than one environment, and 22 QTL and 19 new QTL were detected in combined analysis (E9). A total of 26 QTL hotspot regions were observed on 13 chromosomes and 2 larger linkage groups, and some QTL clusters related to fiber qualities or yield components were also observed. The results obtained in the present study suggested that to map accurate QTL in crops with larger plant types, such as cotton, phenotyping under multiple environments is necessary to effectively apply the obtained results in molecular marker-assisted selection breeding and QTL cloning.

Transcription factors (TFs) play an important role in the regulation of plant growth and development. The study of the structure and function of TFs represents a research frontier in plant molecular biology. The findings of these studies will provide significant information regarding genetic improvement traits in crops. Currently, a large number of TFs have been cloned, and their function has been verified. However, relatively few studies that genetically map TFs in cotton are available. To genetically map TFs in cotton in this study, specific primers were designed for TF genes that were published in the Plant Transcription Factor Database. A total of 977 TF primers were obtained, and 31 TF polymorphic loci were mapped on 15 cotton chromosomes. These polymorphic loci were clearly preferentially distributed on chromosomes 5, 11, 19 and 20; and TFs from the same family mapped to homologous cotton chromosomes. In-silico mapping verified that many mapped TFs were mapped on their corresponding chromosomes or their homologous chromosomes' corresponding chromosomes in the diploid genomes. QTL mapping for fiber quality revealed that TF-Ghi005602-2 mapped on Chr19 was associated with fiber length. Eighty-five TF genes were selected for RT-PCR analysis, and 4 TFs were selected for qRT-PCR analysis, revealing unique expression patterns across different stages of fiber development between the mapping parents. Our data offer an overview of the chromosomal distribution of TFs in cotton, and the comparative expression analysis between Gossypium hirsutum and G. barbadense provides a rough understanding of the regulation of TFs during cotton fiber development.

To provide a resource for cotton genetics and breeding, an interspecific hybridization between Gossypium hirsutum cv. Emian22 and G. barbadense acc. 3-79 was made. A population of 54 BILs (backcross inbred lines, BC1F8) was developed with the aim of transferring G. barbadense genes into G. hirsutum in order to genetically analyze these genes' function in a G. hirsutum background and create new germplasms for breeding. Preliminary investigation of the morphological traits showed that the BILs had diverse variations in plant architecture, seed size, and fuzz color; the related traits of yield and fiber quality evaluated in 4 environments also showed abundant phenotypic variation. In order to explore the molecular diversity of the BIL population, 446 SSR markers selected at an average genetic distance of 10 cM from our interspecific linkage map were used to genotype the BIL population. A total of 393 polymorphic loci accounting for 84.4% MAF (major allele frequency) > 0.05 and 922 allele loci were detected, and the Shannon diversity index (I) was 0.417 per locus. The average introgression segment length was 16.24 cM, and an average of 29.53 segments were introgressed in each BIL line with an average background recovery of 79.8%. QTL mapping revealed 58 QTL associated with fiber quality and yield traits, and 47 favored alleles derived from the donor parent were discovered. This study demonstrated that the interspecific BIL population was enriched with much phenotypic and molecular variation which could be a resource for cotton genetics and breeding.