A new interspecific, Gossypium hirsutum x G. barbadense, RIL population: towards a unified consensus linkage map of tetraploid cotton

Deciphering the dynamic expression network of fiber elongation and the functional role of the GhTUB5 gene for fiber length in cotton based on an introgression population of upland cotton

INTRODUCTION

Interspecific introgression between Gossypium hirsutum and G. barbadense allows breeding cotton with outstanding fiber length (FL). However, the dynamic gene regulatory network of FL-related genes has not been characterized, and the functional mechanism through which the hub gene GhTUB5 mediates fiber elongation has yet to be determined.

METHODS

Coexpression analyses of 277 developing fiber transcriptomes integrated with QTL mapping using 250 introgression lines of different FL phenotypes were conducted to identify genes related to fiber elongation. The function of GhTUB5 was determined by ectopic expression of two TUB5 alleles in Arabidopsis and knockout of GhTUB5 in upland cotton. Yeast two-hybrid, split-luciferase and pull-down assays were conducted to screen for interacting proteins, and upstream genes were identified by yeast one-hybrid, dual-LUC and electrophoretic mobility shift assays.

RESULTS

The 32,612, 30,837 and 30,277 genes expressed at 5, 10 and 15 days postanthesis (dpa) were grouped into 19 distinct coexpression modules, and 988 genes in the MEblack module were enriched in the cell wall process and exhibited significant associations with FL. A total of 20 FL-QTLs were identified, each explaining 3.34-16.04 % of the phenotypic variance in the FL. Furthermore, several FL-QTLs contained 15 genes that were differentially expressed in the MEblack module including the tubulin beta gene (TUB5). Compared with the wild type, the overexpression of GhTUB5 and GbTUB5 in Arabidopsis suppressed root cell length but promoted cellulose synthesis. Knockout of GhTUB5 resulted in longer fiber lines. Protein-based experiments revealed that GhTUB5 interacts with GhZFP6. Additionally, GhTUB5 was directly activated by GhHD-ZIP7, a homeobox-leucine zipper transcription factor, and its paralogous gene was previously reported to mediate fiber elongation.

CONCLUSION

This study opens a new avenue to dissect functional mechanism of cotton fiber elongation. Our findings provide some molecular details on how GhTUB5 mediates the FL phenotype in cotton.

Unravelling the genetic basis and regulation networks related to fibre quality improvement using chromosome segment substitution lines in cotton

The elucidation of genetic architecture and molecular regulatory networks underlying complex traits remains a significant challenge in life science, largely due to the substantial background effects that arise from epistasis and gene-environment interactions. The chromosome segment substitution line (CSSL) is an ideal material for genetic and molecular dissection of complex traits due to its near-isogenic properties; yet a comprehensive analysis, from the basic identification of substitution segments to advanced regulatory network, is still insufficient. Here, we developed two cotton CSSL populations on the Gossypium hirsutum background, representing wide adaptation and high lint yield, with introgression from G. barbadense, representing superior fibre quality. We sequenced 99 CSSLs that demonstrated significant differences from G. hirsutum in fibre, and characterized 836 dynamic fibre transcriptomes in three crucial developmental stages. We developed a workflow for precise resolution of chromosomal substitution segments; the genome sequencing revealed substitutions collectively representing 87.25% of the G. barbadense genome. Together, the genomic and transcriptomic survey identified 18 novel fibre-quality-related quantitative trait loci with high genetic contributions and the comprehensive landscape of fibre development regulation. Furthermore, analysis determined unique cis-expression patterns in CSSLs to be the driving force for fibre quality alteration; building upon this, the co-expression regulatory network revealed biological relationships among the noted pathways and accurately described the molecular interactions of GhHOX3, GhRDL1 and GhEXPA1 during fibre elongation, along with reliable predictions for their interactions with GhTBA8A5. Our study will enhance more strategic employment of CSSL in crop molecular biology and breeding programmes.

Genetic analysis and molecular validation of gene conferring petal spot phenotype in interspecific crosses of cotton

Cotton (Gossypium species) has received considerable interest from the geneticists, cytologists and evolutionary biologists since the last more than a century. Here, we explore the genetics of petal spot in the interspecific derivatives involving tetraploid and diploid cottons; and confirm the location of gene governing petal spot phenotype on chromosome A7 by demonstrating co-segregation of SSR marker NAU 2186 with petal spot phenotype. The presence of petal spot was observed to be dominant over its absence. Petal spot inheritance showed significant deviation from the expected Mendelian ratio in all the segregating populations indicating segregation distortion. The distortion was biased towards the hirsutum parent which has important implications from introgression point of view. We also report a strong association between petal spot and petal margin coloration phenotypes. Extant American cotton varieties generally lack petal spot and margin coloration phenotypes. These petal characteristics can serve as morphological markers during germplasm characterization.

Aspartyl proteases identified as candidate genes of a fiber length QTL, qFLD05, that regulates fiber length in cotton (Gossypium hirsutum L.)

KEY MESSAGE

GhAP genes were identified as the candidates involved in cotton fiber length under the scope of fine mapping a stable fiber length QTL, qFLD05. Moreover, the transcription factor GhWRKY40 positively regulated GhAP3 to decrease fiber length. Fiber length (FL) is an economically important fiber quality trait. Although several genes controlling cotton fiber development have been identified, our understanding of this process remains limited. In this study, an FL QTL (qFLD05) was fine-mapped to a 216.9-kb interval using a secondary F2:3 population derived from the upland hybrid cultivar Ji1518. This mapped genomic segment included 15 coding genes, four of which were annotated as aspartyl proteases (GhAP1-GhAP4). GhAPs were identified as candidates for qFLD05 as the sequence variations in GhAPs were associated with FL deviations in the mapping population, and functional validation of GhAP3 and GhAP4 indicated a longer FL following decreases in their expression levels through virus-induced gene silencing (VIGS). Subsequently, the potential involvement of GhWRKY40 in the regulatory network was revealed: GhWRKY40 positively regulated GhAP3's expression according to transcriptional profiling, VIGS, yeast one-hybrid assays and dual-luciferase experiments. Furthermore, alterations in the expression of the eight previously reported cotton FL-responsive genes from the above three VIGS lines (GhAP3, GhAP4 and GhWRKY40) implied that MYB5_A12 was involved in the GhWRKY40-GhAP network. In short, we unveiled the unprecedented FL regulation roles of GhAPs in cotton, which was possibly further regulated by GhWRKY40. These findings will reveal the genetic basis of FL development associated with qFLD05 and be beneficial for the marker-assisted selection of long-staple cotton.

Transcriptomic datasets of Verticillium wilt resistant and non-resistant Gossypium barbadense varieties during pathogen inoculation

