Rapid, Repeated, and Clustered Loss of Duplicate Genes in Allopolyploid Plant Populations of Independent Origin

Systematic approach of polyploidy as an evolutionary genetic and genomic phenomenon in horticultural crops

Polyploidy is thought to be an evolutionary and systematic mechanism for gene flow and phenotypic advancement in flowering plants. It is a natural phenomenon that promotes diversity by creating new permutations enhancing the prime potentials as compared to progenitors. Two different pathways have been recognized in studying polyploidy in nature; mitotic or somatic chromosome doubling and cytogenetics variation. Secondly, the vital influence of being polyploid is its heritable property (unreduced reproductive cells) formed during first and second-division restitution (FDR & SDR). Different approaches either chemical (Colchicine, Oryzalin, Caffeine, Trifuralin, or phosphoric amides) or gaseous i.e. Nitrous oxide have been deliberated as strong polyploidy causing agents. A wide range of cytogenetic practices like chromosomes study, ploidy, genome analysis, and plant morphology and anatomy have been studied in different plant species. Flow cytometry for ploidy and chromosome analysis through fluorescence and genomic in situ hybridization (FISH & GISH) are the basic methods to evaluate heredity substances sampled from leaves and roots. Many horticultural crops have been developed successfully and released commercially for consumption. Moreover, some deep detailed studies are needed to check the strong relationship between unique morphological features and genetic makeup concerning genes and hormonal expression in a strong approach.

GhGASA14 regulates the flowering time of upland cotton in response to GA3

KEY MESSAGE

The silencing of GhGASA14 and the identification of superior allelic variation in its coding region indicate that GhGASA14 may positively regulate flowering and the response to GA3. Gibberellic acid-stimulated Arabidopsis (GASA), a member of the gibberellin-regulated short amino acid family, has been extensively investigated in several plant species and found to be critical for plant growth and development. However, research on this topic in cotton has been limited. In this study, we identified 38 GhGASAs that were dispersed across 18 chromosomes in upland cotton, and all of these genes had a GASA core domain. Transcriptome expression patterns and qRT-PCR results revealed that GhGASA9 and GhGASA14 exhibited upregulated expression not only in the floral organs but also in the leaves of early-maturing cultivars. The two genes were functionally characterized by virus-induced gene silencing (VIGS), and the budding and flowering times after silencing the target genes were later than those of the control (TRV:00). Compared with that in the water-treated group (MOCK), the flowering period of the different fruiting branches in the GA3-treated group was more concentrated. Interestingly, allelic variation was detected in the coding sequence of GhGASA14 between early-maturing and late-maturing accessions, and the frequency of this favorable allele was greater in high-latitude cotton cultivars than in low-latitude ones. Additionally, a significant linear relationship was observed between the expression level of GhGASA14 and flowering time among the 12 upland cotton accessions. Taken together, these results indicated that GhGASA14 may positively regulate flowering time and respond to GA3. These findings could lead to the use of valuable genetic resources for breeding early-maturing cotton cultivars in the future.

Transgenerational epigenetic inheritance during plant evolution and breeding

Plants can program and reprogram their genomes to create genetic variation and epigenetic modifications, leading to phenotypic plasticity. Although consequences of genetic changes are comprehensible, the basis for transgenerational inheritance of epigenetic variation is elusive. This review addresses contributions of external (environmental) and internal (genomic) factors to the establishment and maintenance of epigenetic memory during plant evolution, crop domestication, and modern breeding. Dynamic and pervasive changes in DNA methylation and chromatin modifications provide a diverse repertoire of epigenetic variation potentially for transgenerational inheritance. Elucidating and harnessing epigenetic inheritance will help us develop innovative breeding strategies and biotechnological tools to improve crop yield and resilience in the face of environmental challenges. Beyond plants, epigenetic principles are shared across sexually reproducing organisms including humans with relevance to medicine and public health.

Genome-wide DNA methylation dynamics following recent polyploidy in the allotetraploid Tragopogon miscellus (Asteraceae)

Polyploidy is an important evolutionary force, yet epigenetic mechanisms, such as DNA methylation, that regulate genome-wide expression of duplicated genes remain largely unknown. Here, we use Tragopogon (Asteraceae) as a model system to discover patterns and temporal dynamics of DNA methylation in recently formed polyploids. The naturally occurring allotetraploid Tragopogon miscellus formed in the last 95-100 yr from parental diploids Tragopogon dubius and T. pratensis. We profiled the DNA methylomes of these three species using whole-genome bisulfite sequencing. Genome-wide methylation levels in T. miscellus were intermediate between its diploid parents. However, nonadditive CG and CHG methylation occurred in transposable elements (TEs), with variation among TE types. Most differentially methylated regions (DMRs) showed parental legacy, but some novel DMRs were detected in the polyploid. Differentially methylated genes (DMGs) were also identified and characterized. This study provides the first assessment of both overall and locus-specific patterns of DNA methylation in a recent natural allopolyploid and shows that novel methylation variants can be generated rapidly after polyploid formation. Together, these results demonstrate that mechanisms to regulate duplicate gene expression may arise soon after allopolyploid formation and that these mechanisms vary among genes.

Whole-genome duplication in the Multicellularity Long Term Evolution Experiment

Whole-genome duplication (WGD) is widespread across eukaryotes and can promote adaptive evolution1-4. However, given the instability of newly-formed polyploid genomes5-7, understanding how WGDs arise in a population, persist, and underpin adaptations remains a challenge. Using our ongoing Multicellularity Long Term Evolution Experiment (MuLTEE)8, we show that diploid snowflake yeast (Saccharomyces cerevisiae) under selection for larger multicellular size rapidly undergo spontaneous WGD. From its origin within the first 50 days of the experiment, tetraploids persist for the next 950 days (nearly 5,000 generations, the current leading edge of our experiment) in ten replicate populations, despite being genomically unstable. Using synthetic reconstruction, biophysical modeling, and counter-selection experiments, we found that tetraploidy evolved because it confers immediate fitness benefits in this environment, by producing larger, longer cells that yield larger clusters. The same selective benefit also maintained tetraploidy over long evolutionary timescales, inhibiting the reversion to diploidy that is typically seen in laboratory evolution experiments. Once established, tetraploidy facilitated novel genetic routes for adaptation, playing a key role in the evolution of macroscopic multicellular size via the origin of evolutionarily conserved aneuploidy. These results provide unique empirical insights into the evolutionary dynamics and impacts of WGD, showing how it can initially arise due to its immediate adaptive benefits, be maintained by selection, and fuel long-term innovations by creating additional dimensions of heritable genetic variation.

Investigation of regulatory divergence between homoeologs in the recently formed allopolyploids, Tragopogon mirus and T. miscellus (Asteraceae)

Polyploidy is an important evolutionary process throughout eukaryotes, particularly in flowering plants. Duplicated gene pairs (homoeologs) in allopolyploids provide additional genetic resources for changes in molecular, biochemical, and physiological mechanisms that result in evolutionary novelty. Therefore, understanding how divergent genomes and their regulatory networks reconcile is vital for unraveling the role of polyploidy in plant evolution. Here, we compared the leaf transcriptomes of recently formed natural allotetraploids (Tragopogon mirus and T. miscellus) and their diploid parents (T. porrifolius X T. dubius and T. pratensis X T. dubius, respectively). Analysis of 35 400 expressed loci showed a significantly higher level of transcriptomic additivity compared to old polyploids; only 22% were non-additively expressed in the polyploids, with 5.9% exhibiting transgressive expression (lower or higher expression in the polyploids than in the diploid parents). Among approximately 7400 common orthologous regions (COREs), most loci in both allopolyploids exhibited expression patterns that were vertically inherited from their diploid parents. However, 18% and 20.3% of the loci showed novel expression bias patterns in T. mirus and T. miscellus, respectively. The expression changes of 1500 COREs were explained by cis-regulatory divergence (the condition in which the two parental subgenomes do not interact) between the diploid parents, whereas only about 423 and 461 of the gene expression changes represent trans-effects (the two parental subgenomes interact) in T. mirus and T. miscellus, respectively. The low degree of both non-additivity and trans-effects on gene expression may present the ongoing evolutionary processes of the newly formed Tragopogon polyploids (~80-90 years).

