Evolutionary Dynamics of Chromatin Structure and Duplicate Gene Expression in Diploid and Allopolyploid Cotton.

Polyploidy is a prominent mechanism of plant speciation and adaptation, yet the mechanistic understandings of duplicated gene regulation remain elusive. Chromatin structure dynamics are suggested to govern gene regulatory control. Here, we characterized genome-wide nucleosome organization and chromatin accessibility in allotetraploid cotton, Gossypium hirsutum (AADD, 2n = 4X = 52), relative to its two diploid parents (AA or DD genome) and their synthetic diploid hybrid (AD), using DNS-seq. The larger A-genome exhibited wider average nucleosome spacing in diploids, and this intergenomic difference diminished in the allopolyploid but not hybrid. Allopolyploidization also exhibited increased accessibility at promoters genome-wide and synchronized cis-regulatory motifs between subgenomes. A prominent cis-acting control was inferred for chromatin dynamics and demonstrated by transposable element removal from promoters. Linking accessibility to gene expression patterns, we found distinct regulatory effects for hybridization and later allopolyploid stages, including nuanced establishment of homoeolog expression bias and expression level dominance. Histone gene expression and nucleosome organization are coordinated through chromatin accessibility. Our study demonstrates the capability to track high-resolution chromatin structure dynamics and reveals their role in the evolution of cis-regulatory landscapes and duplicate gene expression in polyploids, illuminating regulatory ties to subgenomic asymmetry and dominance.

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.

Insights into the Ecology and Evolution of Polyploid Plants through Network Analysis.

Polyploidy is a widespread phenomenon throughout eukaryotes, with important ecological and evolutionary consequences. Although genes operate as components of complex pathways and networks, polyploid changes in genes and gene expression have typically been evaluated as either individual genes or as a part of broad-scale analyses. Network analysis has been fruitful in associating genomic and other 'omic'-based changes with phenotype for many systems. In polyploid species, network analysis has the potential not only to facilitate a better understanding of the complex 'omic' underpinnings of phenotypic and ecological traits common to polyploidy, but also to provide novel insight into the interaction among duplicated genes and genomes. This adds perspective to the global patterns of expression (and other 'omic') change that accompany polyploidy and to the patterns of recruitment and/or loss of genes following polyploidization. While network analysis in polyploid species faces challenges common to other analyses of duplicated genomes, present technologies combined with thoughtful experimental design provide a powerful system to explore polyploid evolution. Here, we demonstrate the utility and potential of network analysis to questions pertaining to polyploidy with an example involving evolution of the transgressively superior cotton fibres found in polyploid Gossypium hirsutum. By combining network analysis with prior knowledge, we provide further insights into the role of profilins in fibre domestication and exemplify the potential for network analysis in polyploid species.

Re-evaluating the phylogeny of allopolyploid Gossypium L.

The formation of allopolyploid cotton precipitated a rapid diversification and colonization of dry coastal American tropical and subtropical regions. Previous phylogenetic analyses, combined with molecular divergence analyses, have offered a temporal framework for this radiation, but provide only weak support for some of the resolved branches. Moreover, these earlier analyses did not include the recently recognized sixth polyploid species, G. ekmanianum Wittmack. Here we use targeted sequence capture of multiple loci in conjunction with both concatenated and Bayesian concordance analyses to reevaluate the phylogeny of allopolyploid cotton species. Although phylogenetic resolution afforded by individual genes is often low, sufficient signal was attained both through the concatenated and concordance analyses to provide robust support for the Gossypium polyploid clade, which is reported here.

Persistence of subgenomes in paleopolyploid cotton after 60 my of evolution.

The importance of whole-genome multiplication (WGM) in plant evolution has long been recognized. In flowering plants, WGM is both ubiquitous and in many lineages cyclical, each round followed by substantial gene loss (fractionation). This process may be biased with respect to duplicated chromosomes, often with overexpression of genes in less fractionated relative to more fractionated regions. This bias is hypothesized to arise through downregulation of gene expression through silencing of local transposable elements (TEs). We assess differences in gene expression between duplicated regions of the paleopolyploid cotton genome and demonstrate that the rate of fractionation is negatively correlated with gene expression. We examine recent hypotheses regarding the source of fractionation bias and show that TE-mediated, positional downregulation is absent in the modern cotton genome, seemingly excluding this phenomenon as the primary driver of biased gene loss. Nevertheless, the paleo subgenomes of diploid cotton are still distinguishable with respect to TE content, targeting of 24-nt-small interfering RNAs and GC content, despite approximately 60 My of evolution. We propose that repeat content per se and differential recombination rates may drive biased fractionation following WGM. These data highlight the likely importance of ancient genomic fractionation biases in shaping modern crop genomes.

A Transcriptome Profile for Developing Seed of Polyploid Cotton.

