Genetic Analysis of the Transition from Wild to Domesticated Cotton (Gossypium hirsutum L.)

Genomic and Cytogenetic Analysis of Synthetic Polyploids between Diploid and Tetraploid Cotton (Gossypium) Species

Cotton (Gossypium spp.) is the most important natural fiber source in the world. The genetic potential of cotton can be successfully and efficiently exploited by identifying and solving the complex fundamental problems of systematics, evolution, and phylogeny, based on interspecific hybridization of cotton. This study describes the results of interspecific hybridization of G. herbaceum L. (A1-genome) and G. mustelinum Miers ex Watt (AD4-genome) species, obtaining fertile hybrids through synthetic polyploidization of otherwise sterile triploid forms with colchicine (C22H25NO6) treatment. The fertile F1C hybrids were produced from five different cross combinations: (1) G. herbaceum subsp. frutescens × G. mustelinum; (2) G. herbaceum subsp. pseudoarboreum × G. mustelinum; (3) G. herbaceum subsp. pseudoarboreum f. harga × G. mustelinum; (4) G. herbaceum subsp. africanum × G. mustelinum; (5) G. herbaceum subsp. euherbaceum (variety A-833) × G. mustelinum. Cytogenetic analysis discovered normal conjugation of bivalent chromosomes in addition to univalent, open, and closed ring-shaped quadrivalent chromosomes at the stage of metaphase I in the F1C and F2C hybrids. The setting of hybrid bolls obtained as a result of these crosses ranged from 13.8-92.2%, the fertility of seeds in hybrid bolls from 9.7-16.3%, and the pollen viability rates from 36.6-63.8%. Two transgressive plants with long fiber of 35.1-37.0 mm and one plant with extra-long fiber of 39.1-41.0 mm were identified in the F2C progeny of G. herbaceum subsp. frutescens × G. mustelinum cross. Phylogenetic analysis with 72 SSR markers that detect genomic changes showed that tetraploid hybrids derived from the G. herbaceum × G. mustelinum were closer to the species G. mustelinum. The G. herbaceum subsp. frutescens was closer to the cultivated form, and its subsp. africanum was closer to the wild form. New knowledge of the interspecific hybridization and synthetic polyploidization was developed for understanding the genetic mechanisms of the evolution of tetraploid cotton during speciation. The synthetic polyploids of cotton obtained in this study would provide beneficial genes for developing new cotton varieties of the G. hirsutum species, with high-quality cotton fiber and strong tolerance to biotic or abiotic stress. In particular, the introduction of these polyploids to conventional and molecular breeding can serve as a bridge of transferring valuable genes related to high-quality fiber and stress tolerance from different cotton species to the new cultivars.

Breeding effects on durum wheat traits detected using GWAS and haplotype block analysis

Introduction

The recent boosting of genomic data in durum wheat (Triticum turgidum subsp. durum) offers the opportunity to better understand the effects of breeding on the genetic structures that regulate the expression of traits of agronomic interest. Furthermore, the identification of DNA markers useful for marker-assisted selection could also improve the reliability of technical protocols used for variety protection and registration.

Methods

Within this motivation context, 123 durum wheat accessions, classified into three groups: landraces (LR), ancient (OC) and modern cultivars (MC), were evaluated in two locations, for 34 agronomic traits, including UPOV descriptors, to assess the impact of changes that occurred during modern breeding.

Results

The association mapping analysis, performed with 4,241 SNP markers and six multi-locus-GWAS models, revealed 28 reliable Quantitative Trait Nucleotides (QTNs) related to plant morphology and kernel-related traits. Some important genes controlling flowering time and plant height were in linkage disequilibrium (LD) decay with QTNs identified in this study. A strong association for yellow berry was found on chromosome 6A (Q.Yb-6A) in a region containing the nadh-ubiquinone oxidoreductase subunit, a gene involved in starch metabolism. The Q.Kcp-2A harbored the PPO locus, with the associated marker (Ku_c13700_1196) in LD decay with Ppo-A1 and Ppo-A2. Interestingly, the Q.FGSGls-2B.1, identified by RAC875_c34512_685 for flag leaf glaucosity, mapped less than 1 Mb from the Epistatic inhibitors of glaucousness (Iw1), thus representing a good candidate for supporting the morphological DUS traits also with molecular markers. LD haplotype block approach revealed a higher diversity, richness and length of haploblocks in MC than OC and LR (580 in LR, 585 in OC and 612 in MC), suggesting a possible effect exerted by breeding programs on genomic regions associated with the agronomic traits.