Transcription factor retention through multiple polyploidization steps in wheat

Whole-genome duplication is widespread in plant evolutionary history and is followed by nonrandom gene loss to return to a diploid state. Across multiple angiosperm species, the retained genes tend to be dosage-sensitive regulatory genes such as transcription factors, yet data for younger polyploid species is sparse. Here, we analyzed the retention, expression, and genetic variation in transcription factors in the recent allohexaploid bread wheat (Triticum aestivum L.). By comparing diploid, tetraploid, and hexaploid wheat, we found that, following each of two hybridization and whole-genome duplication events, the proportion of transcription factors in the genome increased. Transcription factors were preferentially retained over other genes as homoeologous groups in tetraploid and hexaploid wheat. Across cultivars, transcription factor homoeologs contained fewer deleterious missense mutations than nontranscription factors, suggesting that transcription factors are maintained as three functional homoeologs in hexaploid wheat populations. Transcription factor homoeologs were more strongly coexpressed than nontranscription factors, indicating conservation of function between homoeologs. We found that the B3, MADS-M-type, and NAC transcription factor families were less likely to have three homoeologs present than other families, which was associated with low expression levels and high levels of tandem duplication. Together, our results show that transcription factors are preferentially retained in polyploid wheat genomes although there is variation between families. Knocking out one transcription factor homoeolog to alter gene dosage, using TILLING or CRISPR, could generate new phenotypes for wheat breeding.

Global Patterns of Subgenome Evolution in Organelle-Targeted Genes of Six Allotetraploid Angiosperms

Whole-genome duplications (WGDs) are a prominent process of diversification in eukaryotes. The genetic and evolutionary forces that WGD imposes on cytoplasmic genomes are not well understood, despite the central role that cytonuclear interactions play in eukaryotic function and fitness. Cellular respiration and photosynthesis depend on successful interaction between the 3,000+ nuclear-encoded proteins destined for the mitochondria or plastids and the gene products of cytoplasmic genomes in multi-subunit complexes such as OXPHOS, organellar ribosomes, Photosystems I and II, and Rubisco. Allopolyploids are thus faced with the critical task of coordinating interactions between the nuclear and cytoplasmic genes that were inherited from different species. Because the cytoplasmic genomes share a more recent history of common descent with the maternal nuclear subgenome than the paternal subgenome, evolutionary "mismatches" between the paternal subgenome and the cytoplasmic genomes in allopolyploids might lead to the accelerated rates of evolution in the paternal homoeologs of allopolyploids, either through relaxed purifying selection or strong directional selection to rectify these mismatches. We report evidence from six independently formed allotetraploids that the subgenomes exhibit unequal rates of protein-sequence evolution, but we found no evidence that cytonuclear incompatibilities result in altered evolutionary trajectories of the paternal homoeologs of organelle-targeted genes. The analyses of gene content revealed mixed evidence for whether the organelle-targeted genes are lost more rapidly than the non-organelle-targeted genes. Together, these global analyses provide insights into the complex evolutionary dynamics of allopolyploids, showing that the allopolyploid subgenomes have separate evolutionary trajectories despite sharing the same nucleus, generation time, and ecological context.

Tracing the Evolution of the Angiosperm Genome from the Cytogenetic Point of View

Cytogenetics constitutes a branch of genetics that is focused on the cellular components, especially chromosomes, in relation to heredity and genome structure, function and evolution. The use of modern cytogenetic approaches and the latest microscopes with image acquisition and processing systems enables the simultaneous two- or three-dimensional, multicolour visualisation of both single-copy and highly-repetitive sequences in the plant genome. The data that is gathered using the cytogenetic methods in the phylogenetic background enable tracing the evolution of the plant genome that involve changes in: (i) genome sizes; (ii) chromosome numbers and morphology; (iii) the content of repetitive sequences and (iv) ploidy level. Modern cytogenetic approaches such as FISH using chromosome- and genome-specific probes have been widely used in studies of the evolution of diploids and the consequences of polyploidy. Nowadays, modern cytogenetics complements analyses in other fields of cell biology and constitutes the linkage between genetics, molecular biology and genomics.

GC content of plant genes is linked to past gene duplications

