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 incomplete 1-4. Here we report referencegrade 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 genomes 1,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. Cotton represents the largest source of natural textile fibers in the world. Over 90% of annual fiber production comes from allotetraploid cotton (G. hirsutum and G. barbadense), which originated from an allopolyplodization event approximately 1-2 million year ago, followed by millennia of asymmetric subgenome selection 5,6. G. hirsutum is cultivated all over the world because of its high yield and G. barbadense is prized for its superior fiber quality. To cultivate G. hirsutum that produces longer, finer and stronger fibers, one approach is to introduce the superior fiber traits from G. barbadense into G. hirsutum. A genomics-enabled breeding strategy requires a detailed and robust understanding of genomic organization. Genomic feature G. hirsutum G. barbadense

Publisher's copyright statement:Additional information: Use policyThe full-text may be used and/or reproduced, and given to third parties in any format or medium, without prior permission or charge, for personal research or study, educational, or not-for-pro t purposes provided that:• a full bibliographic reference is made to the original source • a link is made to the metadata record in DRO • the full-text is not changed in any way The full-text must not be sold in any format or medium without the formal permission of the copyright holders.Please consult the full DRO policy for further details. Dt, 0.56 × 10 -3 ) ( Fig. 1d and Supplementary Fig. 3). This shows that a large amount is associated with the development of the long fiber trait in cultivated cotton (Fig. 3b). 217Domestication has led to the transformation of cotton fiber from brown to white. 218To understand this phenomenon, we examined two homoeologous gene pairs only 219 subjected to domestication selection in the Dt, 4-COUMARATE:COA LIGASE (4CL) 220 and CHALCONE SYNTHASE (CHS), which encode enzymes involved in the 221 phenylpropanoid metabolic pathway ( Fig. 3c and Supplementary Fig. 6 Fig. 3c). These SNPs display reductions in nucleotide diversity that occurred 225 during domestication (Fig. 3c). Interestingly, we found that the two SNPs in the Fig. 8) 42 . We identified a total of 188,360 DNase I-hypersensitive 248 sites (DHSs) in cotton leaves and fibers, of which ca. 47% are common to both tissues 249 (Fig. 4a). DHSs were preferentially identified in chromosomal arms and 250 approximately half were detected in promoter and intergenic regions ( Fig. 4b and 251 Supplementary Fig. 9). We found DHSs are hypo-methylated, consistent with 252 previous studies 42 (Fig. 4c) H3K4me1 and inactive H3K9me2 (Fig. 4d). Intergenic DHSs were also found to 255 exhibit an enrichment of H3K4me3 and H3K27me3, but depletion of H3K9me2 and 256 no enrichment of H3K4me1 (Fig. 4e). As predicted, the patterns of chromatin 257 modification marks in cotton are different between genic and TE regions 258 ( Supplementary Fig. 10). In addition, genes with promoter DHSs are generally 259 expressed at a higher level in both tissues than those without promoter DHSs (Fig. 4f), 260 and tissue-specific promoter DHSs corresponded to higher levels of gene expression 261 ( Fig. 4g) Hi-C analysis was carried out using the TM-1 accession to characterize global 296 chromatin interactions. We generated 1.1 billion Hi-C paired-end reads, of which ca. possible Hi-C bias, HindIII fragments of less than 2 kb were merged to obtain 299 305,682 chromosomal anchor regions (Fig. 5a). On the basis of a high-quality 300 genome assembly of TM-1 (Supplementary Fig. 11), we used the Hi-C data to 301 characterize the cotton chromatin interactome (Supplementary Fig. 12) and ( Fig. 5b), but many topologically associated domain-like (TAD-like) regions were 305 identified (Fig. 5c, Supplementary Fig. 13 and Supplementary are less frequent at regions marked by H3K9me2 (Fig. 5d). (Fig. 5g). 320We...

