Chromosomal assignment of RFLP linkage groups harboring important QTLs on an intraspecific cotton (Gossypium hirsutum L.) Joinmap

Multiple domestication events explain the origin of Gossypium hirsutum landraces in Mexico

Several Mesoamerican crops constitute wild-to-domesticated complexes generated by multiple initial domestication events, and continuous gene flow among crop populations and between these populations and their wild relatives. It has been suggested that the domestication of cotton (Gossypium hirsutum) started in the northwest of the Yucatán Peninsula, from where it spread to other regions inside and outside of Mexico. We tested this hypothesis by assembling chloroplast genomes of 23 wild, landraces, and breeding lines (transgene-introgressed and conventional). The phylogenetic analysis showed that the evolutionary history of cotton in Mexico involves multiple events of introgression and genetic divergence. From this, we conclude that Mexican landraces arose from multiple wild populations. Our results also revealed that their structural and functional chloroplast organizations had been preserved. However, genetic diversity decreases as a consequence of domestication, mainly in transgene-introgressed (TI) individuals (π = 0.00020, 0.00001, 0.00016, 0, and 0, of wild, TI-wild, landraces, TI-landraces, and breeding lines, respectively). We identified homologous regions that differentiate wild from domesticated plants and indicate a relationship among the samples. A decrease in genetic diversity associated with transgene introgression in cotton was identified for the first time, and our outcomes are therefore relevant to both biosecurity and agrobiodiversity conservation.

Identification of Stable and Multiple Environment Interaction QTLs and Candidate Genes for Fiber Productive Traits Under Irrigated and Water Stress Conditions Using Intraspecific RILs of Gossypium hirsutum var. MCU5 X TCH1218

Cotton productivity under water-stressed conditions is controlled by multiple quantitative trait loci (QTL). Enhancement of these productivity traits under water deficit stress is crucial for the genetic improvement of upland cotton, Gossypium hirsutum. In the present study, we constructed a genetic map with 504 single nucleotide polymorphisms (SNPs) covering a total span length of 4,416 cM with an average inter-marker distance of 8.76 cM. A total of 181 intra-specific recombinant inbred lines (RILs) were derived from a cross between G. hirsutum var. MCU5 and TCH1218 were used. Although 2,457 polymorphic SNPs were detected between the parents using the CottonSNP50K assay, only 504 SNPs were found to be useful for the construction of the genetic map. In the SNP genotyping, a large number of SNPs showed either >20% missing data, duplication, or segregation distortion. However, the mapped SNPs of this study showed collinearity with the physical map of the reference genome (G. hirsutum var.TM-1), indicating that there was no chromosomal rearrangement within the studied mapping population. RILs were evaluated under multi-environments and seasons for which the phenotypic data were acquired. A total of 53 QTL controlling plant height (PH), number of sympodial branches, boll number (BN), and boll weight (BW) were dissected by QTL analysis under irrigated and water stress conditions. Additionally, it was found that nine QTL hot spots not only co-localized for more than one investigated trait but were also stable with major QTL, i.e., with > 10% of phenotypic variation. One QTL hotspot on chromosome 22 flanked by AX-182254626-AX-182264770 with a span length of 89.4 cM co-localized with seven major and stable QTL linked to a number of sympodial branches both under irrigated and water stress conditions. In addition, putative candidate genes associated with water stress in the QTL hotspots were identified. Besides, few QTL from the hotspots were previously reported across various genetic architects in cotton validating the potential applications of these identified QTL for cotton breeding and improvement. Thus, the major and stable QTL identified in the present study would improve the cotton productivity under water-limited environments through marker-assisted selection.

Construction of a high-density genetic map and identification of QTLs related to agronomic and physiological traits in an interspecific (Gossypium hirsutum × Gossypium barbadense) F2 population

Background

Advances in genome sequencing technology, particularly restriction-site associated DNA sequence (RAD-seq) and whole-genome resequencing, have greatly aided the construction of cotton interspecific genetic maps based on single nucleotide polymorphism (SNPs), Indels, and other types of markers. High-density genetic maps can improve accuracy of quantitative trait locus (QTL) mapping, narrow down location intervals, and facilitate identification of the candidate genes.

Result

In this study, 249 individuals from an interspecific F2 population (TM-1 and Hai7124) were re-sequenced, yielding 6303 high-confidence bin markers spanning 5057.13 cM across 26 cotton chromosomes. A total of 3380 recombination hot regions RHRs were identified which unevenly distributed on the 26 chromosomes. Based on this map, 112 QTLs relating to agronomic and physiological traits from seedling to boll opening stage were identified, including 15 loci associated with 14 traits that contained genes harboring nonsynonymous SNPs. We analyzed the sequence and expression of these ten candidate genes and discovered that GhRHD3 (GH_D10G0500) may affect fiber yield while GhGPAT6 (GH_D04G1426) may affect photosynthesis efficiency.

Conclusion

Our research illustrates the efficiency of constructing a genetic map using binmap and QTL mapping on the basis of a certain size of the early-generation population. High-density genetic map features high recombination exchanges in number and distribution. The QTLs and the candidate genes identified based on this high-density genetic map may provide important gene resources for the genetic improvement of cotton.

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.

Identification of small effect quantitative trait loci of plant architectural, flowering, and early maturity traits in reciprocal interspecific introgression population in cotton

Plant architecture, flowering time and maturity traits are important determinants of yield and fiber quality of cotton. Genetic dissection of loci determining these yield and quality components is complicated by numerous loci with alleles conferring small differences. Therefore, mapping populations segregating for smaller numbers and sizes of introgressed segments is expected to facilitate dissection of these complex quantitative traits. At an advanced stage in the development of reciprocal advanced backcross populations from crosses between elite Gossypium hirsutum cultivar 'Acala Maxxa' (GH) and G. barbadense 'Pima S6' (GB), we undertook mapping of plant architectural traits, flowering time and maturity. A total of 284 BC₄F₁ and BC₄F₂ progeny rows, 120 in GH and 164 in GB background, were evaluated for phenotype, with only 4 and 3 (of 7) traits showing significant differences among progenies. Genotyping by sequencing yielded 3,186 and 3,026 SNPs, respectively, that revealed a total of 27 QTLs in GH background and 22 in GB, for plant height, days to flowering, residual flowering at maturity and maturity. More than of 90% QTLs identified in both backgrounds had small effects (%PV < 10), supporting the merit of this population structure to reduce background noise and small effect QTLs. Germplasm developed in this study may serve as potential pre-breeding material to develop improved cotton cultivars.

Quantitative Trait Loci and Transcriptome Analysis Reveal Genetic Basis of Fiber Quality Traits in CCRI70 RIL Population of Gossypium hirsutum

Upland cotton (Gossypium hirsutum) is widely planted around the world for its natural fiber, and producing high-quality fiber is essential for the textile industry. CCRI70 is a hybrid cotton plant harboring superior yield and fiber quality, whose recombinant inbred line (RIL) population was developed from two upland cotton varieties (sGK156 and 901-001) and were used here to investigate the source of high-quality related alleles. Based on the material of the whole population, a high-density genetic map was constructed using specific locus-amplified fragment sequencing (SLAF-seq). It contained 24,425 single nucleotide polymorphism (SNP) markers, spanning a distance of 4,850.47 centimorgans (cM) over 26 chromosomes with an average marker interval of 0.20 cM. In evaluating three fiber quality traits in nine environments to detect multiple environments stable quantitative trait loci (QTLs), we found 289 QTLs, of which 36 of them were stable QTLs and 18 were novel. Based on the transcriptome analysis for two parents and two RILs, 24,941 unique differentially expressed genes (DEGs) were identified, 473 of which were promising genes. For the fiber strength (FS) QTLs, 320 DEGs were identified, suggesting that pectin synthesis, phenylpropanoid biosynthesis, and plant hormone signaling pathways could influence FS, and several transcription factors may regulate fiber development, such as GAE6, C4H, OMT1, AFR18, EIN3, bZIP44, and GAI. Notably, the marker D13_56413025 in qFS-chr18-4 provides a potential basis for enhancing fiber quality of upland cotton via marker-assisted breeding and gene cloning of important fiber quality traits.

Eco-Geography of Feral Cotton: A Missing Piece in the Puzzle of Gene Flow Dynamics Among Members of Gossypium hirsutum Primary Gene Pool

Mexico is the center of origin and genetic diversity of upland cotton (Gossypium hirsutum L.), the most important source of natural fiber in the world. Currently, wild and domesticated populations (including genetically modified [GM] varieties) occur in this country and gene flow among them has shaped the species’ genetic diversity and structure, setting a complex and challenging scenario for its conservation. Moreover, recent gene flow from GM cultivars to wild Mexican cotton populations has been reported since 2011. In situ conservation of G. hirsutum requires knowledge about the extent of its geographic distribution, both wild and domesticated, as well as the possible routes and mechanisms that contribute to gene flow between the members of the species wild-to-domesticated continuum (i.e., the primary gene pool). However, little is known about the distribution of feral populations that could facilitate gene flow by acting as bridges. In this study, we analyzed the potential distribution of feral cotton based on an ecological niche modeling approach and discussed its implications in the light of the distribution of wild and domesticated cotton. Then, we examined the processes that could be leading to the escape of seeds from the cultivated fields. Our results indicate that the climatic suitability of feral plants in the environmental and geographic space is broad and overlaps with areas of wild cotton habitat and crop fields, suggesting a region that could bridge cultivated cotton and its wild relatives by allowing gene flow between them. This study provides information for management efforts focused on the conservation of wild populations, native landraces, and non-GM domesticated cotton at its center of origin and genetic diversity.

