Evolution and expression of homeologous loci in Tragopogon miscellus (Asteraceae), a recent and reciprocally formed allopolyploid

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.

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).

Frequent allopolyploidy with distant progenitors in the moss genera Physcomitrium and Entosthodon (Funariaceae) identified via subgenome phasing of targeted nuclear genes

Allopolyploids represent a new frontier in species discovery among embryophytes. Within mosses, allopolyploid discovery is challenged by low morphological complexity. The rapid expansion of sequencing approaches in addition to computational developments to identifying genome merger and whole-genome duplication using variation among nuclear loci representing homeologs has allowed for increased allopolyploid discovery among mosses. Here, we test a novel approach to phasing homeologs within loci and phasing loci across subgenomes, or subgenome assignment, called Homologizer, in the family Funariaceae. We confirm the intergeneric hybrid nature of Entosthodon hungaricus, and the allopolyploid origin of Physcomitrium eurystomum and one population of Physcomitrium collenchymatum. We also reveal that hybridization gave rise to Physcomitrium immersum, as well as to yet unrecognized lineages sharing the phenotype of Physcomitrium pyriforme and Physcomitrium sphaericum. Our findings demonstrate the utility of our approach when working with polyploid genomes, and its value in identifying progenitor species using target capture data.

Protein nonadditive expression and solubility contribute to heterosis in Arabidopsis hybrids and allotetraploids

Hybrid vigor or heterosis has been widely applied in agriculture and extensively studied using genetic and gene expression approaches. However, the biochemical mechanism underlying heterosis remains elusive. One theory suggests that a decrease in protein aggregation may occur in hybrids due to the presence of protein variants between parental alleles, but it has not been experimentally tested. Here, we report comparative analysis of soluble and insoluble proteomes in Arabidopsis intraspecific and interspecific hybrids or allotetraploids formed between A. thaliana and A. arenosa. Both allotetraploids and intraspecific hybrids displayed nonadditive expression (unequal to the sum of the two parents) of the proteins, most of which were involved in biotic and abiotic stress responses. In the allotetraploids, homoeolog-expression bias was not observed among all proteins examined but accounted for 17-20% of the nonadditively expressed proteins, consistent with the transcriptome results. Among expression-biased homoeologs, there were more A. thaliana-biased than A. arenosa-biased homoeologs. Analysis of the insoluble and soluble proteomes revealed more soluble proteins in the hybrids than their parents but not in the allotetraploids. Most proteins in ribosomal biosynthesis and in the thylakoid lumen, membrane, and stroma were in the soluble fractions, indicating a role of protein stability in photosynthetic activities for promoting growth. Thus, nonadditive expression of stress-responsive proteins and increased solubility of photosynthetic proteins may contribute to heterosis in Arabidopsis hybrids and allotetraploids and possibly hybrid crops.

Patterns, mechanisms, and consequences of homoeologous exchange in allopolyploid angiosperms: a genomic and epigenomic perspective

Allopolyploids result from hybridization between different evolutionary lineages coupled with genome doubling. Homoeologous chromosomes (chromosomes with common shared ancestry) may undergo recombination immediately after allopolyploid formation and continue over successive generations. The outcome of this meiotic pairing behavior is dynamic and complex. Homoeologous exchanges (HEs) may lead to the formation of unbalanced gametes, reduced fertility, and selective disadvantage. By contrast, HEs could act as sources of novel evolutionary substrates, shifting the relative dosage of parental gene copies, generating novel phenotypic diversity, and helping the establishment of neo-allopolyploids. However, HE patterns vary among lineages, across generations, and even within individual genomes and chromosomes. The causes and consequences of this variation are not fully understood, though interest in this evolutionary phenomenon has increased in the last decade. Recent technological advances show promise in uncovering the mechanistic basis of HEs. Here, we describe recent observations of the common patterns among allopolyploid angiosperm lineages, underlying genomic and epigenomic features, and consequences of HEs. We identify critical research gaps and discuss future directions with far-reaching implications in understanding allopolyploid evolution and applying them to the development of important phenotypic traits of polyploid crops.

Meeting in Liguria: hybridisation between Apennine endemic Euphorbia barrelieri and western Mediterranean E. nicaeensis led to the allopolyploid origin of E. ligustica

The Mediterranean Basin is renowned for its extremely rich biota and is considered as one of the 25 Global Biodiversity Hotspots, but its diversity is not homogeneously distributed. Outstanding in the number of (endemic) species are the Ligurian Alps (Italy). At the foot of the Ligurian Alps, little above the Mediterranean Sea, a disjunct occurrence of Italian endemic Euphorbia barrelieri was reported. Using an array of integrative methods ranging from cytogenetic (chromosome number and relative genome size estimation), over phylogenetic approaches (plastid, ITS and RAD sequencing) to multivariate morphometrics we disentangled the origin of these populations that were shown to be tetraploid. We performed phylogenetic analyses of the nuclear ITS and plastid regions of a broad taxonomic sampling of Euphorbia sect. Pithyusa to identify possible species involved in the origin of the tetraploid populations and then applied various analyses of RADseq data to identify the putative parental species. Our results have shown that the Ligurian populations of E. barrelieri are of allotetraploid origin that involved E. barrelieri and western Mediterranean E. nicaeensis as parental species. We thus describe a new species, E. ligustica, and hypothesise that its adaptation to similar environments in which E. barrelieri occurs, triggered development of similar morphology, whereas its genetic composition appears to be closer to that of E. nicaeensis. Our study emphasises the importance of polyploidisation for plant diversification, highlights the value of the Ligurian Alps as a hotspot of biodiversity and endemism and underlines the importance of integrative taxonomic approaches in uncovering cryptic diversity.

Effect of Life-History Traits and Habitat Condition on Genetic Diversity between Invasive and Native Plant Populations

Plant invasions have a huge impact on the health of ecosystems and human well-being. The invasion risk varies with the introduction pathway, the propagule pressure, and the genetic diversity of the founding population. We performed a systematic review and meta-analysis of 30 studies reporting the genetic diversity of 31 plant species in their invasive and native ranges. We evaluated if patterns of genetic diversity differ between ranges and whether these responses are influenced by life-history traits, hybridization, polyploidization, and habitat condition. We found that invasive populations had significantly lower genetic diversity and higher inbreeding than native populations. In fragmented and degraded habitats, the genetic diversity of invaders was lower, but inbreeding was not affected. Polyploid invaders with hybrid capacity also showed lower genetic diversity. Invasive herbs with vegetative propagation were more sensitive to the loss of genetic diversity and had higher levels of inbreeding. Our synthesis showed that the genetic response in the invaded range could result from historical processes, such as founder and bottleneck events. Traits such as selfing are more likely to preserve the signatures of founder events and influence the genetic diversity in invasive populations. Additionally, clonality seems to be the predominant reproduction system in the invaded range.