Cotton ranks among the world's important oilseed crops, yet relative to other oilseeds there are few studies of oil-related biosynthetic and regulatory pathways. We present global transcriptome analyses of cotton seed development using RNA-seq and four developmental time-points. Because Upland cotton (Gossypium hirsutum L.) is an allopolyploid containing two genomes (A/D), we partitioned expression into the individual contributions of each homeologous gene copy. Data were explored with respect to genic and subgenomic patterns of expression, globally and with respect to seed pathways and networks. The most dynamic period of transcriptome change is from 20-30 d postanthesis (DPA), with about 20% of genes showing homeolog expression bias. Co-expression analysis shows largely congruent homeolog networks, but also homeolog-specific divergence. Functional enrichment tests show that flavonoid biosynthesis and lipid related genes were significantly represented early and later in seed development, respectively. An involvement of new features in oil biosynthesis was found, like the contribution of DGAT3 (diacylglycerol acyltransferase) to the total triglyceride expression pool. Also, catechin-based and epicatechin-based proanthocyanidin expression are reciprocally biased with respect to homeolog usage. This study provides the first temporal analysis of duplicated gene expression in cotton seed and a resource for understanding new aspects of oil and flavonoid biosynthetic processes.

CenH3 evolution in diploids and polyploids of three angiosperm genera.

BACKGROUND: Centromeric DNA sequences alone are neither necessary nor sufficient for centromere specification. The centromere specific histone, CenH3, evolves rapidly in many species, perhaps as a coevolutionary response to rapidly evolving centromeric DNA. To gain insight into CenH3 evolution, we characterized patterns of nucleotide and protein diversity among diploids and allopolyploids within three diverse angiosperm genera, Brassica, Oryza, and Gossypium (cotton), with a focus on evidence for diversifying selection in the various domains of the CenH3 gene. In addition, we compare expression profiles and alternative splicing patterns for CenH3 in representatives of each genus. RESULTS: All three genera retain both duplicated CenH3 copies, while Brassica and Gossypium exhibit pronounced homoeologous expression level bias. Comparisons among genera reveal shared and unique aspects of CenH3 evolution, variable levels of diversifying selection in different CenH3 domains, and that alternative splicing contributes significantly to CenH3 diversity. CONCLUSIONS: Since the N terminus is subject to diversifying selection but the DNA binding domains do not appear to be, rapidly evolving centromere sequences are unlikely to be the primary driver of CenH3 sequence diversification. At present, the functional explanation for the diversity generated by both conventional protein evolution in the N terminal domain, as well as alternative splicing, remains unexplained.

An assessment of karyotype restructuring in the neoallotetraploid Tragopogon miscellus (Asteraceae).

Tragopogon miscellus and Tragopogon mirus are two rare examples of allopolyploids that have formed recently in nature. Molecular cytogenetic studies have revealed chromosome copy number variation and intergenomic translocations in both allotetraploids. Due to a lack of interstitial chromosome markers, there remained the possibility of additional karyotype restructuring in these neopolyploids, via intrachromosomal and intragenomic rearrangements. To address this issue, we searched for additional high-copy tandem repeats in genomic sequences of the diploid progenitor species-Tragopogon dubius, Tragopogon pratensis and Tragopogon porrifolius-for application to the chromosomes of the allotetraploids. Eight novel repeats were localised by fluorescence in situ hybridisation (FISH) in the diploids; one of these repeats, TTR3, provided interstitial coverage. TTR3 was included in a cocktail with other previously characterised probes, producing better-resolved karyotypes for the three diploids. The cocktail was then used to test a hypothesis of karyotype restructuring in the recent allotetraploid T. miscellus by comparing repeat distributions to its diploid progenitors, T. dubius and T. pratensis. Five individuals of T. miscellus were selected from across the range of karyotypic variation previously observed in natural populations. FISH signal distributions mostly matched those observed in the diploid progenitors, with the exception of several losses or gains of signal at chromosomal subtermini and previously noted intergenomic translocations. Thus, in T. miscellus, we find most changes restricted to the subterminal regions, and although some were recurrent, none of the changes were common to all individuals analysed. We consider these findings in relation to the gene loss reported previously for T. miscellus.

Extensive chromosomal variation in a recently formed natural allopolyploid species, Tragopogon miscellus (Asteraceae).

Polyploidy, or whole genome duplication, has played a major role in the evolution of many eukaryotic lineages. Although the prevalence of polyploidy in plants is well documented, the molecular and cytological consequences are understood largely from newly formed polyploids (neopolyploids) that have been grown experimentally. Classical cytological and molecular cytogenetic studies both have shown that experimental neoallopolyploids often have meiotic irregularities, producing chromosomally variable gametes and progeny; however, little is known about the extent or duration of chromosomal variation in natural neoallopolyploid populations. We report the results of a molecular cytogenetic study on natural populations of a neoallopolyploid, Tragopogon miscellus, which formed multiple times in the past 80 y. Using genomic and fluorescence in situ hybridization, we uncovered massive and repeated patterns of chromosomal variation in all populations. No population was fixed for a particular karyotype; 76% of the individuals showed intergenomic translocations, and 69% were aneuploid for one or more chromosomes. Importantly, 85% of plants exhibiting aneuploidy still had the expected chromosome number, mostly through reciprocal monosomy-trisomy of homeologous chromosomes (1:3 copies) or nullisomy-tetrasomy (0:4 copies). The extensive chromosomal variation still present after ca. 40 generations in this biennial species suggests that substantial and prolonged chromosomal instability might be common in natural populations after whole genome duplication. A protracted period of genome instability in neoallopolyploids may increase opportunities for alterations to genome structure, losses of coding and noncoding DNA, and changes in gene expression.