Conservation and Divergence in Duplicated Fiber Coexpression Networks Accompanying Domestication of the Polyploid Gossypium hirsutum L

Gossypium hirsutum L. (Upland cotton) has an evolutionary history involving inter-genomic hybridization, polyploidization, and subsequent domestication. We analyzed the developmental dynamics of the cotton fiber transcriptome accompanying domestication using gene coexpression networks for both joint and homoeologous networks. Remarkably, most genes exhibited expression for at least one homoeolog, confirming previous reports of widespread gene usage in cotton fibers. Most coexpression modules comprising the joint network are preserved in each subgenomic network and are enriched for similar biological processes, showing a general preservation of network modular structure for the two co-resident genomes in the polyploid. Interestingly, only one fifth of homoeologs co-occur in the same module when separated, despite similar modular structures between the joint and homoeologous networks. These results suggest that the genome-wide divergence between homoeologous genes is sufficient to separate their co-expression profiles at the intermodular level, despite conservation of intramodular relationships within each subgenome. Most modules exhibit D-homoeolog expression bias, although specific modules do exhibit A-homoeolog bias. Comparisons between wild and domesticated coexpression networks revealed a much tighter and denser network structure in domesticated fiber, as evidenced by its fewer modules, 13-fold increase in the number of development-related module member genes, and the poor preservation of the wild network topology. These results demonstrate the amazing complexity that underlies the domestication of cotton fiber.

Genetic Analysis of the Transition from Wild to Domesticated Cotton (Gossypium hirsutum L.)

The evolution and domestication of cotton is of great interest from both economic and evolutionary standpoints. Although many genetic and genomic resources have been generated for cotton, the genetic underpinnings of the transition from wild to domesticated cotton remain poorly known. Here we generated an intraspecific QTL mapping population specifically targeting domesticated cotton phenotypes. We used 466 F₂ individuals derived from an intraspecific cross between the wild Gossypium hirsutum var. yucatanense (TX2094) and the elite cultivar G. hirsutum cv. Acala Maxxa, in two environments, to identify 120 QTL associated with phenotypic changes under domestication. While the number of QTL recovered in each subpopulation was similar, only 22 QTL were considered coincident (i.e., shared) between the two locations, eight of which shared peak markers. Although approximately half of QTL were located in the A-subgenome, many key fiber QTL were detected in the D-subgenome, which was derived from a species with unspinnable fiber. We found that many QTL are environment-specific, with few shared between the two environments, indicating that QTL associated with G. hirsutum domestication are genomically clustered but environmentally labile. Possible candidate genes were recovered and are discussed in the context of the phenotype. We conclude that the evolutionary forces that shape intraspecific divergence and domestication in cotton are complex, and that phenotypic transformations likely involved multiple interacting and environmentally responsive factors.

Genome-wide quantitative trait loci reveal the genetic basis of cotton fibre quality and yield-related traits in a Gossypium hirsutum recombinant inbred line population

Cotton is widely cultivated globally because it provides natural fibre for the textile industry and human use. To identify quantitative trait loci (QTLs)/genes associated with fibre quality and yield, a recombinant inbred line (RIL) population was developed in upland cotton. A consensus map covering the whole genome was constructed with three types of markers (8295 markers, 5197.17 centimorgans (cM)). Six fibre yield and quality traits were evaluated in 17 environments, and 983 QTLs were identified, 198 of which were stable and mainly distributed on chromosomes 4, 6, 7, 13, 21 and 25. Thirty-seven QTL clusters were identified, in which 92.8% of paired traits with significant medium or high positive correlations had the same QTL additive effect directions, and all of the paired traits with significant medium or high negative correlations had opposite additive effect directions. In total, 1297 genes were discovered in the QTL clusters, 414 of which were expressed in two RNA-Seq data sets. Many genes were discovered, 23 of which were promising candidates. Six important QTL clusters that included both fibre quality and yield traits were identified with opposite additive effect directions, and those on chromosome 13 (qClu-chr13-2) could increase fibre quality but reduce yield; this result was validated in a natural population using three markers. These data could provide information about the genetic basis of cotton fibre quality and yield and help cotton breeders to improve fibre quality and yield simultaneously.

Enhancing Upland cotton for drought resilience, productivity, and fiber quality: comparative evaluation and genetic dissection

To provision the world sustainably, modern society must increase overall crop production, while conserving and preserving natural resources. Producing more with diminishing water resources is an especially daunting endeavor. Toward the goal of genetically improving drought resilience of cultivated Upland cotton (Gossypium hirsutum L.), this study addresses the genetics of differential yield components referred to as productivity and fiber quality traits under regular-water versus low-water (LW) field conditions. We used ten traits to assess water stress deficit, which included six productivity and four fiber quality traits on two recombinant inbred line (RIL) populations from reciprocally crossed cultivars, Phytogen 72 and Stoneville 474. To facilitate genetic inferences, we genotyped RILs with the CottonSNP63K array, assembled high-density linkage maps of over 7000 SNPs and then analyzed quantitative trait variations. Analysis of variance revealed significant differences for all traits (p < 0.05) in these RIL populations. Although the LW irrigation regime significantly reduced all traits, except lint percent, the RILs exhibited a broad phenotypic spectrum of heritable differences across the water regimes. Transgressive segregation occurred among the RILs, suggesting the possibility of genetic gain through phenotypic selection for drought resilience and perhaps through marker-based selection. Analyses revealed more than 150 quantitative trait loci (QTLs) associated with productivity and fiber quality traits (p < 0.005) on different genomic regions of the cotton genome. The multiple-QTL models analysis with LOD > 3.0 detected 21 QTLs associated with productivity and 22 QTLs associated with fiber quality. For fiber traits, strong clustering and QTL associations occurred in c08 and its homolog c24 as well as c10, c14, and c21. Using contemporary genome sequence assemblies and bioinformatically related information, the identification of genomic regions associated with responses to plant stress/drought elevates the possibility of using marker-assisted and omics-based selection to enhance breeding for drought resilient cultivars and identifying candidate genes and networks. RILs with different responses to drought indicated that it is possible to maintain high fiber quality under LW conditions or reduce the of LW impact on quality. The heritable variation among elite bi-parental RILs for productivity and quality under field drought conditions, and their association of QTLs, and thus specific genomic regions, indicate opportunities for breeding-based gains in water resource conservation, i.e., enhancing cotton's agricultural sustainability.

QTL mapping and genetic effect of chromosome segment substitution lines with excellent fiber quality from Gossypium hirsutum × Gossypium barbadense

Chromosome segment substitution lines (CSSLs) are ideal materials for identifying genetic effects. In this study, CSSL MBI7561 with excellent fiber quality that was selected from BC₄F_3:5 of CCRI45 (Gossypium hirsutum) × Hai1 (Gossypium barbadense) was used to construct 3 secondary segregating populations with 2 generations (BC₅F₂ and BC₅F_2:3). Eighty-one polymorphic markers related to 33 chromosome introgressive segments on 18 chromosomes were finally screened using 2292 SSR markers which covered the whole tetraploid cotton genome. A total of 129 quantitative trait loci (QTL) associated with fiber quality (103) and yield-related traits (26) were detected on 17 chromosomes, explaining 0.85-30.35% of the phenotypic variation; 39 were stable (30.2%), 53 were common (41.1%), 76 were new (58.9%), and 86 had favorable effects on the related traits. More QTL were distributed in the Dt subgenome than in the At subgenome. Twenty-five stable QTL clusters (with stable or common QTL) were detected on 22 chromosome introgressed segments. Finally, the 6 important chromosome introgressed segments (Seg-A02-1, Seg-A06-1, Seg-A07-2, Seg-A07-3, Seg-D07-3, and Seg-D06-2) were identified as candidate chromosome regions for fiber quality, which should be given more attention in future QTL fine mapping, gene cloning, and marker-assisted selection (MAS) breeding.

Molecular Evolution and Stress and Phytohormone Responsiveness of SUT Genes in Gossypium hirsutum

