Asymmetric subgenome selection and cis-regulatory divergence during cotton domestication

Comparative population genomics offers an excellent opportunity for unraveling the genetic history of crop domestication. Upland cotton (Gossypium hirsutum) has long been an important economic crop, but a genome-wide and evolutionary understanding of the effects of human selection is lacking. Here, we describe a variation map for 352 wild and domesticated cotton accessions. We scanned 93 domestication sweeps occupying 74 Mb of the A subgenome and 104 Mb of the D subgenome, and identified 19 candidate loci for fiber-quality-related traits through a genome-wide association study. We provide evidence showing asymmetric subgenome domestication for directional selection of long fibers. Global analyses of DNase I–hypersensitive sites and 3D genome architecture, linking functional variants to gene transcription, demonstrate the effects of domestication on cis-regulatory divergence. This study provides new insights into the evolution of gene organization, regulation and adaptation in a major crop, and should serve as a rich resource for genome-based cotton improvement.

[1]  P. Cook.  The organization of replication and transcription. , 1999, Science.

[2]  T. Hattori,et al.  A bZIP factor, TRAB1, interacts with VP1 and mediates abscisic acid-induced transcription. , 1999, Proceedings of the National Academy of Sciences of the United States of America.

[3]  N. Chua,et al.  Characterization of the arabidopsis formin-like protein AFH1 and its interacting protein. , 2000, Plant & cell physiology.

[4]  J. Wendel,et al.  Comparative development of fiber in wild and cultivated cotton , 2001, Evolution & development.

[5]  P. Hussey,et al.  Formins: intermediates in signal-transduction cascades that affect cytoskeletal reorganization. , 2002, Trends in plant science.

[6]  A. Paterson,et al.  Rate variation among nuclear genes and the age of polyploidy in Gossypium. , 2003, Molecular biology and evolution.

[7]  M. Stephens,et al.  Inference of population structure using multilocus genotype data: linked loci and correlated allele frequencies. , 2003, Genetics.

[8]  Jia-Wei Wang,et al.  Control of Plant Trichome Development by a Cotton Fiber MYB Genew⃞ , 2004, The Plant Cell Online.

[9]  H. Matsubara,et al.  SPIRAL1 Encodes a Plant-Specific Microtubule-Localized Protein Required for Directional Control of Rapidly Expanding Arabidopsis Cells , 2004, The Plant Cell Online.

[10]  Alison G. Smith,et al.  Functional identification of AtFao3, a membrane bound long chain alcohol oxidase in Arabidopsis thaliana , 2004, FEBS letters.

[11]  Royston Goodacre,et al.  Identification of Novel Genes in Arabidopsis Involved in Secondary Cell Wall Formation Using Expression Profiling and Reverse Genetics , 2005, The Plant Cell Online.

[12]  Longfu Zhu,et al.  Isolation by suppression-subtractive hybridization of genes preferentially expressed during early and late fiber development stages in cotton , 2006, Molecular Biology.

[13]  Bruce D. Smith,et al.  The Molecular Genetics of Crop Domestication , 2006, Cell.

[14]  D. Reich,et al.  Principal components analysis corrects for stratification in genome-wide association studies , 2006, Nature Genetics.

[15]  A. Paterson,et al.  A rapid method for extraction of cotton (Gossypium spp.) genomic DNA suitable for RFLP or PCR analysis , 1993, Plant Molecular Biology Reporter.

[16]  Edward S. Buckler,et al.  TASSEL: software for association mapping of complex traits in diverse samples , 2007, Bioinform..

[17]  John Z. Yu,et al.  Toward Sequencing Cotton (Gossypium) Genomes , 2007, Plant Physiology.

[18]  Manuel A. R. Ferreira,et al.  PLINK: a tool set for whole-genome association and population-based linkage analyses. , 2007, American journal of human genetics.

[19]  Terrence S. Furey,et al.  F-Seq: a feature density estimator for high-throughput sequence tags , 2008, Bioinform..

[20]  A. Schnittger,et al.  Transcriptional Profiling of Mature Arabidopsis Trichomes Reveals That NOECK Encodes the MIXTA-Like Transcriptional Regulator MYB1061[C][W][OA] , 2008, Plant Physiology.

[21]  Yuxian Zhu,et al.  Molecular cloning, expression profiling, and yeast complementation of 19 β-tubulin cDNAs from developing cotton ovules , 2008, Journal of experimental botany.

[22]  Lex E. Flagel,et al.  The Evolution of Spinnable Cotton Fiber Entailed Prolonged Development and a Novel Metabolism , 2008, PLoS genetics.

[23]  E. Liu,et al.  An Oestrogen Receptor α-bound Human Chromatin Interactome , 2009, Nature.

[24]  I. Amit,et al.  Comprehensive mapping of long range interactions reveals folding principles of the human genome , 2011 .

