Quantitative genomics-enabled selection for simultaneous improvement of lint yield and seed traits in cotton (Gossypium hirsutum L.).

KEY MESSAGE: A Bayesian linkage disequilibrium-based multiple-locus mixed model identified QTLs for fibre, seed and oil traits and predicted breeding worthiness of test lines, enabling their simultaneous improvement in cotton. Improving cotton seed and oil yields has become increasingly important while continuing to breed for higher lint yield. In this study, a novel Bayesian linkage disequilibrium-based multiple-locus mixed model was developed for QTL identification and genomic prediction (GP). A multi-parent population consisting of 256 recombinant inbred lines, derived from four elite cultivars with distinct combinations of traits, was used in the analysis of QTLs for lint percentage, seed index, lint index and seed oil content and their interrelations. All four traits were moderately heritable and correlated but with no large influence of genotype × environment interactions across multiple seasons. Seven to ten major QTLs were identified for each trait with many being adjacent or overlapping for different trait pairs. A fivefold cross-validation of the model indicated prediction accuracies of 0.46-0.62. GP results based on any two-season phenotypes were strongly correlated with phenotypic means of a pooled analysis of three-season experiments (r = 0.83-0.92). When used for selection of improvement in lint, seed and oil yields, GP captured 40-100% of individuals with comparable lint yields of those selected based on the three-season phenotypic results. Thus, this quantitative genomics-enabled approach can not only decipher the genomic variation underlying lint, seed and seed oil traits and their interrelations, but can provide predictions for their simultaneous improvement. We discuss future breeding strategies in cotton that will enhance the entire value of the crop, not just its fibre.

Cotton Breeding in Australia: Meeting the Challenges of the 21st Century.

The Commonwealth Scientific and Industrial Research Organisation (CSIRO) cotton breeding program is the sole breeding effort for cotton in Australia, developing high performing cultivars for the local industry which is worth∼AU$3 billion per annum. The program is supported by Cotton Breeding Australia, a Joint Venture between CSIRO and the program's commercial partner, Cotton Seed Distributors Ltd. (CSD). While the Australian industry is the focus, CSIRO cultivars have global impact in North America, South America, and Europe. The program is unique compared with many other public and commercial breeding programs because it focuses on diverse and integrated research with commercial outcomes. It represents the full research pipeline, supporting extensive long-term fundamental molecular research; native and genetically modified (GM) trait development; germplasm enhancement focused on yield and fiber quality improvements; integration of third-party GM traits; all culminating in the release of new commercial cultivars. This review presents evidence of past breeding successes and outlines current breeding efforts, in the areas of yield and fiber quality improvement, as well as the development of germplasm that is resistant to pests, diseases and abiotic stressors. The success of the program is based on the development of superior germplasm largely through field phenotyping, together with strong commercial partnerships with CSD and Bayer CropScience. These relationships assist in having a shared focus and ensuring commercial impact is maintained, while also providing access to markets, traits, and technology. The historical successes, current foci and future requirements of the CSIRO cotton breeding program have been used to develop a framework designed to augment our breeding system for the future. This will focus on utilizing emerging technologies from the genome to phenome, as well as a panomics approach with data management and integration to develop, test and incorporate new technologies into a breeding program. In addition to streamlining the breeding pipeline for increased genetic gain, this technology will increase the speed of trait and marker identification for use in genome editing, genomic selection and molecular assisted breeding, ultimately producing novel germplasm that will meet the coming challenges of the 21st Century.

Characterization and Genetic Mapping of Black Root Rot Resistance in Gossypium arboreum L.

Black root rot (BRR) is an economically important disease of cotton and other crops, especially in cooler regions with short growing seasons. Symptoms include black discoloration of the roots, reduced number of lateral roots and stunted or slow plant growth. The cultivated tetraploid Gossypium species are susceptible to BRR. Resistance to BRR was identified in G. arboreum accession BM13H and is associated with reduced and restricted hyphal growth and less sporulation. Transcriptome analysis indicates that BM13H responds to infection at early time points 2- and 3-days post-inoculation, but by day 5, few differentially expressed genes are observed between infected and uninfected roots. Inheritance of BM13H resistance to BRR was evaluated in an F6 recombinant inbred population and shows a single semi-dominant locus conferring resistance that was fine mapped to a region on chromosome 1, containing ten genes including five putative resistance-like genes.