Sucrose transporters (SUTs) play key roles in allocating the translocation of assimilates from source to sink tissues. Although the characteristics and biological roles of SUTs have been intensively investigated in higher plants, this gene family has not been functionally characterized in cotton. In this study, we performed a comprehensive analysis of SUT genes in the tetraploid cotton Gossypium hirsutum. A total of 18 G. hirsutum SUT genes were identified and classified into three groups based on their evolutionary relationships. Up to eight SUT genes in G. hirsutum were placed in the dicot-specific SUT1 group, while four and six SUT genes were, respectively, clustered into SUT4 and SUT2 groups together with members from both dicot and monocot species. The G. hirsutum SUT genes within the same group displayed similar exon/intron characteristics, and homologous genes in G. hirsutum At and Dt subgenomes, G. arboreum, and G. raimondii exhibited one-to-one relationships. Additionally, the duplicated genes in the diploid and polyploid cotton species have evolved through purifying selection, suggesting the strong conservation of SUT loci in these species. Expression analysis in different tissues indicated that SUT genes might play significant roles in cotton fiber elongation. Moreover, analyses of cis-acting regulatory elements in promoter regions and expression profiling under different abiotic stress and exogenous phytohormone treatments implied that SUT genes, especially GhSUT6A/D, might participate in plant responses to diverse abiotic stresses and phytohormones. Our findings provide valuable information for future studies on the evolution and function of SUT genes in cotton.

The Mating System of the Wild-to-Domesticated Complex of Gossypium hirsutum L. Is Mixed

The domestication syndrome of many plants includes changes in their mating systems. The evolution of the latter is shaped by ecological and genetic factors that are particular to an area. Thus, the reproductive biology of wild relatives must be studied in their natural distribution to understand the mating system of a crop species as a whole. Gossypium hirsutum (upland cotton) includes both domesticated varieties and wild populations of the same species. Most studies on mating systems describe cultivated cotton as self-pollinated, while studies on pollen dispersal report outcrossing; however, the mating system of upland cotton has not been described as mixed and little is known about its wild relatives. In this study we selected two wild metapopulations for comparison with domesticated plants and one metapopulation with evidence of recent gene flow between wild relatives and the crop to evaluate the mating system of cotton's wild-to-domesticated complex. Using classic reproductive biology methods, our data demonstrate that upland cotton presents a mixed mating system throughout the complex. Given cotton's capacity for outcrossing, differences caused by the domestication process in cultivated individuals can have consequences for its wild relatives. This characterization of the diversity of the wild relatives in their natural distribution, as well as their interactions with the crop, will be useful to design and implement adequate strategies for conservation and biosecurity.

QTL delineation for five fiber quality traits based on an intra-specific Gossypium hirsutum L. recombinant inbred line population

Gossypium hirsutum L. is the most important fiber crop worldwide and contributes to more than 95% of global cotton production. Marker-assisted selection (MAS) is an effective approach for improving fiber quality, and quantitative trait loci (QTL) mapping of fiber quality traits is important for cotton breeding. In this study, a permanent intra-specific recombinant inbred line (RIL) population containing 137 families was used for fiber quality testing. Based on a previously reported high-density genetic map with an average marker distance of 0.63 cM, 186 additive QTLs were obtained for five fiber quality traits over five consecutive years, including 39 for fiber length (FL), 36 for fiber strength (FS), 50 for fiber uniformity (FU), 33 for micronaire (MC) and 28 for fiber elongation (FE). Three stable QTLs, qMC-A4-1, qMC-D2-3 and qFS-D9-1, were detected in four datasets, and another eight stable QTLs, qMC-A4-2, qMC-D11-2, qFU-A9-1, qFU-A10-4, qFS-D11-1, qFL-D9-2, qFL-D11-1 and qFE-A3-2, were detected in three datasets. The annotated genes in these 11 stable QTLs were collected, and these genes included many transcription factors with functions during fiber development. 33 QTL coincidence regions were found, and these involved nearly half of the total QTLs. Four chromosome regions containing at least 6 QTLs were promising for fine mapping. In addition, 41 pairs of epistatic QTLs (e-QTLs) were screened, including 6 for FL, 30 for FS, 2 for FU and 3 for MC. The identification of stable QTLs adds valuable information for further QTL fine mapping and gene positional cloning for fiber quality genetic detection and provides useful markers for further molecular breeding in enhancing fiber quality.

Insights Into Upland Cotton (Gossypium hirsutum L.) Genetic Recombination Based on 3 High-Density Single-Nucleotide Polymorphism and a Consensus Map Developed Independently With Common Parents

High-density linkage maps are vital to supporting the correct placement of scaffolds and gene sequences on chromosomes and fundamental to contemporary organismal research and scientific approaches to genetic improvement, especially in paleopolyploids with exceptionally complex genomes, eg, upland cotton (Gossypium hirsutum L., "2n = 52"). Three independently developed intraspecific upland mapping populations were analyzed to generate 3 high-density genetic linkage single-nucleotide polymorphism (SNP) maps and a consensus map using the CottonSNP63K array. The populations consisted of a previously reported F2, a recombinant inbred line (RIL), and reciprocal RIL population, from "Phytogen 72" and "Stoneville 474" cultivars. The cluster file provided 7417 genotyped SNP markers, resulting in 26 linkage groups corresponding to the 26 chromosomes (c) of the allotetraploid upland cotton (AD)1 arisen from the merging of 2 genomes ("A" Old World and "D" New World). Patterns of chromosome-specific recombination were largely consistent across mapping populations. The high-density genetic consensus map included 7244 SNP markers that spanned 3538 cM and comprised 3824 SNP bins, of which 1783 and 2041 were in the At and Dt subgenomes with 1825 and 1713 cM map lengths, respectively. Subgenome average distances were nearly identical, indicating that subgenomic differences in bin number arose due to the high numbers of SNPs on the Dt subgenome. Examination of expected recombination frequency or crossovers (COs) on the chromosomes within each population of the 2 subgenomes revealed that COs were also not affected by the SNPs or SNP bin number in these subgenomes. Comparative alignment analyses identified historical ancestral At-subgenomic translocations of c02 and c03, as well as of c04 and c05. The consensus map SNP sequences aligned with high congruency to the NBI assembly of Gossypium hirsutum. However, the genomic comparisons revealed evidence of additional unconfirmed possible duplications, inversions and translocations, and unbalance SNP sequence homology or SNP sequence/loci genomic dominance, or homeolog loci bias of the upland tetraploid At and Dt subgenomes. The alignments indicated that 364 SNP-associated previously unintegrated scaffolds can be placed in pseudochromosomes of the NBI G hirsutum assembly. This is the first intraspecific SNP genetic linkage consensus map assembled in G hirsutum with a core of reproducible mendelian SNP markers assayed on different populations and it provides further knowledge of chromosome arrangement of genic and nongenic SNPs. Together, the consensus map and RIL populations provide a synergistically useful platform for localizing and identifying agronomically important loci for improvement of the cotton crop.

Association mapping analysis of fiber yield and quality traits in Upland cotton (Gossypium hirsutum L.)

Fiber yield and quality are the most important traits for Upland cotton (Gossypium hirsutum L.). Identifying high yield and good fiber quality genes are the prime concern of researchers in cotton breeding. Association mapping offers an alternative and powerful method for detecting those complex agronomic traits. In this study, 198 simple sequence repeats (SSRs) were used to screen markers associated with fiber yield and quality traits with 302 elite Upland cotton accessions that were evaluated in 12 locations representing the Yellow River and Yangtze River cotton growing regions of China. Three subpopulations were found after the estimation of population structure. The pair-wise kinship values varied from 0 to 0.867. Only 1.59% of the total marker locus pairs showed significant linkage disequilibrium (LD, p < 0.001). The genome-wide LD decayed within the genetic distance of ~30 to 32 cM at r ² = 0.1, and decreased to ~1 to 2 cM at r ² = 0.2, indicating the potential for association mapping. Analysis based on a mixed linear model detected 57 significant (p < 0.01) marker-trait associations, including seven associations for fiber length, ten for fiber micronaire, nine for fiber strength, eight for fiber elongation, five for fiber uniformity index, five for fiber uniformity ratio, six for boll weight and seven for lint percent, for a total of 35 SSR markers, of which 11 markers were associated with more than one trait. Among marker-trait associations, 24 associations coincided with the previously reported quantitative trait loci (QTLs), the remainder were newly identified QTLs/genes. The QTLs identified in this study will potentially facilitate improvement of fiber yield and quality in the future cotton molecular breeding programs.

A comparative transcriptome analysis of two sets of backcross inbred lines differing in lint-yield derived from a Gossypium hirsutum × Gossypium barbadense population