[25]  Lior Pachter,et al.  Sequence Analysis , 2020, Definitions.

[26]  E. Nielsen,et al.  The Rab GTPase RabA4d Regulates Pollen Tube Tip Growth in Arabidopsis thaliana[W] , 2009, The Plant Cell Online.

[27]  Gonçalo R. Abecasis,et al.  The Sequence Alignment/Map format and SAMtools , 2009, Bioinform..

[28]  Richard Durbin,et al.  Sequence analysis Fast and accurate short read alignment with Burrows – Wheeler transform , 2009 .

[29]  N. Harberd,et al.  High Temperature-Mediated Adaptations in Plant Architecture Require the bHLH Transcription Factor PIF4 , 2009, Current Biology.

[30]  R. Wilson,et al.  BreakDancer: An algorithm for high resolution mapping of genomic structural variation , 2009, Nature Methods.

[31]  Ryan A. Rapp,et al.  Gene expression in developing fibres of Upland cotton (Gossypium hirsutum L.) was massively altered by domestication , 2010, BMC Biology.

[32]  Cole Trapnell,et al.  Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. , 2010, Nature biotechnology.

[33]  M. DePristo,et al.  The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. , 2010, Genome research.

[34]  C. Glass,et al.  Simple combinations of lineage-determining transcription factors prime cis-regulatory elements required for macrophage and B cell identities. , 2010, Molecular cell.

[35]  Qifa Zhang,et al.  Genome-wide association studies of 14 agronomic traits in rice landraces , 2010, Nature Genetics.

[36]  David Reich,et al.  Population differentiation as a test for selective sweeps. , 2010, Genome research.

[37]  Aaron R. Quinlan,et al.  Bioinformatics Applications Note Genome Analysis Bedtools: a Flexible Suite of Utilities for Comparing Genomic Features , 2022 .

[38]  J. Stewart,et al.  Physiology of cotton , 2010 .

[39]  D. Delmer,et al.  TRICHOME BIREFRINGENCE and Its Homolog AT5G01360 Encode Plant-Specific DUF231 Proteins Required for Cellulose Biosynthesis in Arabidopsis1[W][OA] , 2010, Plant Physiology.

[40]  K. Olsen,et al.  Genetic perspectives on crop domestication. , 2010, Trends in plant science.

[41]  Felix Krueger,et al.  Bismark: a flexible aligner and methylation caller for Bisulfite-Seq applications , 2011, Bioinform..

[42]  Thomas Zichner,et al.  DELLY: structural variant discovery by integrated paired-end and split-read analysis , 2012, Bioinform..

[43]  James C. Schnable,et al.  High-resolution mapping of open chromatin in the rice genome. , 2012, Genome research.

[44]  Shane J. Neph,et al.  An expansive human regulatory lexicon encoded in transcription factor footprints , 2012, Nature.

[45]  Steven L Salzberg,et al.  Fast gapped-read alignment with Bowtie 2 , 2012, Nature Methods.

[46]  Jeremy D. DeBarry,et al.  MCScanX: a toolkit for detection and evolutionary analysis of gene synteny and collinearity , 2012, Nucleic acids research.

[47]  Yong Zhang,et al.  Identifying ChIP-seq enrichment using MACS , 2012, Nature Protocols.

[48]  Xun Xu,et al.  Comparative population genomics of maize domestication and improvement , 2012, Nature Genetics.

[49]  Raymond K. Auerbach,et al.  Extensive Promoter-Centered Chromatin Interactions Provide a Topological Basis for Transcription Regulation , 2012, Cell.

[50]  Pablo Cingolani,et al.  © 2012 Landes Bioscience. Do not distribute. , 2022 .

[51]  Tao Zhang,et al.  Genome-Wide Identification of Regulatory DNA Elements and Protein-Binding Footprints Using Signatures of Open Chromatin in Arabidopsis[C][W][OA] , 2012, Plant Cell.

[52]  P. Wittkopp,et al.  Cis-regulatory elements: molecular mechanisms and evolutionary processes underlying divergence , 2011, Nature Reviews Genetics.

[53]  A. Fujiyama,et al.  A map of rice genome variation reveals the origin of cultivated rice , 2012, Nature.

[54]  W. Sung,et al.  Chromatin connectivity maps reveal dynamic promoter–enhancer long-range associations , 2013, Nature.

[55]  Y. Pei,et al.  The Dual Functions of WLIM1a in Cell Elongation and Secondary Wall Formation in Developing Cotton Fibers[C][W] , 2013, Plant Cell.

[56]  Lili Tu,et al.  A Genetic and Metabolic Analysis Revealed that Cotton Fiber Cell Development Was Retarded by Flavonoid Naringenin1[W][OA] , 2013, Plant Physiology.

