Search

Allotetraploid cotton is an economically important natural-fiber-producing crop worldwide. After polyploidization, Gossypium hirsutum L. evolved to produce a higher fiber yield and to better survive harsh environments than Gossypium barbadense, which produces superior-quality fibers. The global genetic and molecular bases for these interspecies divergences were unknown. Here we report high-quality de novo-assembled genomes for these two cultivated allotetraploid species with pronounced improvement in repetitive-DNA-enriched centromeric regions. Whole-genome comparative analyses revealed that species-specific alterations in gene expression, structural variations and expanded gene families were responsible for speciation and the evolutionary history of these species. These findings help to elucidate the evolution of cotton genomes and their domestication history. The information generated not only should enable breeders to improve fiber quality and resilience to ever-changing environmental conditions but also can be translated to other crops for better understanding of their domestication history and use in improvement.

Cotton, with cellulose-enriched mature fibers, is the largest source of natural textiles. Through a map-based cloning strategy, we isolated an industrially important lint fiber development gene (Li 3 ) that encodes an MYB-MIXTA-like transcription factor (MML) on chromosome D12 (GhMML4_D12). Virus-induced gene silencing or decreasing the expression of the GhMML4_D12 gene in n 2 NSM plants resulted in a significant reduction in epidermal cell prominence and lint fiber production. GhMML4_D12 is arranged in tandem with GhMML3, another MIXTA gene responsible for fuzz fiber development. These two very closely related MIXTA genes direct fiber initiation production in two specialized cell forms: lint and fuzz fibers. They may control the same metabolic pathways in different cell types. The MIXTAs expanded in Malvaceae during their evolution and produced a Malvaceae-specific family that regulates epidermal cell differentiation, different from the gene family that regulates leaf hair trichome development. Cotton has developed a unique transcriptional regulatory network for fiber development. Characterization of target genes regulating fiber production has provided insights into the molecular mechanisms underlying cotton fiber development and has allowed the use of genetic engineering to increase lint yield by inducing more epidermal cells to develop into lint rather than fuzz fibers., (© 2017 The Authors. New Phytologist © 2017 New Phytologist Trust.)

Map-based gene cloning is a vital strategy for the identification of the quantitative trait loci or genes underlying important agronomic traits. The conventional map-based cloning method is powerful but generally time-consuming and labor-intensive. In this context, we introduce an improved bulked segregant analysis method in combination with a virus-induced gene silencing (VIGS) strategy for rapid and reliable gene mapping, identification and functional verification. This method was applied to a multiple recessive marker line of upland cotton, Texas 582 (T582), and identified unique genomic positions harboring mutant loci, showing the reliability and efficacy of this method. The v1 locus was further fine-mapped. Only one gene, GhCHLI, which encodes one of the subunits of Mg chelatase, was differentially down-regulated in T582 compared with TM-1. A point mutation occurred in the AAA+ conserved region of GhCHLI and led to an amino acid substitution. Suppression of its expression by VIGS in TM-1 resulted in a yellow blade phenotype that was similar to T582. This integrated approach provides a paradigm for the rapid mapping and identification of the candidate genes underlying the genetic traits in plants with large and complex genomes in the future., (© The Author 2017. Published by Oxford University Press on behalf of the Society for Experimental Biology.)

Background: Cotton has been cultivated and used to make fabrics for at least 7000 years. Two allotetraploid species of great commercial importance, Gossypium hirsutum and Gossypium barbadense, were domesticated after polyploidization and are cultivated worldwide. Although the overall genetic diversity between these two cultivated species has been studied with limited accessions, their population structure and genetic variations remain largely unknown., Results: We resequence the genomes of 147 cotton accessions, including diverse wild relatives, landraces, and modern cultivars, and construct a comprehensive variation map to provide genomic insights into the divergence and dual domestication of these two important cultivated tetraploid cotton species. Phylogenetic analysis shows two divergent groups for G. hirsutum and G. barbadense, suggesting a dual domestication processes in tetraploid cottons. In spite of the strong genetic divergence, a small number of interspecific reciprocal introgression events are found between these species and the introgression pattern is significantly biased towards the gene flow from G. hirsutum into G. barbadense. We identify selective sweeps, some of which are associated with relatively highly expressed genes for fiber development and seed germination., Conclusions: We report a comprehensive analysis of the evolution and domestication history of allotetraploid cottons based on the whole genomic variation between G. hirsutum and G. barbadense and between wild accessions and modern cultivars. These results provide genomic bases for improving cotton production and for further evolution analysis of polyploid crops.

The morphology of cotton leaves varies considerably. Phenotypes, including okra, sea-island, super-okra, and broad leaf, are controlled by a multiple allele locus, L 2 Okra leaf (L 2 °) is an incomplete mutation that alters leaf shape by increasing the length of lobes with deeper sinuses. Using a map-based cloning strategy, we cloned the L 2 locus gene, which encodes a LATE MERISTEM IDENTITY 1 (LMI1)-like transcription factor (GhOKRA). Silencing GhOKRA leads to a change in phenotype from okra to broad leaf. Overexpression of GhOKRA in Arabidopsis thaliana greatly increases the degree of the leaf lobes and changes the leaf shape. Premature termination of translation in GhOKRA results in the production of broad leaves. The sequences of OKRA from diploid progenitor D-genome species, and wild races and domesticated allotetraploid cottons in Gossypium hirsutum show that a premature termination mutation occurred before and after the formation of tetraploid cotton, respectively. This study provides genomic insights into the two interspecific hybridization events: one produced the present broad leaf and another formed okra leaf phenotype with complete OKRA, that occurred during allotetraploid cotton formation., (Copyright © 2016 by the Genetics Society of America.)