Photosynthetic Efficiency and Glyco-Metabolism Changes in Artificial Triploid Loquats Contribute to Heterosis Manifestation

Previous studies indicated that extensive genetic variations could be generated due to polyploidy, which is considered to be closely associated with the manifestation of polyploid heterosis. Our previous studies confirmed that triploid loquats demonstrated significant heterosis, other than the ploidy effect, but the underlying mechanisms are largely unknown. This study aimed to overcome the narrow genetic distance of loquats, increase the genetic variation level of triploid loquats, and systematically illuminate the heterosis mechanisms of triploid loquats derived from two cross combinations. Here, inter-simple sequence repeats (ISSRs) and simple sequence repeats (SSRs) were adopted for evaluating the genetic diversity, and transcriptome sequencing (RNA-Seq) was performed to investigate gene expression as well as pathway changes in the triploids. We found that extensive genetic variations were produced during the formation of triploid loquats. The polymorphism ratios of ISSRs and SSRs were 43.75% and 19.32%, respectively, and almost all their markers had a PIC value higher than 0.5, suggesting that both ISSRs and SSRs could work well in loquat assisted breeding. Furthermore, our results revealed that by broadening the genetic distance between the parents, genetic variations in triploids could be promoted. Additionally, RNA-Seq results suggested that numerous genes differentially expressed between the triploids and parents were screened out. Moreover, KEGG analyses revealed that “photosynthetic efficiency” and “glyco-metabolism” were significantly changed in triploid loquats compared with the parents, which was consistent with the results of physiological indicator analyses, leaf micro-structure observations, and qRT-PCR validation. Collectively, our results suggested that extensive genetic variations occurred in the triploids and that the changes in the “photosynthetic efficiency” as well as “glyco-metabolism” of triploids might have further resulted in heterosis manifestation in the triploid loquats.

Expression Level Dominance and Homeolog Expression Bias Upon Cold Stress in the F1 Hybrid Between the Invasive Sphagneticola trilobata and the Native S. calendulacea in South China, and Implications for Its Invasiveness

The role of hybridization is significant in biological invasion, and thermotolerance is a trait critical to range expansions. The South American Sphagneticola trilobata is now widespread in South China, threatening the native S. calendulacea by competition and hybridization. Furthermore, upon formation, their F1 hybrid can quickly replace both parents. In this study, the three taxa were used as a model to investigate the consequences of hybridization on cold tolerance, particularly the effect of subgenome dominance in the hybrid. Upon chilling treatments, physiological responses and transcriptome profiles were compared across different temperature points to understand their differential responses to cold. While both parents showed divergent responses, the hybrid’s responses showed an overall resemblance to S. calendulacea, but the contribution of homeolog expression bias to cold stress was not readily evident in the F1 hybrid possibly due to inherent bias that comes with the sampling location. Our findings provided insights into the role of gene expression in differential cold tolerance, and further contribute to predicting the invasive potential of other hybrids between S. trilobata and its congeners around the world.

Rapid and Synchronous Breeding of Cytoplasmic Male Sterile and Maintainer Line Through Mitochondrial DNA Rearrangement Using Doubled Haploid Inducer in Brassica napus

When homozygously fertile plants were induced using doubled haploid (DH) induction lines Y3380 and Y3560, the morphology of the induced F1 generation was basically consistent with the female parent, but the fertility was separated, showing characteristics similar to cytoplasmic male sterile (CMS) and maintainer lines. In this study, the morphology, fertility, ploidy, and cytoplasm genotype of the induced progeny were identified, and the results showed that the sterile progeny was polima cytoplasm sterile (pol CMS) and the fertile progeny was nap cytoplasm. The molecular marker and test-cross experimental results showed that the fertile progeny did not carry the restorer gene of pol CMS and the genetic distance between the female parent and the offspring was 0.002. This suggested that those inductions which produced sterile and fertile progeny were coordinated to CMS and maintainer lines. Through the co-linearity analysis of the mitochondrial DNA (mtDNA), it was found that the rearrangement of mtDNA by DH induction was the key factor that caused the transformation of fertility (nap) into sterility (pol). Also, when heterozygous females were induced with DH induction lines, the induction F2 generation also showed the segregation of fertile and sterile lines, and the genetic distance between sterile and fertile lines was approximately 0.075. Therefore, the induction line can induce different types of female parents, and the breeding of the sterile line and the maintainer line can be achieved through the rapid synchronization of sister crosses and self-crosses. The induction of DH inducer in B. napus can provide a new model for the innovation of germplasm resources and open up a new way for its application.

An Eruption of LTR Retrotransposons in the Autopolyploid Genomes of Chrysanthemum nankingense (Asteraceae)

Whole genome duplication, associated with the induction of widespread genetic changes, has played an important role in the evolution of many plant taxa. All extant angiosperm species have undergone at least one polyploidization event, forming either an auto- or allopolyploid organism. Compared with allopolyploidization, however, few studies have examined autopolyploidization, and few studies have focused on the response of genetic changes to autopolyploidy. In the present study, newly synthesized C. nankingense autotetraploids (Asteraceae) were employed to characterize the genome shock following autopolyploidization. Available evidence suggested that the genetic changes primarily involved the loss of old fragments and the gain of novel fragments, and some novel sequences were potential long terminal repeat (LTR) retrotransposons. As Ty1-copia and Ty3-gypsy elements represent the two main superfamilies of LTR retrotransposons, the dynamics of Ty1-copia and Ty3-gypsy were evaluated using RT-PCR, transcriptome sequencing, and LTR retrotransposon-based molecular marker techniques. Additionally, fluorescence in situ hybridization(FISH)results suggest that autopolyploidization might also be accompanied by perturbations of LTR retrotransposons, and emergence retrotransposon insertions might show more rapid divergence, resulting in diploid-like behaviour, potentially accelerating the evolutionary process among progenies. Our results strongly suggest a need to expand the current evolutionary framework to include a genetic dimension when seeking to understand genomic shock following autopolyploidization in Asteraceae.

Genomic and Meiotic Changes Accompanying Polyploidization

Hybridization and polyploidy have been considered as significant evolutionary forces in adaptation and speciation, especially among plants. Interspecific gene flow generates novel genetic variants adaptable to different environments, but it is also a gene introgression mechanism in crops to increase their agronomical yield. An estimate of 9% of interspecific hybridization has been reported although the frequency varies among taxa. Homoploid hybrid speciation is rare compared to allopolyploidy. Chromosome doubling after hybridization is the result of cellular defects produced mainly during meiosis. Unreduced gametes, which are formed at an average frequency of 2.52% across species, are the result of altered spindle organization or orientation, disturbed kinetochore functioning, abnormal cytokinesis, or loss of any meiotic division. Meiotic changes and their genetic basis, leading to the cytological diploidization of allopolyploids, are just beginning to be understood especially in wheat. However, the nature and mode of action of homoeologous recombination suppressor genes are poorly understood in other allopolyploids. The merger of two independent genomes causes a deep modification of their architecture, gene expression, and molecular interactions leading to the phenotype. We provide an overview of genomic changes and transcriptomic modifications that particularly occur at the early stages of allopolyploid formation.