Summary Gossypium hirsutum is an allotetraploid with a complex genome. Most genes have multiple copies that belong to At and Dt subgenomes. Sequence similarity is also very high between gene homologues. To efficiently achieve site/gene‐specific mutation is quite needed. Due to its high efficiency and robustness, the CRISPR (clustered regularly interspaced short palindromic repeats)/Cas9 system has exerted broad site‐specific genome editing from prokaryotes to eukaryotes. In this study, we utilized a CRISPR/Cas9 system to generate two sgRNAs in a single vector to conduct multiple sites genome editing in allotetraploid cotton. An exogenously transformed gene Discosoma red fluorescent protein2(DsRed2) and an endogenous gene GhCLA1 were chosen as targets. The DsRed2‐edited plants in T0 generation reverted its traits to wild type, with vanished red fluorescence the whole plants. Besides, the mutated phenotype and genotype were inherited to their T1 progenies. For the endogenous gene GhCLA1, 75% of regenerated plants exhibited albino phenotype with obvious nucleotides and DNA fragments deletion. The efficiency of gene editing at each target site is 66.7–100%. The mutation genotype was checked for both genes with Sanger sequencing. Barcode‐based high‐throughput sequencing, which could be highly efficient for genotyping to a population of mutants, was conducted in GhCLA1‐edited T0 plants and it matched well with Sanger sequencing results. No off‐target editing was detected at the potential off‐target sites. These results prove that the CRISPR/Cas9 system is highly efficient and reliable for allotetraploid cotton genome editing.

Plants are constantly challenged by a multitude of pathogens and pests, which causes massive yield and quality losses annually. A promising approach to reduce such losses is to enhance the immune system of plants through genetic engineering. Previous work has shown that laccases (p-diphenol:dioxygen oxidoreductase, EC 1.10.3.2) function as lignin polymerization enzymes. Here we demonstrate that transgenic manipulation of the expression of the laccase gene GhLac1 in cotton (Gossypium hirsutum) can confer an enhanced defense response to both pathogens and pests. Overexpression of GhLac1 leads to increased lignification, associated with increased tolerance to the fungal pathogen Verticillium dahliae and to the insect pests cotton bollworm (Helicoverpa armigera) and cotton aphid (Aphis gosypii). Suppression of GhLac1 expression leads to a redirection of metabolic flux in the phenylpropanoid pathway, causing the accumulation of JA and secondary metabolites that confer resistance to V. dahliae and cotton bollworm; it also leads to increased susceptibility to cotton aphid. Plant laccases therefore provide a new molecular tool to engineer pest and pathogen resistance in crops.

Male sterility caused by long-term high-temperature (HT) stress occurs widely in crops. MicroRNAs (miRNAs), a class of endogenous non-coding small RNAs, play an important role in the plant response to various abiotic stresses. To dissect the working principle of miRNAs in male sterility under HT stress in cotton, a total of 112 known miRNAs, 270 novel miRNAs and 347 target genes were identified from anthers of HT-insensitive (84021) and HT-sensitive (H05) cotton cultivars under normal-temperature and HT conditions through small RNA and degradome sequencing. Quantitative reverse transcriptase-polymerase chain reaction and 5'-RNA ligase-mediated rapid amplification of cDNA ends experiments were used to validate the sequencing data. The results show that miR156 was suppressed by HT stress in both 84021 and H05; miR160 was suppressed in 84021 but induced in H05. Correspondingly, SPLs (target genes of miR156) were induced both in 84021 and H05; ARF10 and ARF17 (target genes of miR160) were induced in 84021 but suppressed in H05. Overexpressing miR160 increased cotton sensitivity to HT stress seen as anther indehiscence, associated with the suppression of ARF10 and ARF17 expression, thereby activating the auxin response that leads to anther indehiscence. Supporting this role for auxin, exogenous Indole-3-acetic acid (IAA) leads to a stronger male sterility phenotype both in 84021 and H05 under HT stress. Cotton plants overexpressing miR157 suppressed the auxin signal, and also showed enhanced sensitivity to HT stress, with microspore abortion and anther indehiscence. Thus, we propose that the auxin signal, mediated by miRNAs, is essential for cotton anther fertility under HT stress.