Correction to: Improved lint yield under field conditions in cotton over-expressing transcription factors regulating fibre initiation.

Due to an unfortunate misunderstanding, an extra middle initial erroneously appeared in the original publication and the full name of the first author should read Shi Ming Liu.

Improved lint yield under field conditions in cotton over-expressing transcription factors regulating fibre initiation.

Only a few transcription factors (TFs) regulating which cells of the ovule epidermis differentiate into lint fibres have been identified in cotton (Gossypium hirsutum L.). In this study, the effect on lint yield and fibre quality of over-expressing three TFs in cotton, GhHD-1, GhMYB25 and GhMYB25Like, and their double and triple combinations, were evaluated in field experiments over two seasons. The expression of single or stacked TFs were all driven either by an ovule-specific promoter, FBP 7, or a constitutive promoter, Stunt 7, in a Coker 315 background. TF type, either singly or in combination, was found to be the most significant factor affecting lint yield. Among 64 transgenic lines tested, seven were higher yielding than null segregant lines in one or both seasons and were all from the sets with single and double over-expressed TF combinations. A reduced yield was associated with the set of triple combinations. The two most stable high yielding lines across the seasons recorded 12-22% higher yields than the nulls, although were not competitive to locally adapted commercial controls. Over-expression of TFs singly or in combination did not significantly alter fibre length and strength, but sometimes increased fibre micronaire. There were positive relationships between lint yield and lint percentage and lint yield and fibre density amongst the transgenic lines. Our preliminary results suggest that manipulating TF expression, either singly or in pairs, can increase the density of fibres initiated on developing seeds and fibre yields under field conditions while maintaining overall fibre quality.

Expansion of MIR482/2118 by a class-II transposable element in cotton.

Some plant microRNA (miRNA) families contain multiple members generating identical or highly similar mature miRNA variants. Mechanisms underlying the expansion of miRNA families remain elusive, although tandem and/or segmental duplications have been proposed. In this study of two tetraploid cottons, Gossypium hirsutum and Gossypium barbadense, and their extant diploid progenitors, Gossypium arboreum and Gossypium raimondii, we investigated the gain and loss of members of the miR482/2118 superfamily, which modulates the expression of nucleotide-binding site leucine-rich repeat (NBS-LRR) disease resistance genes. We found significant expansion of MIR482/2118d in G. barbadense, G. hirsutum and G. raimondii, but not in G. arboreum. Several newly expanded MIR482/2118d loci have mutated to produce different miR482/2118 variants with altered target-gene specificity. Based on detailed analysis of sequences flanking these MIR482/2118 loci, we found that this expansion of MIR482/2118d and its derivatives resulted from an initial capture of an MIR482/2118d by a class-II DNA transposable element (TE) in G. raimondii prior to the tetraploidization event, followed by transposition to new genomic locations in G. barbadense, G. hirsutum and G. raimondii. The 'GosTE' involved in the capture and proliferation of MIR482/2118d and its derivatives belongs to the PIF/Harbinger superfamily, generating a 3-bp target site duplication upon insertion at new locations. All orthologous MIR482/2118 loci in the two diploids were retained in the two tetraploids, but mutation(s) in miR482/2118 were observed across all four species as well as in different cultivars of both G. barbadense and G. hirsutum, suggesting a dynamic co-evolution of miR482/2118 and its NBS-LRR targets. Our results provide fresh insights into the mechanisms contributing to MIRNA proliferation and enrich our knowledge on TEs.

Trichomes at the Base of the Petal Are Regulated by the Same Transcription Factors as Cotton Seed Fibers.