Natural antisense transcripts (NATs) are commonly observed in eukaryotic genomes, but only a limited number of such genes have been identified as being involved in gene regulation in plants. In this research, we investigated the function of small RNA derived from a NAT in fiber cell development. Using a map-based cloning strategy for the first time in tetraploid cotton, we cloned a naked seed mutant gene (N1 ) encoding a MYBMIXTA-like transcription factor 3 (MML3)/GhMYB25-like in chromosome A12, GhMML3_A12, that is associated with fuzz fiber development. The extremely low expression of GhMML3_A12 in N1 is associated with NAT production, driven by its 3' antisense promoter, as indicated by the promoter-driven histochemical staining assay. In addition, small RNA deep sequencing analysis suggested that the bidirectional transcriptions of GhMML3_A12 form double-stranded RNAs and generate 21-22 nt small RNAs. Therefore, in a fiber-specific manner, small RNA derived from the GhMML3_A12 locus can mediate GhMML3_A12 mRNA self-cleavage and result in the production of naked seeds followed by lint fiber inhibition in N1 plants. The present research reports the first observation of gene-mediated NATs and siRNA directly controlling fiber development in cotton., (© 2016 The Authors. New Phytologist © 2016 New Phytologist Trust.)

Background: SNPs are the most abundant polymorphism type, and have been explored in many crop genomic studies, including rice and maize. SNP discovery in allotetraploid cotton genomes has lagged behind that of other crops due to their complexity and polyploidy. In this study, genome-wide SNPs are detected systematically using next-generation sequencing and efficient SNP genotyping methods, and used to construct a linkage map and characterize the structural variations in polyploid cotton genomes., Results: We construct an ultra-dense inter-specific genetic map comprising 4,999,048 SNP loci distributed unevenly in 26 allotetraploid cotton linkage groups and covering 4,042 cM. The map is used to order tetraploid cotton genome scaffolds for accurate assembly of G. hirsutum acc. TM-1. Recombination rates and hotspots are identified across the cotton genome by comparing the assembled draft sequence and the genetic map. Using this map, genome rearrangements and centromeric regions are identified in tetraploid cotton by combining information from the publicly-available G. raimondii genome with fluorescent in situ hybridization analysis., Conclusions: We report the genotype-by-sequencing method used to identify millions of SNPs between G. hirsutum and G. barbadense. We construct and use an ultra-dense SNP map to correct sequence mis-assemblies, merge scaffolds into pseudomolecules corresponding to chromosomes, detect genome rearrangements, and identify centromeric regions in allotetraploid cottons. We find that the centromeric retro-element sequence of tetraploid cotton derived from the D subgenome progenitor might have invaded the A subgenome centromeres after allotetrapolyploid formation. This study serves as a valuable genomic resource for genetic research and breeding of cotton.

Upland cotton is a model for polyploid crop domestication and transgenic improvement. Here we sequenced the allotetraploid Gossypium hirsutum L. acc. TM-1 genome by integrating whole-genome shotgun reads, bacterial artificial chromosome (BAC)-end sequences and genotype-by-sequencing genetic maps. We assembled and annotated 32,032 A-subgenome genes and 34,402 D-subgenome genes. Structural rearrangements, gene loss, disrupted genes and sequence divergence were more common in the A subgenome than in the D subgenome, suggesting asymmetric evolution. However, no genome-wide expression dominance was found between the subgenomes. Genomic signatures of selection and domestication are associated with positively selected genes (PSGs) for fiber improvement in the A subgenome and for stress tolerance in the D subgenome. This draft genome sequence provides a resource for engineering superior cotton lines.

The cotton fiber, as a single-celled trichome, is a biological model system for studying cell differentiation and elongation. However, the complexity of gene expression and regulation in the fiber complicates genetic research. In this study, we investigated the genome-wide transcriptome profiling in Texas Marker-1 (TM-1) and five naked seed or fuzzless mutants (three dominant and two recessive) during the fuzz initial development stage. More than three million clean tags were generated from each sample representing the expression data for 27,325 genes, which account for 72.8% of the annotated Gossypium raimondii primary transcript genes. Thousands of differentially expressed genes (DEGs) were identified between TM-1 and the mutants. Based on functional enrichment analysis, the DEGs downregulated in the mutants were enriched in protein synthesis-related genes and transcription factors, while DEGs upregulated in the mutants were enriched in DNA/chromatin structure-related genes and transcription factors. Pathway analysis showed that ATP synthesis, and sugar and lipid metabolism-related pathways play important roles in fuzz initial development. Also, we identified a large number of transcription factors such as MYB, bHLH, HB, WRKY, AP2/EREBP, bZIP and C2H2 zinc finger families that were differently expressed between TM-1 and the mutants, and were also related to trichome development in Arabidopsis.

Cotton is the source of the most important, renewable natural textile fiber and oil in the world. MicroRNAs (miRNAs) are endogenous, non-coding, approximately 18-24 nucleotides long RNAs and function in the negative regulation of their target genes. Two mostly overlapping libraries of small RNA molecules were constructed and sequenced, and served as repetition sets of data to identify miRNAs involved in fiber initiation and seed development. The D genome sequence of Gossypium raimondii was used in conjunction with EST sequences to predict miRNA precursors. Overall, 93 new miRNA precursors were identified, of which 28 belonged to 10 known families and the other 65 were considered to be novel miRNAs. Seven hundred EST sequences were proposed to be candidate target genes which involved in the regulation of a diverse group of genes with diverse functions and transcription factors. Some of the novel miRNAs and candidate target genes were validated by the Northern blot and rapid amplification of 5' cDNA ends (5' RACE).