Upland cotton (Gossypium hirsutum L.) is the most important fiber crop, and its lint-yield improvement is impeded due to its narrow genetic base and the lack of understanding of the genetic basis of yield. Backcross inbred lines (BILs) or near-isogenic lines (NILs) in the same genetic background differing in lint yield, developed through advanced backcrossing, provide an important genomic resource to study the molecular genetic basis of lint yield. In the present study, a high-yield (HY) group and a low-yield (LY) group each with three BILs were selected from a BIL population between G. hirsutum and G. barbadense. Using a microarray-based comparative transcriptome analysis on developing fibers at 10 days post-anthesis (DPA) between the two groups, 1486 differentially expressed genes (DEGs) were identified. A total of 212 DEGs were further mapped in the regions of 24 yield QTL and 11 yield trait QTL hotspots as reported previously, and 81 DEGs mapped with the 7 lint-yield QTL identified in the BIL population from which the two sets of BILs were selected. Gene Ontology annotations and Blast-Mapping-Annotation-KEGG analysis via Blast2GO revealed that more DEGs were associated with catalytic activity and binding, followed by transporters, nucleic acid binding transcription factors, structural molecules and molecular transducer activities. Six DEGs were chosen for a quantitative RT-PCR assay, and the results were consistent with the microarray analysis. The development of DEGs-based markers revealed that 7 single strand conformation polymorphism-based single nucleotide polymorphic (SSCP-SNP) markers were associated with yield traits, and 3 markers with lint yield. In the present study, we identified a number of yield and yield component QTL-co-localizing DEGs and developed several DEG-based SSCP-SNP markers for the traits, thereby providing a set of candidate genes for molecular breeding and genetic manipulation of lint yield in cotton.

Identification of stable quantitative trait loci (QTLs) for fiber quality traits across multiple environments in Gossypium hirsutum recombinant inbred line population

BACKGROUND

The identification of quantitative trait loci (QTLs) that are stable and consistent across multiple environments and populations plays an essential role in marker-assisted selection (MAS). In the present study, we used 28,861 simple sequence repeat (SSR) markers, which included 12,560 Gossypium raimondii (D genome) sequence-based SSR markers to identify polymorphism between two upland cotton strains 0-153 and sGK9708. A total of 851 polymorphic primers were finally selected and used to genotype 196 recombinant inbred lines (RIL) derived from a cross between 0 and 153 and sGK9708 and used to construct a linkage map. The RIL population was evaluated for fiber quality traits in six locations in China for five years. Stable QTLs identified in this intraspecific cross could be used in future cotton breeding program and with fewer obstacles.

RESULTS

The map covered a distance of 4,110 cM, which represents about 93.2 % of the upland cotton genome, and with an average distance of 5.2 cM between adjacent markers. We identified 165 QTLs for fiber quality traits, of which 47 QTLs were determined to be stable across multiple environments. Most of these QTLs aggregated into clusters with two or more traits. A total of 30 QTL clusters were identified which consisted of 103 QTLs. Sixteen clusters in the At sub-genome comprised 44 QTLs, whereas 14 clusters in the Dt sub-genome that included 59 QTLs for fiber quality were identified. Four chromosomes, including chromosome 4 (c4), c7, c14, and c25 were rich in clusters harboring 5, 4, 5, and 6 clusters respectively. A meta-analysis was performed using Biomercator V4.2 to integrate QTLs from 11 environmental datasets on the RIL populations of the above mentioned parents and previous QTL reports. Among the 165 identified QTLs, 90 were identified as common QTLs, whereas the remaining 75 QTLs were determined to be novel QTLs. The broad sense heritability estimates of fiber quality traits were high for fiber length (0.93), fiber strength (0.92), fiber micronaire (0.85), and fiber uniformity (0.80), but low for fiber elongation (0.27). Meta-clusters on c4, c7, c14 and c25 were identified as stable QTL clusters and were considered more valuable in MAS for the improvement of fiber quality of upland cotton.

CONCLUSION

Multiple environmental evaluations of an intraspecific RIL population were conducted to identify stable QTLs. Meta-QTL analyses identified a common chromosomal region that plays an important role in fiber development. Therefore, QTLs identified in the present study are an ideal candidate for MAS in cotton breeding programs to improve fiber quality.

Identification of stable quantitative trait loci (QTLs) for fiber quality traits across multiple environments in Gossypium hirsutum recombinant inbred line population

BACKGROUND

The identification of quantitative trait loci (QTLs) that are stable and consistent across multiple environments and populations plays an essential role in marker-assisted selection (MAS). In the present study, we used 28,861 simple sequence repeat (SSR) markers, which included 12,560 Gossypium raimondii (D genome) sequence-based SSR markers to identify polymorphism between two upland cotton strains 0-153 and sGK9708. A total of 851 polymorphic primers were finally selected and used to genotype 196 recombinant inbred lines (RIL) derived from a cross between 0 and 153 and sGK9708 and used to construct a linkage map. The RIL population was evaluated for fiber quality traits in six locations in China for five years. Stable QTLs identified in this intraspecific cross could be used in future cotton breeding program and with fewer obstacles.

RESULTS

The map covered a distance of 4,110 cM, which represents about 93.2 % of the upland cotton genome, and with an average distance of 5.2 cM between adjacent markers. We identified 165 QTLs for fiber quality traits, of which 47 QTLs were determined to be stable across multiple environments. Most of these QTLs aggregated into clusters with two or more traits. A total of 30 QTL clusters were identified which consisted of 103 QTLs. Sixteen clusters in the At sub-genome comprised 44 QTLs, whereas 14 clusters in the Dt sub-genome that included 59 QTLs for fiber quality were identified. Four chromosomes, including chromosome 4 (c4), c7, c14, and c25 were rich in clusters harboring 5, 4, 5, and 6 clusters respectively. A meta-analysis was performed using Biomercator V4.2 to integrate QTLs from 11 environmental datasets on the RIL populations of the above mentioned parents and previous QTL reports. Among the 165 identified QTLs, 90 were identified as common QTLs, whereas the remaining 75 QTLs were determined to be novel QTLs. The broad sense heritability estimates of fiber quality traits were high for fiber length (0.93), fiber strength (0.92), fiber micronaire (0.85), and fiber uniformity (0.80), but low for fiber elongation (0.27). Meta-clusters on c4, c7, c14 and c25 were identified as stable QTL clusters and were considered more valuable in MAS for the improvement of fiber quality of upland cotton.

CONCLUSION

Multiple environmental evaluations of an intraspecific RIL population were conducted to identify stable QTLs. Meta-QTL analyses identified a common chromosomal region that plays an important role in fiber development. Therefore, QTLs identified in the present study are an ideal candidate for MAS in cotton breeding programs to improve fiber quality.

Construction of an EST-SSR-based interspecific transcriptome linkage map of fibre development in cotton

Quantitative trait locus (QTL) mapping is an important method in marker-assisted selection breeding. Many studies on the QTLs focus on cotton fibre yield and quality; however, most are conducted at the DNA level, which may reveal null QTLs. Hence, QTL mapping based on transcriptome maps at the cDNA level is often more reliable. In this study, an interspecific transcriptome map of allotetraploid cotton was developed based on an F2 population (Emian22 x 3-79) by amplifying cDNA using EST-SSRs. The map was constructed using cDNA obtained from developing fibres at five days post anthesis (DPA). A total of 1270 EST-SSRs were screened for polymorphisms between the mapping parents. The resulting transcriptome linkage map contained 242 markers that were distributed in 32 linkage groups (26 chromosomes). The full length of this map is 1938.72 cM with a mean marker distance of 8.01 cM. The functions of some ESTs have been annotated by exploring homologous sequences. Some markers were related to the differentiation and elongation of cotton fibre, while most were related to the basic metabolism. This study demonstrates that constructing a transcriptome linkage map by amplifying cDNAs using EST-SSRs is a simple and practical method as well as a powerful tool to map eQTLs for fibre quality and other traits in cotton.

High Resolution Consensus Mapping of Quantitative Trait Loci for Fiber Strength, Length and Micronaire on Chromosome 25 of the Upland Cotton (Gossypium hirsutum L.)

Cotton (Gossypium hirsutum L.) is an important agricultural crop that provides renewable natural fiber resources for the global textile industry. Technological developments in the textile industry and improvements in human living standards have increased the requirement for supplies and better quality cotton. Upland cotton 0-153 is an elite cultivar harboring strong fiber strength genes. To conduct quantitative trait locus (QTL) mapping for fiber quality in 0-153, we developed a population of 196 recombinant inbred lines (RILs) from a cross between 0-153 and sGK9708. The fiber quality traits in 11 environments were measured and a genetic linkage map of chromosome 25 comprising 210 loci was constructed using this RIL population, mainly using simple sequence repeat markers and single nucleotide polymorphism markers. QTLs were identified across diverse environments using the composite interval mapping method. A total of 37 QTLs for fiber quality traits were identified on chromosome 25, of which 17 were stably expressed in at least in two environments. A stable fiber strength QTL, qFS-chr25-4, which was detected in seven environments and was located in the marker interval between CRI-SNP120491 and BNL2572, could explain 6.53%-11.83% of the observed phenotypic variations. Meta-analysis also confirmed the above QTLs with previous reports. Application of these QTLs could contribute to improving fiber quality and provide information for marker-assisted selection.

The Negative Correlation between Fiber Color and Quality Traits Revealed by QTL Analysis