Cotton is a significant cash crop and the primary source of natural fiber globally. Among the numerous diseases encountered in cotton production, Verticillium wilt is one of the most serious, caused by the pathogen Verticillium dahliae (V. dahliae). Unfortunately, there are no effective targeted methods to combat this disease. Genomic resources for Verticillium wilt resistance primarily exist in Gossypium barbadense (G. barbadense). Regrettably, there have been limited transcriptomic comparisons between V. dahliae-resistant and -susceptible varieties of G. barbadense due to the scarcity of susceptible resources. In this study, we conducted a transcriptome analysis on both V. dahliae-resistant and -susceptible varieties of G. barbadense at the 0, 12, 24 and 48 hours after V. dahliae inoculation. This comparative transcriptome analysis yielded high-quality data and offered new insights into the molecular mechanisms underlying cotton's resistance against this destructive pathogen.

Unraveling genomic regions and candidate genes for multiple disease resistance in upland cotton using meta-QTL analysis

Verticillium wilt (VW), Fusarium wilt (FW) and Root-knot nematode (RKN) are the main diseases affecting cotton production. However, many reported quantitative trait loci (QTLs) for cotton resistance have not been used for agricultural practices because of inconsistencies in the cotton genetic background. The integration of existing cotton genetic resources can facilitate the discovery of important genomic regions and candidate genes involved in disease resistance. Here, an improved and comprehensive meta-QTL analysis was conducted on 487 disease resistant QTLs from 31 studies in the last two decades. A consensus linkage map with genetic overall length of 3006.59 cM containing 8650 markers was constructed. A total of 28 Meta-QTLs (MQTLs) were discovered, among which nine MQTLs were identified as related to resistance to multiple diseases. Candidate genes were predicted based on public transcriptome data and enriched in pathways related to disease resistance. This study used a method based on the integration of Meta-QTL, known genes and transcriptomics to reveal major genomic regions and putative candidate genes for resistance to multiple diseases, providing a new basis for marker-assisted selection of high disease resistance in cotton breeding.

Expressed genes and their new alleles identification during fibre elongation reveal the genetic factors underlying improvements of fibre length in cotton

Interspecific breeding in cotton takes advantage of genetic recombination among desirable genes from different parental lines. However, the expression new alleles (ENAs) from crossovers within genic regions and their significance in fibre length (FL) improvement are currently not understood. Here, we generated resequencing genomes of 191 interspecific backcross inbred lines derived from CRI36 (Gossypium hirsutum) × Hai7124 (Gossypium barbadense) and 277 dynamic fibre transcriptomes to identify the ENAs and extremely expressed genes (eGenes) potentially influencing FL, and uncovered the dynamic regulatory network of fibre elongation. Of 35 420 eGenes in developing fibres, 10 366 ENAs were identified and preferentially distributed in chromosomes subtelomeric regions. In total, 1056-1255 ENAs showed transgressive expression in fibres at 5-15 dpa (days post-anthesis) of some BILs, 520 of which were located in FL-quantitative trait locus (QTLs) and GhFLA9 (recombination allele) was identified with a larger effect for FL than GhFLA9 of CRI36 allele. Using ENAs as a type of markers, we identified three novel FL-QTLs. Additionally, 456 extremely eGenes were identified that were preferentially distributed in recombination hotspots. Importantly, 34 of them were significantly associated with FL. Gene expression quantitative trait locus analysis identified 1286, 1089 and 1059 eGenes that were colocalized with the FL trait at 5, 10 and 15 dpa, respectively. Finally, we verified the Ghir_D10G011050 gene linked to fibre elongation by the CRISPR-cas9 system. This study provides the first glimpse into the occurrence, distribution and expression of the developing fibres genes (especially ENAs) in an introgression population, and their possible biological significance in FL.

Linkage and association analyses reveal that hub genes in energy-flow and lipid biosynthesis pathways form a cluster in upland cotton

Upland cotton is an important allotetraploid crop that provides both natural fiber for the textile industry and edible vegetable oil for the food or feed industry. To better understand the genetic mechanism that regulates the biosynthesis of storage oil in cottonseed, we identified the genes harbored in the major quantitative trait loci/nucleotides (QTLs/QTNs) of kernel oil content (KOC) in cottonseed via both multiple linkage analyses and genome-wide association studies (GWAS). In ‘CCRI70′ RILs, six stable QTLs were simultaneously identified by linkage analysis of CHIP and SLAF-seq strategies. In ‘0-153′ RILs, eight stable QTLs were detected by consensus linkage analysis integrating multiple strategies. In the natural panel, thirteen and eight loci were associated across multiple environments with two algorithms of GWAS. Within the confidence interval of a major common QTL on chromosome 3, six genes were identified as participating in the interaction network highly correlated with cottonseed KOC. Further observations of gene differential expression showed that four of the genes, LtnD, PGK, LPLAT1, and PAH2, formed hub genes and two of them, FER and RAV1, formed the key genes in the interaction network. Sequence variations in the coding regions of LtnD, FER, PGK, LPLAT1, and PAH2 genes may support their regulatory effects on oil accumulation in mature cottonseed. Taken together, clustering of the hub genes in the lipid biosynthesis interaction network provides new insights to understanding the mechanism of fatty acid biosynthesis and TAG assembly and to further genetic improvement projects for the KOC in cottonseeds.

Features of Chromosome Introgression from Gossypium barbadense L. into G. hirsutum L. during the Development of Alien Substitution Lines

The creation of G. barbadense L./G. hirsutum L. chromosome-substitution lines is an important method to transfer agronomically valuable traits from G. barbadense into G. hirsutum. In this study, 30 monosomic lines of G. hirsutum from the Cytogenetic Collection of Uzbekistan, created in the genotypic background of line L-458, were used in crosses with the G. barbadense line Pima 3-79 to create substitution lines. In the course of this work, new monosomic lines were identified for chromosome 12 and monotelodisome 6 of the A_tsubgenome and for chromosomes 17, 21, and 22 of the D_tsubgenome using chromosome-specific SSR markers and a well-defined tester set of cotton translocation lines (USA). Compared to those in the F₁ hybrids, a strong decrease in the crossing and setting rates was found in the BC₁F₁ backcross lines, with the substitution of chromosomes 2, 4, 6, 7, and 12 of the A_tsubgenome and 17, 18, 21, and 22 of the D_tsubgenome. The F₁ and BC₁F₁ offspring from interspecific crosses differed in their transmission of univalents. Despite the regular pairing of chromosomes and the high meiotic index, interspecific aneuploid hybrids were characterized by a decrease in pollen fertility, which may indicate hidden structural variability in these genomes that did not affect meiotic division. The identification of chromosomes using chromosome-specific SSR markers in the early stages of plant development has greatly accelerated the detection of monosomic plants. The analysis of morphobiological traits revealed that monosomic F₁ hybrids were more similar to the donor line, while BC₁F₁ hybrids were more similar to the recurrent parent but also showed previously undetected traits.

Examining two sets of introgression lines across multiple environments reveals background-independent and stably expressed quantitative trait loci of fiber quality in cotton