The Role of Interspecific Hybridisation in Adaptation and Speciation: Insights From Studies in Senecio

Hybridisation is well documented in many species, especially plants. Although hybrid populations might be short-lived and do not evolve into new lineages, hybridisaiton could lead to evolutionary novelty, promoting adaptation and speciation. The genus Senecio (Asteraceae) has been actively used to unravel the role of hybridisation in adaptation and speciation. In this article, we first briefly describe the process of hybridisation and the state of hybridisation research over the years. We then discuss various roles of hybridisation in plant adaptation and speciation illustrated with examples from different Senecio species, but also mention other groups of organisms whenever necessary. In particular, we focus on the genomic and transcriptomic consequences of hybridisation, as well as the ecological and physiological aspects from the hybrids' point of view. Overall, this article aims to showcase the roles of hybridisation in speciation and adaptation, and the research potential of Senecio, which is part of the ecologically and economically important family, Asteraceae.

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.

Untangling structural factors driving genome stabilization in nascent Brassica napus allopolyploids

Allopolyploids have globally higher fitness than their diploid progenitors; however, by comparison, most resynthesized allopolyploids have poor fertility and highly unstable genome. Elucidating the evolutionary processes promoting genome stabilization and fertility is thus essential to comprehend allopolyploid success. Using the Brassica model, we mimicked the speciation process of a nascent allopolyploid species by resynthesizing allotetraploid Brassica napus and systematically selecting for euploid individuals over eight generations in four independent allopolyploidization events with contrasted genetic backgrounds, cytoplasmic donors, and polyploid formation type. We evaluated the evolution of meiotic behavior and fertility and identified rearrangements in S₁ to S₉ lineages to explore the positive consequences of euploid selection on B. napus genome stability. Recurrent selection of euploid plants for eight generations drastically reduced the percentage of aneuploid progenies as early as the fourth generation, concomitantly with a decrease in number of newly fixed homoeologous rearrangements. The consequences of homoeologous rearrangements on meiotic behavior and seed number depended strongly on the genetic background and cytoplasm donor. The combined use of both self-fertilization and recurrent euploid selection allowed identification of genomic regions associated with fertility and meiotic behavior, providing complementary evidence to explain B. napus speciation success.

Imprints of independent allopolyploid formations on patterns of gene expression in two sibling yarrow species (Achillea, Asteraceae)

BACKGROUND

Polyploid species often originate recurrently. While this is well known, there is little information on the extent to which distinct allotetraploid species formed from the same parent species differ in gene expression. The tetraploid yarrow species Achillea alpina and A. wilsoniana arose independently from allopolyploidization between diploid A. acuminata and A. asiatica. The genetics and geography of these origins are clear from previous studies, providing a solid basis for comparing gene expression patterns of sibling allopolyploid species that arose independently.

RESULTS

We conducted comparative RNA-sequencing analyses on the two Achillea tetraploid species and their diploid progenitors to evaluate: 1) species-specific gene expression and coexpression across the four species; 2) patterns of inheritance of parental gene expression; 3) parental contributions to gene expression in the allotetraploid species, and homeolog expression bias. Diploid A. asiatica showed a higher contribution than diploid A. acuminata to the transcriptomes of both tetraploids and also greater homeolog bias in these transcriptomes, possibly reflecting a maternal effect. Comparing expressed genes in the two allotetraploids, we found expression of ca. 30% genes were species-specific in each, which were most enriched for GO terms pertaining to "defense response". Despite species-specific and differentially expressed genes between the two allotetraploids, they display similar transcriptome changes in comparison to their diploid progenitors.

CONCLUSION

Two independently originated Achillea allotetraploid species exhibited difference in gene expression, some of which must be related to differential adaptation during their post-speciation evolution. On the other hand, they showed similar expression profiles when compared to their progenitors. This similarity might be expected when pairs of merged diploid genomes in tetraploids are similar, as is the case in these two particular allotetraploids.

Evolution of Homeologous Gene Expression in Polyploid Wheat

Polyploidization has played a prominent role in the evolutionary history of plants. Two recent and sequential allopolyploidization events have resulted in the formation of wheat species with different ploidies, and which provide a model to study the effects of polyploidization on the evolution of gene expression. In this study, we identified differentially expressed genes (DEGs) between four BBAA tetraploid wheats of three different ploidy backgrounds. DEGs were found to be unevenly distributed among functional categories and duplication modes. We observed more DEGs in the extracted tetraploid wheat (ETW) than in natural tetraploid wheats (TD and TTR13) as compared to a synthetic tetraploid (AT2). Furthermore, DEGs showed higher Ka/Ks ratios than those that did not show expression changes (non-DEGs) between genotypes, indicating DEGs and non-DEGs experienced different selection pressures. For A-B homeolog pairs with DEGs, most of them had only one differentially expressed copy, however, when both copies of a homeolog pair were DEGs, the A and B copies were more likely to be regulated to the same direction. Our results suggest that both cis- and inter-subgenome trans-regulatory changes are important drivers in the evolution of homeologous gene expression in polyploid wheat, with ploidy playing a significant role in the process.

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.

The mitogen activated protein kinase (MAPK) gene family functions as a cohort during the Glycine max defense response to Heterodera glycines

Mitogen activated protein kinases (MAPKs) play important signal transduction roles. However, little is known regarding how they influence the gene expression of other family members and the relationship to a biological process, including the Glycine max defense response to Heterodera glycines. Transcriptomics have identified MAPK gene expression occurring within root cells undergoing a defense response to a pathogenic event initiated by H. glycines in the allotetraploid Glycine max. Functional analyses are presented for its 32 MAPKs revealing 9 have a defense role, including homologs of Arabidopsis thaliana MAPK (MPK) MPK2, MPK3, MPK4, MPK5, MPK6, MPK13, MPK16 and MPK20. Defense signaling occurring through pathogen activated molecular pattern (PAMP) triggered immunity (PTI) and effector triggered immunity (ETI) have been determined in relation to these MAPKs. Five different types of gene expression relate to MAPK expression, influencing PTI and ETI gene expression and proven defense genes including an ABC-G transporter, 20S membrane fusion particle components, glycoside biosynthesis, carbon metabolism, hemicellulose modification, transcription and secretion. The experiments show MAPKs broadly influence defense MAPK gene expression, including the co-regulation of parologous MAPKs and reveal its relationship to proven defense genes. The experiments reveal each defense MAPK induces the expression of a G. max homolog of a PATHOGENESIS RELATED1 (PR1), itself shown to function in defense in the studied pathosystem.