Discussion

Our findings pave new ways to support the phenotypic characterization necessary for variety registration by using a panel of cost-effectiveness SNP markers associated also to the UPOV descriptors. Moreover, the panel of associated SNPs might represent a reservoir of favourable alleles to use in durum wheat breeding and genetics.

Domestication over Speciation in Allopolyploid Cotton Species: A Stronger Transcriptomic Pull

Cotton has been domesticated independently four times for its fiber, but the genomic targets of selection during each domestication event are mostly unknown. Comparative analysis of the transcriptome during cotton fiber development in wild and cultivated materials holds promise for revealing how independent domestications led to the superficially similar modern cotton fiber phenotype in upland (G. hirsutum) and Pima (G. barbadense) cotton cultivars. Here we examined the fiber transcriptomes of both wild and domesticated G. hirsutum and G. barbadense to compare the effects of speciation versus domestication, performing differential gene expression analysis and coexpression network analysis at four developmental timepoints (5, 10, 15, or 20 days after flowering) spanning primary and secondary wall synthesis. These analyses revealed extensive differential expression between species, timepoints, domestication states, and particularly the intersection of domestication and species. Differential expression was higher when comparing domesticated accessions of the two species than between the wild, indicating that domestication had a greater impact on the transcriptome than speciation. Network analysis showed significant interspecific differences in coexpression network topology, module membership, and connectivity. Despite these differences, some modules or module functions were subject to parallel domestication in both species. Taken together, these results indicate that independent domestication led G. hirsutum and G. barbadense down unique pathways but that it also leveraged similar modules of coexpression to arrive at similar domesticated phenotypes.

Multiple domestication events explain the origin of Gossypium hirsutum landraces in Mexico

Several Mesoamerican crops constitute wild-to-domesticated complexes generated by multiple initial domestication events, and continuous gene flow among crop populations and between these populations and their wild relatives. It has been suggested that the domestication of cotton (Gossypium hirsutum) started in the northwest of the Yucatán Peninsula, from where it spread to other regions inside and outside of Mexico. We tested this hypothesis by assembling chloroplast genomes of 23 wild, landraces, and breeding lines (transgene-introgressed and conventional). The phylogenetic analysis showed that the evolutionary history of cotton in Mexico involves multiple events of introgression and genetic divergence. From this, we conclude that Mexican landraces arose from multiple wild populations. Our results also revealed that their structural and functional chloroplast organizations had been preserved. However, genetic diversity decreases as a consequence of domestication, mainly in transgene-introgressed (TI) individuals (π = 0.00020, 0.00001, 0.00016, 0, and 0, of wild, TI-wild, landraces, TI-landraces, and breeding lines, respectively). We identified homologous regions that differentiate wild from domesticated plants and indicate a relationship among the samples. A decrease in genetic diversity associated with transgene introgression in cotton was identified for the first time, and our outcomes are therefore relevant to both biosecurity and agrobiodiversity conservation.

The role of non-additive gene action on gene expression variation in plant domestication

BACKGROUND

Plant domestication is a remarkable example of rapid phenotypic transformation of polygenic traits, such as organ size. Evidence from a handful of study cases suggests this transformation is due to gene regulatory changes that result in non-additive phenotypes. Employing data from published genetic crosses, we estimated the role of non-additive gene action in the modulation of transcriptional landscapes in three domesticated plants: maize, sunflower, and chili pepper. Using A. thaliana, we assessed the correlation between gene regulatory network (GRN) connectivity properties, transcript abundance variation, and gene action. Finally, we investigated the propagation of non-additive gene action in GRNs.

RESULTS