Many polypetalous plants have a constriction at the base of the petal that leaves a small gap that can provide entry into the young flower bud before the reproductive organs are fully developed. In cotton (Gossypium hirsutum L.), this gap is occluded by tufts of short unicellular trichomes superficially resembling the fibers found on cotton seeds. We are just beginning to understand the developmental regulation of the seed fibers and have previously characterized several MIXTA-like MYB transcription factors (TFs) that are critical for correct seed fiber development but know little about the molecular regulation of other types of cotton trichomes. Here, using RNAi or dominant suppression transgenic cotton lines and natural fiber mutants, we investigated the development and regulation of the petal base trichomes. Petal base trichomes and seed trichomes were also examined across several different species within and outside of the Malvoideae. We found that the petal base trichomes are regulated by the same MYB TFs as cotton seed fibers and, since they are more widely distributed across different taxa than the seed fibers, could have preceded them in the evolution of these important textile fibers produced by some cotton species.

Engineering Trienoic Fatty Acids into Cottonseed Oil Improves Low-Temperature Seed Germination, Plant Photosynthesis and Cotton Fiber Quality.

Alpha-linolenic acid (ALA, 18:3Δ9,12,15) and γ-linolenic acid \ (GLA, 18:3Δ6,9,12) are important trienoic fatty acids, which are beneficial for human health in their own right, or as precursors for the biosynthesis of long-chain polyunsaturated fatty acids. ALA and GLA in seed oil are synthesized from linoleic acid (LA, 18:2Δ9,12) by the microsomal ω-3 fatty acid desaturase (FAD3) and Δ6 desaturase (D6D), respectively. Cotton (Gossypium hirsutum L.) seed oil composition was modified by transforming with an FAD3 gene from Brassica napus and a D6D gene from Echium plantagineum, resulting in approximately 30% ALA and 20% GLA, respectively. The total oil content in transgenic seeds remained unaltered relative to parental seeds. Despite the use of a seed-specific promoter for transgene expression, low levels of GLA and increased levels of ALA were found in non-seed cotton tissues. At low temperature, the germinating cottonseeds containing the linolenic acid isomers elongated faster than the untransformed controls. ALA-producing lines also showed higher photosynthetic rates at cooler temperature and better fiber quality compared to both untransformed controls and GLA-producing lines. The oxidative stability of the novel cottonseed oils was assessed, providing guidance for potential food, pharmaceutical and industrial applications of these oils.

Correction to: Tissue and cell-specific transcriptomes in cotton reveal the subtleties of gene regulation underlying the diversity of plant secondary cell walls.

Upon publication of the original article [1], the authors had flagged that Fig. 1 had been published twice, as both Fig. 1 and Additional file 3.

Historical Datasets Support Genomic Selection Models for the Prediction of Cotton Fiber Quality Phenotypes Across Multiple Environments.

Genomic selection (GS) has successfully been used in plant breeding to improve selection efficiency and reduce breeding time and cost. However, there has not been a study to evaluate GS prediction models that may be used for predicting cotton breeding lines across multiple environments. In this study, we evaluated the performance of Bayes Ridge Regression, BayesA, BayesB, BayesC and Reproducing Kernel Hilbert Spaces regression models. We then extended the single-site GS model to accommodate genotype × environment interaction (G×E) in order to assess the merits of multi- over single-environment models in a practical breeding and selection context in cotton, a crop for which this has not previously been evaluated. Our study was based on a population of 215 upland cotton (Gossypium hirsutum) breeding lines which were evaluated for fiber length and strength at multiple locations in Australia and genotyped with 13,330 single nucleotide polymorphic (SNP) markers. BayesB, which assumes unique variance for each marker and a proportion of markers to have large effects, while most other markers have zero effect, was the preferred model. GS accuracy for fiber length based on a single-site model varied across sites, ranging from 0.27 to 0.77 (mean = 0.38), while that of fiber strength ranged from 0.19 to 0.58 (mean = 0.35) using randomly selected sub-populations as the training population. Prediction accuracies from the M×E model were higher than those for single-site and across-site models, with an average accuracy of 0.71 and 0.59 for fiber length and strength, respectively. The use of the M×E model could therefore identify which breeding lines have effects that are stable across environments and which ones are responsible for G×E and so reduce the amount of phenotypic screening required in cotton breeding programs to identify adaptable genotypes.