Differential expression of HPG-axis genes in autotetraploids derived from red crucian carp Carassius auratus red var., ♀ × blunt snout bream Megalobrama amblycephala, ♂

Autotetraploid fish (4n = 200, abbreviated as 4nRR), which reach sexual maturity at 1 year of age, were derived from the whole genome duplication of red crucian carp Carassius auratus red var. (RCC; 2n = 100) and possess four sets of chromosomes from RCC. The histological features of the gonads showed that the RCC and 4nRR both possessed normal gonadal structure and could arrive at maturation. To understand the expression characteristics of genes related to reproductive development in the autotetraploid fish, we analysed the nucleotide sequence and expression characteristics of the gnrh2, gthb and gthr genes, which are the pivotal genes of the hypothalamic-pituitary-gonadal (HPG) axis. We found that the gnrh2, gthb and gthr genes in 4nRR share remarkable homology with RCC, but there were obvious differences in expression levels between 4nRR and RCC. These results demonstrate that autotetraploidization can lead to gene expression changes. This study provides insights into the molecular mechanism underlying the reproductive development of autotetraploid fish and is expected to be of great significance for subsequent research on polyploidization.

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.

Extensive genetic and DNA methylation variation contribute to heterosis in triploid loquat hybrids

We aim to overcome the unclear origin of the loquat and elucidate the heterosis mechanism of the triploid loquat. Here we investigated the genetic and epigenetic variations between the triploid plant and its parental lines using amplified fragment length polymorphism (AFLP) and methylation-sensitive amplified fragment length polymorphism (MSAP) analyses. We show that in addition to genetic variations, extensive DNA methylation variation occurred during the formation process of triploid loquat, with the triploid hybrid having increased DNA methylation compared to the parents. Furthermore, a correlation existed between genetic variation and DNA methylation remodeling, suggesting that genome instability may lead to DNA methylation variation or vice versa. Sequence analysis of the MSAP bands revealed that over 53% of them overlap with protein-coding genes, which may indicate a functional role of the differential DNA methylation in gene regulation and hence heterosis phenotypes. Consistent with this, the genetic and epigenetic alterations were associated closely to the heterosis phenotypes of triploid loquat, and this association varied for different traits. Our results suggested that the formation of triploid is accompanied by extensive genetic and DNA methylation variation, and these changes contribute to the heterosis phenotypes of the triploid loquats from the two cross lines.

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.

Homoeologous exchanges cause extensive dosage-dependent gene expression changes in an allopolyploid crop

Structural variation is a major source of genetic diversity and an important substrate for selection. In allopolyploids, homoeologous exchanges (i.e. between the constituent subgenomes) are a very frequent type of structural variant. However, their direct impact on gene content and gene expression had not been determined. Here, we used a tissue-specific mRNA-Seq dataset to measure the consequences of homoeologous exchanges (HE) on gene expression in Brassica napus, a representative allotetraploid crop. We demonstrate that expression changes are proportional to the change in gene copy number triggered by the HEs. Thus, when homoeologous gene pairs have unbalanced transcriptional contributions before the HE, duplication of one copy does not accurately compensate for loss of the other and combined homoeologue expression also changes. These effects are, however, mitigated over time. This study sheds light on the origins, timing and functional consequences of homeologous exchanges in allopolyploids. It demonstrates that the interplay between new structural variation and the resulting impacts on gene expression, influences allopolyploid genome evolution.

Karyotypic variation and pollen stainability in resynthesized allopolyploids Tragopogon miscellus and T. mirus

PREMISE OF THE STUDY

Polyploidy has extensively shaped the evolution of plants, but the early stages of polyploidy are still poorly understood. The neoallopolyploid species Tragopogon mirus and T. miscellus are both characterized by widespread karyotypic variation, including frequent aneuploidy and intergenomic translocations. Our study illuminates the origins and early impacts of this variation by addressing two questions: How quickly does karyotypic variation accumulate in Tragopogon allopolyploids following whole-genome duplication (WGD), and how does the fertility of resynthesized Tragopogon allopolyploids evolve shortly after WGD?

METHODS

We used genomic in situ hybridization and lactophenol-cotton blue staining to estimate the karyotypic variation and pollen stainability, respectively, of resynthesized T. mirus and T. miscellus during the first five generations after WGD.

KEY RESULTS

Widespread karyotypic variation developed quickly in synthetics and resembled that of naturally occurring T. mirus and T. miscellus by generation S₄ . Pollen stainability in resynthesized allopolyploids was consistently lower than that of natural T. mirus and T. miscellus, as well as their respective diploid progenitor species. Logistic regression showed that mean pollen stainability increased slightly over four generations in resynthesized T. mirus but remained at equivalent levels in T. miscellus.

CONCLUSIONS

Our results clarify some of the changes that occur in T. mirus and T. miscellus immediately following their origin, most notably the rapid onset of karyotypic variation within these species immediately following WGD.

Evolution of Gene Expression Balance Among Homeologs of Natural Polyploids

Polyploidy is a major evolutionary process in eukaryotes, yet the expression balance of homeologs in natural polyploids is largely unknown. To study this expression balance, the expression patterns of 2180 structurally well-characterized genes of wheat were studied, of which 813 had the expected three copies and 375 had less than three. Copy numbers of the remaining 992 ranged from 4 to 14, including homeologs, orthologs, and paralogs. Of the genes with three structural copies corresponding to homeologs, 55% expressed from all three, 38% from two, and the remaining 7% expressed from only one of the three copies. Homeologs of 76-87% of the genes showed differential expression patterns in different tissues, thus have evolved different gene expression controls, possibly resulting in novel functions. Homeologs of 55% of the genes showed tissue-specific expression, with the largest percentage (14%) in the anthers and the smallest (7%) in the pistils. The highest number (1.72/3) of homeologs/gene expression was in the roots and the lowest (1.03/3) in the anthers. As the expression of homeologs changed with changes in structural copy number, about 30% of the genes showed dosage dependence. Chromosomal location also impacted expression pattern as a significantly higher proportion of genes in the proximal regions showed expression from all three copies compared to that present in the distal regions.

Correlation analysis of the mRNA and miRNA expression profiles in the nascent synthetic allotetraploid Raphanobrassica