Background-independent (BI) and stably expressed (SE) quantitative trait loci (QTLs) were identified using two sets of introgression lines across multiple environments. Genetic background more greatly affected fiber quality traits than environmental factors. Sixty-one SE-QTLs, including two BI-QTLs, were novel and 48 SE-QTLs, including seven BI-QTLs, were previously reported. Cotton fiber quality traits are controlled by QTLs and are susceptible to environmental influence. Fiber quality improvement is an essential goal in cotton breeding but is hindered by limited knowledge of the genetic basis of fiber quality traits. In this study, two sets of introgression lines of Gossypium hirsutum × G. barbadense were used to dissect the QTL stability of three fiber quality traits (fiber length, strength and micronaire) across environments using 551 simple sequence repeat markers selected from our high-density genetic map. A total of 76 and 120 QTLs were detected in the CCRI36 and CCRI45 backgrounds, respectively. Nine BI-QTLs were found, and 78 (41.71%) of the detected QTLs were reported previously. Thirty-nine and 79 QTLs were SE-QTLs in at least two environments in the CCRI36 and CCRI45 backgrounds, respectively. Forty-eight SE-QTLs, including seven BI-QTLs, were confirmed in previous reports, and 61 SE-QTLs, including two BI-QTLs, were considered novel. These results indicate that genetic background more strongly impacts on fiber quality traits than environmental factors. Twenty-three clusters with BI- and/or SE-QTLs were identified, 19 of which harbored favorable alleles from G. barbadense for two or three fiber quality traits. This study is the first report using two sets of introgression lines to identify fiber quality QTLs across environments in cotton, providing insights into the effect of genetic backgrounds and environments on the QTL expression of fiber quality and important information for the genetic basis underlying fiber quality traits toward QTL cloning and molecular breeding.

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.

GWAS Analysis and QTL Identification of Fiber Quality Traits and Yield Components in Upland Cotton Using Enriched High-Density SNP Markers

It is of great importance to identify quantitative trait loci (QTL) controlling fiber quality traits and yield components for future marker-assisted selection (MAS) and candidate gene function identifications. In this study, two kinds of traits in 231 F6:8 recombinant inbred lines (RILs), derived from an intraspecific cross between Xinluzao24, a cultivar with elite fiber quality, and Lumianyan28, a cultivar with wide adaptability and high yield potential, were measured in nine environments. This RIL population was genotyped by 122 SSR and 4729 SNP markers, which were also used to construct the genetic map. The map covered 2477.99 cM of hirsutum genome, with an average marker interval of 0.51 cM between adjacent markers. As a result, a total of 134 QTLs for fiber quality traits and 122 QTLs for yield components were detected, with 2.18-24.45 and 1.68-28.27% proportions of the phenotypic variance explained by each QTL, respectively. Among these QTLs, 57 were detected in at least two environments, named stable QTLs. A total of 209 and 139 quantitative trait nucleotides (QTNs) were associated with fiber quality traits and yield components by four multilocus genome-wide association studies methods, respectively. Among these QTNs, 74 were detected by at least two algorithms or in two environments. The candidate genes harbored by 57 stable QTLs were compared with the ones associated with QTN, and 35 common candidate genes were found. Among these common candidate genes, four were possibly "pleiotropic." This study provided important information for MAS and candidate gene functional studies.

Identification of Introgressed Alleles Conferring High Fiber Quality Derived From Gossypium barbadense L. in Secondary Mapping Populations of G. hirsutum L

The improvement of fiber quality is an essential goal in cotton breeding. In our previous studies, several quantitative trait loci (QTLs) contributing to improved fiber quality were identified in different introgressed chromosomal regions from Sea Island cotton (Gossypium barbadense L.) in a primary introgression population (Pop. A) of upland cotton (G. hirsutum L.). In the present study, to finely map introgressed major QTLs and accurately dissect the genetic contribution of the target introgressed chromosomal segments, we backcrossed two selected recombinant inbred lines (RILs) that presented desirable high fiber quality with their high lint-yielding recurrent parent to ultimately develop two secondary mapping populations (Pop. B and Pop. C). Totals of 20 and 27 QTLs for fiber quality were detected in Pop. B and Pop. C, respectively, including four and five for fiber length, four and eight for fiber micronaire, two and four for fiber uniformity, five and four for fiber elongation, and six and four for fiber strength, respectively. Two QTLs for lint percentage were detected only in Pop. C. In addition, seven stable QTLs were identified, including two for both fiber length and fiber strength and three for fiber elongation. Five QTL clusters for fiber quality were identified in the introgressed chromosomal regions, and negative effects of these chromosomal regions on lint percentage (a major lint yield parameter) were not observed. Candidate genes with a QTL-cluster associated with fiber strength and fiber length in the introgressed region of Chr.7 were further identified. The results may be helpful for revealing the genetic basis of superior fiber quality contributed by introgressed alleles from G. barbadense. Possible strategies involving marker-assisted selection (MAS) for simultaneously improving upland cotton fiber quality and lint yield in breeding programs was also discussed.

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.

Chromosome structural variation of two cultivated tetraploid cottons and their ancestral diploid species based on a new high-density genetic map

A high-resolution genetic map is a useful tool for assaying genomic structural variation and clarifying the evolution of polyploid cotton. A total of 36956 SSRs, including 11289 released in previous studies and 25567 which were newly developed based on the genome sequences of G. arboreum and G. raimondii, were utilized to construct a new genetic map. The new high-density genetic map includes 6009 loci and spanned 3863.97 cM with an average distance of 0.64 cM between consecutive markers. Four inversions (one between Chr08 and Chr24, one between Chr09 and Chr23 and two between Chr10 and Chr20) were identified by homology analysis. Comparative genomic analysis between genetic map and two diploid cottons showed that structural variations between the A genome and At subgenome are more extensive than between D genome and Dt subgenome. A total of 17 inversions, seven simple translocations and two reciprocal translocations were identified between genetic map and G. raimondii. Good colinearity was revealed between the corresponding chromosomes of tetraploid G. hirsutum and G. barbadense genomes, but a total of 16 inversions were detected between them. These results will accelerate the process of evolution analysis of Gossipium genus.

A meta-analysis of quantitative trait loci for abiotic and biotic stress resistance in tetraploid cotton

The number and location of mapped quantitative trait loci (QTL) depend on genetic populations and testing environments. The identification of consistent QTL across genetic backgrounds and environments is a pre-requisite to marker-assisted selection. This study analyzed a total of 661 abiotic and biotic stress resistance QTL based on our previous work and other publications using the meta-analysis software Biomercator. It identified chromosomal regions containing QTL clusters for different resistance traits and hotspots for a particular resistance trait in cotton from 98 QTL for drought tolerance under greenhouse (DT) and 150 QTL in field conditions (FDT), 80 QTL for salt tolerance in the greenhouse conditions (ST), 201 QTL for resistance to Verticillium wilt (VW, Verticillium dahliae), 47 QTL for resistance to Fusarium wilt (FW, Fusarium oxysporum f. sp. vasinfectum), and 85 QTL for resistance to root-knot nematodes (RKN, Meloiodogyne incognita) and reniform nematodes (RN, Rotylenchulus reniformis). The traits used in QTL mapping for abiotic stress tolerance included morphological traits-plant height and fresh and dry shoot and root weights, physiological traits-chlorophyll content, osmotic potential, carbon isotope ratio, stomatal conductance, photosynthetic rate, transpiration, canopy temperature, and leaf area index, agronomic traits-seedcotton yield, lint yield, boll weight, and lint percent, and fiber quality traits-fiber length, uniformity, strength, elongation, and micronaire. The results showed that resistance QTL are not uniformly distributed across the cotton genome; some chromosomes carried disproportionally more QTL, QTL clusters, or hotspots. Twenty-three QTL clusters were found on 15 chromosomes (c3, c4, c5, c6, c7, c11, c14, c15, c16, c19, c20, c23, c24, c25, and c26). Moreover, 28 QTL hotshots were associated with different resistance traits including one hotspot on c4 for Verticillium wilt resistance, two QTL hotspots on c24 for chlorophyll content measured under both drought and salt stress conditions, and three other hotspots on c19 for the resistance to Verticillium wilt and Fusarium wilt, and micronaire under drought stress conditions. This meta-analysis of stress tolerance QTL provides an important foundation for cotton breeding and further studies on the genetic mechanisms of abiotic and biotic stress resistance in cotton.

