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.

GRASSY TILLERS1 (GT1) and SIX-ROWED SPIKE1 (VRS1) homologs share conserved roles in growth repression.

Crop engineering and de novo domestication using gene editing are new frontiers in agriculture. However, outside of well-studied crops and model systems, prioritizing engineering targets remains challenging. Evolution can guide us, revealing genes with deeply conserved roles that have repeatedly been selected in the evolution of plant form. Homologs of the transcription factor genes GRASSY TILLERS1 (GT1) and SIX-ROWED SPIKE1 (VRS1) have repeatedly been targets of selection in domestication and evolution, where they repress growth in many developmental contexts. This suggests a conserved role for these genes in regulating growth repression. To test this, we determined the roles of GT1 and VRS1 homologs in maize (Zea mays) and the distantly related grass brachypodium (Brachypodium distachyon) using gene editing and mutant analysis. In maize, gt1; vrs1-like1 (vrl1) mutants have derepressed growth of floral organs. In addition, gt1; vrl1 mutants bore more ears and more branches, indicating broad roles in growth repression. In brachypodium, Bdgt1; Bdvrl1 mutants have more branches, spikelets, and flowers than wild-type plants, indicating conserved roles for GT1 and VRS1 homologs in growth suppression over ca. 59 My of grass evolution. Importantly, many of these traits influence crop productivity. Notably, maize GT1 can suppress growth in arabidopsis (Arabidopsis thaliana) floral organs, despite ca. 160 My of evolution separating the grasses and arabidopsis. Thus, GT1 and VRS1 maintain their potency as growth regulators across vast timescales and in distinct developmental contexts. This work highlights the power of evolution to inform gene editing in crop improvement.

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.

Structural evolution drives diversification of the large LRR-RLK gene family.

●Cells are continuously exposed to chemical signals that they must discriminate between and respond to appropriately. In embryophytes, the leucine-rich repeat receptor-like kinases (LRR-RLKs) are signal receptors critical in development and defense. LRR-RLKs have diversified to hundreds of genes in many plant genomes. Although intensively studied, a well-resolved LRR-RLK gene tree has remained elusive. ●To resolve the LRR-RLK gene tree, we developed an improved gene discovery method based on iterative hidden Markov model searching and phylogenetic inference. We used this method to infer complete gene trees for each of the LRR-RLK subclades and reconstructed the deepest nodes of the full gene family. ●We discovered that the LRR-RLK gene family is even larger than previously thought, and that protein domain gains and losses are prevalent. These structural modifications, some of which likely predate embryophyte diversification, led to misclassification of some LRR-RLK variants as members of other gene families. Our work corrects this misclassification. ●Our results reveal ongoing structural evolution generating novel LRR-RLK genes. These new genes are raw material for the diversification of signaling in development and defense. Our methods also enable phylogenetic reconstruction in any large gene family.

Reticulate Evolution Helps Explain Apparent Homoplasy in Floral Biology and Pollination in Baobabs (Adansonia; Bombacoideae; Malvaceae).

Baobabs (Adansonia) are a cohesive group of tropical trees with a disjunct distribution in Australia, Madagascar, and continental Africa, and diverse flowers associated with two pollination modes. We used custom-targeted sequence capture in conjunction with new and existing phylogenetic comparative methods to explore the evolution of floral traits and pollination systems while allowing for reticulate evolution. Our analyses suggest that relationships in Adansonia are confounded by reticulation, with network inference methods supporting at least one reticulation event. The best supported hypothesis involves introgression between Adansonia rubrostipa and core Longitubae, both of which are hawkmoth pollinated with yellow/red flowers, but there is also some support for introgression between the African lineage and Malagasy Brevitubae, which are both mammal-pollinated with white flowers. New comparative methods for phylogenetic networks were developed that allow maximum-likelihood inference of ancestral states and were applied to study the apparent homoplasy in floral biology and pollination mode seen in Adansonia. This analysis supports a role for introgressive hybridization in morphological evolution even in a clade with highly divergent and geographically widespread species. Our new comparative methods for discrete traits on species networks are implemented in the software PhyloNetworks. [Comparative methods; Hyb-Seq; introgression; network inference; population trees; reticulate evolution; species tree inference; targeted sequence capture.].