The frequency of G and C nucleotides in genomes varies from species to species, and sometimes even between different genes in the same genome. The monocot grasses have a bimodal distribution of genic GC content absent in dicots. We categorized plant genes from 5 dicots and 4 monocot grasses by synteny to related species and determined that syntenic genes have significantly higher GC content than non-syntenic genes at their 5`-end in the third position within codons for all 9 species. Lower GC content is correlated with gene duplication, as lack of synteny to distantly related genomes is associated with past interspersed gene duplications. Two mutation types can account for biased GC content, mutation of methylated C to T and gene conversion from A to G. Gene conversion involves non-reciprocal exchanges between homologous alleles and is not detectable when the alleles are identical or heterozygous for presence-absence variation, both likely situations for genes duplicated to new loci. Gene duplication can cause production of siRNA which can induce targeted methylation, elevating mC→T mutations. Recently duplicated plant genes are more frequently methylated and less likely to undergo gene conversion, each of these factors synergistically creating a mutational environment favoring AT nucleotides. The syntenic genes with high GC content in the grasses compose a subset that have undergone few duplications, or for which duplicate copies were purged by selection. We propose a “biased gene duplication / biased mutation” (BDBM) model that may explain the origin and trajectory of the observed link between duplication and genic GC bias. The BDBM model is supported by empirical data based on joint analyses of 9 angiosperm species with their genes categorized by duplication status, GC content, methylation levels and functional classes.

MicroRNAs as Indicators into the Causes and Consequences of Whole-Genome Duplication Events

Whole-genome duplications (WGDs) have long been considered the causal mechanism underlying dramatic increases to morphological complexity due to the neo-functionalization of paralogs generated during these events. Nonetheless, an alternative hypothesis suggests that behind the retention of most paralogs is not neo-functionalization, but instead the degree of the inter-connectivity of the intended gene product, as well as the mode of the WGD itself. Here, we explore both the causes and consequences of WGD by examining the distribution, expression, and molecular evolution of microRNAs (miRNAs) in both gnathostome vertebrates as well as chelicerate arthropods. We find that although the number of miRNA paralogs tracks the number of WGDs experienced within the lineage, few of these paralogs experienced changes to the seed sequence, and thus are functionally equivalent relative to their mRNA targets. Nonetheless, in gnathostomes, although the retention of paralogs following the 1R autotetraploidization event is similar across the two subgenomes, the paralogs generated by the gnathostome 2R allotetraploidization event are retained in higher numbers on one subgenome relative to the second, with the miRNAs found on the preferred subgenome showing both higher expression of mature miRNA transcripts and slower molecular evolution of the precursor miRNA sequences. Importantly, WGDs do not result in the creation of miRNA novelty, nor do WGDs correlate to increases in complexity. Instead, it is the number of miRNA seed sequences in the genome itself that not only better correlate to instances in complexification, but also mechanistically explain why complexity increases when new miRNA families are established.

Impact of polyploidy on plant tolerance to abiotic and biotic stresses

Polyploidy, defined as the coexistence of three or more complete sets of chromosomes in an organism's cells, is considered as a pivotal moving force in the evolutionary history of vascular plants and has played a major role in the domestication of several crops. In the last decades, improved cultivars of economically important species have been developed artificially by inducing autopolyploidy with chemical agents. Studies on diverse species have shown that the anatomical and physiological changes generated by either natural or artificial polyploidization can increase tolerance to abiotic and biotic stresses as well as disease resistance, which may positively impact on plant growth and net production. The aim of this work is to review the current literature regarding the link between plant ploidy level and tolerance to abiotic and biotic stressors, with an emphasis on the physiological and molecular mechanisms responsible for these effects, as well as their impact on the growth and development of both natural and artificially generated polyploids, during exposure to adverse environmental conditions. We focused on the analysis of those types of stressors in which more progress has been made in the knowledge of the putative morpho-physiological and/or molecular mechanisms involved, revealing both the factors in common, as well as those that need to be addressed in future research.

Repeatome-Based Phylogenetics in Pelargonium Section Ciconium (Sweet) Harvey

The repetitive part of the genome (the repeatome) contains a wealth of often overlooked information that can be used to resolve phylogenetic relationships and test evolutionary hypotheses for clades of related plant species such as Pelargonium. We have generated genome skimming data for 18 accessions of Pelargonium section Ciconium and one outgroup. We analyzed repeat abundancy and repeat similarity in order to construct repeat profiles and then used these for phylogenetic analyses. We found that phylogenetic trees based on read similarity were largely congruent with previous work based on morphological and chloroplast sequence data. For example, results agreed in identifying a “Core Ciconium” group which evolved after the split with P. elongatum. We found that this group was characterized by a unique set of repeats, which confirmed currently accepted phylogenetic hypotheses. We also found four species groups within P. sect. Ciconium that reinforce previous plastome-based reconstructions. A second repeat expansion was identified in a subclade which contained species that are considered to have dispersed from Southern Africa into Eastern Africa and the Arabian Peninsula. We speculate that the Core Ciconium repeat set correlates with a possible WGD event leading to this branch.

Meiosis in Polyploids and Implications for Genetic Mapping: A Review

Plant cytogenetic studies have provided essential knowledge on chromosome behavior during meiosis, contributing to our understanding of this complex process. In this review, we describe in detail the meiotic process in auto- and allopolyploids from the onset of prophase I through pairing, recombination, and bivalent formation, highlighting recent findings on the genetic control and mode of action of specific proteins that lead to diploid-like meiosis behavior in polyploid species. During the meiosis of newly formed polyploids, related chromosomes (homologous in autopolyploids; homologous and homoeologous in allopolyploids) can combine in complex structures called multivalents. These structures occur when multiple chromosomes simultaneously pair, synapse, and recombine. We discuss the effectiveness of crossover frequency in preventing multivalent formation and favoring regular meiosis. Homoeologous recombination in particular can generate new gene (locus) combinations and phenotypes, but it may destabilize the karyotype and lead to aberrant meiotic behavior, reducing fertility. In crop species, understanding the factors that control pairing and recombination has the potential to provide plant breeders with resources to make fuller use of available chromosome variations in number and structure. We focused on wheat and oilseed rape, since there is an abundance of elucidating studies on this subject, including the molecular characterization of the Ph1 (wheat) and PrBn (oilseed rape) loci, which are known to play a crucial role in regulating meiosis. Finally, we exploited the consequences of chromosome pairing and recombination for genetic map construction in polyploids, highlighting two case studies of complex genomes: (i) modern sugarcane, which has a man-made genome harboring two subgenomes with some recombinant chromosomes; and (ii) hexaploid sweet potato, a naturally occurring polyploid. The recent inclusion of allelic dosage information has improved linkage estimation in polyploids, allowing multilocus genetic maps to be constructed.

The role of genome duplication in big sagebrush growth and fecundity

PREMISE

Adaptive traits can be dramatically altered by genome duplication. The study of interactions among traits, ploidy, and the environment are necessary to develop an understanding of how polyploidy affects niche differentiation and to develop restoration strategies for resilient native ecosystems.

METHODS

Growth and fecundity were measured in common gardens for 39 populations of big sagebrush (Artemisia tridentata) containing two subspecies and two ploidy levels. General linear mixed-effect models assessed how much of the trait variation could be attributed to genetics (i.e., ploidy and climatic adaptation), environment, and gene-environment interactions.

RESULTS

Growth and fecundity variation were explained well by the mixed models (80% and 91%, respectively). Much of the trait variation was attributed to environment, and 15% of variation in growth and 34% of variation in seed yield were attributed to genetics. Genetic trait variation was mostly attributable to ploidy, with much higher growth and seed production in diploids, even in a warm-dry environment typically dominated by tetraploids. Population-level genetic variation was also evident and was related to the climate of each population's origin.

CONCLUSIONS

Ploidy is a strong predictor growth and seed yield, regardless of common-garden environment. The superior growth and fecundity of diploids across environments raises the question as to how tetraploids can be more prevalent than diploids, especially in warm-dry environments. Two hypotheses that may explain the abundance of tetraploids on the landscape include selection for drought resistance at the seedling stage, and greater competitive ability in water uptake in the upper soil horizon.

Characterization and expression analysis of wall-associated kinase (WAK) and WAK-like family in cotton

The wall-associated kinases (WAKs) and WAK-like kinases (WAKLs) form a group of receptor-like kinases (RLKs) with extracellular domains tightly linked to the cell wall. The WAKs/WAKLs have been known to be involved in plant growth, development, and stress responses. However, the functions of WAKs/WAKLs are less well known in cotton. In this study, 58, 66, and 99 WAK/WAKL genes were identified in Gossypium arboreum, G. raimondii, and G. hirsutum, respectively. Phylogenetic analysis showed they were classified into five groups, with two groups specific to cotton. Collinearity analysis revealed that segmental and tandem duplications resulted in expansion of the WAK/WAKL gene family in cotton. Moreover, the Ka/Ks ratios indicated this family was exposed to purifying selection pressure during evolution. The structures of the GhWAK/WAKL genes and encoded proteins suggested the functions of WAKs/WAKLs in cotton were conserved. Transient expression of four WAK/WAKL-GFP fusion constructs in Arabidopsis protoplasts indicated that they were localized on the plasma membrane. The cis-elements in the GhWAK/WAKL promoters were responsive to multiple phytohormones and abiotic stresses. Expression profiling showed that GhWAK/WAKL genes were induced by various abiotic stresses. This study provides insights into the evolution of WAK/WAKL genes and presents fundamental information for further analysis in cotton.

Trajectories of Homoeolog-Specific Expression in Allotetraploid Tragopogon castellanus Populations of Independent Origins

Polyploidization can have a significant ecological and evolutionary impact by providing substantially more genetic material that may result in novel phenotypes upon which selection may act. While the effects of polyploidization are broadly reviewed across the plant tree of life, the reproducibility of these effects within naturally occurring, independently formed polyploids is poorly characterized. The flowering plant genus Tragopogon (Asteraceae) offers a rare glimpse into the intricacies of repeated allopolyploid formation with both nascent (< 90 years old) and more ancient (mesopolyploids) formations. Neo- and mesopolyploids in Tragopogon have formed repeatedly and have extant diploid progenitors that facilitate the comparison of genome evolution after polyploidization across a broad span of evolutionary time. Here, we examine four independently formed lineages of the mesopolyploid Tragopogon castellanus for homoeolog expression changes and fractionation after polyploidization. We show that expression changes are remarkably similar among these independently formed polyploid populations with large convergence among expressed loci, moderate convergence among loci lost, and stochastic silencing. We further compare and contrast these results for T. castellanus with two nascent Tragopogon allopolyploids. While homoeolog expression bias was balanced in both nascent polyploids and T. castellanus, the degree of additive expression was significantly different, with the mesopolyploid populations demonstrating more non-additive expression. We suggest that gene dosage and expression noise minimization may play a prominent role in regulating gene expression patterns immediately after allopolyploidization as well as deeper into time, and these patterns are conserved across independent polyploid lineages.

Patterns and Processes of Diploidization in Land Plants

Most land plants are now known to be ancient polyploids that have rediploidized. Diploidization involves many changes in genome organization that ultimately restore bivalent chromosome pairing and disomic inheritance, and resolve dosage and other issues caused by genome duplication. In this review, we discuss the nature of polyploidy and its impact on chromosome pairing behavior. We also provide an overview of two major and largely independent processes of diploidization: cytological diploidization and genic diploidization/fractionation. Finally, we compare variation in gene fractionation across land plants and highlight the differences in diploidization between plants and animals. Altogether, we demonstrate recent advancements in our understanding of variation in the patterns and processes of diploidization in land plants and provide a road map for future research to unlock the mysteries of diploidization and eukaryotic genome evolution.

Whole-Genome Sequence of Synthesized Allopolyploids in Cucumis Reveals Insights into the Genome Evolution of Allopolyploidization

The importance of allopolyploidy in plant evolution has been widely recognized. The genetic changes triggered by allopolyploidy, however, are not yet fully understood due to inconsistent phenomena reported across diverse species. The construction of synthetic polyploids offers a controlled approach to systematically reveal genomic changes that occur during the process of polyploidy. This study reports the first fully sequenced synthetic allopolyploid constructed from a cross between Cucumis sativus and C. hystrix, with high-quality assembly. The two subgenomes are confidently partitioned and the C. sativus-originated subgenome predominates over the C. hystrix-originated subgenome, retaining more sequences and showing higher homeologous gene expression. Most of the genomic changes emerge immediately after interspecific hybridization. Analysis of a series of genome sequences from several generations (S0, S4-S13) of C. ×hytivus confirms that genomic changes occurred in the very first generations, subsequently slowing down as the process of diploidization is initiated. The duplicated genome of the allopolyploid with double genes from both parents broadens the genetic base of C. ×hytivus, resulting in enhanced phenotypic plasticity. This study provides novel insights into plant polyploid genome evolution and demonstrates a promising strategy for the development of a wide array of novel plant species and varieties through artificial polyploidization.

Resolving phylogeny and polyploid parentage using genus-wide genome-wide sequence data from birch trees

Numerous plant genera have a history including frequent hybridisation and polyploidisation (allopolyploidisation), which means that their phylogeny is a network of reticulate evolution that cannot be accurately depicted as a bifurcating tree with a single tip per species. The genus Betula, which contains many ecologically important tree species, is a case in point. We generated genome-wide sequence reads for 27 diploid and 36 polyploid Betula species or subspecies using restriction site associated DNA (RAD) sequences. These reads were assembled into contigs with a mean length of 675 bp. We reconstructed the evolutionary relationships among diploid Betula species using both supermatrix (concatenation) and species tree methods. We identified the closest diploid relatives of the polyploids according to the relative rates at which reads from polyploids mapped to contigs from different diploid species within a concatenated reference sequence. By mapping reads from allopolyploids to their different putative diploid relatives we assembled contigs from the putative sub-genomes of allopolyploid taxa. We used these to build new phylogenies that included allopolyploid sub-genomes as separate tips. This approach yielded a highly evidenced phylogenetic hypothesis for the genus Betula, including the complex reticulate origins of the majority of its polyploid taxa. Our phylogeny divides the genus into two well supported clades, which, interestingly, differ in their seed-wing morphology. We therefore propose to split Betula into two subgenera.

Complete chloroplast genome sequencing and comparative analysis reveals changes to the chloroplast genome after allopolyploidization in Cucumis

Allopolyploids undergo "genomic shock" leading to significant genetic and epigenetic modifications. Previous studies have mainly focused on nuclear changes, while little is known about the inheritance and changes of organelle genome in allopolyploidization. The synthetic allotetraploid Cucumis ×hytivus, which is generated via hybridization between C. hystrix and C. sativus, is a useful model system for studying cytonuclear variation. Here, we report the chloroplast genome of allotetraploid C. ×hytivus and its diploid parents via sequencing and comparative analysis. The size of the obtained chloroplast genomes ranged from 154 673 to 155 760 bp, while their gene contents, gene orders, and GC contents were similar to each other. Comparative genome analysis supports chloroplast maternal inheritance. However, we identified 51 indels and 292 SNP genetic variants in the chloroplast genome of the allopolyploid C. ×hytivus relative to its female parent C. hystrix. Nine intergenic regions with rich variation were identified through comparative analysis of the chloroplast genomes within the subgenus Cucumis. The phylogenetic network based on the chloroplast genome sequences clarified the evolution and taxonomic position of the synthetic allotetraploid C. ×hytivus. The results of this study provide us with an insight into the changes of organelle genome after allopolyploidization, and a new understanding of the cytonuclear evolution.

Genetic and Epigenetic Changes Are Rapid Responses of the Genome to the Newly Synthesized Autotetraploid Carassius auratus

Whole genome duplication events have occurred frequently during the course of vertebrate evolution. To better understand the influence of polyploidization on the fish genome, we herein used the autotetraploid Carassius auratus (4n = 200, RRRR) (4nRR) resulting from the whole genome duplication of Carassius auratus (2n = 100, RR) (RCC) to explore the genomic and epigenetic alterations after polyploidization. We subsequently performed analyses of full-length transcriptome dataset, amplified fragment length polymorphism (AFLP) and methylation sensitive amplification polymorphism (MSAP) on 4nRR and RCC. By matching the results of 4nRR and RCC isoforms with reference genome in full-length transcriptome dataset, 649 and 1,971 novel genes were found in the RCC and 4nRR full-length geneset, respectively. Compared to Carassius auratus and Megalobrama amblycephala, 4nRR presented 3,661 unexpressed genes and 2,743 expressed genes. Furthermore, GO enrichment analysis of expressed genes in 4nRR revealed that they were enriched in meiosis I, whereas KEGG enrichment analysis displayed that they were mainly enriched in proteasome. Using AFLP analysis, we noted that 32.61% of RCC fragments had disappeared, while 32.79% of new bands were uncovered in 4nRR. Concerning DNA methylation, 4nRR exhibited a lower level of global DNA methylation than RCC. Additionally, 60.31% of methylation patterns in 4nRR were altered compared to RCC. These observations indicated that transcriptome alterations, genomic changes and regulation of DNA methylation levels and patterns had occurred in the newly established autotetraploid genomes, suggesting that genetic and epigenetic alterations were influenced by autotetraploidization. In summary, this study provides valuable novel insights into vertebrate genome evolution and generates relevant information for fish breeding.

Transcriptome Dynamics of the Inflorescence in Reciprocally Formed Allopolyploid Tragopogon miscellus (Asteraceae)

Polyploidy is an important evolutionary mechanism and is prevalent among land plants. Most polyploid species examined have multiple origins, which provide genetic diversity and may enhance the success of polyploids. In some polyploids, recurrent origins can result from reciprocal crosses between the same diploid progenitors. Although great progress has been made in understanding the genetic consequences of polyploidy, the genetic implications of reciprocal polyploidization remain poorly understood, especially in natural polyploids. Tragopogon (Asteraceae) has become an evolutionary model system for studies of recent and recurrent polyploidy. Allotetraploid T. miscellus has formed reciprocally in nature with resultant distinctive floral and inflorescence morphologies (i.e., short- vs. long-liguled forms). In this study, we performed comparative inflorescence transcriptome analyses of reciprocally formed T. miscellus and its diploid parents, T. dubius and T. pratensis. In both forms of T. miscellus, homeolog expression of ∼70% of the loci showed vertical transmission of the parental expression patterns (i.e., parental legacy), and ∼20% of the loci showed biased homeolog expression, which was unbalanced toward T. pratensis. However, 17.9% of orthologous pairs showed different homeolog expression patterns between the two forms of T. miscellus. No clear effect of cytonuclear interaction on biased expression of the maternal homeolog was found. In terms of the total expression level of the homeologs studied, 22.6% and 16.2% of the loci displayed non-additive expression in short- and long-liguled T. miscellus, respectively. Unbalanced expression level dominance toward T. pratensis was observed in both forms of T. miscellus. Significantly, genes annotated as being involved in pectin catabolic processes were highly expressed in long-liguled T. miscellus relative to the short-liguled form, and the majority of these differentially expressed genes were transgressively down-regulated in short-liguled T. miscellus. Given the known role of these genes in cell expansion, they may play a role in the differing floral and inflorescence morphologies of the two forms. In summary, the overall inflorescence transcriptome profiles are highly similar between reciprocal origins of T. miscellus. However, the dynamic homeolog-specific expression and non-additive expression patterns observed in T. miscellus emphasize the importance of reciprocal origins in promoting the genetic diversity of polyploids.

Recurrent genome duplication events likely contributed to both the ancient and recent rise of ferns

Ferns, the second largest group of vascular plants, originated ~400 million years ago (Mya). They became dominant in the ancient Earth landscape before the angiosperms and are still important in current ecosystems. Many ferns have exceptionally high chromosome numbers, possibly resulting from whole-genome duplications (WGDs). However, WGDs have not been investigated molecularly across fern diversity. Here we detected and dated fern WGDs using a phylogenomic approach and by calculating synonymous substitution rates (Ks). We also investigated a possible correlation between proposed WGDs and shifts in species diversification rates. We identified 19 WGDs: three ancient events along the fern phylogenetic backbone that are shared by 66%-97% of extant ferns, with additional lineage-specific WGDs for eight orders, providing strong evidence for recurring genome duplications across fern evolutionary history. We also observed similar Ks peak values for more than half of these WGDs, with multiple WGDs occurring close to the Cretaceous (~145-66 Mya). Despite the repeated WGD events, the biodiversity of ferns declined during the Cretaceous, implying that other factors probably contributed to the floristic turnover from ferns to angiosperms. This study provides molecular evidence for recurring WGDs in ferns and offers important clues to the genomic evolutionary history of ferns.

Molecular Evolution and Expansion of the KUP Family in the Allopolyploid Cotton Species Gossypium hirsutum and Gossypium barbadense

The comprehensive analysis of gene family evolution will elucidate the origin and evolution of gene families. The K+ uptake (KUP) gene family plays important roles in K+ uptake and transport, plant growth and development, and abiotic stress responses. However, the current understanding of the KUP family in cotton is limited. In this study, 51 and 53 KUPs were identified in Gossypium barbadense and Gossypium hirsutum, respectively. These KUPs were divided into five KUP subfamilies, with subfamily 2 containing three groups. Different subfamilies had different member numbers, conserved motifs, gene structures, regulatory elements, and gene expansion and loss rates. A paleohexaploidization event caused the expansion of GhKUP and GbKUP in cotton, and duplication events in G. hirsutum and G. barbadense have happened in a common ancestor of Gossypium. Meanwhile, the KUP members of the two allopolyploid subgenomes of G. hirsutum and G. barbadense exhibited unequal gene proportions, gene structural diversity, uneven chromosomal distributions, asymmetric expansion rates, and biased gene loss rates. In addition, the KUP families of G. hirsutum and G. barbadense displayed evolutionary conservation and divergence. Taken together, these results illustrated the molecular evolution and expansion of the KUP family in allopolyploid cotton species.

Genome-wide analysis and characterization of F-box gene family in Gossypium hirsutum L

Background

F-box proteins are substrate-recognition components of the Skp1-Rbx1-Cul1-F-box protein (SCF) ubiquitin ligases. By selectively targeting the key regulatory proteins or enzymes for ubiquitination and 26S proteasome mediated degradation, F-box proteins play diverse roles in plant growth/development and in the responses of plants to both environmental and endogenous signals. Studies of F-box proteins from the model plant Arabidopsis and from many additional plant species have demonstrated that they belong to a super gene family, and function across almost all aspects of the plant life cycle. However, systematic exploration of F-box family genes in the important fiber crop cotton (Gossypium hirsutum) has not been previously performed. The genome-wide analysis of the cotton F-box gene family is now possible thanks to the completion of several cotton genome sequencing projects.

Results

In current study, we first conducted a genome-wide investigation of cotton F-box family genes by reference to the published F-box protein sequences from other plant species. 592 F-box protein encoding genes were identified in the Gossypium hirsutume acc.TM-1 genome and, subsequently, we were able to present their gene structures, chromosomal locations, syntenic relationships with their parent species. In addition, duplication modes analysis showed that cotton F-box genes were distributed to 26 chromosomes, with the maximum number of genes being detected on chromosome 5. Although the WGD (whole-genome duplication) mode seems play a dominant role during cotton F-box gene expansion process, other duplication modes including TD (tandem duplication), PD (proximal duplication), and TRD (transposed duplication) also contribute significantly to the evolutionary expansion of cotton F-box genes. Collectively, these bioinformatic analysis suggest possible evolutionary forces underlying F-box gene diversification. Additionally, we also conducted analyses of gene ontology, and expression profiles in silico, allowing identification of F-box gene members potentially involved in hormone signal transduction.

Conclusion

The results of this study provide first insights into the Gossypium hirsutum F-box gene family, which lays the foundation for future studies of functionality, particularly those involving F-box protein family members that play a role in hormone signal transduction.

Nuclear-Cytoplasmic Coevolution Analysis of RuBisCO in Synthesized Cucumis Allopolyploid

Allopolyploids are often faced with the challenge of maintaining well-coordination between nuclear and cytoplasmic genes inherited from different species. The synthetic allotetraploid Cucumis × hytivus is a useful model to explore cytonuclear coevolution. In this study, the sequences and expression of cytonuclear enzyme complex RuBisCO as well as its content and activity in C. × hytivus were compared to its parents to explore plastid–nuclear coevolution. The plastome-coded rbcL gene sequence was confirmed to be stable maternal inheritance, and parental copy of nuclear rbcS genes were both preserved in C. × hytivus. Thus, the maternal plastid may interact with the biparentally inherited rbcS alleles. The expression of the rbcS gene of C-homoeologs (paternal) was significantly higher than that of H-homoeologs (maternal) in C. × hytivus (HHCC). Protein interaction prediction analysis showed that the rbcL protein has stronger binding affinity to the paternal copy of rbcS protein than that of maternal copy in C. × hytivus, which might explain the transcriptional bias of the rbcS homoeologs. Moreover, both the activity and content of RuBisCO in C. × hytivus showed mid-parent heterosis. In summary, our results indicate a paternal transcriptional bias of the rbcS genes in C. × hytivus, and we found new nuclear–cytoplasmic combination may be one of the reasons for allopolyploids heterosis.

Alignment of Rutaceae Genomes Reveals Lower Genome Fractionation Level Than Eudicot Genomes Affected by Extra Polyploidization

Owing to their nutritional and commercial values, the genomes of several citrus plants have been sequenced, and the genome of one close relative in the Rutaceae family, atalantia (Atalantia buxifolia), has also been sequenced. Here, we show a family-level comparative analysis of Rutaceae genomes. By using grape as the outgroup and checking cross-genome gene collinearity, we systematically performed a hierarchical and event-related alignment of Rutaceae genomes, and produced a gene list defining homologous regions based on ancestral polyploidization or speciation. We characterized genome fractionation resulting from gene loss or relocation, and found that erosion of gene collinearity could largely be described by a geometric distribution. Moreover, we found that well-assembled Rutaceae genomes retained significantly more genes (65-82%) than other eudicots affected by recursive polyploidization. Additionally, we showed divergent evolutionary rates among Rutaceae plants, with sweet orange evolving faster than others, and by performing evolutionary rate correction, re-dated major evolutionary events during their evolution. We deduced that the divergence between the Rutaceae family and grape occurred about 81.15-91.74 million years ago (mya), while the split between citrus and atalantia plants occurred <10 mya. In addition, we showed that polyploidization led to a copy number expansion of key gene families contributing to the biosynthesis of vitamin C. Overall, the present effort provides an important comparative genomics resource and lays a foundation to understand the evolution and functional innovation of Rutaceae genomes.

Transcriptome Profile Analysis on Ovarian Tissues of Autotetraploid Fish and Diploid Red Crucian Carp

Polyploidization can significantly alter the size of animal gametes. Autotetraploid fish (RRRR, 4nRR = 200) (4nRR) possessing four sets of chromosomes were derived from whole-genome duplication in red crucian carp (RR, 2n = 100) (RCC). The diploid eggs of the 4nRR fish were significantly larger than the eggs of RCC. To explore the differences between the ovaries of these two ploidies of fishes at the molecular level, we compared the ovary transcriptome profiles of 4nRR fish and RCC and identified differentially expressed genes (DEGs). A total of 19,015 unigenes were differentially expressed between 4nRR fish and RCC, including 12,591 upregulated and 6,424 downregulated unigenes in 4nRR fish. Functional analyses revealed that eight genes (CDKL1, AHCY, ARHGEF3, TGFβ, WNT11, CYP27A, GDF7, and CKB) were involved in the regulation of cell proliferation, cell division, gene transcription, ovary development and energy metabolism, suggesting that these eight genes were related to egg size in 4nRR fish and RCC. We validated the expression levels of these eight DEGs in 4nRR fish and RCC using quantitative PCR. The study results provided insights into the regulatory mechanisms underlying the differences in crucian carp egg sizes.

Parental legacy, demography, and admixture influenced the evolution of the two subgenomes of the tetraploid Capsella bursa-pastoris (Brassicaceae)

Allopolyploidy is generally perceived as a major source of evolutionary novelties and as an instantaneous way to create isolation barriers. However, we do not have a clear understanding of how two subgenomes evolve and interact once they have fused in an allopolyploid species nor how isolated they are from their relatives. Here, we address these questions by analyzing genomic and transcriptomic data of allotetraploid Capsella bursa-pastoris in three differentiated populations, Asia, Europe, and the Middle East. We phased the two subgenomes, one descended from the outcrossing and highly diverse Capsella grandiflora (CbpCg) and the other one from the selfing and genetically depauperate Capsella orientalis (CbpCo). For each subgenome, we assessed its relationship with the diploid relatives, temporal changes of effective population size (Ne), signatures of positive and negative selection, and gene expression patterns. In all three regions, Ne of the two subgenomes decreased gradually over time and the CbpCo subgenome accumulated more deleterious changes than CbpCg. There were signs of widespread admixture between C. bursa-pastoris and its diploid relatives. The two subgenomes were impacted differentially depending on geographic region suggesting either strong interploidy gene flow or multiple origins of C. bursa-pastoris. Selective sweeps were more common on the CbpCg subgenome in Europe and the Middle East, and on the CbpCo subgenome in Asia. In contrast, differences in expression were limited with the CbpCg subgenome slightly more expressed than CbpCo in Europe and the Middle-East. In summary, after more than 100,000 generations of co-existence, the two subgenomes of C. bursa-pastoris still retained a strong signature of parental legacy but their evolutionary trajectory strongly varied across geographic regions.

Divergent Expression Patterns and Function of Two cxcr4 Paralogs in Hermaphroditic Epinephelus coioides

Chemokine receptor Cxcr4 evolved two paralogs in the teleost lineage. However, cxcr4a and cxcr4b have been characterized only in a few species. In this study, we identified two cxcr4 paralogs from the orange-spotted grouper, Epinephelus coioides. The phylogenetic relationship and gene structure and synteny suggest that the duplicated cxcr4a/b should result from the teleost-specific genome duplication (Ts3R). The teleost cxcr4 gene clusters in two paralogous chromosomes exhibit a complementary gene loss/retention pattern. Ec_cxcr4a and Ec_cxcr4b show differential and biased expression patterns in grouper adult tissue, gonads, and embryos at different stages. During embryogenesis, Ec_cxcr4a/b are abundantly transcribed from the neurula stage and mainly expressed in the neural plate and sensory organs, indicating their roles in neurogenesis. Ec_Cxcr4a and Ec_Cxcr4b possess different chemotactic migratory abilities from the human SDF-1α, Ec_Cxcl12a, and Ec_Cxcl12b. Moreover, we uncovered the N-terminus and TM5 domain as the key elements for specific ligand⁻receptor recognition of Ec_Cxcr4a-Ec_Cxcl12b and Ec_Cxcr4b-Ec_Cxcl12a. Based on the biased and divergent expression patterns of Eccxcr4a/b, and specific ligand⁻receptor recognition of Ec_Cxcl12a/b⁻Ec_Cxcr4b/a, the current study provides a paradigm of sub-functionalization of two teleost paralogs after Ts3R.

A Robust Methodology for Assessing Differential Homeolog Contributions to the Transcriptomes of Allopolyploids

Polyploidy has played a pivotal and recurring role in angiosperm evolution. Allotetraploids arise from hybridization between species and possess duplicated gene copies (homeologs) that serve redundant roles immediately after polyploidization. Although polyploidization is a major contributor to plant evolution, it remains poorly understood. We describe an analytical approach for assessing homeolog-specific expression that begins with de novo assembly of parental transcriptomes and effectively (i) reduces redundancy in de novo assemblies, (ii) identifies putative orthologs, (iii) isolates common regions between orthologs, and (iv) assesses homeolog-specific expression using a robust Bayesian Poisson-Gamma model to account for sequence bias when mapping polyploid reads back to parental references. Using this novel methodology, we examine differential homeolog contributions to the transcriptome in the recently formed allopolyploids Tragopogon mirus and T. miscellus (Compositae). Notably, we assess a larger Tragopogon gene set than previous studies of this system. Using carefully identified orthologous regions and filtering biased orthologs, we find in both allopolyploids largely balanced expression with no strong parental bias. These new methods can be used to examine homeolog expression in any tetrapolyploid system without requiring a reference genome.

Application of CRISPR/Cas9 to Tragopogon (Asteraceae), an evolutionary model for the study of polyploidy

Tragopogon (Asteraceae) is an excellent natural system for studies of recent polyploidy. Development of an efficient CRISPR/Cas9-based genome editing platform in Tragopogon will facilitate novel studies of the genetic consequences of polyploidy. Here, we report our initial results of developing CRISPR/Cas9 in Tragopogon. We have established a feasible tissue culture and transformation protocol for Tragopogon. Through protoplast transient assays, use of the TragCRISPR system (i.e. the CRISPR/Cas9 system adapted for Tragopogon) was capable of introducing site-specific mutations in Tragopogon protoplasts. Agrobacterium-mediated transformation with Cas9-sgRNA constructs targeting the phytoene desaturase gene (TraPDS) was implemented in this model polyploid system. Sequencing of PCR amplicons from the target regions indicated simultaneous mutations of two alleles and four alleles of TraPDS in albino shoots from Tragopogon porrifolius (2x) and Tragopogon mirus (4x), respectively. The average proportions of successfully transformed calli with the albino phenotype were 87% and 78% in the diploid and polyploid, respectively. This appears to be the first demonstration of CRISPR/Cas9-based genome editing in any naturally formed neopolyploid system. Although a more efficient tissue culture system should be developed in Tragopogon, application of a robust CRISPR/Cas9 system will permit unique studies of biased fractionation, the gene-balance hypothesis and cytonuclear interactions in polyploids. In addition, the CRISPR/Cas9 platform enables investigations of those genes involved in phenotypic changes in polyploids and will also facilitate novel functional biology studies in Asteraceae. Our workflow provides a guide for applying CRISPR/Cas9 to other nongenetic model plant systems.

Whole Genome Incorporation and Epigenetic Stability in a Newly Synthetic Allopolyploid of Gynogenetic Gibel Carp

Allopolyploidization plays an important role in speciation, and some natural or synthetic allopolyploid fishes have been extensively applied to aquaculture. Although genetic and epigenetic inheritance and variation associated with plant allopolyploids have been well documented, the relative research in allopolyploid animals is scarce. In this study, the genome constitution and DNA methylation inheritance in a newly synthetic allopolyploid of gynogenetic gibel carp were analyzed. The incorporation of a whole genome of paternal common carp sperm in the allopolyploid was confirmed by genomic in situ hybridization, chromosome localization of 45S rDNAs, and sequence comparison. Pooled sample-based methylation sensitive amplified polymorphism (MSAP) revealed that an overwhelming majority (98.82%) of cytosine methylation patterns in the allopolyploid were inherited from its parents of hexaploid gibel carp clone D and common carp. Compared to its parents, 11 DNA fragments in the allopolyploid were proved to be caused by interindividual variation, recombination, deletion, and mutation through individual sample-based MSAP and sequencing. Contrast to the rapid and remarkable epigenetic changes in most of analyzed neopolyploids, no cytosine methylation variation was detected in the gynogenetic allopolyploid. Therefore, the newly synthetic allopolyploid of gynogenetic gibel carp combined genomes from its parents and maintained genetic and epigenetic stability after its formation and subsequently seven successive gynogenetic generations. Our current results provide a paradigm for recurrent polyploidy consequences in the gynogenetic allopolyploid animals.

Evolutionary Model for the Statistical Divergence of Paralogous and Orthologous Gene Pairs Generated by Whole Genome Duplication and Speciation

We outline a principled approach to the analysis of duplicate gene similarity distributions, based on a model integrating sequence divergence and the process of fractionation of duplicate genes resulting from whole genome duplication (WGD). This model allows us to predict duplicate gene similarity distributions for a series of two or three WGD, for whole genome triplication followed by a WGD, and for triplication, followed by speciation, followed by WGD. We calculate the probabilities of all possible fates of a gene pair as its two members proliferate or are lost, predicting the number of surviving pairs from each event. We discuss how to calculate maximum likelihood estimators for the parameters of these models, illustrating with an analysis of the distribution of paralog similarities in the poplar genome.

Epigenetic perspectives on the evolution and domestication of polyploid plant and crops

Polyploidy or whole genome duplication (WGD) is a prominent feature for genome evolution of some animals and all flowering plants, including many important crops such as wheat, cotton, and canola. In autopolyploids, genome duplication often perturbs dosage regulation on biological networks. In allopolyploids, interspecific hybridization could induce genetic and epigenetic changes, the effects of which could be amplified by genome doubling (ploidy changes). Albeit the importance of genetic changes, some epigenetic changes can be stabilized and transmitted as epialleles into the progeny, which are subject to natural selection, adaptation, and domestication. Here we review recent advances for general and specific roles of epigenetic changes in the evolution of flowering plants and domestication of agricultural crops.

Functional Divergence between Subgenomes and Gene Pairs after Whole Genome Duplications

Gene loss following whole genome duplication (WGD) is often biased, with one subgenome retaining more ancestral genes and the other sustaining more gene deletions. While bias toward the greater expression of gene copies on one subgenome can explain bias in gene loss, this raises the question to what drives differences in gene expression levels between subgenomes. Differences in chromatin modifications and epigenetic markers between subgenomes in several model species are now being identified, providing an explanation for bias in gene expression between subgenomes. WGDs can be classified into duplications with higher, biased gene loss and bias in gene expression between subgenomes versus those with lower, unbiased rates of gene loss and an absence of detectable bias between subgenomes; however, the originally proposed link between these two classes and whether WGD results from an allo- or autopolyploid event is inconsistent with recent data from the allopolyploid Capsella bursa-pastoris. The gene balance hypothesis can explain bias in the functional categories of genes retained following WGD, the difference in gene loss rates between unbiased and biased WGDs, and how plant genomes have avoided being overrun with genes encoding dose-sensitive subunits of multiprotein complexes. Comparisons of gene expression patterns between retained transcription factor pairs in maize suggest the high degree of retention for WGD-derived pairs of transcription factors may instead be explained by the older duplication-degeneration-complementation model.

Asymmetrical changes of gene expression, small RNAs and chromatin in two resynthesized wheat allotetraploids

Polyploidy occurs in some animals and all flowering plants, including important crops such as wheat. The consequences of polyploidy in crops remain elusive, partly because their progenitors are unknown. Using two resynthesized wheat allotetraploids S^l S^l AA and AADD with known diploid progenitors, we analyzed mRNA and small RNA transcriptomes in the endosperm, compared transcriptomes between endosperm and root in AADD, and examined chromatin changes in the allotetraploids. In the endosperm, there were more non-additively expressed genes in S^l S^l AA than in AADD. In AADD, non-additively expressed genes were developmentally regulated, and the majority (62-70%) were repressed. The repressed genes in AADD included a group of histone methyltransferase gene homologs, which correlated with reduced histone H3K9me2 levels and activation of various transposable elements in AADD. In S^l S^l AA, there was a tendency for expression dominance of S^l over A homoeologs, but the histone methyltransferase gene homologs were additively expressed, correlating with insignificant changes in histone H3K9me2 levels. Moreover, more 24-nucleotide small inferring RNAs (siRNAs) in the A subgenome were disrupted in AADD than in S^l S^l AA, which were associated with expression changes of siRNA-associated genes. Our results indicate that asymmetrical changes in siRNAs, chromatin modifications, transposons and gene expression coincide with unstable AADD genomes and stable S^l S^l AA genomes, which could help explain the evolutionary trajectories of wheat allotetraploids formed by different progenitors.

Decomposing Additive Genetic Variance Revealed Novel Insights into Trait Evolution in Synthetic Hexaploid Wheat

Whole genome duplication (WGD) is an evolutionary phenomenon, which causes significant changes to genomic structure and trait architecture. In recent years, a number of studies decomposed the additive genetic variance explained by different sets of variants. However, they investigated diploid populations only and none of the studies examined any polyploid organism. In this research, we extended the application of this approach to polyploids, to differentiate the additive variance explained by the three subgenomes and seven sets of homoeologous chromosomes in synthetic allohexaploid wheat (SHW) to gain a better understanding of trait evolution after WGD. Our SHW population was generated by crossing improved durum parents (Triticum turgidum; 2n = 4x = 28, AABB subgenomes) with the progenitor species Aegilops tauschii (syn Ae. squarrosa, T. tauschii; 2n = 2x = 14, DD subgenome). The population was phenotyped for 10 fungal/nematode resistance traits as well as two abiotic stresses. We showed that the wild D subgenome dominated the additive effect and this dominance affected the A more than the B subgenome. We provide evidence that this dominance was not inflated by population structure, relatedness among individuals or by longer linkage disequilibrium blocks observed in the D subgenome within the population used for this study. The cumulative size of the three homoeologs of the seven chromosomal groups showed a weak but significant positive correlation with their cumulative explained additive variance. Furthermore, an average of 69% for each chromosomal group's cumulative additive variance came from one homoeolog that had the highest explained variance within the group across all 12 traits. We hypothesize that structural and functional changes during diploidization may explain chromosomal group relations as allopolyploids keep balanced dosage for many genes. Our results contribute to a better understanding of trait evolution mechanisms in polyploidy, which will facilitate the effective utilization of wheat wild relatives in breeding.

Asymmetric Evolution and Expansion of the NAC Transcription Factor in Polyploidized Cotton

Polyploidy in Gossypium hirsutum conferred different properties from its diploid ancestors under the regulation of transcription factors. The NAC transcription factor is a plant-specific family that can be related to plant growth and development. So far, little is known about the NAC family in cotton. This study identified 495 NAC genes in three cotton species and investigated the evolution and expansion of different genome-derived NAC genes in cotton. We revealed 15 distinct NAC subfamilies in cotton. Different subfamilies had different gene proportions, expansion rate, gene loss rate, and orthologous exchange rate. Paleohexaploidization (35%) and cotton-specific decaploidy (32%) might have primarily led to the expansion of the NAC family in cotton. Half of duplication events in G. hirsutum were inherited from its diploid ancestor, and others might have occurred after interspecific hybridization. In addition, NAC genes in the At and Dt subgenomes displayed asymmetric molecular evolution, as evidenced by their different gene loss rates, orthologous exchange, evolutionary rates, and expression levels. The dominant duplication event was different during the cotton evolutionary history. Different genome-derived NACs might have interacted with each other, which ultimately resulted in morphogenetic evolution. This study delineated the expansion and evolutionary history of the NAC family in cotton and illustrated the different fates of NAC genes during polyploidization.

The role of reproductive isolation in allopolyploid speciation patterns: empirical insights from the progenitors of common wheat

The ability to cause reproductive isolation often varies among individuals within a plant species. We addressed whether such polymorphism influenced speciation of the allopolyploid common wheat (Triticum aestivum L., AABBDD genome) by evaluating the expression of pre-pollination (outcrossing potential) and post-pollination (crossability) barriers in Aegilops tauschii Coss. (the D genome progenitor). In total, 201 Ae. tauschii accessions representing the entire natural habitat range of the species were used for anther length measurement and artificial crosses with a Triticum turgidum L. (the AB genome progenitor) tester. Intraspecific comparisons showed that both barriers were more strongly expressed in the TauL1 lineage than in the TauL2 lineage. The ability of Ae. tauschii to cause reproductive isolation in the hybridisation with T. turgidum might have markedly influenced common wheat’s speciation by inducing lineage-associated patterns of gene flow. The TauL2 accessions with high potential for natural hybridisation with T. turgidum clustered in the southern coastal Caspian region. This provided phenotypic support for the derivation of the D genome of common wheat from southern Caspian populations. The present study underscored the importance of approaches that incorporate the genealogical and geographic structure of the parental species’ reproductive isolation in understanding the mechanism of plant allopolyploid speciation.

Cytonuclear Variation of Rubisco in Synthesized Rice Hybrids and Allotetraploids

The allopolyploid speciation process faces the genomic challenge of stoichiometric disruption caused by merging biparental nuclear genomes with only one (usually maternal) of the two sets of progenitor cytoplasmic genomes. The photosynthetic protein 1,5-bisphosphate carboxylase/oxygenase (Rubisco) is composed of nuclear-encoded small subunits (SSUs) and plastome-encoded large subunits (LSUs), making it an ideal enzyme to explore the evolution process of cytonuclear accommodation. We investigated the variation of SSUs and their encoding genes in synthetic nascent rice ( L.) allotetraploid lineages, formed from the parental subspecies and of Asian rice. The evolution of genes in rice subspecies involves both mutation and concerted homogenization. Within reciprocal rice hybrids and allopolyploids, there was no consistent pattern of biased expression of alleles or homeologs, nor was there biased gene conversion favoring the maternal gene copies. Instead, we observed an apparently stochastic pattern of intergenomic gene conversions and biased expression of homeologs. We conclude that in young rice allopolyploids, cytonuclear coordination either is not selectively favored because of high parental sequence similarity or because there has been insufficient time for subtle selective effects to become observable.

Segmental allotetraploidy generates extensive homoeologous expression rewiring and phenotypic diversity at the population level in rice

Allopolyploidization, that is, concomitant merging and doubling of two or more divergent genomes in a common nucleus/cytoplasm, is known to instantly alter genomewide transcriptome dynamics, a phenomenon referred to as "transcriptomic shock." However, the immediate effects of transcriptomic alteration in generating phenotypic diversity at the population level remain underinvestigated. Here, we employed the MassARRAY-based Sequenom platform to assess and compare orthologous, allelic and homoeologous gene expression status in two tissues (leaf and root) of a set of randomly chosen individuals from populations of parental rice subspecies (indica and japonica), in vitro "hybrids" (parental mixes), reciprocal F1 hybrids and reciprocal tetraploids at the 5th-selfed generation (S5). We show that hybridization and whole genome duplication (WGD) have opposing effects on allelic and homoeologous expression in the F1 hybrids and tetraploids, respectively. Whereas hybridization exerts strong attenuating effects on allelic expression differences in diploid hybrids, WGD augments the intrinsic parental differences and generates extensive and variable homoeolog content which triggers diversification in expression patterning among the tetraploid plants. Coupled with the vast phenotypic diversity observed among the tetraploid individuals, our results provide experimental evidence in support of the notion that allopolyploidy catalyses rapid phenotypic diversification in higher plants. Our data further suggest that largely stochastic homoeolog content reshuffling rather than alteration in total expression level may be an important feature of evolution in young segmental allopolyploids, which underlies rapid expression diversity at the population level.

Diverse genome organization following 13 independent mesopolyploid events in Brassicaceae contrasts with convergent patterns of gene retention

Hybridization and polyploidy followed by genome-wide diploidization had a significant impact on the diversification of land plants. The ancient At-α whole-genome duplication (WGD) preceded the diversification of crucifers (Brassicaceae). Some genera and tribes also experienced younger, mesopolyploid WGDs concealed by subsequent genome diploidization. Here we tested if multiple base chromosome numbers originated due to genome diploidization after independent mesopolyploid WGDs and how diploidization affected post-polyploid gene retention. Sixteen species representing 10 Brassicaceae tribes were analyzed by comparative chromosome painting and/or whole-transcriptome analysis of gene age distributions and phylogenetic analyses of gene duplications. Overall, we found evidence for at least 13 independent mesopolyploidies followed by different degrees of diploidization across the Brassicaceae. New mesotetraploid events were uncovered for the tribes Anastaticeae, Iberideae and Schizopetaleae, and mesohexaploid WGDs for Cochlearieae and Physarieae. In contrast, we found convergent patterns of gene retention and loss among these independent WGDs. Our combined analyses of genomic data for Brassicaceae indicate that extant chromosome number variation in many plant groups, and especially monophyletic taxa with multiple base chromosome numbers, can result from clade-specific genome duplications followed by diploidization. Our observation of parallel gene retention and loss across multiple independent WGDs provides one of the first multi-species tests of the predictability of patterns of post-polyploid genome evolution.

The abundance of homoeologue transcripts is disrupted by hybridization and is partially restored by genome doubling in synthetic hexaploid wheat

Background

The formation of an allopolyploid is a two step process, comprising an initial wide hybridization event, which is later followed by a whole genome doubling. Both processes can affect the transcription of homoeologues. Here, RNA-Seq was used to obtain the genome-wide leaf transcriptome of two independent Triticum turgidum × Aegilops tauschii allotriploids (F1), along with their spontaneous allohexaploids (S1) and their parental lines. The resulting sequence data were then used to characterize variation in homoeologue transcript abundance.

Results

The hybridization event strongly down-regulated D-subgenome homoeologues, but this effect was in many cases reversed by whole genome doubling. The suppression of D-subgenome homoeologue transcription resulted in a marked frequency of parental transcription level dominance, especially with respect to genes encoding proteins involved in photosynthesis. Singletons (genes where no homoeologues were present) were frequently transcribed at both the allotriploid and allohexaploid plants.

Conclusions

The implication is that whole genome doubling helps to overcome the phenotypic weakness of the allotriploid, restoring a more favourable gene dosage in genes experiencing transcription level dominance in hexaploid wheat.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3558-0) contains supplementary material, which is available to authorized users.

Gene conversion events and variable degree of homogenization of rDNA loci in cultivars of Brassica napus

BACKGROUND AND AIMS

Brassica napus (AACC, 2n = 38, oilseed rape) is a relatively recent allotetraploid species derived from the putative progenitor diploid species Brassica rapa (AA, 2n = 20) and Brassica oleracea (CC, 2n = 18). To determine the influence of intensive breeding conditions on the evolution of its genome, we analysed structure and copy number of rDNA in 21 cultivars of B. napus, representative of genetic diversity.

METHODS

We used next-generation sequencing genomic approaches, Southern blot hybridization, expression analysis and fluorescence in situ hybridization (FISH). Subgenome-specific sequences derived from rDNA intergenic spacers (IGS) were used as probes for identification of loci composition on chromosomes.

KEY RESULTS

Most B. napus cultivars (18/21, 86 %) had more A-genome than C-genome rDNA copies. Three cultivars analysed by FISH ('Darmor', 'Yudal' and 'Asparagus kale') harboured the same number (12 per diploid set) of loci. In B. napus 'Darmor', the A-genome-specific rDNA probe hybridized to all 12 rDNA loci (eight on the A-genome and four on the C-genome) while the C-genome-specific probe showed weak signals on the C-genome loci only. Deep sequencing revealed high homogeneity of arrays suggesting that the C-genome genes were largely overwritten by the A-genome variants in B. napus 'Darmor'. In contrast, B. napus 'Yudal' showed a lack of gene conversion evidenced by additive inheritance of progenitor rDNA variants and highly localized hybridization signals of subgenome-specific probes on chromosomes. Brassica napus 'Asparagus kale' showed an intermediate pattern to 'Darmor' and 'Yudal'. At the expression level, most cultivars (95 %) exhibited stable A-genome nucleolar dominance while one cultivar ('Norin 9') showed co-dominance.

CONCLUSIONS

The B. napus cultivars differ in the degree and direction of rDNA homogenization. The prevalent direction of gene conversion (towards the A-genome) correlates with the direction of expression dominance indicating that gene activity may be needed for interlocus gene conversion.