Fine mapping and candidate gene analysis of qFL-chr1, a fiber length QTL in cotton

A fiber length QTL, qFL-chr1, was fine mapped to a 0.9 cM interval of cotton chromosome 1. Two positional candidate genes showed positive correlation between gene expression level and fiber length. Prior analysis of a backcross-self mapping population derived from a cross between Gossypium hirsutum L. and G. barbadense L. revealed a QTL on chromosome 1 associated with increased fiber length (qFL-chr1), which was confirmed in three independent populations of near-isogenic introgression lines (NIILs). Here, a single NIIL, R01-40-08, was used to develop a large population segregating for the target region. Twenty-two PCR-based polymorphic markers used to genotype 1672 BC₄F₂ plants identified 432 recombinants containing breakpoints in the target region. Substitution mapping using 141 informative recombinants narrowed the position of qFL-chr1 to a 1.0-cM interval between SSR markers MUSS084 and CIR018. To exclude possible effects of non-target introgressions on fiber length, different heterozygous BC₄F₃ plants introgressed between SSR markers NAU3384 and CGR5144 were selected to develop sub-NILs. The qFL-chr1 was further mapped at 0.9-cM interval between MUSS422 and CIR018 by comparisons of sub-NIL phenotype, and increased fiber length by ~1 mm. The 2.38-Mb region between MUSS422 and CIR018 in G. barbadense contained 19 annotated genes. Expression levels of two of these genes, GOBAR07705 (encoding 1-aminocyclopropane-1-carboxylate synthase) and GOBAR25992 (encoding amino acid permease), were positively correlated with fiber length in a small F₂ population, supporting these genes as candidates for qFL-chr1.

High-density linkage map construction and QTL analysis for earliness-related traits in Gossypium hirsutum L

Background

Gossypium hirsutum L., or upland cotton, is an important renewable resource for textile fiber. To enhance understanding of the genetic basis of cotton earliness, we constructed an intra-specific recombinant inbred line population (RIL) containing 137 lines, and performed linkage map construction and quantitative trait locus (QTL) mapping.

Results

Using restriction-site associated DNA sequencing, a genetic map composed of 6,434 loci, including 6,295 single nucleotide polymorphisms and 139 simple sequence repeat loci, was developed from RIL population. This map spanned 4,071.98 cM, with an average distance of 0.63 cM between adjacent markers. A total of 247 QTLs for six earliness-related traits were detected in 6 consecutive years. In addition, 55 QTL coincidence regions representing more than 60 % of total QTLs were found on 22 chromosomes, which indicated that several earliness-related traits might be simultaneously improved. Fine-mapping of a 2-Mb region on chromosome D3 associated with five stable QTLs between Marker25958 and Marker25963 revealed that lines containing alleles derived from CCRI36 in this region exhibited smaller phenotypes and earlier maturity. One candidate gene (EMF2) was predicted and validated by quantitative real-time PCR in early-, medium- and late-maturing cultivars from 3- to 6-leaf stages, with highest expression level in early-maturing cultivar, CCRI74, lowest expression level in late-maturing cultivar, Bomian1.

Conclusions

We developed an SNP-based genetic map, and this map is the first high-density genetic map for short-season cotton and has the potential to provide deeper insights into earliness. Cotton earliness-related QTLs and QTL coincidence regions will provide useful materials for QTL fine mapping, gene positional cloning and MAS. And the gene, EMF2, is promising for further study.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3269-y) contains supplementary material, which is available to authorized users.

Mapping of fiber quality QTLs reveals useful variation and footprints of cotton domestication using introgression lines

Fiber quality improvement is a driving force for further cotton domestication and breeding. Here, QTLs for fiber quality were mapped in 115 introgression lines (ILs) first developed from two intraspecific populations of cultivated and feral cotton landraces. A total of 60 QTLs were found, which explained 2.03–16.85% of the phenotypic variance found in fiber quality traits. A total of 36 markers were associated with five fiber traits, 33 of which were found to be associated with QTLs in multiple environments. In addition, nine pairs of common QTLs were identified; namely, one pair of QTLs for fiber elongation, three pairs for fiber length, three pairs for fiber strength and two pairs for micronaire (qMICs). All common QTLs had additive effects in the same direction in both IL populations. We also found five QTL clusters, allowing cotton breeders to focus their efforts on regions of QTLs with the highest percentages of phenotypic variance. Our results also reveal footprints of domestication; for example, fourteen QTLs with positive effects were found to have remained in modern cultivars during domestication, and two negative qMICs that had never been reported before were found, suggesting that the qMICs regions may be eliminated during artificial selection.

A Genetic Map Between Gossypium hirsutum and the Brazilian Endemic G. mustelinum and Its Application to QTL Mapping

Among the seven tetraploid cotton species, little is known about transmission genetics and genome organization in Gossypium mustelinum, the species most distant from the source of most cultivated cotton, G. hirsutum. In this research, an F₂ population was developed from an interspecific cross between G. hirsutum and G. mustelinum (HM). A genetic linkage map was constructed mainly using simple sequence repeat (SSRs) and restriction fragment length polymorphism (RFLP) DNA markers. The arrangements of most genetic loci along the HM chromosomes were identical to those of other tetraploid cotton species. However, both major and minor structural rearrangements were also observed, for which we propose a parsimony-based model for structural divergence of tetraploid cottons from common ancestors. Sequences of mapped markers were used for alignment with the 26 scaffolds of the G. hirsutum draft genome, and showed high consistency. Quantitative trait locus (QTL) mapping of fiber elongation in advanced backcross populations derived from the same parents demonstrated the value of the HM map. The HM map will serve as a valuable resource for QTL mapping and introgression of G. mustelinum alleles into G. hirsutum, and help clarify evolutionary relationships between the tetraploid cotton genomes.

Genome Wide SSR High Density Genetic Map Construction from an Interspecific Cross of Gossypium hirsutum × Gossypium tomentosum