Naturally existing colored cotton was far from perfection due to having genetic factors for lower yield, poor fiber quality and monotonous color. These factors posed a challenge to colored cotton breeding and innovation. To identify novel quantitative trait loci (QTL) for fiber color along with understanding of correlation between fiber color and quality in colored cotton, a RIL and two F2 populations were generated from crosses among Zong128 (Brown fiber cotton) and two white fiber cotton lines which were then analyzed in four environments. Two stable and major QTLs (qLC-7-1, qFC-7-1) for fiber lint and fuzz color were detected accounting for 16.01%-59.85% of the phenotypic variation across multiple generations and environments. Meanwhile, some minor QTLs were also identified on chromosomes 5, 14, 21 and 24 providing low phenotypic variation (<5%) from only F2 populations, not from the RILs population. Especially, a multiple-effect locus for fiber color and quality has been detected between flanking markers NAU1043 and NAU3654 on chromosome 7 (A genome) over multiple environments. Of which, qLC-7-1, qFC-7-1 were responsible for positive effects and improved fiber color in offsprings. Meanwhile, the QTLs (qFL-7-1, qFU-7-1, qFF-7-1, qFE-7-1, and qFS-7-1) for fiber quality had negative effects and explained 2.19%-8.78% of the phenotypic variation. This multiple-effect locus for fiber color and quality may reveal the negative correlation between the two types of above traits, so paving the way towards cotton genetic improvement.

Mapping genomic loci for cotton plant architecture, yield components, and fiber properties in an interspecific (Gossypium hirsutum L. × G. barbadense L.) RIL population

A quantitative trait locus (QTL) mapping was conducted to better understand the genetic control of plant architecture (PA), yield components (YC), and fiber properties (FP) in the two cultivated tetraploid species of cotton (Gossypium hirsutum L. and G. barbadense L.). One hundred and fifty-nine genomic regions were identified on a saturated genetic map of more than 2,500 SSR and SNP markers, constructed with an interspecific recombinant inbred line (RIL) population derived from the genetic standards of the respective cotton species (G. hirsutum acc. TM-1 × G. barbadense acc. 3-79). Using the single nonparametric and MQM QTL model mapping procedures, we detected 428 putative loci in the 159 genomic regions that confer 24 cotton traits in three diverse production environments [College Station F&B Road (FB), TX; Brazos Bottom (BB), TX; and Shafter (SH), CA]. These putative QTL loci included 25 loci for PA, 60 for YC, and 343 for FP, of which 3, 12, and 60, respectively, were strongly associated with the traits (LOD score ≥ 3.0). Approximately 17.7 % of the PA putative QTL, 32.9 % of the YC QTL, and 48.3 % of the FP QTL had trait associations under multiple environments. The At subgenome (chromosomes 1-13) contributed 72.7 % of loci for PA, 46.2 % for YC, and 50.4 % for FP while the Dt subgenome (chromosomes 14-26) contributed 27.3 % of loci for PA, 53.8 % for YC, and 49.6 % for FP. The data obtained from this study augment prior evidence of QTL clusters or gene islands for specific traits or biological functions existing in several non-homoeologous cotton chromosomes. DNA markers identified in the 159 genomic regions will facilitate further dissection of genetic factors underlying these important traits and marker-assisted selection in cotton.

The Li2 mutation results in reduced subgenome expression bias in elongating fibers of allotetraploid cotton (Gossypium hirsutum L.)

Next generation sequencing (RNA-seq) technology was used to evaluate the effects of the Ligon lintless-2 (Li₂) short fiber mutation on transcriptomes of both subgenomes of allotetraploid cotton (Gossypium hirsutum L.) as compared to its near-isogenic wild type. Sequencing was performed on 4 libraries from developing fibers of Li₂ mutant and wild type near-isogenic lines at the peak of elongation followed by mapping and PolyCat categorization of RNA-seq data to the reference D₅ genome (G. raimondii) for homeologous gene expression analysis. The majority of homeologous genes, 83.6% according to the reference genome, were expressed during fiber elongation. Our results revealed: 1) approximately two times more genes were induced in the A_T subgenome comparing to the D_T subgenome in wild type and mutant fiber; 2) the subgenome expression bias was significantly reduced in the Li₂ fiber transcriptome; 3) Li₂ had a significantly greater effect on the D_T than on the A_T subgenome. Transcriptional regulators and cell wall homeologous genes significantly affected by the Li₂ mutation were reviewed in detail. This is the first report to explore the effects of a single mutation on homeologous gene expression in allotetraploid cotton. These results provide deeper insights into the evolution of allotetraploid cotton gene expression and cotton fiber development.

Characterization of expressed sequence tags from developing fibers of Gossypium barbadense and evaluation of insertion-deletion variation in tetraploid cultivated cotton species

Background

Cotton is the leading fiber crop worldwide. Gossypium barbadense is an important species of cotton because of its extra-long staple fibers with superior luster and silkiness. However, a systematic analysis and utilization of cDNA sequences from G. barbadense fiber development remains understudied.

Results

A total of 21,079 high quality sequences were generated from two non-normalized cDNA libraries prepared by using a mixture of G. barbadense Hai7124 fibers and ovules. After assembly processing, a set of 8,653 unigenes were obtained. Of those, 7,786 were matched to known proteins and 7,316 were assigned to functional categories. The molecular functions of these unigenes were mostly related to binding and catalytic activity, and carbohydrate, amino acid, and energy metabolisms were major contributors among the subsets of metabolism. Sequences comparison between G. barbadense and G. hirsutum revealed that 8,245 unigenes from G. barbadense were detected the similarity with those released publicly in G. hirsutum, however, the remaining 408 sequences had no hits against G. hirsutum unigenes database. Furthermore, 13,275 putative ESTs InDels loci involved in the orthologous and/or homoeologous differences between/within G. barbadense and G. hirsutum were discovered by in silico analyses, and 2,160 InDel markers were developed by ESTs with more than five insertions or deletions. By gel electrophoresis combined with sequencing verification, 71.11% candidate InDel loci were reconfirmed orthologous and/or homoeologous loci polymorphisms using G. hirsutum acc TM-1 and G. barbadense cv Hai7124. Blastx result showed among 2,160 InDel loci, 81 with significant function similarity with known genes associated with secondary wall synthesis process, indicating the important roles in fiber quality in tetraploid cultivated cotton species.

Conclusion

Sequence comparisons and InDel markers development will lay the groundwork for promoting the identification of genes related to superior agronomic traits, genetic differentiation and comparative genomic studies between G. hirsutum and G. barbadense.

Construction of genetic linkage map with chromosomal assigment and quantitative trait loci associated with some important agronomic traits in cotton

Cotton is the world's leading natural fiber and second most important oilseed crop and has been a focus of genetic, systematic and breeding research. The genetic and physiological bases of some important agronomic traits in cotton were investigated by QTL mapping through constructing of genetic map with chromosomal assignment. A segregating F2 population derived from an interspecific cross (G. barbadense x G. hirsutum) between two genotypes, cvs. "Giza 83" and "Deltapine" was used in this study. Different molecular markers including SSR, EST, EST-SSR, AFLP and RAPD were employed to identify markers that reveal differences between the parents. In total 42 new markers were merged with 140 previously mapped markers to produce a new map with 182 loci covering a total length of 2370.5 cM. Among these new markers, some of them were used to assign chromosomes to the produced 26 linkage groups. The LG2, LG3, LG11 and LG26 were assigned to chromosomes 1, 6, 5 and 20 respectively. Single point analysis was used to identify genomic regions controlling traits for plant height, number of nodes at flowering time, bolling date, days to flowering and number of bolls. In total 40 significant QTL were identified for the five traits on 11 linkage groups (1, 2, 3, 4, 5, 10, 11, 12, 18, 19 and 23). This work represents an improvement of the previously constructed genetic map in addition to chromosomal assignment and detection of new significant QTL for the five traits in Egyptian cotton. The Significant QTLs detected in this study can be employed in marker assisted selection for molecular breeding programs aiming at developing cotton cultivars with improved agronomic traits.

Identification of quantitative trait loci controlling fibre length and lignin content in Arabidopsis thaliana stems

Fibre properties and the biochemical composition of cell walls are important traits in many applications. For example, the lengths of fibres define the strength and quality of paper, and lignin content is a critical parameter for the use of biomass in biofuel production. Identifying genes controlling these traits is comparatively difficult in woody species, because of long generation times and limited amenability to high-resolution genetic mapping. To address this problem, this study mapped quantitative trait loci (QTLs) defining fibre length and lignin content in the Arabidopsis recombinant inbred line population Col-4 × Ler-0. Adapting high-throughput phenotyping techniques for both traits for measurements in Arabidopsis inflorescence stems identified significant QTLs for fibre length on chromosomes 2 and 5, as well as one significant QTL affecting lignin content on chromosome 2. For fibre length, total variation within the population was 208% higher than between parental lines and the identified QTLs explained 50.58% of the observed variation. For lignin content, the values were 261 and 26.51%, respectively. Bioinformatics analysis of the associated intervals identified a number of candidate genes for fibre length and lignin content. This study demonstrates that molecular mapping of QTLs pertaining to wood and fibre properties is possible in Arabidopsis, which substantially broadens the use of Arabidopsis as a model species for the functional characterization of plant genes.

QTL analysis for transgressive resistance to root-knot nematode in interspecific cotton (Gossypium spp.) progeny derived from susceptible parents