Raphanobrassica is an allopolyploid species derived from inter-generic hybridization that combines the R genome from R. sativus and the C genome from B. oleracea var. alboglabra. In the present study, we used a high-throughput sequencing method to identify the mRNA and miRNA profiles in Raphanobrassica and its parents. A total of 33,561 mRNAs and 283 miRNAs were detected, 9,209 mRNAs and 134 miRNAs were differentially expressed respectively, 7,633 mRNAs and 39 miRNAs showed ELD expression, 5,219 mRNAs and 57 miRNAs were non-additively expressed in Raphanobrassica. Remarkably, differentially expressed genes (DEGs) were up-regulated and maternal bias was detected in Raphanobrassica. In addition, a miRNA-mRNA interaction network was constructed based on reverse regulated miRNA-mRNAs, which included 75 miRNAs and 178 mRNAs, 31 miRNAs were non-additively expressed target by 13 miRNAs. The related target genes were significantly enriched in the GO term 'metabolic processes'. Non-additive related target genes regulation is involved in a range of biological pathways, like providing a driving force for variation and adaption in this allopolyploid. The integrative analysis of mRNA and miRNA profiling provides more information to elucidate gene expression mechanism and may supply a comprehensive and corresponding method to study genetic and transcription variation of allopolyploid.

Rapid genomic changes in allopolyploids of Carassius auratus red var. (♀) × Megalobrama amblycephala (♂)

To better understand genomic changes in the early generations after polyploidisation, we examined the chromosomal consequences of genomic merger in allotetraploid hybrids (4 nF₁) (AABB, 4n = 148) of Carassius auratus red var. (RCC) (AA, 2n = 100) (♀) × Megalobrama amblycephala (BSB) (BB, 2n = 48) (♂). Complete loss of the paternal 5S rDNA sequence and the expected number of maternal chromosomal loci were found in 4 nF₁, suggesting directional genomic changes occurred in the first generations after polyploidisation. Recent studies have reported instability of newly established allotetraploid genomes. To assess this in the newly formed 4 nF₁ genome, we performed fluorescence in situ hybridisation on an allotetraploid gynogenetic hybrid (4 nG) (AABB, 4n = 148) and an allopentaploid hybrid (5 nH) (AABBB, 5n = 172) from 4 nF₁ (♀) × BSB (♂) with 5S rDNA gene and centromere probes from RCC, the original diploid parent. The expected numbers of maternal chromosomal loci were found in 4 nG, while chromosomal locus deletions and chromosome recombinations were detected in 5 nH. These observations suggest that abnormal meiosis did not lead to obvious genomic changes in the newly established allotetraploid genomes, but hybridisation with the original diploid parent resulted in obvious genomic changes in the newly established allotetraploid genomes, as was found for the maternal genome.

Critical factors in the establishment of allopolyploids

PREMISE OF THE STUDY

The growth and spread of new polyploid populations have been explained in terms of fitness advantages over their diploid progenitors. However, a fitness advantage is not sufficient to insure the establishment of a polyploid; it must also overcome the obstacles of demographic stochasticity and minority disadvantage. Several studies have addressed the population dynamics of autopolyploids, but the present study is the first to consider allopolyploids, which are affected by more factors than autopolyploids.

METHODS

We constructed a population dynamic model of four types of plants (two parent species, hybrids, allopolyploids) that also included an explicit breeding system.

KEY RESULTS

The numbers of plants of each type were the most important factors determining whether the new allopolyploid would become established. More polyploid plants greatly increased the likelihood of polyploid persistence. More plants of the parent species and more hybrids resulted in more polyploids being produced. The model parameters with the most effect on polyploid establishment were potential population size (K), individual plant fecundity, and niche separation (α). The most important breeding system parameters were selfing rates, which mitigated minority disadvantage imposed by pollen limitation.

CONCLUSIONS

The importance of population sizes, and the parameters that controlled them, in overcoming demographic stochasticity parallels the well-recognized role of propagule pressure in determining the success of invasive species. We modeled the establishment of a new allopolyploid; analogous considerations would affect the establishment of a new autopolyploid. The critical role of population sizes in polyploid establishment should be more widely recognized.

Comparative structural analysis of Bru1 region homeologs in Saccharum spontaneum and S. officinarum

BACKGROUND

Sugarcane is a major sugar and biofuel crop, but genomic research and molecular breeding have lagged behind other major crops due to the complexity of auto-allopolyploid genomes. Sugarcane cultivars are frequently aneuploid with chromosome number ranging from 100 to 130, consisting of 70-80 % S. officinarum, 10-20 % S. spontaneum, and 10 % recombinants between these two species. Analysis of a genomic region in the progenitor autoploid genomes of sugarcane hybrid cultivars will reveal the nature and divergence of homologous chromosomes.

RESULTS

To investigate the origin and evolution of haplotypes in the Bru1 genomic regions in sugarcane cultivars, we identified two BAC clones from S. spontaneum and four from S. officinarum and compared to seven haplotype sequences from sugarcane hybrid R570. The results clarified the origin of seven homologous haplotypes in R570, four haplotypes originated from S. officinarum, two from S. spontaneum and one recombinant.. Retrotransposon insertions and sequences variations among the homologous haplotypes sequence divergence ranged from 18.2 % to 60.5 % with an average of 33.7 %. Gene content and gene structure were relatively well conserved among the homologous haplotypes. Exon splitting occurred in haplotypes of the hybrid genome but not in its progenitor genomes. Tajima's D analysis revealed that S. spontaneum hapotypes in the Bru1 genomic regions were under strong directional selection. Numerous inversions, deletions, insertions and translocations were found between haplotypes within each genome.

CONCLUSIONS

This is the first comparison among haplotypes of a modern sugarcane hybrid and its two progenitors. Tajima's D results emphasized the crucial role of this fungal disease resistance gene for enhancing the fitness of this species and indicating that the brown rust resistance gene in R570 is from S. spontaneum. Species-specific InDel, sequences similarity and phylogenetic analysis of homologous genes can be used for identifying the origin of S. spontaneum and S. officinarum haplotype in Saccharum hybrids. Comparison of exon splitting among the homologous haplotypes suggested that the genome rearrangements in Saccharum hybrids after hybridization. The combined minimum difference at 19.5 % among homologous chromosomes in S. officinarum would be sufficient for proper genome assembly of this autopolyploid genome. Retrotransposon insertions and sequences variations among the homologous haplotypes sequence divergence may allow sequencing and assembling the autopolyploid Saccharum genomes and the auto-allopolyploid hybrid genomes using whole genome shotgun sequencing.

The origin of exon 3 skipping of paternal GLOBOSA pre-mRNA in some Nicotiana tabacum lines correlates with a point mutation of the very last nucleotide of the exon