A high density genetic map was constructed using F2 population derived from an interspecific cross of G. hirsutum × G. tomentosum. The map consisted of 3093 marker loci distributed across all the 26 chromosomes and covered 4365.3 cM of cotton genome with an average inter-marker distance of 1.48 cM. The maximum length of chromosome was 218.38 cM and the minimum was 122.09 cM with an average length of 167.90 cM. A sub-genome covers more genetic distance (2189.01 cM) with an average inter loci distance of 1.53 cM than D sub-genome which covers a length of 2176.29 cM with an average distance of 1.43 cM. There were 716 distorted loci in the map accounting for 23.14% and most distorted loci were distributed on D sub-genome (25.06%), which were more than on A sub-genome (21.23%). In our map 49 segregation hotspots (SDR) were distributed across the genome with more on D sub-genome as compared to A genome. Two post-polyploidization reciprocal translocations of "A2/A3 and A4/A5" were suggested by seven pairs of duplicate loci. The map constructed through these studies is one of the three densest genetic maps in cotton however; this is the first dense genome wide SSR interspecific genetic map between G. hirsutum and G. tomentosum.

Natural Allelic Variations in Highly Polyploidy Saccharum Complex

Sugarcane (Saccharum spp.) is an important sugar and biofuel crop with high polyploid and complex genomes. The Saccharum complex, comprised of Saccharum genus and a few related genera, are important genetic resources for sugarcane breeding. A large amount of natural variation exists within the Saccharum complex. Though understanding their allelic variation has been challenging, it is critical to dissect allelic structure and to identify the alleles controlling important traits in sugarcane. To characterize natural variations in Saccharum complex, a target enrichment sequencing approach was used to assay 12 representative germplasm accessions. In total, 55,946 highly efficient probes were designed based on the sorghum genome and sugarcane unigene set targeting a total of 6 Mb of the sugarcane genome. A pipeline specifically tailored for polyploid sequence variants and genotype calling was established. BWA-mem and sorghum genome approved to be an acceptable aligner and reference for sugarcane target enrichment sequence analysis, respectively. Genetic variations including 1,166,066 non-redundant SNPs, 150,421 InDels, 919 gene copy number variations, and 1,257 gene presence/absence variations were detected. SNPs from three different callers (Samtools, Freebayes, and GATK) were compared and the validation rates were nearly 90%. Based on the SNP loci of each accession and their ploidy levels, 999,258 single dosage SNPs were identified and most loci were estimated as largely homozygotes. An average of 34,397 haplotype blocks for each accession was inferred. The highest divergence time among the Saccharum spp. was estimated as 1.2 million years ago (MYA). Saccharum spp. diverged from Erianthus and Sorghum approximately 5 and 6 MYA, respectively. The target enrichment sequencing approach provided an effective way to discover and catalog natural allelic variation in highly polyploid or heterozygous genomes.

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.

A high-density SSR genetic map constructed from a F2 population of Gossypium hirsutum and Gossypium darwinii

The cultivated allotetraploid species Gossypium hirsutum, accounts for 90% of the world cotton production, has narrow genetic basis that's why its yield, quality or stress resistance breeding is stagnant. It is therefore, essential to explore desirable genes from Gossypium darwinii which has enviable traits such as high fiber fineness, drought tolerance, fusarium and verticillium resistance. We used G. darwinii as primary plant materials in this study not only to enrich the genetic diversity of exiting germplasm but also to better understand its genome structure. An interspecific high density linkage map of allotetraploid cotton was constructed using F2 population (G. hirsutum×G. darwinii). The map was based entirely on genome-wide simple sequence repeat (SSR) markers. A total of 2763 markers were mapped in 26 linkage groups (chromosomes) covering a genome length of 4176.7cM with an average inter-locus distance of 1.5cM. The length of the chromosomes ranged from 84.7 to 238.5cM with an average length of 160.6cM. At subgenome length was 2160.7cM with an average distance of 1.6cM, where as Dt genome length was 2016cM with an average distance of 1.4cM. There were 601 distorted SSR loci. Less number of segregation distortion loci were located in At subgenome than in Dt subgenome. Two post-polyploidization reciprocal translocations of "A2/A3 and A4/A5" were suggested by 44 pairs of duplicate loci.

Genetic analysis of Verticillium wilt resistance in a backcross inbred line population and a meta-analysis of quantitative trait loci for disease resistance in cotton

Background

Verticillium wilt (VW) and Fusarium wilt (FW), caused by the soil-borne fungi Verticillium dahliae and Fusarium oxysporum f. sp. vasinfectum, respectively, are two most destructive diseases in cotton production worldwide. Root-knot nematodes (Meloidogyne incognita, RKN) and reniform nematodes (Rotylenchulus reniformis, RN) cause the highest yield loss in the U.S. Planting disease resistant cultivars is the most cost effective control method. Numerous studies have reported mapping of quantitative trait loci (QTLs) for disease resistance in cotton; however, very few reliable QTLs were identified for use in genomic research and breeding.

Results

This study first performed a 4-year replicated test of a backcross inbred line (BIL) population for VW resistance, and 10 resistance QTLs were mapped based on a 2895 cM linkage map with 392 SSR markers. The 10 VW QTLs were then placed to a consensus linkage map with other 182 VW QTLs, 75 RKN QTLs, 27 FW QTLs, and 7 RN QTLs reported from 32 publications. A meta-analysis of QTLs identified 28 QTL clusters including 13, 8 and 3 QTL hotspots for resistance to VW, RKN and FW, respectively. The number of QTLs and QTL clusters on chromosomes especially in the A-subgenome was significantly correlated with the number of nucleotide-binding site (NBS) genes, and the distribution of QTLs between homeologous A- and D- subgenome chromosomes was also significantly correlated.

Conclusions

Ten VW resistance QTL identified in a 4-year replicated study have added useful information to the understanding of the genetic basis of VW resistance in cotton. Twenty-eight disease resistance QTL clusters and 24 hotspots identified from a total of 306 QTLs and linked SSR markers provide important information for marker-assisted selection and high resolution mapping of resistance QTLs and genes. The non-overlapping of most resistance QTL hotspots for different diseases indicates that their resistances are controlled by different genes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1682-2) contains supplementary material, which is available to authorized users.

Inheritance and QTL Mapping of Leaf Nutrient Concentration in a Cotton Inter-Specific Derived RIL Population

Developing and deploying cotton cultivars with high nutrient uptake, use efficiency and tolerance to nutrient related soil stresses is desirable to assist sustainable soil management. Genetic variation, heritability, selection response and quantitative trait loci (QTLs) were investigated for five macronutrients (P, K, Ca, Mg, S) and five micronutrients (Fe, Mn, B, Zn, and Cu) in a recombinant inbred line (RIL) population from an inter-specific cross between Gossypium hirsutum cv. Guazuncho 2, and G. barbadense accession VH8-4602. Na and K/Na ratio were also studied as the imbalance between Na and other nutrients is detrimental to cotton growth and development. The concentrations of nutrients were measured for different plant parts of the two parents and for leaf samples of the whole population collected at early to peak flowering in field experiments over two years in a sodic Vertosol soil. Parental contrast was large for most nutrient concentrations in leaves when compared with other plant parts. Segregation for leaf nutrient concentration was observed within the population with transgression for P, K, K/Na ratio and all micronutrients. Genotypic difference was the major factor behind within-population variation for most nutrients, while narrow sense heritability was moderate (0.27 for Mn and Cu, and 0.43 for B). At least one significant QTL was identified for each nutrient except K and more than half of those QTLs were clustered on chromosomes 14, 18 and 22. Selection response was predicted to be low for P and all micronutrients except B, high for K, Na and B, and very high for K/Na ratio. Correlations were more common between macronutrients, Na and K/Na ratio where the nature and strength of the relations varied (r=-0.69 to 0.76). We conclude that there is sufficient genetic diversity between these two tetraploid cotton species that could be exploited to improve cotton nutrient status by introgressing species-unique favourable alleles.

