Emmanuel Szadkowski

Repeated polyploidization of Gossypium genomes and the evolution of spinnable cotton fibres

Polyploidy often confers emergent properties, such as the higher fibre productivity and quality of tetraploid cottons than diploid cottons bred for the same environments. Here we show that an abrupt five- to sixfold ploidy increase approximately 60 million years (Myr) ago, and allopolyploidy reuniting divergent Gossypium genomes approximately 1–2 Myr ago, conferred about 30–36-fold duplication of ancestral angiosperm (flowering plant) genes in elite cottons (Gossypium hirsutum and Gossypium barbadense), genetic complexity equalled only by Brassica among sequenced angiosperms. Nascent fibre evolution, before allopolyploidy, is elucidated by comparison of spinnable-fibred Gossypium herbaceum A and non-spinnable Gossypium longicalyx F genomes to one another and the outgroup D genome of non-spinnable Gossypium raimondii. The sequence of a G. hirsutum AtDt (in which ‘t’ indicates tetraploid) cultivar reveals many non-reciprocal DNA exchanges between subgenomes that may have contributed to phenotypic innovation and/or other emergent properties such as ecological adaptation by polyploids. Most DNA-level novelty in G. hirsutum recombines alleles from the D-genome progenitor native to its New World habitat and the Old World A-genome progenitor in which spinnable fibre evolved. Coordinated expression changes in proximal groups of functionally distinct genes, including a nuclear mitochondrial DNA block, may account for clusters of cotton-fibre quantitative trait loci affecting diverse traits. Opportunities abound for dissecting emergent properties of other polyploids, particularly angiosperms, by comparison to diploid progenitors and outgroups.

Homoeolog expression bias and expression level dominance in allopolyploid cotton

Allopolyploidy is an evolutionary and mechanistically intriguing process, in that it entails the reconciliation of two or more sets of diverged genomes and regulatory interactions. In this study, we explored gene expression patterns in interspecific hybrid F1, and synthetic and natural allopolyploid cotton using RNA-Seq reads from leaf transcriptomes. We determined how the extent and direction of expression level dominance (total level of expression for both homoeologs) and homoeolog expression bias (relative contribution of homoeologs to the transcriptome) changed from hybridization through evolution at the polyploid level and following cotton domestication. Genome-wide expression level dominance was biased toward the A-genome in the diploid hybrid and natural allopolyploids, whereas the direction was reversed in the synthetic allopolyploid. This biased expression level dominance was mainly caused by up- or downregulation of the homoeolog from the ‘non-dominant’ parent. Extensive alterations in homoeolog expression bias and expression level dominance accompany the initial merger of two diverged diploid genomes, suggesting a combination of regulatory (cis or trans) and epigenetic interactions that may arise and propagate through the transcriptome network. The extent of homoeolog expression bias and expression level dominance increases over time, from genome merger through evolution at the polyploid level. Higher rates of transgressive and novel gene expression patterns as well as homoeolog silencing were observed in natural allopolyploids than in F1 hybrid and synthetic allopolyploid cottons. These observations suggest that natural selection reconciles the regulatory mismatches caused by initial genomic merger, while new gene expression conditions are generated for evaluation by selection.

The first meiosis of resynthesized Brassica napus, a genome blender

Polyploidy promotes the restructuring of merged genomes within initial generations of resynthesized Brassica napus, possibly caused by homoeologous recombination at meiosis. However, little is known about the impact of the first confrontation of two genomes at the first meiosis which could lead to genome exchanges in progeny. Here, we assessed the role of the first meiosis in the genome instability of synthetic B. napus. We used three different newly resynthesized B. napus plants and established meiotic pairing frequencies for the A and C genomes. We genotyped the three corresponding progenies in a cross to a natural B. napus on the two homoeologous A1 and C1 chromosomes. Pairing at meiosis in a set of progenies with various rearrangements was scored. Here, we confirmed that the very first meiosis of resynthesized plants of B. napus acts as a genome blender, with many of the meiotic-driven genetic changes transmitted to the progenies, in proportions that depend significantly on the cytoplasm background inherited from the progenitors. We conclude that the first meiosis generates rearrangements on both genomes and promotes subsequent restructuring in further generations. Our study advances the knowledge on the timing of genetic changes and the mechanisms that may bias their transmission.

Homoeolog expression bias and expression level dominance in allopolyploids