Genetic dissection of the fuzzless seed trait in Gossypium barbadense.

Cotton fibres are single-celled trichomes arising from the epidermal cells of the seed coat and may be either long (lint) or very short (fuzz). The dominant fuzzless N1 of Gossypium hirsutum is a defective allele of the At-subgenome homoeolog of MYB25-like, but the genetic components underlying the recessive fuzzless trait from G. barbadense (Gb) are unknown. We have identified five genetic loci, including a major contributing locus containing MYB25-like_Dt, associated with Gb fuzzless seeds based on genotyping of fuzzy and fuzzless near isogenic lines (NILs) from an interspecies cross (G. barbadense × G. hirsutum). At 3 d post-anthesis when fuzz fibres are initiating, expression of MYB25-like_Dt was significantly lower in fuzzless NILs than in fuzzy seeded NILs, while higher MYB25-like_Dt expression was associated with more seed fuzz across different cotton genotypes. Phenotypic and genotypic analysis of MYB25-like homoeoalleles in cottons showing different fibre phenotypes and their crossing progeny indicated that both MYB25-like_At and MYB25-like_Dt are associated with lint development, and that fuzz development is mainly determined by the expression level of MYB25-like_Dt at ~3 d post-anthesis. Expression of Gb fuzzless seeds depends on genetic background and interactions amongst the multiple loci identified. MYB25-like_Dt is one of the best candidates for N2.

Tissue and cell-specific transcriptomes in cotton reveal the subtleties of gene regulation underlying the diversity of plant secondary cell walls.

BACKGROUND: Knowledge of plant secondary cell wall (SCW) regulation and deposition is mainly based on the Arabidopsis model of a 'typical' lignocellulosic SCW. However, SCWs in other plants can vary from this. The SCW of mature cotton seed fibres is highly cellulosic and lacks lignification whereas xylem SCWs are lignocellulosic. We used cotton as a model to study different SCWs and the expression of the genes involved in their formation via RNA deep sequencing and chemical analysis of stem and seed fibre. RESULTS: Transcriptome comparisons from cotton xylem and pith as well as from a developmental series of seed fibres revealed tissue-specific and developmentally regulated expression of several NAC transcription factors some of which are likely to be important as top tier regulators of SCW formation in xylem and/or seed fibre. A so far undescribed hierarchy was identified between the top tier NAC transcription factors SND1-like and NST1/2 in cotton. Key SCW MYB transcription factors, homologs of Arabidopsis MYB46/83, were practically absent in cotton stem xylem. Lack of expression of other lignin-specific MYBs in seed fibre relative to xylem could account for the lack of lignin deposition in seed fibre. Expression of a MYB103 homolog correlated with temporal expression of SCW CesAs and cellulose synthesis in seed fibres. FLAs were highly expressed and may be important structural components of seed fibre SCWs. Finally, we made the unexpected observation that cell walls in the pith of cotton stems contained lignin and had a higher S:G ratio than in xylem, despite that tissue's lacking many of the gene transcripts normally associated with lignin biosynthesis. CONCLUSIONS: Our study in cotton confirmed some features of the currently accepted gene regulatory cascade for 'typical' plant SCWs, but also revealed substantial differences, especially with key downstream NACs and MYBs. The lignocellulosic SCW of cotton xylem appears to be achieved differently from that in Arabidopsis. Pith cell walls in cotton stems are compositionally very different from that reported for other plant species, including Arabidopsis. The current definition of a 'typical' primary or secondary cell wall might not be applicable to all cell types in all plant species.

Diversity analysis of cotton (Gossypium hirsutum L.) germplasm using the CottonSNP63K Array.