The CLAVATA receptor FASCIATED EAR2 responds to distinct CLE peptides by signaling through two downstream effectors.

Meristems contain groups of indeterminate stem cells, which are maintained by a feedback loop between CLAVATA (CLV) and WUSCHEL (WUS) signaling. CLV signaling involves the secretion of the CLV3 peptide and its perception by a number of Leucine-Rich-Repeat (LRR) receptors, including the receptor-like kinase CLV1 and the receptor-like protein CLV2 coupled with the CORYNE (CRN) pseudokinase. CLV2, and its maize ortholog FASCIATED EAR2 (FEA2) appear to function in signaling by CLV3 and several related CLV3/EMBRYO-SURROUNDING REGION (CLE) peptide ligands. Nevertheless, how signaling specificity is achieved remains unknown. Here we show that FEA2 transmits signaling from two distinct CLE peptides, the maize CLV3 ortholog ZmCLE7 and ZmFON2-LIKE CLE PROTEIN1 (ZmFCP1) through two different candidate downstream effectors, the alpha subunit of the maize heterotrimeric G protein COMPACT PLANT2 (CT2), and ZmCRN. Our data provide a novel framework to understand how diverse signaling peptides can activate different downstream pathways through common receptor proteins.

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.

Candidate Gene Identification of Flowering Time Genes in Cotton.

Flowering time control is critically important to all sexually reproducing angiosperms in both natural ecological and agronomic settings. Accordingly, there is much interest in defining the genes involved in the complex flowering-time network and how these respond to natural and artificial selection, the latter often entailing transitions in day-length responses. Here we describe a candidate gene analysis in the cotton genus Gossypium, which uses homologs from the well-described Arabidopsis flowering network to bioinformatically and phylogenetically identify orthologs in the published genome sequence from G. raimondii Ulbr., one of the two model diploid progenitors of the commercially important allopolyploid cottons, G. hirsutum L. and G. barbadense L. Presence and patterns of expression were evaluated from 13 aboveground tissues related to flowering for each of the candidate genes using allopolyploid G. hirsutum as a model. Furthermore, we use a comparative context to determine copy number variability of each key gene family across 10 published angiosperm genomes. Data suggest a pattern of repeated loss of duplicates following ancient whole-genome doubling events in diverse lineages. The data presented here provide a foundation for understanding both the parallel evolution of day-length neutrality in domesticated cottons and the flowering-time network, in general, in this important crop plant.

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.

Ancient gene duplicates in Gossypium (cotton) exhibit near-complete expression divergence.

Whole genome duplication (WGD) is widespread in flowering plants and is a driving force in angiosperm diversification. The redundancy introduced by WGD allows the evolution of novel gene interactions and functions, although the patterns and processes of diversification are poorly understood. We identified ∼ 2,000 pairs of paralogous genes in Gossypium raimondii (cotton) resulting from an approximately 60 My old 5- to 6-fold ploidy increase. Gene expression analyses revealed that, in G. raimondii, 99.4% of the gene pairs exhibit differential expression in at least one of the three tissues (petal, leaf, and seed), with 93% to 94% exhibiting differential expression on a per-tissue basis. For 1,666 (85%) pairs, differential expression was observed in all tissues. These observations were mirrored in a time series of G. raimondii seed, and separately in leaf, petal, and seed of G. arboreum, indicating expression level diversification before species divergence. A generalized linear model revealed 92.4% of the paralog pairs exhibited expression divergence, with most exhibiting significant gene and tissue interactions indicating complementary expression patterns in different tissues. These data indicate massive, near-complete expression level neo- and/or subfunctionalization among ancient gene duplicates, suggesting these processes are essential in their maintenance over ∼ 60 Ma.

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.