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.

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

The incompatible pathosystem between resistant cotton (Gossypium barbadense cv. 7124) and Verticillium dahliae strain V991 was used to study the cotton transcriptome changes after pathogen inoculation by RNA-Seq. Of 32 774 genes detected by mapping the tags to assembly cotton contigs, 3442 defence-responsive genes were identified. Gene cluster analyses and functional assignments of differentially expressed genes indicated a significant transcriptional complexity. Quantitative real-time PCR (qPCR) was performed on selected genes with different expression levels and functional assignments to demonstrate the utility of RNA-Seq for gene expression profiles during the cotton defence response. Detailed elucidation of responses of leucine-rich repeat receptor-like kinases (LRR-RLKs), phytohormone signalling-related genes, and transcription factors described the interplay of signals that allowed the plant to fine-tune defence responses. On the basis of global gene regulation of phenylpropanoid metabolism-related genes, phenylpropanoid metabolism was deduced to be involved in the cotton defence response. A closer look at the expression of these genes, enzyme activity, and lignin levels revealed differences between resistant and susceptible cotton plants. Both types of plants showed an increased level of expression of lignin synthesis-related genes and increased phenylalanine-ammonia lyase (PAL) and peroxidase (POD) enzyme activity after inoculation with V. dahliae, but the increase was greater and faster in the resistant line. Histochemical analysis of lignin revealed that the resistant cotton not only retains its vascular structure, but also accumulates high levels of lignin. Furthermore, quantitative analysis demonstrated increased lignification and cross-linking of lignin in resistant cotton stems. Overall, a critical role for lignin was believed to contribute to the resistance of cotton to disease.

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...

Verticillium wilt causes massive annual losses of cotton yield, but the mechanism of cotton resistance to Verticillium dahliae is complex and poorly understood. In this study, a comparative proteomic analysis was performed in resistant cotton (Gossypium barbadense cv7124) on infection with V. dahliae. A total of 188 differentially expressed proteins were identified by mass spectrometry (MALDI-TOF/TOF) analysis and could be classified into 17 biological processes based on Gene Ontology annotation. Most of these proteins were implicated in stimulus response, cellular processes and metabolic processes. Based on the proteomic analysis, several genes involved in secondary metabolism, reactive oxygen burst and phytohormone signaling pathways were identified for further physiological and molecular analysis. The roles of the corresponding genes were further characterized by employing virus-induced gene silencing (VIGS). Based on the results, we suggest that the production of gossypol is sufficient to affect the cotton resistance to V. dahliae. Silencing of GbCAD1, a key enzyme involving in gossypol biosynthesis, compromised cotton resistance to V. dahliae. Reactive oxygen species and salicylic acid signaling may be also implicated as regulators in cotton responsive to V. dahliae according to the analysis of GbSSI2, an important regulator in the crosstalk between salicylic acid and jasmonic acid signal pathways. Moreover, brassinosteroids and jasmonic acid signaling may play essential roles in the cotton disease resistance to V. dahliae. The brassinosteroids signaling was activated in cotton on inoculation with V. dahliae and the disease resistance of cotton was enhanced after exogenous application of brassinolide. Meanwhile, jasmonic acid signaling was also activated in cotton after inoculation with V. dahliae and brassinolide application. These data provide highlights in the molecular basis of cotton resistance to V. dahliae. Molecular & Cellular