Constructing a high-density linkage map for Gossypium hirsutum × Gossypium barbadense and identifying QTLs for lint percentage

To introgress the good fiber quality and yield from Gossypium barbadense into a commercial Upland cotton variety, a high-density simple sequence repeat (SSR) genetic linkage map was developed from a BC1 F1 population of Gossypium hirsutum × Gossypium barbadense. The map comprised 2,292 loci and covered 5115.16 centiMorgan (cM) of the cotton AD genome, with an average marker interval of 2.23 cM. Of the marker order for 1,577 common loci on this new map, 90.36% agrees well with the marker order on the D genome sequence genetic map. Compared with five published high-density SSR genetic maps, 53.14% of marker loci were newly discovered in this map. Twenty-six quantitative trait loci (QTLs) for lint percentage (LP) were identified on nine chromosomes. Nine stable or common QTLs could be used for marker-assisted selection. Fifty percent of the QTLs were from G. barbadense and increased LP by 1.07%-2.41%. These results indicated that the map could be used for screening chromosome substitution segments from G. barbadense in the Upland cotton background, identifying QTLs or genes from G. barbadense, and further developing the gene pyramiding effect for improving fiber yield and quality.

Distribution and evolution of cotton fiber development genes in the fibreless Gossypium raimondii genome

Cotton fiber represents the largest single cell in plants and they serve as models to study cell development. This study investigated the distribution and evolution of fiber Unigenes anchored to recombination hotspots between tetraploid cotton (Gossypium hirsutum) At and Dt subgenomes, and within a parental diploid cotton (Gossypium raimondii) D genome. Comparative analysis of At vs D and Dt vs D showed that 1) the D genome provides many fiber genes after its merger with another parental diploid cotton (Gossypium arboreum) A genome although the D genome itself does not produce any spinnable fiber; 2) similarity of fiber genes is higher between At vs D than between Dt vs D genomic hotspots. This is the first report that fiber genes have higher similarity between At and D than between Dt and D. The finding provides new insights into cotton genomic regions that would facilitate genetic improvement of natural fiber properties.

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

KEY MESSAGE

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

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.

Genetic structure, linkage disequilibrium and association mapping of Verticillium wilt resistance in elite cotton (Gossypium hirsutum L.) germplasm population

Understanding the population structure and linkage disequilibrium in an association panel can effectively avoid spurious associations and improve the accuracy in association mapping. In this study, one hundred and fifty eight elite cotton (Gossypium hirsutum L.) germplasm from all over the world, which were genotyped with 212 whole genome-wide marker loci and phenotyped with an disease nursery and greenhouse screening method, were assayed for population structure, linkage disequilibrium, and association mapping of Verticillium wilt resistance. A total of 480 alleles ranging from 2 to 4 per locus were identified from all collections. Model-based analysis identified two groups (G1 and G2) and seven subgroups (G1a–c, G2a–d), and differentiation analysis showed that subgroup having a single origin or pedigree was apt to differentiate with those having a mixed origin. Only 8.12% linked marker pairs showed significant LD (P<0.001) in this association panel. The LD level for linked markers is significantly higher than that for unlinked markers, suggesting that physical linkage strongly influences LD in this panel, and LD level was elevated when the panel was classified into groups and subgroups. The LD decay analysis for several chromosomes showed that different chromosomes showed a notable change in LD decay distances for the same gene pool. Based on the disease nursery and greenhouse environment, 42 marker loci associated with Verticillium wilt resistance were identified through association mapping, which widely were distributed among 15 chromosomes. Among which 10 marker loci were found to be consistent with previously identified QTLs and 32 were new unreported marker loci, and QTL clusters for Verticillium wilt resistanc on Chr.16 were also proved in our study, which was consistent with the strong linkage in this chromosome. Our results would contribute to association mapping and supply the marker candidates for marker-assisted selection of Verticillium wilt resistance in cotton.

A whole-genome DNA marker map for cotton based on the D-genome sequence of Gossypium raimondii L

We constructed a very-high-density, whole-genome marker map (WGMM) for cotton by using 18,597 DNA markers corresponding to 48,958 loci that were aligned to both a consensus genetic map and a reference genome sequence. The WGMM has a density of one locus per 15.6 kb, or an average of 1.3 loci per gene. The WGMM was anchored by the use of colinear markers to a detailed genetic map, providing recombinational information. Mapped markers occurred at relatively greater physical densities in distal chromosomal regions and lower physical densities in the central regions, with all 1 Mb bins having at least nine markers. Hotspots for quantitative trait loci and resistance gene analog clusters were aligned to the map and DNA markers identified for targeting of these regions of high practical importance. Based on the cotton D genome reference sequence, the locations of chromosome structural rearrangements plotted on the map facilitate its translation to other Gossypium genome types. The WGMM is a versatile genetic map for marker assisted breeding, fine mapping and cloning of genes and quantitative trait loci, developing new genetic markers and maps, genome-wide association mapping, and genome evolution studies.

Transmission ratio distortion of molecular markers in a doubled haploid population originated from a natural hybrid between Osmunda japonica and O. lancea

In ferns, intra-gametophytic selfing occurs as a mode of reproduction where two gametes from the same gametophyte form a completely homozygous sporophyte. Intra-gametophytic selfing is considered to be prevented by lethal or deleterious recessive genes in several diploid species. In order to investigate the modes and tempo of selection acting different developmental stages, doubled haploids obtained from intra-gametophytic selfing within isolated gametophytes of a putative F1 hybrid between Osmunda japonica and O. lancea were analyzed with EST_derived molecular markers, and the distribution pattern of transmission ratio distortion (TRD) along linkage map was clarified. As the results, the markers with skewness were clustered in two linkage groups. For the two highly distorted regions, gametophytes and F2 population were also examined. The markers skewed towards O. japonica on a linkage group (LG_2) showed skewness also in gametophytes, and the TRD was generated in the process of spore formation or growth of gametophytes. Also, selection appeared to be operating in the gametophytic stage. The markers on other linkage group (LG_11) showed highest skewness towards O. lancea in doubled haploids, and it was suggested that the segregation of LG_11 were influenced by zygotic lethality or genotypic evaluation and that some deleterious recessive genes exist in LG_11 and reduce the viability of homozygotes with O. japonica alleles. It is very likely that a region of LG_11were responsible for the low frequencies of intra-gametophytic selfing in O. japonica.

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

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