We compared crosses between domesticated plants and their wild relatives to a set of control crosses that included a pair of subspecies evolving under natural selection and a set of inbred lines evolving under domestication. We found abundance differences on a higher portion of transcripts in crosses between domesticated-wild plants relative to the control crosses. These transcripts showed non-additive gene action more often in crosses of domesticated-wild plants than in our control crosses. This pattern was strong for genes associated with cell cycle and cell fate determination, which control organ size. We found weak but significant negative correlations between the number of targets of trans-acting genes (Out-degree) and both the magnitude of transcript abundance difference a well as the absolute degree of dominance. Likewise, we found that the number of regulators that control a gene's expression (In-degree) is weakly but negatively correlated with the magnitude of transcript abundance differences. We observed that dominant-recessive gene action is highly propagable through GRNs. Finally, we found that transgressive gene action is driven by trans-acting regulators showing additive gene action.

CONCLUSIONS

Our study highlights the role of non-additive gene action on modulating domestication-related traits, such as organ size via regulatory divergence. We propose that GRNs are shaped by regulatory changes at genes with modest connectivity, which reduces the effects of antagonistic pleiotropy. Finally, we provide empirical evidence of the propagation of non-additive gene action in GRNs, which suggests a transcriptional epistatic model for the control of polygenic traits, such as organ size.

Dual Domestication, Diversity, and Differential Introgression in Old World Cotton Diploids

Domestication in the cotton genus is remarkable in that it has occurred independently four different times at two different ploidy levels. Relatively little is known about genome evolution and domestication in the cultivated diploid species Gossypium herbaceum and Gossypium arboreum, due to the absence of wild representatives for the latter species, their ancient domestication, and their joint history of human-mediated dispersal and interspecific gene flow. Using in-depth resequencing of a broad sampling from both species, we provide support for their independent domestication, as opposed to a progenitor-derivative relationship, showing that diversity (mean π = 6 × 10-3) within species is similar, and that divergence between species is modest (FST = 0.413). Individual accessions were homozygous for ancestral single-nucleotide polymorphisms at over half of variable sites, while fixed, derived sites were at modest frequencies. Notably, two chromosomes with a paucity of fixed, derived sites (i.e., chromosomes 7 and 10) were also strongly implicated as having experienced high levels of introgression. Collectively, these data demonstrate variable permeability to introgression among chromosomes, which we propose is due to divergent selection under domestication and/or the phenomenon of F2 breakdown in interspecific crosses. Our analyses provide insight into the evolutionary forces that shape diversity and divergence in the diploid cultivated species and establish a foundation for understanding the contribution of introgression and/or strong parallel selection to the extensive morphological similarities shared between species.

Status and prospects of genome-wide association studies in cotton

Over the last two decades, the use of high-density SNP arrays and DNA sequencing have allowed scientists to uncover the majority of the genotypic space for various crops, including cotton. Genome-wide association study (GWAS) links the dots between a phenotype and its underlying genetics across the genomes of populations. It was first developed and applied in the field of human disease genetics. Many areas of crop research have incorporated GWAS in plants and considerable literature has been published in the recent decade. Here we will provide a comprehensive review of GWAS studies in cotton crop, which includes case studies on biotic resistance, abiotic tolerance, fiber yield and quality traits, current status, prospects, bottlenecks of GWAS and finally, thought-provoking question. This review will serve as a catalog of GWAS in cotton and suggest new frontiers of the cotton crop to be studied with this important tool.

The Current Progresses in the Genes and Networks Regulating Cotton Plant Architecture

Cotton is the most important source of natural fiber in the world as well as a key source of edible oil. The plant architecture and flowering time in cotton are crucial factors affecting cotton yield and the efficiency of mechanized harvest. In the model plant arabidopsis, the functions of genes related to plant height, inflorescence structure, and flowering time have been well studied. In the model crops, such as tomato and rice, the similar genetic explorations have greatly strengthened the economic benefits of these crops. Plants of the Gossypium genus have the characteristics of perennials with indeterminate growth and the cultivated allotetraploid cottons, G. hirsutum (Upland cotton), and G. barbadense (Sea-island cotton), have complex branching patterns. In this paper, we review the current progresses in the identification of genes affecting cotton architecture and flowering time in the cotton genome and the elucidation of their functional mechanisms associated with branching patterns, branching angle, fruit branch length, and plant height. This review focuses on the following aspects: (i) plant hormone signal transduction pathway; (ii) identification of cotton plant architecture QTLs and PEBP gene family members; (iii) functions of FT/SFT and SP genes; (iv) florigen and anti-florigen systems. We highlight areas that require further research, and should lay the groundwork for the targeted bioengineering of improved cotton cultivars with flowering times, plant architecture, growth habits and yields better suited for modern, mechanized cultivation.