[57]  V. Brunaud,et al.  ARABIDOPSIS THALIANA HOMEOBOX25 Uncovers a Role for Gibberellins in Seed Longevity1[C][W] , 2013, Plant Physiology.

[58]  J Crossa,et al.  Genomic prediction in CIMMYT maize and wheat breeding programs , 2013, Heredity.

[59]  Y. Pei,et al.  Transcriptome and Biochemical Analyses Revealed a Detailed Proanthocyanidin Biosynthesis Pathway in Brown Cotton Fiber , 2014, PloS one.

[60]  Jie Sun,et al.  Molecular analysis of proanthocyanidins related to pigmentation in brown cotton fibre (Gossypium hirsutum L.). , 2014, Journal of experimental botany.

[61]  William Stafford Noble,et al.  Statistical confidence estimation for Hi-C data reveals regulatory chromatin contacts , 2014, Genome research.

[62]  Zhongxu Lin,et al.  A comparative meta-analysis of QTL between intraspecific Gossypium hirsutum and interspecific G. hirsutum × G. barbadense populations , 2015, Molecular Genetics and Genomics.

[63]  R. Terauchi,et al.  Harvesting the Promising Fruits of Genomics: Applying Genome Sequencing Technologies to Crop Breeding , 2014, PLoS biology.

[64]  Jonathan F. Wendel,et al.  Ecology, Evolution and Organismal Biology Publications Ecology, Evolution and Organismal Biology Comparative Evolutionary and Developmental Dynamics of the Cotton (gossypium Hirsutum) Fiber Transcriptome Comparative Evolutionary and Developmental Dynamics of the Cotton (gossypium Hirsutum) Fiber Tra , 2022 .

[65]  Björn Usadel,et al.  Trimmomatic: a flexible trimmer for Illumina sequence data , 2014, Bioinform..

[66]  Kui Lin,et al.  Genomic analyses provide insights into the history of tomato breeding , 2014, Nature Genetics.

[67]  Hui Xiang,et al.  Resequencing 302 wild and cultivated accessions identifies genes related to domestication and improvement in soybean , 2015, Nature Biotechnology.

[68]  Z. Weng,et al.  Local sequence assembly reveals a high-resolution profile of somatic structural variations in 97 cancer genomes , 2015, Nucleic acids research.

[69]  He Zhang,et al.  Genome sequence of cultivated Upland cotton (Gossypium hirsutum TM-1) provides insights into genome evolution , 2015, Nature Biotechnology.

[70]  Xue-Bao Li,et al.  Cotton (Gossypium hirsutum) 14-3-3 proteins participate in regulation of fibre initiation and elongation by modulating brassinosteroid signalling. , 2015, Plant biotechnology journal.

[71]  Jean-Philippe Vert,et al.  HiC-Pro: an optimized and flexible pipeline for Hi-C data processing , 2015, Genome Biology.

[72]  J. T. Erichsen,et al.  Enhancer Evolution across 20 Mammalian Species , 2015, Cell.

[73]  Lei Fang,et al.  Sequencing of allotetraploid cotton (Gossypium hirsutum L. acc. TM-1) provides a resource for fiber improvement , 2015, Nature Biotechnology.

[74]  Damian Roqueiro,et al.  Genome-wide analysis of local chromatin packing in Arabidopsis thaliana , 2015, Genome research.

[75]  Jiming Jiang,et al.  Genome-wide mapping of DNase I hypersensitive sites in plants. , 2015, Methods in molecular biology.

[76]  Eric Haugen,et al.  Large-scale identification of sequence variants impacting human transcription factor occupancy in vivo , 2015, Nature Genetics.

[77]  Michael Freeling,et al.  Advances in understanding cis regulation of the plant gene with an emphasis on comparative genomics. , 2015, Current opinion in plant biology.

[78]  Qing-Yong Yang,et al.  De novo plant genome assembly based on chromatin interactions: a case study of Arabidopsis thaliana. , 2015, Molecular plant.

[79]  X. Zhou,et al.  TopDom: an efficient and deterministic method for identifying topological domains in genomes , 2015, Nucleic acids research.

[80]  Dongyuan Liu,et al.  Subgenome parallel selection is associated with morphotype diversification and convergent crop domestication in Brassica rapa and Brassica oleracea , 2016, Nature Genetics.

[81]  Maojun Wang,et al.  Multi-omics maps of cotton fibre reveal epigenetic basis for staged single-cell differentiation , 2016, Nucleic acids research.

[82]  Maojun Wang,et al.  Fibre elongation requires normal redox homeostasis modulated by cytosolic ascorbate peroxidase in cotton (Gossypium hirsutum) , 2016, Journal of experimental botany.

[83]  Maojun Wang,et al.  Genome-wide SSR-based association mapping for fiber quality in nation-wide upland cotton inbreed cultivars in China , 2016, BMC Genomics.