The southern root-knot nematode (RKN, Meloidogyne incognita) is a major soil-inhabiting plant parasite that causes significant yield losses in cotton (Gossypium spp.). Progeny from crosses between cotton genotypes susceptible to RKN produced segregants in subsequent populations which were highly resistant to this parasite. A recombinant inbred line (RIL) population of 138 lines developed from a cross between Upland cotton TM-1 (G. hirsutum L.) and Pima 3–79 (G. barbadense L.), both susceptible to RKN, was used to identify quantitative trait loci (QTLs) determining responses to RKN in greenhouse infection assays with simple sequence repeat (SSR) markers. Compared to both parents, 53.6% and 52.1% of RILs showed less (P<0.05) root-galling index (GI) and had lower (P<0.05) nematode egg production (eggs per gram root, EGR). Highly resistant lines (transgressive segregants) were identified in this RIL population for GI and/or EGR in two greenhouse experiments. QTLs were identified using the single-marker analysis nonparametric mapping Kruskal-Wallis test. Four major QTLs located on chromosomes 3, 4, 11, and 17 were identified to account for 8.0 to 12.3% of the phenotypic variance (R²) in root-galling. Two major QTLs accounting for 9.7% and 10.6% of EGR variance were identified on chromosomes 14 and 23 (P<0.005), respectively. In addition, 19 putative QTLs (P<0.05) accounted for 4.5–7.7% of phenotypic variance (R²) in GI, and 15 QTLs accounted for 4.2–7.3% of phenotypic variance in EGR. In lines with alleles positive for resistance contributed by both parents in combinations of two to four QTLs, dramatic reductions of >50% in both GI and EGR were observed. The transgressive segregants with epistatic effects derived from susceptible parents indicate that high levels of nematode resistance in cotton may be attained by pyramiding positive alleles using a QTL mapping approach.

A high-density simple sequence repeat and single nucleotide polymorphism genetic map of the tetraploid cotton genome

Genetic linkage maps play fundamental roles in understanding genome structure, explaining genome formation events during evolution, and discovering the genetic bases of important traits. A high-density cotton (Gossypium spp.) genetic map was developed using representative sets of simple sequence repeat (SSR) and the first public set of single nucleotide polymorphism (SNP) markers to genotype 186 recombinant inbred lines (RILs) derived from an interspecific cross between Gossypium hirsutum L. (TM-1) and G. barbadense L. (3-79). The genetic map comprised 2072 loci (1825 SSRs and 247 SNPs) and covered 3380 centiMorgan (cM) of the cotton genome (AD) with an average marker interval of 1.63 cM. The allotetraploid cotton genome produced equivalent recombination frequencies in its two subgenomes (At and Dt). Of the 2072 loci, 1138 (54.9%) were mapped to 13 At-subgenome chromosomes, covering 1726.8 cM (51.1%), and 934 (45.1%) mapped to 13 Dt-subgenome chromosomes, covering 1653.1 cM (48.9%). The genetically smallest homeologous chromosome pair was Chr. 04 (A04) and 22 (D04), and the largest was Chr. 05 (A05) and 19 (D05). Duplicate loci between and within homeologous chromosomes were identified that facilitate investigations of chromosome translocations. The map augments evidence of reciprocal rearrangement between ancestral forms of Chr. 02 and 03 versus segmental homeologs 14 and 17 as centromeric regions show homeologous between Chr. 02 (A02) and 17 (D02), as well as between Chr. 03 (A03) and 14 (D03). This research represents an important foundation for studies on polyploid cottons, including germplasm characterization, gene discovery, and genome sequence assembly.

Differential gene expression and associated QTL mapping for cotton yield based on a cDNA-AFLP transcriptome map in an immortalized F2

cDNA-AFLP techniques have found new applications in recent years. Currently, the methodology is used to establish differential gene expression and construct linkage maps. In the present study, a transcriptome map based on cDNA-AFLP techniques was constructed using an immortalized F(2) (IF(2)) population of 171 lines. The lines were derived from intercrosses between 180 recombinant inbred lines (RILs) of the cotton hybrid Xiangzamian 2 (Gossypium. hirsutum L.). A total of 302 transcriptome-derived fragments (TDFs) were mapped onto 26 linkage groups that covered 2,477.06 cM in length with an average distance of 8.23 cM between two markers. Seventy-one QTL for yield and yield component traits were detected by CIM procedures based on four environments, with 13 QTL identified in at least two environments. Some TDFs co-located with yield QTL were subsequently sequenced and analyzed by online homology searches. Potential candidate genes for yield and yield component traits were found to encode proteins involved in DNA replication, transcription, translation, and biosynthesis regulation. Furthermore, genes regulating metabolic processes signal transduction, transport, and structural components of organelles were identified. Correlation analysis between expression patterns of TDFs and trait performance detected six TDFs positively correlated to both yield and yield heterosis: six TDFs positively correlated to yield, and seven TDFs to yield heterosis. These TDFs have potential for cloning the functional genes responsible for each corresponding trait and have future value in marker-assisted selection.

Mapping Fusarium wilt race 1 resistance genes in cotton by inheritance, QTL and sequencing composition

Knowledge of the inheritance of disease resistance and genomic regions housing resistance (R) genes is essential to prevent expanding pathogen threats such as Fusarium wilt [Fusarium oxysporum f.sp. vasinfectum (FOV) Atk. Sny & Hans] in cotton (Gossypium spp.). We conducted a comprehensive study combining conventional inheritance, genetic and quantitative trait loci (QTL) mapping, QTL marker-sequence composition, and genome sequencing to examine the distribution, structure and organization of disease R genes to race 1 of FOV in the cotton genome. Molecular markers were applied to F(2) and recombinant inbred line (RIL) interspecific mapping populations from the crosses Pima-S7 (G. barbadense L.) × 'Acala NemX' (G. hirsutum L.) and Upland TM-1 (G. hirsutum) × Pima 3-79 (G. barbadense), respectively. Three greenhouse tests and one field test were used to obtain sequential estimates of severity index (DSI) of leaves, and vascular stem and root staining (VRS). A single resistance gene model was observed for the F(2) population based on inheritance of phenotypes. However, additional inheritance analyses and QTL mapping indicated gene interactions and inheritance from nine cotton chromosomes, with major QTLs detected on five chromosomes [Fov1-C06, Fov1-C08, (Fov1-C11 ( 1 ) and Fov1-C11 ( 2)) , Fov1-C16 and Fov1-C19 loci], explaining 8-31% of the DSI or VRS variation. The Fov1-C16 QTL locus identified in the F(2) and in the RIL populations had a significant role in conferring FOV race 1 resistance in different cotton backgrounds. Identified molecular markers may have important potential for breeding effective FOV race 1 resistance into elite cultivars by marker-assisted selection. Reconciliation between genetic and physical mapping of gene annotations from marker-DNA and new DNA sequences of BAC clones tagged with the resistance-associated QTLs revealed defenses genes induced upon pathogen infection and gene regions rich in disease-response elements, respectively. These offer candidate gene targets for Fusarium wilt resistance response in cotton and other host plants.

Marker-assisted breeding as next-generation strategy for genetic improvement of productivity and quality: can it be realized in cotton?

The dawdling development in genetic improvement of cotton with conventional breeding program is chiefly due to lack of complete knowledge on and precise manipulation of fiber productivity and quality. Naturally available cotton continues to be a resource for the upcoming breeding program, and contemporary technologies to exploit the available natural variation are outlined in this paper for further improvement of fiber. Particularly emphasis is given to application, obstacles, and perspectives of marker-assisted breeding since it appears to be more promising in manipulating novel genes that are available in the cotton germplasm. Deployment of system quantitative genetics in marker-assisted breeding program would be essential to realize its role in cotton. At the same time, role of genetic engineering and in vitro mutagenesis cannot be ruled out in genetic improvement of cotton.

Association of AFLP and SSR markers with agronomic and fibre quality traits in Gossypium hirsutum L

Molecular markers linked to QTL contributing to agronomic and fibre quality traits would be useful for cotton improvement. We have attempted to tag yield and fibre quality traits with AFLP and SSR markers using F(2) and F(3) populations of a cross between two Gossypium hirsutum varieties, PS56-4 and RS2013. Out of 50 AFLP primer combinations and 177 SSR primer pairs tested, 32 AFLP and four SSR primers were chosen for genotyping F(2) individuals. Marker-trait associations were studied for eight agronomic and five fibre quality traits through simple and multiple regression analysis (MRA) using a set of 92 AFLP polymorphic loci and four SSR markers. Simple linear regression analysis (SLRA) identified 23 markers for eight different traits whereas multiple regression analysis identified 30 markers for at least one of the 13 traits. SSR marker BNL 3502 was consistently identified to be associated with fibre strength. While all the markers identified in SLRA were also detected in MRA, as many as 16 of the 30 markers were identified to be associated with respective traits in both F2 and F3 generations. The markers explained up to 41 per cent of phenotypic variation for individual traits. A number of markers were found to be associated with multiple traits suggesting clustering of QTLs for fibre quality traits in cotton.