Parallel and Intertwining Threads of Domestication in Allopolyploid Cotton

The two cultivated allopolyploid cottons, Gossypium hirsutum and Gossypium barbadense, represent a remarkable example of parallel independent domestication, both involving dramatic morphological transformations under selection from wild perennial plants to annualized row crops. Deep resequencing of 643 newly sampled accessions spanning the wild-to-domesticated continuum of both species, and their allopolyploid relatives, are combined with existing data to resolve species relationships and elucidate multiple aspects of their parallel domestication. It is confirmed that wild G. hirsutum and G. barbadense were initially domesticated in the Yucatan Peninsula and NW South America, respectively, and subsequently spread under domestication over 4000-8000 years to encompass most of the American tropics. A robust phylogenomic analysis of infraspecific relationships in each species is presented, quantify genetic diversity in both, and describe genetic bottlenecks associated with domestication and subsequent diffusion. As these species became sympatric over the last several millennia, pervasive genome-wide bidirectional introgression occurred, often with striking asymmetries involving the two co-resident genomes of these allopolyploids. Diversity scans revealed genomic regions and genes unknowingly targeted during domestication and additional subgenomic asymmetries. These analyses provide a comprehensive depiction of the origin, divergence, and adaptation of cotton, and serve as a rich resource for cotton improvement.

Cotton architecture: examining the roles of SINGLE FLOWER TRUSS and SELF-PRUNING in regulating growth habits of a woody perennial crop

By specifying patterns of determinate and indeterminate growth, FLOWERING LOCUS T/SINGLE FLOWER TRUSS (SFT) and TERMINAL FLOWER 1/SELF-PRUNING (SP) regulate plant architecture. Though well characterized in Arabidopsis, the impacts of these genes on the architectures of diverse crops cultivated in different environments, and their potential to enhance crop productivity and management, are less well addressed. Cotton (Gossypium spp.) is naturally a short-day photoperiodic perennial that is now grown primarily as a day-neutral, annual row crop. Different environments and cultivation practices favor specific growth habits to optimize yield, and in cotton, especially in regions that rely heavily on mechanized harvest, the trend has been to more determinate varieties. Identifying and functionally characterizing SFT and SP homologs in cotton uncovered new aspects of how ratios of indeterminate and determinate growth are balanced, and unraveling their genetic networks emphasized how broadly these gene products affect cotton growth habits.

Conservation and Divergence in Duplicated Fiber Coexpression Networks Accompanying Domestication of the Polyploid Gossypium hirsutum L

Gossypium hirsutum L. (Upland cotton) has an evolutionary history involving inter-genomic hybridization, polyploidization, and subsequent domestication. We analyzed the developmental dynamics of the cotton fiber transcriptome accompanying domestication using gene coexpression networks for both joint and homoeologous networks. Remarkably, most genes exhibited expression for at least one homoeolog, confirming previous reports of widespread gene usage in cotton fibers. Most coexpression modules comprising the joint network are preserved in each subgenomic network and are enriched for similar biological processes, showing a general preservation of network modular structure for the two co-resident genomes in the polyploid. Interestingly, only one fifth of homoeologs co-occur in the same module when separated, despite similar modular structures between the joint and homoeologous networks. These results suggest that the genome-wide divergence between homoeologous genes is sufficient to separate their co-expression profiles at the intermodular level, despite conservation of intramodular relationships within each subgenome. Most modules exhibit D-homoeolog expression bias, although specific modules do exhibit A-homoeolog bias. Comparisons between wild and domesticated coexpression networks revealed a much tighter and denser network structure in domesticated fiber, as evidenced by its fewer modules, 13-fold increase in the number of development-related module member genes, and the poor preservation of the wild network topology. These results demonstrate the amazing complexity that underlies the domestication of cotton fiber.