Allotetraploid cotton species (Gossypium hirsutum and Gossypium barbadense) have long been cultivated worldwide for natural renewable textile fibers. The draft genome sequences of both species are available but they are highly fragmented and incomplete1-4. Here we report reference-grade genome assemblies and annotations for G. hirsutum accession Texas Marker-1 (TM-1) and G. barbadense accession 3-79 by integrating single-molecule real-time sequencing, BioNano optical mapping and high-throughput chromosome conformation capture techniques. Compared with previous assembled draft genomes1,3, these genome sequences show considerable improvements in contiguity and completeness for regions with high content of repeats such as centromeres. Comparative genomics analyses identify extensive structural variations that probably occurred after polyploidization, highlighted by large paracentric/pericentric inversions in 14 chromosomes. We constructed an introgression line population to introduce favorable chromosome segments from G. barbadense to G. hirsutum, allowing us to identify 13 quantitative trait loci associated with superior fiber quality. These resources will accelerate evolutionary and functional genomic studies in cotton and inform future breeding programs for fiber improvement.

Summary: The CRISPR/Cas9 system has been extensively applied for crop improvement. However, our understanding of Cas9 specificity is very limited in Cas9‐edited plants. To identify on‐ and off‐target mutation in an edited crop, we described whole genome sequencing (WGS) of 14 Cas9‐edited cotton plants targeted to three genes, and three negative (Ne) control and three wild‐type (WT) plants. In total, 4188–6404 unique single‐nucleotide polymorphisms (SNPs) and 312–745 insertions/deletions (indels) were detected in 14 Cas9‐edited plants compared to WT, negative and cotton reference genome sequences. Since the majority of these variations lack a protospacer‐adjacent motif (PAM), we demonstrated that the most variations following Cas9‐edited are due either to somaclonal variation or/and pre‐existing/inherent variation from maternal plants, but not off‐target effects. Of a total of 4413 potential off‐target sites (allowing ≤5 mismatches within the 20‐bp sgRNA and 3‐bp PAM sequences), the WGS data revealed that only four are bona fide off‐target indel mutations, validated by Sanger sequencing. Moreover, inherent genetic variation of WT can generate novel off‐target sites and destroy PAMs, which suggested great care should be taken to design sgRNA for the minimizing of off‐target effect. These findings suggested that CRISPR/Cas9 system is highly specific for cotton plants.

Background: Millennia of directional human selection has reshaped the genomic architecture of cultivated cotton relative to wild counterparts, but we have limited understanding of the selective retention and fractionation of genomic components.Results: We construct a comprehensive genomic variome based on 1961 cottons and identify 456 Mb and 357 Mb of sequence with domestication and improvement selection signals and 162 loci, 84 of which are novel, including 47 loci associated with 16 agronomic traits. Using pan-genome analyses, we identify 32,569 and 8851 non-reference genes lost fromGossypium hirsutumandGossypium barbadensereference genomes respectively, of which 38.2% (39,278) and 14.2% (11,359) of genes exhibit presence/absence variation (PAV). We document the landscape of PAV selection accompanied by asymmetric gene gain and loss and identify 124 PAVs linked to favorable fiber quality and yield loci.Conclusions: This variation repertoire points to genomic divergence during cotton domestication and improvement, which informs the characterization of favorable gene alleles for improved breeding practice using a pan-genome-based approach.