In plants, genome duplication followed by genome diversification and selection is recognized as a major evolutionary process. Rapid epigenetic and genetic changes that affect the transcription of parental genes are frequently observed after polyploidization. The pattern of alternative splicing is also frequently altered, yet the related molecular processes remain largely unresolved. Here, we study the inheritance and expression of parental variants of three floral organ identity genes in allotetraploid tobacco. DEFICIENS and GLOBOSA are B-class genes, and AGAMOUS is a C-class gene. Parental variants of these genes were found to be maintained in the tobacco genome, and the respective mRNAs were present in flower buds in comparable amounts. However, among five tobacco cultivars, we identified two in which the majority of paternal GLOBOSA pre-mRNA transcripts undergo exon 3 skipping, producing an mRNA with a premature termination codon. At the DNA level, we identified a G-A transition at the very last position of exon 3 in both cultivars. Although alternative splicing resulted in a dramatic decrease in full-length paternal GLOBOSA mRNA, no phenotypic effect was observed. Our finding likely serves as an example of the initiation of homoeolog diversification in a relatively young polyploid genome.

Structure, phylogeny, allelic haplotypes and expression of sucrose transporter gene families in Saccharum

BACKGROUND

Sugarcane is an economically important crop contributing to about 80% of the world sugar production. Increasing efforts in molecular biological studies have been performed for improving the sugar yield and other relevant important agronomic traits. However, due to sugarcane's complicated genomes, it is still challenging to study the genetic basis of traits, such as sucrose accumulation. Sucrose transporters (SUTs) are critical for both phloem loading in source tissue and sucrose uptaking in sink tissue, and are considered to be the control points for regulating sucrose storage. However, no genomic study for sugarcane sucrose transporter (SsSUT) families has been reported up to date.

RESULTS

By using comparative genomics and bacterial artificial chromosomes (BACs), six SUT genes were identified and characterized in S. spontaenum. Phylogenetic analyses revealed that the two pairs SsSUTs (SsSUT1/SsSUT3 and SsSUT5/SsSUT6) could be clustered together into two separate monocot specific SUT groups, while SsSUT2 and SsSUT4 were separated into the other two groups, with members from both dicot and monocot species. Gene structure comparison demonstrated that the number and position of exons/introns in SUTs were highly conserved among the close orthologs; in contrast, there were variations among the paralogous SUTs in Sacchuarm. Though with the high polyploidy level, gene allelic haplotype comparative analysis showed that the examined four SsSUT members exhibited conservations of gene structures and amino acid sequences among the allelic haplotypes accompanied by variations of intron sizes. Gene expression analyses were performed for tissues from seedlings under drought stress and mature plants of three Saccharum species (S.officinarnum, S.spotaneum and S.robustum). Both SUT1 and SUT4 expressed abundantly at different conditions. SUT2 had similar expression level in all of the examined tissues, but SUT3 was undetectable. Both of SUT5 and SUT6 had lower expression level than other gene member, and expressed stronger in source leaves and are likely to play roles in phloem loading. In the seeding plant leave under water stress, four genes SUT1, SUT2, SUT4 and SUT5 were detectable. In these detectable genes, SUT1 and SUT4 were down regulated, while, SUT2 and SUT5 were up regulated.

CONCLUSIONS

In this study, we presented the first comprehensive genomic study for a whole gene family, the SUT family, in Saccharum. We speculated that there were six SUT members in the S. spotaneum genome. Out of the six members, SsSUTs, SsSUT5 and SsSUT6 were recent duplication genes accompanied by rapid evolution, while, SsSUT2 and SsSUT4 were the ancient members in the families. Despite the high polypoidy genome, functional redundancy may not exist among the SUTs allelic haplotypes supported by the evidence of strong purifying selection of the gene allele. SUT3 could be a low active member in the family because it is undetectable in our study, but it might not be a pseudogene because it harbored integrated gene structure. SUT1 and SUT4 were the main members for the sucrose transporter, while, these SUTs had sub-functional divergence in response to sucrose accumulation and plant development in Saccharum.

Transcriptome asymmetry in synthetic and natural allotetraploid wheats, revealed by RNA-sequencing

Allopolyploidization has occurred frequently within the Triticum-Aegilops complex which provides a suitable system to investigate how allopolyploidization shapes the expression patterns of duplicated homeologs. We have conducted transcriptome-profiling of leaves and young inflorescences in wild and domesticated tetraploid wheats, Triticum turgidum ssp. dicoccoides (BBAA) and ssp. durum (BBAA), an extracted tetraploid (BBAA), and a synthetic tetraploid (S(l) S(l) AA) wheat together with its diploid parents, Aegilops longissima (S(l) S(l) ) and Triticum urartu (AA). The two diploid species showed tissue-specific differences in genome-wide ortholog expression, which plays an important role in transcriptome shock-mediated homeolog expression rewiring and hence transcriptome asymmetry in the synthetic tetraploid. Further changes of homeolog expression apparently occurred in natural tetraploid wheats, which led to novel transcriptome asymmetry between the two subgenomes. In particular, our results showed that extremely biased homeolog expression can occur rapidly upon the allotetraploidzation and this trend is further enhanced in the course of domestication and evolution of polyploid wheats. Our results suggest that allopolyploidization is accompanied by distinct phases of homeolog expression changes, with parental legacy playing major roles in the immediate rewiring of homeolog expression upon allopolyploidization, while evolution and domestication under allotetraploidy drive further homeolog-expression changes toward re-established subgenome expression asymmetry.

Cytonuclear Coordination Is Not Immediate upon Allopolyploid Formation in Tragopogon miscellus (Asteraceae) Allopolyploids

Allopolyploids, formed by hybridization and chromosome doubling, face the immediate challenge of having duplicated nuclear genomes that interact with the haploid and maternally inherited cytoplasmic (plastid and mitochondrial) genomes. Most of our knowledge of the genomic consequences of allopolyploidy has focused on the fate of the duplicated nuclear genes without regard to their potential interactions with cytoplasmic genomes. As a step toward understanding the fates of nuclear-encoded subunits that are plastid-targeted, here we examine the retention and expression of the gene encoding the small subunit of Ribulose-1, 5-bisphosphate carboxylase/oxygenase (Rubisco; rbcS) in multiple populations of allotetraploid Tragopogon miscellus (Asteraceae). These polyploids formed recently (~80 years ago) and repeatedly from T. dubius and T. pratensis in the northwestern United States. Examination of 79 T. miscellus individuals from 10 natural populations, as well as 25 synthetic allotetraploids, including reciprocally formed plants, revealed a low percentage of naturally occurring individuals that show a bias in either gene (homeolog) loss (12%) or expression (16%), usually toward maintaining the maternal nuclear copy of rbcS. For individuals showing loss, seven retained the maternally derived rbcS homeolog only, while three had the paternally derived copy. All of the synthetic polyploid individuals examined (S0 and S1 generations) retained and expressed both parental homeologs. These results demonstrate that cytonuclear coordination does not happen immediately upon polyploid formation in Tragopogon miscellus.

250 years of hybridization between two biennial herb species without speciation