BACKGROUND: Cotton germplasm resources contain beneficial alleles that can be exploited to develop germplasm adapted to emerging environmental and climate conditions. Accessions and lines have traditionally been characterized based on phenotypes, but phenotypic profiles are limited by the cost, time, and space required to make visual observations and measurements. With advances in molecular genetic methods, genotypic profiles are increasingly able to identify differences among accessions due to the larger number of genetic markers that can be measured. A combination of both methods would greatly enhance our ability to characterize germplasm resources. Recent efforts have culminated in the identification of sufficient SNP markers to establish high-throughput genotyping systems, such as the CottonSNP63K array, which enables a researcher to efficiently analyze large numbers of SNP markers and obtain highly repeatable results. In the current investigation, we have utilized the SNP array for analyzing genetic diversity primarily among cotton cultivars, making comparisons to SSR-based phylogenetic analyses, and identifying loci associated with seed nutritional traits. RESULTS: The SNP markers distinctly separated G. hirsutum from other Gossypium species and distinguished the wild from cultivated types of G. hirsutum. The markers also efficiently discerned differences among cultivars, which was the primary goal when designing the CottonSNP63K array. Population structure within the genus compared favorably with previous results obtained using SSR markers, and an association study identified loci linked to factors that affect cottonseed protein content. CONCLUSIONS: Our results provide a large genome-wide variation data set for primarily cultivated cotton. Thousands of SNPs in representative cotton genotypes provide an opportunity to finely discriminate among cultivated cotton from around the world. The SNPs will be relevant as dense markers of genome variation for association mapping approaches aimed at correlating molecular polymorphisms with variation in phenotypic traits, as well as for molecular breeding approaches in cotton.

Trichomes control flower bud shape by linking together young petals.

Trichomes are widespread in plants and develop from surface cells on different tissues(1). They have many forms and functions, from defensive spines to physical barriers that trap layers of air to insulate against desiccation, but there is growing evidence that trichomes can also have developmental roles in regulating flower structure(2,3). We report here that the trichomes on petals of cotton, Gossypium hirsutum L., are essential for correct flower bud shape through a mechanical entanglement of the trichomes on adjacent petals that anchor the edges to counter the opposing force generated by asymmetric expansion of overlapping petals. Silencing a master regulator of petal trichomes, GhMYB-MIXTA-Like10 (GhMYBML10), by RNA interference (RNAi) suppressed petal trichome growth and resulted in flower buds forming into abnormal corkscrew shapes that exposed developing anthers and stigmas to desiccation damage. Artificially gluing petal edges together could partially restore correct bud shape and fertility. Such petal 'Velcro' is present in other Malvaceae and perhaps more broadly in other plant families, although it is not ubiquitous. This mechanism for physical association between separate organs to regulate flower shape and function is different from the usual organ shape control(4) exerted through cell-to-cell communication and differential cell expansion within floral tissues(5,6).

Integrated mapping and characterization of the gene underlying the okra leaf trait in Gossypium hirsutum L.

Diverse leaf morphology has been observed among accessions of Gossypium hirsutum, including okra leaf, which has advantages and disadvantages in cotton production. The okra leaf locus has been mapped to chromosome 15 of the Dt subgenome, but the underlying gene has yet to be identified. In this study, we used a combination of targeted association analysis, F2 population-based fine mapping, and comparative sequencing of orthologues to identify a candidate gene underlying the okra leaf trait in G. hirsutum. The okra leaf gene identified, GhOKRA, encoded a homeodomain leucine-zipper class I protein, whose closely related genes in several other plant species have been shown to be involved in regulating leaf morphology. The transcript levels of GhOKRA in shoot apices were positively correlated with the phenotypic expression of the okra leaf trait. Of the multiple sequence variations observed in the coding region among GrOKRA of Gossypium raimondii and GhOKRA-Dt of normal and okra/superokra leaf G. hirsutum accessions, a non-synonymous substitution near the N terminus and the variable protein sequences at the C terminus may be related to the leaf shape difference. Our results suggest that both transcription and protein activity of GhOKRA may be involved in regulating leaf shape. Furthermore, we found that non-reciprocal homoeologous recombination, or gene conversion, may have played a role in the origin of the okra leaf allele. Our results provided tools for further investigating and understanding the fundamental biological processes that are responsible for the cotton leaf shape variation and will help in the design of cotton plants with an ideal leaf shape for enhanced cotton production.

GbEXPATR, a species-specific expansin, enhances cotton fibre elongation through cell wall restructuring.