Polyploidy is now recognized as a characteristic feature of all angiosperm genomes (Jiao et al., 2011), and remains an important speciation process today (Wendel, 2000; Comai, 2005; Doyle et al., 2008; Leitch & Leitch, 2008; Soltis & Soltis, 2009; Soltis et al., 2010). In allopolyploids, genomic merger and doubling are associated with myriad non-Mendelian interactions and processes, including sequence elimination (Shaked et al., 2001; Ozkan et al., 2003; Han et al., 2005; Skalicka et al., 2005; Anssour et al., 2009; Tate et al., 2009; Jackson & Chen, 2010), alterations of epigenetic marks (Shaked et al., 2001;Madlung et al., 2002; Rapp&Wendel, 2005; Chen, 2007; Doyle et al., 2008; Kovarik et al., 2008b; Soltis & Soltis, 2009; Soltis et al., 2010), activation of genes and retroelements (O’Neill et al., 1998; Kashkush et al., 2003; Kraitshtein et al., 2010) and several kinds of homoeologous interactions and exchanges (Gaeta et al., 2007; Kovarik et al., 2008a; Salmon et al., 2010; Szadkowski et al., 2010). Changes in duplicate gene expression are no less diverse, spanning the spectrum from expression conservation, relative to that of the diploid progenitors, to silencing of one homoeolog, to novel patterns of upand down-regulation (transgressive expression). Each of these transcriptomic responses varies in magnitude among allopolyploid species and individuals, among tissues and organ types within any one system, and with respect to the time since polyploid formation (Flagel et al., 2008; Flagel & Wendel, 2010). The phenotypic consequences of alterations in gene expression associated with hybridization and polyploidy are many and varied (Ni et al., 2009; Swanson-Wagner et al., 2009), underscoring the importance of understanding the expression level consequences of genomemerger and doubling. The advent and subsequent widespread utilization ofmicroarray and next-generation sequencing technologies has led to a rapid increase in explorations of gene expression in a variety of polyploid plants. These many efforts have generated a sufficient body of empirical data that generalizations are beginning to emerge concerning transcriptome changes in allopolyploids. For example, in every allopolyploid examined to date, some fraction of the duplicate gene pairs will be expressed unequally, and this suite of unequally expressed genes may itself favor one of the co-resident genomes, leading to a transcriptome that is unequally expressed with respect to the component genomes. While these generalizations are broadly applicable, much remains to be learned regarding the mechanistic underpinnings of duplicate gene expression change, the proximate and ultimate causes of inter-taxon and inter-organ variation in the response dynamics to polyploidy, and the functional, ecological, and evolutionary significance of duplicate gene expression modification. In addition to unequal expression of two homoeologs, other phenomena have been described which are even more poorly understood and for which fewer examples have yet been published. One of these is the concept of genome dominance (or genome expression dominance), which describes the expression condition in an allopolyploid where, for a given gene, the total expression of homoeologs is statistically the same as only one of the polyploid parents. This phenomenon was originally described for cotton allopolyploids by Rapp et al. (2009), confirmed and extended by Flagel & Wendel (2010), and subsequently described for both Spartina (Chelaifa et al., 2010) and Coffea (Bardil et al., 2011). This phenomenon is distinct from homoeolog expression bias (sometimes referred to as transcriptome dominance on a genome-wide basis), which describes the relative expression of homoeologs. Moreover, similar words are being used for rather different phenomena. Schnable et al. (2011), for example, invoked the term genomic dominance in maize, in a paper in which they demonstrated that the two subgenomes derived from the most recent polyploidy event in maize have experienced differential gene loss, with an accompanying gene expression bias favoring the more conserved subgenome (Schnable et al., 2011). By other accounts (Chen, 2007; Flagel & Wendel, 2010), this would be considered homoeolog expression bias (or transcriptome dominance) of ancient homoeologs. This inconsistency of conceptual application of the term genomic dominance also applies to the preferential expression of one subgenome of wheat (Akhunova et al., 2010), and to the patterns of biased expression in the fractionated subgenomes of paleohexaploid Brassica rapa (Cheng et al., 2012; Tang et al., 2012). This semantic and conceptual confusion appears to be gaining foothold in the literature; the phenomenon of preferential expression of one parental genome relative to the other in a polyploid species is termed genomic dominance in two recent reviews (Freeling et al., 2012; Schnable et al., 2012), citing both Schnable et al. (2011) and Flagel & Wendel (2010), and the term has also been applied to genomic modifications (Nicolas et al., 2012). Further complicating matters is the classical genetic concept associated with the term ‘dominance’, which conveys the relative expression hierarchy among a set of alleles. Against this backdrop of terminological and conceptual inconsistency, we thought it might be useful to briefly review the primary phenotypes of gene expression modification associated with allopolyploidy. Toward that end we describe and distinguish expression pattern changes observed in hybrid and polyploid species, and suggest a terminology (homoeolog expression bias and expression level dominance; Table 1; Fig. 1) that we hope will increase clarity of communication.

Polyploid formation pathways have an impact on genetic rearrangements in resynthesized Brassica napus

