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.

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.

Gossypium hirsutumcontributes the most production of cotton fibre, butG. barbadenseis valued for its better comprehensive resistance and superior fibre properties. However, the allotetraploid genome ofG. barbadensehas not been comprehensively analysed. Here we present a high-quality assembly of the 2.57 gigabase genome ofG. barbadense, including 80,876 protein-coding genes. The double-sized genome of the A (or At) (1.50 Gb) against D (or Dt) (853 Mb) primarily resulted from the expansion of Gypsy elements, including Peabody and Retrosat2 subclades in the Del clade and the Athila subclade in the Athila/Tat clade. Substantial gene expansion and contraction were observed and rich homoeologous gene pairs with biased expression patterns were identified, suggesting abundant gene sub-functionalization occurred by allopolyploidization. More specifically, theCesAgene family has adapted differentially temporal expression patterns, suggesting an integrated regulatory mechanism ofCesAgenes from At and Dt subgenomes for the primary and secondary cellulose biosynthesis of cotton fibre in a “relay race”-like fashion. We anticipate that theG. barbadensegenome sequence will advance our understanding the mechanism of genome polyploidization and underpin genome-wide comparison research in this genus.

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.

Summary: The Arabidopsis gene AVP1 encodes a vacuolar pyrophosphatase that functions as a proton pump on the vacuolar membrane. Overexpression of AVP1 in Arabidopsis, tomato and rice enhances plant performance under salt and drought stress conditions, because up‐regulation of the type I H⁺‐PPase from Arabidopsis may result in a higher proton electrochemical gradient, which facilitates enhanced sequestering of ions and sugars into the vacuole, reducing water potential and resulting in increased drought‐ and salt tolerance when compared to wild‐type plants. Furthermore, overexpression of AVP1 stimulates auxin transport in the root system and leads to larger root systems, which helps transgenic plants absorb water more efficiently under drought conditions. Using the same approach, AVP1‐expressing cotton plants were created and tested for their performance under high‐salt and reduced irrigation conditions. The AVP1‐expressing cotton plants showed more vigorous growth than wild‐type plants in the presence of 200 mm NaCl under hydroponic growth conditions. The soil‐grown AVP1‐expressing cotton plants also displayed significantly improved tolerance to both drought and salt stresses in greenhouse conditions. Furthermore, the fibre yield of AVP1‐expressing cotton plants is at least 20% higher than that of wild‐type plants under dry‐land conditions in the field. This research indicates that AVP1 has the potential to be used for improving crop’s drought‐ and salt tolerance in areas where water and salinity are limiting factors for agricultural productivity.