Polyploidization altered gene functions in cotton (Gossypium spp.)

Cotton (Gossypium spp.) is an important crop plant that is widely grown to produce both natural textile fibers and cottonseed oil. Cotton fibers, the economically more important product of the cotton plant, are seed trichomes derived from individual cells of the epidermal layer of the seed coat. It has been known for a long time that large numbers of genes determine the development of cotton fiber, and more recently it has been determined that these genes are distributed across At and Dt subgenomes of tetraploid AD cottons. In the present study, the organization and evolution of the fiber development genes were investigated through the construction of an integrated genetic and physical map of fiber development genes whose functions have been verified and confirmed. A total of 535 cotton fiber development genes, including 103 fiber transcription factors, 259 fiber development genes, and 173 SSR-contained fiber ESTs, were analyzed at the subgenome level. A total of 499 fiber related contigs were selected and assembled. Together these contigs covered about 151 Mb in physical length, or about 6.7% of the tetraploid cotton genome. Among the 499 contigs, 397 were anchored onto individual chromosomes. Results from our studies on the distribution patterns of the fiber development genes and transcription factors between the At and Dt subgenomes showed that more transcription factors were from Dt subgenome than At, whereas more fiber development genes were from At subgenome than Dt. Combining our mapping results with previous reports that more fiber QTLs were mapped in Dt subgenome than At subgenome, the results suggested a new functional hypothesis for tetraploid cotton. After the merging of the two diploid Gossypium genomes, the At subgenome has provided most of the genes for fiber development, because it continues to function similar to its fiber producing diploid A genome ancestor. On the other hand, the Dt subgenome, with its non-fiber producing D genome ancestor, provides more transcription factors that regulate the expression of the fiber genes in the At subgenome. This hypothesis would explain previously published mapping results. At the same time, this integrated map of fiber development genes would provide a framework to clone individual full-length fiber genes, to elucidate the physiological mechanisms of the fiber differentiation, elongation, and maturation, and to systematically study the functional network of these genes that interact during the process of fiber development in the tetraploid cottons.

A draft physical map of a D-genome cotton species (Gossypium raimondii)

BACKGROUND

Genetically anchored physical maps of large eukaryotic genomes have proven useful both for their intrinsic merit and as an adjunct to genome sequencing. Cultivated tetraploid cottons, Gossypium hirsutum and G. barbadense, share a common ancestor formed by a merger of the A and D genomes about 1-2 million years ago. Toward the long-term goal of characterizing the spectrum of diversity among cotton genomes, the worldwide cotton community has prioritized the D genome progenitor Gossypium raimondii for complete sequencing.

RESULTS

A whole genome physical map of G. raimondii, the putative D genome ancestral species of tetraploid cottons was assembled, integrating genetically-anchored overgo hybridization probes, agarose based fingerprints and 'high information content fingerprinting' (HICF). A total of 13,662 BAC-end sequences and 2,828 DNA probes were used in genetically anchoring 1585 contigs to a cotton consensus genetic map, and 370 and 438 contigs, respectively to Arabidopsis thaliana (AT) and Vitis vinifera (VV) whole genome sequences.

CONCLUSION

Several lines of evidence suggest that the G. raimondii genome is comprised of two qualitatively different components. Much of the gene rich component is aligned to the Arabidopsis and Vitis vinifera genomes and shows promise for utilizing translational genomic approaches in understanding this important genome and its resident genes. The integrated genetic-physical map is of value both in assembling and validating a planned reference sequence.

Genetic dissection of chromosome substitution lines of cotton to discover novel Gossypium barbadense L. alleles for improvement of agronomic traits

We recently released a set of 17 chromosome substitution (CS-B) lines (2n = 52) that contain Gossypium barbadense L. doubled-haploid line '3-79' germplasm systematically introgressed into the Upland inbred 'TM-1' of G. hirsutum (L.). TM-1 yields much more than 3-79, but cotton from the latter has superior fiber properties. To explore the use of these quasi-isogenic lines in studying gene interactions, we created a partial diallel among six CS-B lines and the inbred TM-1, and characterized their descendents for lint percentage, boll weight, seedcotton yield and lint yield across four environments. Phenotypic data on the traits were analyzed according to the ADAA genetic model to detect significant additive, dominance, and additive-by-additive epistasis effects at the chromosome and chromosome-by-chromosome levels of CS-B lines. For example, line 3-79 had the lowest boll weight, seedcotton yield and lint yield, but CS-B22Lo homozygous dominance genetic effects on seedcotton and lint yield were nearly four times those of TM-1, and its hybrids with TM-1 had the highest additive-by-additive epistatic effects on seedcotton and lint yield. CS-B14sh, 17, 22Lo and 25 produced positive homozygous dominance effects on lint yield, whereas doubly heterozygous combinations of CS-B14sh with CS-B17, 22Lo and 25 produced negative dominance effects, suggesting that epistatic effects between genes in these chromosomes strongly affect lint yield. The results underscore the opportunities to systematically identify genomic regions harboring genes that impart agronomically significant effects via epistatic interactions. The chromosome-by-chromosome approach significantly complements other strategies to detect and quantify epistatic interaction effects, and the quasi-isogenic nature of families and lines from CS-B intermatings will facilitate high-resolution localization, development of markers for selection and map-assisted identification of genes involved in strong epistatic effects.

The construction of a tetraploid cotton genome wide comprehensive reference map