• Polyploids can be produced by the union of unreduced gametes or through somatic doubling of F(1) interspecific hybrids. The first route is suspected to produce allopolyploid species under natural conditions, whereas experimental data have only been thoroughly gathered for the latter. • We analyzed the meiotic behavior of an F(1) interspecific hybrid (by crossing Brassica oleracea and B.rapa, progenitors of B.napus) and the extent to which recombined homoeologous chromosomes were transmitted to its progeny. These results were then compared with results obtained for a plant generated by somatic doubling of this F₁ hybrid (CD.S₀) and an amphidiploid (UG.S₀) formed via a pathway involving unreduced gametes; we studied the impact of this method of polyploid formation on subsequent generations. • This study revealed that meiosis of the F₁ interspecific hybrid generated more gametes with recombined chromosomes than did meiosis of the plant produced by somatic doubling, although the size of these translocations was smaller. In the progeny of the UG.S₀ plant, there was an unexpected increase in the frequency at which the C1 chromosome was replaced by the A1 chromosome. • We conclude that polyploid formation pathways differ in their genetic outcome. Our study opens up perspectives for the understanding of polyploid origins.

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.

Molecular confirmation of species status for the allopolyploid cotton species, Gossypium ekmanianum Wittmack

Understanding the relationship between domesticated crop species and their wild relatives is paramount to germplasm maintenance and the utilization of wild relatives in breeding programs. Recently, Gossypium ekmanianum was resurrected as an independent species based on morphological analysis of specimens obtained from the Dominican Republic, where the original type specimen was collected. The molecular data presented here support the recognition of G. ekmanianum Wittmack as a distinct species that is phylogenetically close to G. hirsutum L. Analyses of chloroplast DNA data reveal species-specific, indel polymorphisms that unambiguously distinguish G. ekmanianum samples from other polyploid congeners. Furthermore, analysis of accessions that originated from the Dominican Republic demonstrate the cryptic inclusion of this sister taxon within the US National Plant Germplasm System, a germplasm collection maintained for diversity preservation and future breeding resources. The data presented here indicate that “wild” G. hirsutum accessions may include the closely related G. ekmanianum, and provide a method to easily distinguish the two.

Crossover rate between homologous chromosomes and interference are regulated by the addition of specific unpaired chromosomes in Brassica

Recombination is a major mechanism generating genetic diversity, but the control of the crossover rate remains a key question. In Brassica napus (AACC, 2n = 38), we can increase the homologous recombination between A genomes in AAC hybrids. Hypotheses for this effect include the number of C univalent chromosomes, the ratio between univalents and bivalents and, finally, which of the chromosomes are univalents. To test these hypotheses, we produced AA hybrids with zero, one, three, six or nine additional C chromosomes and four different hybrids carrying 2n = 32 and 2n = 35 chromosomes. The genetic map lengths for each hybrid were established to compare their recombination rates. The rates were 1.4 and 2.7 times higher in the hybrids having C6 or C9 alone than in the control (0C). This enhancement reached 3.1 and 4.1 times in hybrids carrying six and nine C chromosomes, and it was also higher for each pair of hybrids carrying 2n = 32 or 2n = 35 chromosomes, with a dependence on which chromosomes remained as univalents. We have shown, for the first time, that the presence of one chromosome, C9 , affects significantly the recombination rate and reduces crossover interference. This result will have fundamental implications on the regulation of crossover frequency.

Nucleotide diversity in the two co-resident genomes of allopolyploid cotton

Genetic diversity within and among populations lies at the heart of evolution. Unraveling the extent to which each intrinsic or extrinsic factor determines levels of diversity among genes, populations, and species is challenging, given the difficulty of isolating any single potentially important variable from all others. Allopolyploid species provide an opportunity to disentangle external and intrinsic factors, as the two (or more) homoeologous genomes co-occur in the same nucleus, often exhibiting high collinearity along homoeologous chromosomes. Here we evaluate the pace of molecular evolution and intraspecific, intragenomic diversity in two species of allopolyploid Gossypium, G. hirsutum and G. barbadense, using several hundred genes sequenced from multiple accessions of each species. Genic diversity in both species is low, having been influenced both by the polyploid bottleneck and a domestication bottleneck (for cultivated accessions), but with a directional bias in homoeolog diversity favoring the same genome in both allopolyploids. Total diversity is remarkably similar for the two homoeologous genomes overall, but the two copies of many gene pairs have accumulated statistically different diversity levels, and in a biased fashion with respect to genome. Domesticated accessions show reduced diversity in both genomes, as expected, but with a much greater reduction in one of the two homoeologous genomes. Furthermore, this biased reduction affects opposite homoeologous genomes in the two species. Interspecific introgression has played a role in shaping diversity within each species. Introgression was only detected for certain accessions, and only from G. barbadense into G. hirsutum in one of the two co-resident genomes.