Richard G. Percy

The draft genome of a diploid cotton Gossypium raimondii

We have sequenced and assembled a draft genome of G. raimondii, whose progenitor is the putative contributor of the D subgenome to the economically important fiber-producing cotton species Gossypium hirsutum and Gossypium barbadense. Over 73% of the assembled sequences were anchored on 13 G. raimondii chromosomes. The genome contains 40,976 protein-coding genes, with 92.2% of these further confirmed by transcriptome data. Evidence of the hexaploidization event shared by the eudicots as well as of a cotton-specific whole-genome duplication approximately 13–20 million years ago was observed. We identified 2,355 syntenic blocks in the G. raimondii genome, and we found that approximately 40% of the paralogous genes were present in more than 1 block, which suggests that this genome has undergone substantial chromosome rearrangement during its evolution. Cotton, and probably Theobroma cacao, are the only sequenced plant species that possess an authentic CDN1 gene family for gossypol biosynthesis, as revealed by phylogenetic analysis.

Genome sequence of the cultivated cotton Gossypium arboreum

The complex allotetraploid nature of the cotton genome (AADD; 2n = 52) makes genetic, genomic and functional analyses extremely challenging. Here we sequenced and assembled the Gossypium arboreum (AA; 2n = 26) genome, a putative contributor of the A subgenome. A total of 193.6 Gb of clean sequence covering the genome by 112.6-fold was obtained by paired-end sequencing. We further anchored and oriented 90.4% of the assembly on 13 pseudochromosomes and found that 68.5% of the genome is occupied by repetitive DNA sequences. We predicted 41,330 protein-coding genes in G. arboreum. Two whole-genome duplications were shared by G. arboreum and Gossypium raimondii before speciation. Insertions of long terminal repeats in the past 5 million years are responsible for the twofold difference in the sizes of these genomes. Comparative transcriptome studies showed the key role of the nucleotide binding site (NBS)-encoding gene family in resistance to Verticillium dahliae and the involvement of ethylene in the development of cotton fiber cells.

Genome sequence of cultivated Upland cotton ( Gossypium hirsutum TM-1) provides insights into genome evolution

Gossypium hirsutum has proven difficult to sequence owing to its complex allotetraploid (AtDt) genome. Here we produce a draft genome using 181-fold paired-end sequences assisted by fivefold BAC-to-BAC sequences and a high-resolution genetic map. In our assembly 88.5% of the 2,173-Mb scaffolds, which cover 89.6%∼96.7% of the AtDt genome, are anchored and oriented to 26 pseudochromosomes. Comparison of this G. hirsutum AtDt genome with the already sequenced diploid Gossypium arboreum (AA) and Gossypium raimondii (DD) genomes revealed conserved gene order. Repeated sequences account for 67.2% of the AtDt genome, and transposable elements (TEs) originating from Dt seem more active than from At. Reduction in the AtDt genome size occurred after allopolyploidization. The A or At genome may have undergone positive selection for fiber traits. Concerted evolution of different regulatory mechanisms for Cellulose synthase (CesA) and 1-Aminocyclopropane-1-carboxylic acid oxidase1 and 3 (ACO1,3) may be important for enhanced fiber production in G. hirsutum.

Toward Sequencing Cotton (Gossypium) Genomes

Despite rapidly decreasing costs and innovative technologies, sequencing of angiosperm genomes is not yet undertaken lightly. Generating larger amounts of sequence data more quickly does not address the difficulties of sequencing and assembling complex genomes de novo. The cotton ( Gossypium spp.)

Allozyme evidence for the origin and diversification of Gossypium barbadense L.

SummaryGossypium barbadense L. is a commercially important cotton species of tropical South American origin presently grownin many regions of the world. The species is morphologically diverse, consisting of a wide range of wild (or feral), commensal, landrace, and highly improvedcommercial forms. We performed allozyme analysis on 153 accessions representing the spectrum of G. barbadense diversityto ascertain the geographic origin of the species, its patterns of diffusion subsequent to domestication, and to reveal infraspecific relationships. Levels ofgenetic variation in G. barbadense are moderate. Of 59 loci scored, 24 were polymorphic, with a mean number of alleles perlocus of 1.69 and an average panmictic heterozygosity of 0.062. Principal component analysis revealed geographic clustering of accessions into six relativelydiscrete regions. Gene frequencies at many loci are significantly heterogeneous among these regions, with an average GSTof 0.272. Northwestern South America contains the greatest genetic variability; we suggest that this region is the ancestral home of the species. The data indicate separate diffusion pathways from this region into Argentina-Paraguay and into eastern and northern South America east of the Andes. Caribbean Island and Central American forms appear to be derived from the latter. These diffusion pathways are in accordance with morphological evidence and historical record. In contrast to expectations based on geographic proximity, Pacific Island forms have their closest affinity to accessions from eastern South America. Advanced cultivated stocks seem largely derived from western Andean material, but also contain introgressed G. hirsutum germ plasm. Introgression was relatively high (22%–50% of accessions) in commercial stocks and in forms from Argentina-Paraguay and various Pacific Islands, but was conspicuously low or absent in material from Central America and the Caribbean, where commensal and commercial forms of both species are sympatric.

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.