Molecular markers associated with the immature fiber (im) gene affecting the degree of fiber cell wall thickening in cotton (Gossypium hirsutum L.)

Cotton fiber fineness and maturity measured indirectly as micronaire (MIC) are important properties of determining fiber grades in the textile market. To understand the genetic control and molecular mechanisms of fiber fineness and maturity, we studied two near isogenic lines, Gossypium hirsutum, Texas Marker-1 wild type (TM-1) and immature fiber (im) mutant showing a significant difference in MIC values. The fibers from im mutant plants were finer and less mature with lower MIC values than those from the recurrent parent, TM-1. A comprehensive fiber property analysis of TM-1 and im mutant showed that the lower MIC of fibers in im mutant was due to the lower degree of fiber cell wall thickening as compared to the TM-1 fibers. Using an F(2) population comprising 366 progenies derived from a cross between TM-1 and im mutant, we confirmed that the immature fiber phenotype present in a mutant plant was controlled by one single recessive gene im. Furthermore, we identified 13 simple sequence repeat markers that were closely linked to the im gene located on chromosome 3. Molecular markers associated with the im gene will lay the foundation to further investigate genetic information required for improving cotton fiber fineness and maturity.

Mapping quantitative trait loci for lint yield and fiber quality across environments in a Gossypium hirsutum × Gossypium barbadense backcross inbred line population

Identification of stable quantitative trait loci (QTLs) across different environments and mapping populations is a prerequisite for marker-assisted selection (MAS) for cotton yield and fiber quality. To construct a genetic linkage map and to identify QTLs for fiber quality and yield traits, a backcross inbred line (BIL) population of 146 lines was developed from a cross between Upland cotton (Gossypium hirsutum) and Egyptian cotton (Gossypium barbadense) through two generations of backcrossing using Upland cotton as the recurrent parent followed by four generations of self pollination. The BIL population together with its two parents was tested in five environments representing three major cotton production regions in China. The genetic map spanned a total genetic distance of 2,895 cM and contained 392 polymorphic SSR loci with an average genetic distance of 7.4 cM per marker. A total of 67 QTLs including 28 for fiber quality and 39 for yield and its components were detected on 23 chromosomes, each of which explained 6.65-25.27% of the phenotypic variation. Twenty-nine QTLs were located on the At subgenome originated from a cultivated diploid cotton, while 38 were on the Dt subgenome from an ancestor that does not produce spinnable fibers. Of the eight common QTLs (12%) detected in more than two environments, two were for fiber quality traits including one for fiber strength and one for uniformity, and six for yield and its components including three for lint yield, one for seedcotton yield, one for lint percentage and one for boll weight. QTL clusters for the same traits or different traits were also identified. This research represents one of the first reports using a permanent advanced backcross inbred population of an interspecific hybrid population to identify QTLs for fiber quality and yield traits in cotton across diverse environments. It provides useful information for transferring desirable genes from G. barbadense to G. hirsutum using MAS.

Integrated consensus genetic and physical maps of flax (Linum usitatissimum L.)

Three linkage maps of flax (Linum usitatissimum L.) were constructed from populations CDC Bethune/Macbeth, E1747/Viking and SP2047/UGG5-5 containing between 385 and 469 mapped markers each. The first consensus map of flax was constructed incorporating 770 markers based on 371 shared markers including 114 that were shared by all three populations and 257 shared between any two populations. The 15 linkage group map corresponds to the haploid number of chromosomes of this species. The marker order of the consensus map was largely collinear in all three individual maps but a few local inversions and marker rearrangements spanning short intervals were observed. Segregation distortion was present in all linkage groups which contained 1-52 markers displaying non-Mendelian segregation. The total length of the consensus genetic map is 1,551 cM with a mean marker density of 2.0 cM. A total of 670 markers were anchored to 204 of the 416 fingerprinted contigs of the physical map corresponding to ~274 Mb or 74 % of the estimated flax genome size of 370 Mb. This high resolution consensus map will be a resource for comparative genomics, genome organization, evolution studies and anchoring of the whole genome shotgun sequence.

Deep sequencing reveals differences in the transcriptional landscapes of fibers from two cultivated species of cotton

Cotton (Gossypium) fiber is the most prevalent natural product used in the textile industry. The two major cultivated species, G. hirsutum (Gh) and G. barbadense (Gb), are allotetraploids with contrasting fiber quality properties. To better understand the molecular basis for their fiber differences, EST pyrosequencing was used to document the fiber transcriptomes at two key development stages, 10 days post anthesis (dpa), representing the peak of fiber elongation, and 22 dpa, representing the transition to secondary cell wall synthesis. The 617,000 high quality reads (89% of the total 692,000 reads) from 4 libraries were assembled into 46,072 unigenes, comprising 38,297 contigs and 7,775 singletons. Functional annotation of the unigenes together with comparative digital gene expression (DGE) revealed a diverse set of functions and processes that were partly linked to specific fiber stages. Globally, 2,770 contigs (7%) showed differential expression (>2-fold) between 10 and 22 dpa (irrespective of genotype), with 70% more highly expressed at 10 dpa, while 2,248 (6%) were differentially expressed between the genotypes (irrespective of stage). The most significant genes with differential DGE at 10 dpa included expansins and lipid transfer proteins (higher in Gb), while at 22 dpa tubulins, cellulose, and sucrose synthases showed higher expression in Gb. DGE was compared with expression data of 10 dpa-old fibers from Affymetrix microarrays. Among 543 contigs showing differential expression on both platforms, 74% were consistent in being either over-expressed in Gh (242 genes) or in Gb (161 genes). Furthermore, the unigene set served to identify 339 new SSRs and close to 21,000 inter-genotypic SNPs. Subsets of 88 SSRs and 48 SNPs were validated through mapping and added 65 new loci to a RIL genetic map. The new set of fiber ESTs and the gene-based markers complement existing available resources useful in basic and applied research for crop improvement in cotton.

A high density consensus genetic map of tetraploid cotton that integrates multiple component maps through molecular marker redundancy check

A consensus genetic map of tetraploid cotton was constructed using six high-density maps and after the integration of a sequence-based marker redundancy check. Public cotton SSR libraries (17,343 markers) were curated for sequence redundancy using 90% as a similarity cutoff. As a result, 20% of the markers (3,410) could be considered as redundant with some other markers. The marker redundancy information had been a crucial part of the map integration process, in which the six most informative interspecific Gossypium hirsutum×G. barbadense genetic maps were used for assembling a high density consensus (HDC) map for tetraploid cotton. With redundant markers being removed, the HDC map could be constructed thanks to the sufficient number of collinear non-redundant markers in common between the component maps. The HDC map consists of 8,254 loci, originating from 6,669 markers, and spans 4,070 cM, with an average of 2 loci per cM. The HDC map presents a high rate of locus duplications, as 1,292 markers among the 6,669 were mapped in more than one locus. Two thirds of the duplications are bridging homoeologous A(T) and D(T) chromosomes constitutive of allopolyploid cotton genome, with an average of 64 duplications per A(T)/D(T) chromosome pair. Sequences of 4,744 mapped markers were used for a mutual blast alignment (BBMH) with the 13 major scaffolds of the recently released Gossypium raimondii genome indicating high level of homology between the diploid D genome and the tetraploid cotton genetic map, with only a few minor possible structural rearrangements. Overall, the HDC map will serve as a valuable resource for trait QTL comparative mapping, map-based cloning of important genes, and better understanding of the genome structure and evolution of tetraploid cotton.