Hybridization between plant species can generate novel morphological diversity and lead to speciation at homoploid or polyploid levels. Hybrids between biennial herbs Tragopogon pratensis and T. porrifolius have been studied in experimental and natural populations for over 250 years. Here we examine their current status in natural populations in southeast England. All hybrids found were diploid; they tended to grow taller and with more buds than their parental species; many showed partial fertility; a few showed evidence of backcrossing. However, we found no evidence to suggest that the hybrids are establishing as a new species, nor can we find literature documenting speciation of these hybrids elsewhere. This lack of speciation despite at least 250 years of hybridization contrasts with the fact that both parental species have formed new allopolyploid species through hybridization with another diploid, T. dubius. Understanding why hybrids often do not speciate, despite repeated opportunities, would enhance our understanding of both the evolutionary process and risk assessments of invasive species.

Origin and Reticulate Evolutionary Process of Wheatgrass Elymus trachycaulus (Triticeae: Poaceae)

To study origin and evolutionary dynamics of tetraploid Elymus trachycaulus that has been cytologically defined as containing StH genomes, thirteen accessions of E. trachycaulus were analyzed using two low-copy nuclear gene Pepc (phosphoenolpyruvate carboxylase) and Rpb2 (the second largest subunit of RNA polymerase II), and one chloroplast region trnL–trnF (spacer between the tRNA Leu (UAA) gene and the tRNA-Phe (GAA) gene). Our chloroplast data indicated that Pseudoroegneria (St genome) was the maternal donor of E. trachycaulus. Rpb2 data indicated that the St genome in E. trachycaulus was originated from either P. strigosa, P. stipifolia, P. spicata or P. geniculate. The Hordeum (H genome)-like sequences of E. trachycaulus are polyphyletic in the Pepc tree, suggesting that the H genome in E. trachycaulus was contributed by multiple sources, whether due to multiple origins or introgression resulting from subsequent hybridization. Failure to recovering St copy of Pepc sequence in most accessions of E. trachycaulus might be caused by genome convergent evolution in allopolyploids. Multiple copies of H-like Pepc sequence from each accession with relative large deletions and insertions might be caused by either instability of Pepc sequence in H- genome or incomplete concerted evolution. Our results highlighted complex evolutionary history of E. trachycaulus.

Genetic and epigenetic changes in somatic hybrid introgression lines between wheat and tall wheatgrass

Broad phenotypic variations were induced in derivatives of an asymmetric somatic hybridization of bread wheat (Triticum aestivum) and tall wheatgrass (Thinopyrum ponticum Podp); however, how these variations occurred was unknown. We explored the nature of these variations by cytogenetic assays and DNA profiling techniques to characterize six genetically stable somatic introgression lines. Karyotyping results show the six lines similar to their wheat parent, but GISH analysis identified the presence of a number of short introgressed tall wheatgrass chromatin segments. DNA profiling revealed many genetic and epigenetic differences, including sequences deletions, altered regulation of gene expression, changed patterns of cytosine methylation, and the reactivation of retrotransposons. Phenotypic variations appear to result from altered repetitive sequences combined with the epigenetic regulation of gene expression and/or retrotransposon transposition. The extent of genetic and epigenetic variation due to the maintenance of parent wheat cells in tissue culture was assessed and shown to be considerably lower than had been induced in the introgression lines. Asymmetric somatic hybridization provides appropriate material to explore the nature of the genetic and epigenetic variations induced by genomic shock.

Patterns of chromosomal variation in natural populations of the neoallotetraploid Tragopogon mirus (Asteraceae)