Integration of multiple genomic maps provides a higher density of markers and greater genome coverage, which not only facilitates the identification and positioning of QTLs and candidate genes, but it also provides a basic structure for the genome sequence assembly. However, the diversity in markers and populations used in individual mapping studies limits the ability to fully integrate the available data. By concentrating on marker orders rather than marker distances, published map data could be used to produce a comprehensive reference map (CRM) that includes a majority of known markers with optimally estimated order of those markers across the genome. In this study, a tetraploid cotton genome-wide CRM was constructed from 28 public cotton genetic maps. The initial CRM contained 7,424 markers and represented over 93% of the combined mapping information from the 28 individual maps. The current output is stored and displayed through CottonDB (http://www.cottondb.org), the public cotton genome database.

Molecular mapping of Verticillium wilt resistance QTL clustered on chromosomes D7 and D9 in upland cotton

Verticillium wilt is a destructive disease with international consequences for cotton production. Breeding broad-spectrum resistant cultivars is considered to be one of the most effective means for reducing crop losses. A resistant cotton cultivar, 60182, was crossed with a susceptible cultivar, Junmian 1, to identify markers for Verticillium resistance genes and validate the mode of its inheritance. Genetic segregation analysis for Verticillium wilt resistance was evaluated based upon infected leaf percentage in the seedling stage using major gene-polygene mixed inheritance models and joint analysis of P(1), P(2), F(1), B(1), B(2) and F(2) populations obtained from the cultivar cross. We found that resistance of upland cotton cultivar 60182 to isolates BP2, VD8 and T9, and their isoconcentration mixture was controlled by two major genes with additive-dominance-epistatic effects, and the inheritance of the major gene was dominant. Furthermore, a genetic linkage map was constructed using F(2) segregating population and resistance phenotypic data were obtained using F(2:3) families inoculated with different isolates and detected in different developmental stages. The genetic linkage map with 139 loci was comprised of 31 linkage groups covering 1165 cM, with an average distance of 8.38 cM between two markers, or 25.89% of the cotton genome length. From 60182, we found 4 QTL on chromosome D7 and 4 QTL on D9 for BP2, 5 QTL on D7 and 9 QTL on D9 for VD8, 4 QTL on D7 and 5 QTL on D9 for T9 and 3 QTL on D7 and 7 QTL on D7 for mixed pathogens. The QTL mapping results revealed that QTL clusters with high contribution rates were screened simultaneously on chromosomes D9 and D7 by multiple interval mapping (CIM), whether from resistance phenotypic data from different developmental stages or for different isolates. The result is consistent with the genetic model of two major genes in 60182 and suggests broad-spectrum resistance to both defoliating isolates of V. dahliae and nondefoliating isolates. The markers associated with resistance QTL may facilitate the use of Verticillium wilt resistance genes in improving breeding programs for cotton.

Using three overlapped RILs to dissect genetically clustered QTL for fiber strength on Chro.D8 in Upland cotton

Fiber strength is an important trait among cotton fiber qualities due to ongoing changes in spinning technology. Major quantitative trait loci (QTL) for fiber quality enable molecular marker-assisted selection (MAS) to effectively improve fiber quality of cotton cultivars. We previously identified a major QTL for fiber strength derived from 7235 in Upland cotton. In the present study, in order to fine-map fiber strength QTL, we chose three recombinant inbred lines (RIL), 7TR-133, 7TR-132, and 7TR-214, developed from a cross between 7235 and TM-1 for backcrossing to TM-1 to develop three large mapping populations. Phenotypic data for fiber strength traits were collected in Nanjing (JES/NAU) and Xinjiang (BES/XJ) in 2006 and 2007. Three simple sequence repeat (SSR) genetic linkage maps on Chro.24(D8) were constructed using these three backcrossed populations. The SSR genetic maps were constructed using 907 individuals in (7TR-133 x TM-1)F(2) (Pop A), 670 in (7TR-132 x TM-1)F(2) (Pop B), and 940 in (7TR-214 x TM-1)F(2) (Pop C). The average distance between SSR loci was 0.62, 1.7, and 0.56 cM for the three maps. MapQTL 5 software detected five-clustered QTL (2.5 < LOD < 29.8) on Chro.D8 for fiber strength following analysis of three RIL backcrossed F(2)/F(2:3) progenies at JES/NAU and BES/XJ over 2 years. Five QTL for fiber strength exhibited a total phenotypic variance (PV) of 28.8-59.6%.

QTL mapping in A-genome diploid Asiatic cotton and their congruence analysis with AD-genome tetraploid cotton in genus Gossypium

Asiatic cotton (Gossypium arboreum L.) is an Old World cultivated cotton species. The sinense race was planted extensively in China. Due to the advances in spinning technology during the last century, the species was replaced by the New World allotetraploid cotton G. hirsutum L. Gossypium arboreum is still grown in India and Pakistan and also used as an elite in current cotton breeding programs. In addition, G. arboreum serves as a model for genomic research in Gossypium. In the present study, we generated an A-genome diploid cotton intraspecific genetic map including 264 SSR loci with three morphological markers mapped to 13 linkage groups. The map spans 2,508.71 cM with an average distance of 9.4 cM between adjacent loci. A population containing 176 F(2:3) families was used to perform quantitative trait loci (QTL) mapping for 17 phenotypes using Multiple QTL Model (MQM) of MapQTL ver 5.0. Overall, 108 QTLs were detected on 13 chromosomes. Thirty-one QTLs for yield and its components were detected in the F2 population. Forty-one QTLs for yield and its components were detected in the F(2:3) families with a total of 43 QTLs for fiber qualities. Two QTLs for seed cotton weight/plant and lint index and three QTLs for seed index were consistently detected both in F2 and F(2:3). Most QTLs for fiber qualities and yields were located at the same interval or neighboring intervals. These results indicated that the negative correlation between fiber qualities and yield traits may result from either pleiotropic effect of one gene or linkage effects of multiple closely linked genes.

QTL mapping of yield and yield components for elite hybrid derived-RILs in upland cotton

A population of 180 recombinant inbred lines (RILs) was developed by single seed descended from the cross of high yield Upland cotton (Gossypium hirsutum L.) varieties Zhongmiansuo12 (ZMS12) and 8891, the two parents of Xiangzamian2 (XZM2). A genetic linkage map consisting of 132 loci and covering 865.20 cM was constructed using the RIL population chiefly with SSR markers. Yield and yield components were investigated for RILs in three environments in China. The purpose of the present research was to analyze the relationship between yield and its components and to map QTL for yield and yield components in cotton. QTL were tagged with data sets from single environment (separate analysis) and a set of data from means of the three environments (joint analysis). A total of 34 QTL for yield and yield components were independently detected in three environments, whereas fifteen QTL were found in joint analysis. Notably, a stable lint percentage QTL qLP-A10-1 was detected both in joint analysis and in two environments of separate analysis, which might be of special value for marker-assisted selection. The QTL detected in the present study provide new information on improving yield and yield components. Results of path analysis showed that bolls/plant had the largest contribution to lint yield, which is consistent with the mid-parent heterosis value in F(1). Accordingly, in cotton breeding, bolls/plant can be considered first and other yield components measured as a whole to implement variety enhancement and hybrid selection of cotton.

QTL mapping of yield and fiber traits based on a four-way cross population in Gossypium hirsutum L

Four-way cross (4WC) involving four different inbred lines frequently appears in the cotton breeding programs. However, linkage analysis and quantitative trait loci (QTL) mapping with molecular markers in cotton has largely been applied to populations derived from a cross between two inbred lines, and few results of QTL dissection were conducted in a 4WC population. In this study, an attempt was made to construct a linkage map and identify QTL for yield and fiber quality traits in 4WC derived from four different inbred lines in Gossypium hirsutum L. A linkage map was constructed with 285 SSR loci and one morphological locus, covering 2113.3 cM, approximately 42% of the total recombination length of the cotton genome. A total of 31 QTL with 5.1-25.8% of the total phenotypic variance explained were detected. Twenty-four common QTL across environments showed high stability, and six QTL were environment-specific. Several genomic segments affecting multiple traits were identified. The advantage of QTL mapping using a 4WC were discussed. This study presents the first example of QTL mapping using a 4WC population in upland cotton. The results presented here will enhance the understanding of the genetic basis of yield and fiber quality traits and enable further marker-assisted selection in cultivar populations in upland cotton.

Gene-rich islands for fiber development in the cotton genome

Cotton fiber is an economically important seed trichome and the world's leading natural fiber used in the manufacture of textiles. As a step toward elucidating the genomic organization and distribution of gene networks responsible for cotton fiber development, we investigated the distribution of fiber genes in the cotton genome. Results revealed the presence of gene-rich islands for fiber genes with a biased distribution in the tetraploid cotton (Gossypium hirsutum L.) genome that was also linked to discrete fiber developmental stages based on expression profiles. There were 3 fiber gene-rich islands associated with fiber initiation on chromosome 5, 3 islands for the early to middle elongation stage on chromosome 10, 3 islands for the middle to late elongation stage on chromosome 14, and 1 island on chromosome 15 for secondary cell wall deposition, for a total of 10 fiber gene-rich islands. Clustering of functionally related gene clusters in the cotton genome displaying similar transcriptional regulation indicates an organizational hierarchy with significant implications for the genetic enhancement of particular fiber quality traits. The relationship between gene-island distribution and functional expression profiling suggests for the first time the existence of functional coupling gene clusters in the cotton genome.

A preliminary analysis of genome structure and composition in Gossypium hirsutum

Background

Upland cotton has the highest yield, and accounts for > 95% of world cotton production. Decoding upland cotton genomes will undoubtedly provide the ultimate reference and resource for structural, functional, and evolutionary studies of the species. Here, we employed GeneTrek and BAC tagging information approaches to predict the general composition and structure of the allotetraploid cotton genome.

Results

142 BAC sequences from Gossypium hirsutum cv. Maxxa were downloaded and confirmed. These BAC sequence analysis revealed that the tetraploid cotton genome contains over 70,000 candidate genes with duplicated gene copies in homoeologous A- and D-subgenome regions. Gene distribution is uneven, with gene-rich and gene-free regions of the genome. Twenty-one percent of the 142 BACs lacked genes. BAC gene density ranged from 0 to 33.2 per 100 kb, whereas most gene islands contained only one gene with an average of 1.5 genes per island. Retro-elements were found to be a major component, first an enriched LTR/gypsy and second LTR/copia. Most LTR retrotransposons were truncated and in nested structures. In addition, 166 polymorphic loci amplified with SSRs developed from 70 BAC clones were tagged on our backbone genetic map. Seventy-five percent (125/166) of the polymorphic loci were tagged on the D-subgenome. By comprehensively analyzing the molecular size of amplified products among tetraploid G. hirsutum cv. Maxxa, acc. TM-1, and G. barbadense cv. Hai7124, and diploid G. herbaceum var. africanum and G. raimondii, 37 BACs, 12 from the A- and 25 from the D-subgenome, were further anchored to their corresponding subgenome chromosomes. After a large amount of genes sequence comparison from different subgenome BACs, the result showed that introns might have no contribution to different subgenome size in Gossypium.

Conclusion

This study provides us with the first glimpse of cotton genome complexity and serves as a foundation for tetraploid cotton whole genomesequencing in the future.

Simple sequence repeat genetic linkage maps of A-genome diploid cotton (Gossypium arboreum)

This study introduces the construction of the first intraspecific genetic linkage map of the A-genome diploid cotton with newly developed simple sequence repeat (SSR) markers using 189 F(2) plants derived from the cross of two Asiatic cotton cultivars (Gossypium arboreum L.) Jianglingzhongmian x Zhejiangxiaoshanlüshu. Polymorphisms between the two parents were detected using 6 092 pairs of SSR primers. Two-hundred and sixty-eight pairs of SSR primers with better polymorphisms were picked out to analyze the F(2) population. In total, 320 polymorphic bands were generated and used to construct a linkage map with JoinMap3.0. Two-hundred and sixty-seven loci, including three phenotypic traits were mapped at a logarithms of odds ratio (LOD) > or = 3.0 on 13 linkage groups. The total length of the map was 2 508.71 cM, and the average distance between adjacent markers was 9.40 cM. Chromosome assignments were according to the association of linkages with our backbone tetraploid specific map using the 89 similar SSR loci. Comparisons among the 13 suites of orthologous linkage groups revealed that the A-genome chromosomes are largely collinear with the A(t) and D(t) sub-genome chromosomes. Chromosomes associated with inversions suggested that allopolyploidization was accompanied by homologous chromosomal rearrangement. The inter-chromosomal duplicated loci supply molecular evidence that the A-genome diploid Asiatic cotton is paleopolyploid.