Cotton provides us the most important natural fibre. High fibre quality is the major goal of cotton breeding, and introducing genes conferring longer, finer and stronger fibre from Gossypium barbadense to Gossypium hirsutum is an important breeding strategy. We previously analysed the G. barbadense fibre development mechanism by gene expression profiling and found two homoeologous fibre-specific α-expansins from G. barbadense, GbEXPA2 and GbEXPATR. GbEXPA2 (from the DT genome) is a classical α-expansin, while its homoeolog, GbEXPATR (AT genome), encodes a truncated protein lacking the normal C-terminal polysaccharide-binding domain of other α-expansins and is specifically expressed in G. barbadense. Silencing EXPA in G. hirsutum induced shorter fibres with thicker cell walls. GbEXPA2 overexpression in G. hirsutum had no effect on mature fibre length, but produced fibres with a slightly thicker wall and increased crystalline cellulose content. Interestingly, GbEXPATR overexpression resulted in longer, finer and stronger fibres coupled with significantly thinner cell walls. The longer and thinner fibre was associated with lower expression of a number of secondary wall-associated genes, especially chitinase-like genes, and walls with lower cellulose levels but higher noncellulosic polysaccharides which advocated that a delay in the transition to secondary wall synthesis might be responsible for better fibre. In conclusion, we propose that α-expansins play a critical role in fibre development by loosening the cell wall; furthermore, a truncated form, GbEXPATR, has a more dramatic effect through reorganizing secondary wall synthesis and metabolism and should be a candidate gene for developing G. hirsutum cultivars with superior fibre quality.

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.

Sequencing of allotetraploid cotton (Gossypium hirsutum L. acc. TM-1) provides a resource for fiber improvement.

Upland cotton is a model for polyploid crop domestication and transgenic improvement. Here we sequenced the allotetraploid Gossypium hirsutum L. acc. TM-1 genome by integrating whole-genome shotgun reads, bacterial artificial chromosome (BAC)-end sequences and genotype-by-sequencing genetic maps. We assembled and annotated 32,032 A-subgenome genes and 34,402 D-subgenome genes. Structural rearrangements, gene loss, disrupted genes and sequence divergence were more common in the A subgenome than in the D subgenome, suggesting asymmetric evolution. However, no genome-wide expression dominance was found between the subgenomes. Genomic signatures of selection and domestication are associated with positively selected genes (PSGs) for fiber improvement in the A subgenome and for stress tolerance in the D subgenome. This draft genome sequence provides a resource for engineering superior cotton lines.

Development of a 63K SNP Array for Cotton and High-Density Mapping of Intraspecific and Interspecific Populations of Gossypium spp.

High-throughput genotyping arrays provide a standardized resource for plant breeding communities that are useful for a breadth of applications including high-density genetic mapping, genome-wide association studies (GWAS), genomic selection (GS), complex trait dissection, and studying patterns of genomic diversity among cultivars and wild accessions. We have developed the CottonSNP63K, an Illumina Infinium array containing assays for 45,104 putative intraspecific single nucleotide polymorphism (SNP) markers for use within the cultivated cotton species Gossypium hirsutum L. and 17,954 putative interspecific SNP markers for use with crosses of other cotton species with G. hirsutum. The SNPs on the array were developed from 13 different discovery sets that represent a diverse range of G. hirsutum germplasm and five other species: G. barbadense L., G. tomentosum Nuttal × Seemann, G. mustelinum Miers × Watt, G. armourianum Kearny, and G. longicalyx J.B. Hutchinson and Lee. The array was validated with 1,156 samples to generate cluster positions to facilitate automated analysis of 38,822 polymorphic markers. Two high-density genetic maps containing a total of 22,829 SNPs were generated for two F2 mapping populations, one intraspecific and one interspecific, and 3,533 SNP markers were co-occurring in both maps. The produced intraspecific genetic map is the first saturated map that associates into 26 linkage groups corresponding to the number of cotton chromosomes for a cross between two G. hirsutum lines. The linkage maps were shown to have high levels of collinearity to the JGI G. raimondii Ulbrich reference genome sequence. The CottonSNP63K array, cluster file and associated marker sequences constitute a major new resource for the global cotton research community.