Allozyme diversity and introgression in the Galapagos Islands endemic Gossypium darwinii and its relationship to continental G. barbadense

Abstract Gossypium darwinii Watt is a tetraploid cotton endemic to the Galapagos Islands. Opinion has been divided as to whether or not it deserves recognition at the specific rank, with some considering it a variety of its presumed progenitor, the widely distributed South American species G. barbadense L. A previous hypothesis states that much of the perceived intergradation between the two taxa arose as a consequence of introgression from G. barbadense following its introduction to the archipelago during the past several hundred years. We performed allozyme analysis on 58 accessions of G. darwinii from six islands, using 17 enzymes collectively encoded by 59 loci. Levels of variation were high for an island endemic, with a mean number of alleles per locus of 1.34 and an average panmictic heterozygosity of 0.062. Principal component analysis revealed clustering of accessions according to their island of origin, and a spatial pattern of island-clusters that approximates geographical relationships among islands. Genetic relationships of G. darwinii with G. barbadense and G. hirsutum L. were studied using previously generated allozyme data. Significant introgression of G. hirsutum alleles was detected; however morphological considerations support the hypothesis that much of G. darwiniis diversity stems from interspecific gene flow from G. barbadense , Evidence is presented suggesting that the occurrence of G. hirsutum alleles in G. darwinii derives not from direct hybridization, but from a mediated transfer through introduced, G. hirsutum -introgressed; G. barbadense . Gossypium darwinii and G. barbadense are nearly fixed for different alleles at four loci and each contains a large number of unique alleles. Notwithstanding the high interspecific Neis genetic identity (0.949), the allozyme data support geographical and morphological evidence in suggesting that a specific rank for G. darwinii is warranted.

Effects of Chromosome-Specific Introgression in Upland Cotton on Fiber and Agronomic Traits

Interspecific chromosome substitution is among the most powerful means of introgression and steps toward quantitative trait locus (QTL) identification. By reducing the genetic “noise” from other chromosomes, it greatly empowers the detection of genetic effects by specific chromosomes on quantitative traits. Here, we report on such results for 14 cotton lines (CS-B) with specific chromosomes or chromosome arms from G. barbadense L. substituted into G. hirsutum and chromosome-specific F2 families. Boll size, lint percentage, micronaire, 2.5% span length, elongation, strength, and yield were measured by replicated field experiments in five diverse environments and analyzed under an additive–dominance (AD) genetic model with genotype and environment interaction. Additive effects were significant for all traits and dominance effects were significant for all traits except 2.5% span length. CS-B25 had additive effects increasing fiber strength and fiber length and decreasing micronaire. CS-B16 and CS-B18 had additive effects related to reduced yields. The results point toward specific chromosomes of G. barbadense 3-79 as the probable locations of the genes that significantly affect quantitative traits of importance. Our results provided a scope to analyze individual chromosomes of the genome in homozygous and heterozygous conditions and thus detected novel effects of alleles controlling important QTL.

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.

Linkage Map Construction and Quantitative Trait Locus Analysis of Agronomic and Fiber Quality Traits in Cotton

The superior fiber properties of Gossypium barbadense L. serve as a source of novel variation for improving fiber quality in Upland cotton (G. hirsutum L.), but introgression from G. barbadense has been largely unsuccessful due to hybrid breakdown and a lack of genetic and genomic resources. In an effort to overcome these limitations, we constructed a linkage map and conducted a quantitative trait locus (QTL) analysis of 10 agronomic and fiber quality traits in a recombinant inbred mapping population derived from a cross between TM‐1, an Upland cotton line, and NM24016, an elite G. hirsutum line with stabilized introgression from G. barbadense. The linkage map consisted of 429 simple‐sequence repeat (SSR) and 412 genotyping‐by‐sequencing (GBS)‐based single‐nucleotide polymorphism (SNP) marker loci that covered half of the tetraploid cotton genome. Notably, the 841 marker loci were unevenly distributed among the 26 chromosomes of tetraploid cotton. The 10 traits evaluated on the TM‐1 × NM24016 population in a multienvironment trial were highly heritable, and most of the fiber traits showed considerable transgressive variation. Through the QTL analysis, we identified a total of 28 QTLs associated with the 10 traits. Our study provides a novel resource that can be used by breeders and geneticists for the genetic improvement of agronomic and fiber quality traits in Upland cotton.