Development and mapping of SNP assays in allotetraploid cotton

A narrow germplasm base and a complex allotetraploid genome have made the discovery of single nucleotide polymorphism (SNP) markers difficult in cotton (Gossypium hirsutum). To generate sequence for SNP discovery, we conducted a genome reduction experiment (EcoRI, BafI double digest, followed by adapter ligation, biotin-streptavidin purification, and agarose gel separation) on two accessions of G. hirsutum and two accessions of G. barbadense. From the genome reduction experiment, a total of 2.04 million genomic sequence reads were assembled into contigs with an N(50) of 508 bp and analyzed for SNPs. A previously generated assembly of expressed sequence tags (ESTs) provided an additional source for SNP discovery. Using highly conservative parameters (minimum coverage of 8× at each SNP and 20% minor allele frequency), a total of 11,834 and 1,679 non-genic SNPs were identified between accessions of G. hirsutum and G. barbadense in genome reduction assemblies, respectively. An additional 4,327 genic SNPs were also identified between accessions of G. hirsutum in the EST assembly. KBioscience KASPar assays were designed for a portion of the intra-specific G. hirsutum SNPs. From 704 non-genic and 348 genic markers developed, a total of 367 (267 non-genic, 100 genic) mapped in a segregating F(2) population (Acala Maxxa × TX2094) using the Fluidigm EP1 system. A G. hirsutum genetic linkage map of 1,688 cM was constructed based entirely on these new SNP markers. Of the genic-based SNPs, we were able to identify within which genome ('A' or 'D') each SNP resided using diploid species sequence data. Genetic maps generated by these newly identified markers are being used to locate quantitative, economically important regions within the cotton genome.

cDNA-AFLP-based genetical genomics in cotton fibers

Genetical genomics, or genetic analysis applied to gene expression data, has not been widely used in plants. We used quantitative cDNA-AFLP to monitor the variation in the expression level of cotton fiber transcripts among a population of inter-specific Gossypium hirsutum × G. barbadense recombinant inbred lines (RILs). Two key fiber developmental stages, elongation (10 days post anthesis, dpa), and secondary cell wall thickening (22 dpa), were studied. Normalized intensity ratios of 3,263 and 1,201 transcript-derived fragments (TDFs) segregating over 88 RILs were analyzed for quantitative trait loci (QTL) mapping for the 10 and 22 dpa fibers, respectively. Two-thirds of all TDFs mapped between 1 and 6 eQTLs (LOD > 3.5). Chromosome 21 had a higher density of eQTLs than other chromosomes in both data sets and, within chromosomes, hotspots of presumably trans-acting eQTLs were identified. The eQTL hotspots were compared to the location of phenotypic QTLs for fiber characteristics among the RILs, and several cases of co-localization were detected. Quantitative RT-PCR for 15 sequenced TDFs showed that 3 TDFs had at least one eQTL at a similar location to those identified by cDNA-AFLP, while 3 other TDFs mapped an eQTL at a similar location but with opposite additive effect. In conclusion, cDNA-AFLP proved to be a cost-effective and highly transferable platform for genome-wide and population-wide gene expression profiling. Because TDFs are anonymous, further validation and interpretation (in silico analysis, qPCR gene profiling) of the eQTL and eQTL hotspots will be facilitated by the increasing availability of cDNA and genomic sequence resources in cotton.

Development and mapping of SNP assays in allotetraploid cotton

A narrow germplasm base and a complex allotetraploid genome have made the discovery of single nucleotide polymorphism (SNP) markers difficult in cotton (Gossypium hirsutum). To generate sequence for SNP discovery, we conducted a genome reduction experiment (EcoRI, BafI double digest, followed by adapter ligation, biotin-streptavidin purification, and agarose gel separation) on two accessions of G. hirsutum and two accessions of G. barbadense. From the genome reduction experiment, a total of 2.04 million genomic sequence reads were assembled into contigs with an N(50) of 508 bp and analyzed for SNPs. A previously generated assembly of expressed sequence tags (ESTs) provided an additional source for SNP discovery. Using highly conservative parameters (minimum coverage of 8× at each SNP and 20% minor allele frequency), a total of 11,834 and 1,679 non-genic SNPs were identified between accessions of G. hirsutum and G. barbadense in genome reduction assemblies, respectively. An additional 4,327 genic SNPs were also identified between accessions of G. hirsutum in the EST assembly. KBioscience KASPar assays were designed for a portion of the intra-specific G. hirsutum SNPs. From 704 non-genic and 348 genic markers developed, a total of 367 (267 non-genic, 100 genic) mapped in a segregating F(2) population (Acala Maxxa × TX2094) using the Fluidigm EP1 system. A G. hirsutum genetic linkage map of 1,688 cM was constructed based entirely on these new SNP markers. Of the genic-based SNPs, we were able to identify within which genome ('A' or 'D') each SNP resided using diploid species sequence data. Genetic maps generated by these newly identified markers are being used to locate quantitative, economically important regions within the cotton genome.

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.

Saturation of an intra-gene pool linkage map: towards a unified consensus linkage map for fine mapping and synteny analysis in common bean

Map-based cloning and fine mapping to find genes of interest and marker assisted selection (MAS) requires good genetic maps with reproducible markers. In this study, we saturated the linkage map of the intra-gene pool population of common bean DOR364 × BAT477 (DB) by evaluating 2,706 molecular markers including SSR, SNP, and gene-based markers. On average the polymorphism rate was 7.7% due to the narrow genetic base between the parents. The DB linkage map consisted of 291 markers with a total map length of 1,788 cM. A consensus map was built using the core mapping populations derived from inter-gene pool crosses: DOR364 × G19833 (DG) and BAT93 × JALO EEP558 (BJ). The consensus map consisted of a total of 1,010 markers mapped, with a total map length of 2,041 cM across 11 linkage groups. On average, each linkage group on the consensus map contained 91 markers of which 83% were single copy markers. Finally, a synteny analysis was carried out using our highly saturated consensus maps compared with the soybean pseudo-chromosome assembly. A total of 772 marker sequences were compared with the soybean genome. A total of 44 syntenic blocks were identified. The linkage group Pv6 presented the most diverse pattern of synteny with seven syntenic blocks, and Pv9 showed the most consistent relations with soybean with just two syntenic blocks. Additionally, a co-linear analysis using common bean transcript map information against soybean coding sequences (CDS) revealed the relationship with 787 soybean genes. The common bean consensus map has allowed us to map a larger number of markers, to obtain a more complete coverage of the common bean genome. Our results, combined with synteny relationships provide tools to increase marker density in selected genomic regions to identify closely linked polymorphic markers for indirect selection, fine mapping or for positional cloning.

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.