Cytological studies have shown many newly formed allopolyploids (neoallopolyploids) exhibit chromosomal variation as a result of meiotic irregularities, but few naturally occurring neoallopolyploids have been examined. Little is known about how long chromosomal variation may persist and how it might influence the establishment and evolution of allopolyploids in nature. In this study we assess chromosomal composition in a natural neoallotetraploid, Tragopogon mirus, and compare it with T. miscellus, which is an allotetraploid of similar age (~40 generations old). We also assess whether parental gene losses in T. mirus correlate with entire or partial chromosome losses. Of 37 T. mirus individuals that were karyotyped, 23 (62%) were chromosomally additive of the parents, whereas the remaining 14 individuals (38%) had aneuploid compositions. The proportion of additive versus aneuploid individuals differed from that found previously in T. miscellus, in which aneuploidy was more common (69%; Fisher's exact test, P=0.0033). Deviations from parental chromosome additivity within T. mirus individuals also did not reach the levels observed in T. miscellus, but similar compensated changes were observed. The loss of T. dubius-derived genes in two T. mirus individuals did not correlate with any chromosomal changes, indicating a role for smaller-scale genetic alterations. Overall, these data for T. mirus provide a second example of prolonged chromosomal instability in natural neoallopolyploid populations.

The legacy of diploid progenitors in allopolyploid gene expression patterns

Allopolyploidization (hybridization and whole-genome duplication) is a common phenomenon in plant evolution with immediate saltational effects on genome structure and gene expression. New technologies have allowed rapid progress over the past decade in our understanding of the consequences of allopolyploidy. A major question, raised by early pioneer of this field Leslie Gottlieb, concerned the extent to which gene expression differences among duplicate genes present in an allopolyploid are a legacy of expression differences that were already present in the progenitor diploid species. Addressing this question necessitates phylogenetically well-understood natural study systems, appropriate technology, availability of genomic resources and a suitable analytical framework, including a sufficiently detailed and generally accepted terminology. Here, we review these requirements and illustrate their application to a natural study system that Gottlieb worked on and recommended for this purpose: recent allopolyploids of Tragopogon (Asteraceae). We reanalyse recent data from this system within the conceptual framework of parental legacies on duplicate gene expression in allopolyploids. On a broader level, we highlight the intellectual connection between Gottlieb's phrasing of this issue and the more contemporary framework of cis- versus trans-regulation of duplicate gene expression in allopolyploid plants.

Silenced rRNA genes are activated and substitute for partially eliminated active homeologs in the recently formed allotetraploid, Tragopogon mirus (Asteraceae)

To study the relationship between uniparental rDNA (encoding 18S, 5.8S and 26S ribosomal RNA) silencing (nucleolar dominance) and rRNA gene dosage, we studied a recently emerged (within the last 80 years) allotetraploid Tragopogon mirus (2n=24), formed from the diploid progenitors T. dubius (2n=12, D-genome donor) and T. porrifolius (2n=12, P-genome donor). Here, we used molecular, cytogenetic and genomic approaches to analyse rRNA gene activity in two sibling T. mirus plants (33A and 33B) with widely different rRNA gene dosages. Plant 33B had ~400 rRNA genes at the D-genome locus, which is typical for T. mirus, accounting for ~25% of total rDNA. We observed characteristic expression dominance of T. dubius-origin genes in all organs. Its sister plant 33A harboured a homozygous macrodeletion that reduced the number of T. dubius-origin genes to about 70 copies (~4% of total rDNA). It showed biparental rDNA expression in root, flower and callus, but not in leaf where D-genome rDNA dominance was maintained. There was upregulation of minor rDNA variants in some tissues. The RNA polymerase I promoters of reactivated T. porrifolius-origin rRNA genes showed reduced DNA methylation, mainly at symmetrical CG and CHG nucleotide motifs. We hypothesise that active, decondensed rDNA units are most likely to be deleted via recombination. The silenced homeologs could be used as a 'first reserve' to ameliorate mutational damage and contribute to evolutionary success of polyploids. Deletion and reactivation cycles may lead to bidirectional homogenisation of rRNA arrays in the long term.

A complex recombination pattern in the genome of allotetraploid Brassica napus as revealed by a high-density genetic map

Polyploidy plays a crucial role in plant evolution. Brassica napus (2n = 38, AACC), the most important oil crop in the Brassica genus, is an allotetraploid that originated through natural doubling of chromosomes after the hybridization of its progenitor species, B. rapa (2n = 20, AA) and B. oleracea (2n = 18, CC). A better understanding of the evolutionary relationship between B. napus and B. rapa, B. oleracea, as well as Arabidopsis, which has a common ancestor with these three species, will provide valuable information about the generation and evolution of allopolyploidy. Based on a high-density genetic map with single nucleotide polymorphism (SNP) and simple sequence repeat (SSR) markers, we performed a comparative genomic analysis of B. napus with Arabidopsis and its progenitor species B. rapa and B. oleracea. Based on the collinear relationship of B. rapa and B. oleracea in the B. napus genetic map, the B. napus genome was found to consist of 70.1% of the skeleton components of the chromosomes of B. rapa and B. oleracea, with 17.7% of sequences derived from reciprocal translocation between homoeologous chromosomes between the A- and C-genome and 3.6% of sequences derived from reciprocal translocation between non-homologous chromosomes at both intra- and inter-genomic levels. The current study thus provides insights into the formation and evolution of the allotetraploid B. napus genome, which will allow for more accurate transfer of genomic information from B. rapa, B. oleracea and Arabidopsis to B. napus.

Gene silencing via DNA methylation in naturally occurring Tragopogon miscellus (Asteraceae) allopolyploids

Background

Hybridization coupled with whole-genome duplication (allopolyploidy) leads to a variety of genetic and epigenetic modifications in the resultant merged genomes. In particular, gene loss and gene silencing are commonly observed post-polyploidization. Here, we investigated DNA methylation as a potential mechanism for gene silencing in Tragopogon miscellus (Asteraceae), a recent and recurrently formed allopolyploid. This species, which also exhibits extensive gene loss, was formed from the diploids T. dubius and T. pratensis.

Results

Comparative bisulfite sequencing revealed CG methylation of parental homeologs for three loci (S2, S18 and TDF-44) that were previously identified as silenced in T. miscellus individuals relative to the diploid progenitors. One other locus (S3) examined did not show methylation, indicating that other transcriptional and post-transcriptional mechanisms are likely responsible for silencing that homeologous locus.

Conclusions

These results indicate that Tragopogon miscellus allopolyploids employ diverse mechanisms, including DNA methylation, to respond to the potential shock of genome merger and doubling.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-701) contains supplementary material, which is available to authorized users.

Rapid genetic and epigenetic alterations under intergeneric genomic shock in newly synthesized Chrysanthemum morifolium x Leucanthemum paludosum hybrids (Asteraceae)

The Asteraceae family is at the forefront of the evolution due to frequent hybridization. Hybridization is associated with the induction of widespread genetic and epigenetic changes and has played an important role in the evolution of many plant taxa. We attempted the intergeneric cross Chrysanthemum morifolium × Leucanthemum paludosum. To obtain the success in cross, we have to turn to ovule rescue. DNA profiling of the amphihaploid and amphidiploid was investigated using amplified fragment length polymorphism, sequence-related amplified polymorphism, start codon targeted polymorphism, and methylation-sensitive amplification polymorphism (MSAP). Hybridization induced rapid changes at the genetic and the epigenetic levels. The genetic changes mainly involved loss of parental fragments and gaining of novel fragments, and some eliminated sequences possibly from the noncoding region of L. paludosum. The MSAP analysis indicated that the level of DNA methylation was lower in the amphiploid (∼45%) than in the parental lines (51.5-50.6%), whereas it increased after amphidiploid formation. Events associated with intergeneric genomic shock were a feature of C. morifolium × L. paludosum hybrid, given that the genetic relationship between the parental species is relatively distant. Our results provide genetic and epigenetic evidence for understanding genomic shock in wide crosses between species in Asteraceae and suggest a need to expand our current evolutionary framework to encompass a genetic/epigenetic dimension when seeking to understand wide crosses.

Polyploidy and novelty: Gottlieb's legacy

Nearly four decades ago, Roose & Gottlieb (Roose & Gottlieb 1976 Evolution 30, 818-830. (doi:10.2307/2407821)) showed that the recently derived allotetraploids Tragopogon mirus and T. miscellus combined the allozyme profiles of their diploid parents (T. dubius and T. porrifolius, and T. dubius and T. pratensis, respectively). This classic paper addressed the link between genotype and biochemical phenotype and documented enzyme additivity in allopolyploids. Perhaps more important than their model of additivity, however, was their demonstration of novelty at the biochemical level. Enzyme multiplicity-the production of novel enzyme forms in the allopolyploids-can provide an extensive array of polymorphism for a polyploid individual and may explain, for example, the expanded ranges of polyploids relative to their diploid progenitors. In this paper, we extend the concept of evolutionary novelty in allopolyploids to a range of genetic and ecological features. We observe that the dynamic nature of polyploid genomes-with alterations in gene content, gene number, gene arrangement, gene expression and transposon activity-may generate sufficient novelty that every individual in a polyploid population or species may be unique. Whereas certain combinations of these features will undoubtedly be maladaptive, some unique combinations of newly generated variation may provide tremendous evolutionary potential and adaptive capabilities.

The place of asymmetric somatic hybridization in wheat breeding

Since its first development some 40 years ago, the application of the somatic hybridization technique has generated a body of hybrid plant material involving a wide combination of parental species. Until the late 1990s, the technique was ineffective in wheat, as regeneration from protoplasts was proving difficult to achieve. Since this time, however, a successful somatic hybridization protocol for wheat has been established and used to generate a substantial number of both symmetric and asymmetric somatic hybrids and derived materials, especially involving the parental combination bread wheat and tall wheatgrass (Thinopyrum ponticum). This review describes the current state of the art for somatic hybridization in wheat and focuses on its potential application for wheat improvement.