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

Dosage effect genes modulate grain development in synthesized Triticum durum-Haynaldia villosa allohexaploid

Polyploidization in plants often leads to increased cell size and grain size, which may be affected by the increased genome dosage and transcription abundance. The synthesized Triticum durum (AABB)-Haynaldia villosa (VV) amphiploid (AABBVV) has significantly increased grain size, especially grain length, than the tetraploid and diploid parents. To investigate how polyploidization affects grain development at the transcriptional level, we perform transcriptome analysis using the immature seeds of T. durum, H. villosa, and the amphiploid. The dosage effect genes are contributed more by differentially expressed genes from genome V of H. villosa. The dosage effect genes overrepresent grain development-related genes. Interestingly, the vernalization gene TaVRN1 is among the positive dosage effect genes in the T. durum‒H. villosa and T. turgidum‒Ae. tauschii amphiploids. The expression levels of TaVRN1 homologs are positively correlated with the grain size and weight. The TaVRN1-B1 or TaVRN1-D1 mutation shows delayed florescence, decreased cell size, grain size, and grain yield. These data indicate that dosage effect genes could be one of the important explanations for increased grain size by regulating grain development. The identification and functional validation of dosage effect genes may facilitate the finding of valuable genes for improving wheat yield.

Assessing small RNA profiles in potato diploid hybrid and its resynthesized allopolyploid reveals conserved abundance with distinct genomic distribution

KEY MESSAGE

The shock produced by the allopolyploidization process on a potato interspecific diploid hybrid displays a non-random remobilization of the small RNAs profile on a variety of genomic features. Allopolyploidy, a complex process involving interspecific hybridization and whole genome duplication, significantly impacts plant evolution, leading to the emergence of novel phenotypes. Polyploids often present phenotypic nuances that enhance adaptability, enabling them to compete better and occasionally to colonize new habitats. Whole-genome duplication represents a genomic "shock" that can trigger genetic and epigenetic changes that yield novel expression patterns. In this work, we investigate the polyploidization effect on a diploid interspecific hybrid obtained through the cross between the cultivated potato Solanum tuberosum and the wild potato Solanum kurtzianum, by assessing the small RNAs (sRNAs) profile of the parental diploid hybrid and its derived allopolyploid. Small RNAs are key components of the epigenetic mechanisms involved in silencing by RNA-directed DNA Methylation (RdDM). A sRNA sequencing (sRNA-Seq) analysis was performed to individually profile the 21 to 22 nucleotide (21 to 22-nt) and 24-nt sRNA size classes due to their unique mechanism of biogenesis and mode of function. The composition and distribution of different genomic features and differentially accumulated (DA) sRNAs were evaluated throughout the potato genome. We selected a subset of genes associated with DA sRNAs for messenger RNA (mRNA) expression analysis to assess potential impacts on the transcriptome. Interestingly, we noted that 24-nt DA sRNAs that exclusively mapped to exons were correlated with differentially expressed mRNAs between genotypes, while this behavior was not observed when 24-nt DA sRNAs were mapped to intronic regions. These findings collectively emphasize the nonstochastic nature of sRNA remobilization in response to the genomic shock induced by allopolyploidization.

Modeling transposable elements dynamics during polyploidization in plants

In this work we study the proliferation of transposable elements (TEs) and the epigenetic response of plants during the process of polyploidization. Through a deterministic model, expanding on our previous work on TE proliferation under epigenetic regulation, we study the long-term TE distribution and TE stability in the subgenomes of both autopolyploids and allopolyploids. We also explore different small-interfering RNA (siRNA) action modes on the subgenomes, including a model where siRNAs are not directed to specific genomes and one where siRNAs are directed - i.e. more active - in subgenomes with higher TE loads. In the autopolyploid case, we find long-term stable equilbria that tend to equilibrate the number of active TEs between subgenomes. In the allopolyploid case, directed siRNA action is fundamental to avoid a "winner takes all" outcome of the competition between the TEs in the different subgenomes. We also show that decaying oscillations in the number of TEs occur naturally in all cases, perhaps explaining some of the observed features of 'genomic shock' after hybridization events, and that the balance in the dynamics of the different types of siRNA is determinant for the synchronization of these oscillations.

The allotetraploid horseradish genome provides insights into subgenome diversification and formation of critical traits

Polyploidization can provide a wealth of genetic variation for adaptive evolution and speciation, but understanding the mechanisms of subgenome evolution as well as its dynamics and ultimate consequences remains elusive. Here, we report the telomere-to-telomere (T2T) gap-free reference genome of allotetraploid horseradish (Armoracia rusticana) sequenced using a comprehensive strategy. The (epi)genomic architecture and 3D chromatin structure of the A and B subgenomes differ significantly, suggesting that both the dynamics of the dominant long terminal repeat retrotransposons and DNA methylation have played critical roles in subgenome diversification. Investigation of the genetic basis of biosynthesis of glucosinolates (GSLs) and horseradish peroxidases reveals both the important role of polyploidization and subgenome differentiation in shaping the key traits. Continuous duplication and divergence of essential genes of GSL biosynthesis (e.g., FMO_GS-OX, IGMT, and GH1 gene family) contribute to the broad GSL profile in horseradish. Overall, the T2T assembly of the allotetraploid horseradish genome expands our understanding of polyploid genome evolution and provides a fundamental genetic resource for breeding and genetic improvement of horseradish.

Epigenetic Regulation of Subgenomic Gene Expression in Allotetraploid Brassica napus

The allotetraploid Brasscia napus has now been extensively utilized to reveal the genetic processes involved in hybridization and polyploidization. Here, transcriptome, WGBS, and Chip-Seq sequencing data were obtained to explore the regulatory consequences of DNA methylation and histone modifications on gene expression in B. napus. When compared with diploid parents, the expression levels of 14,266 (about 32%) and 17,054 (about 30%) genes were altered in the A_n and C_n subgenomes, respectively, and a total of 4982 DEGs were identified in B. napus. Genes with high or no expression in diploid parents often shifted to medium or low expression in B. napus. The number of genes with elevated methylation levels in gene promoters and gene body regions has increased in A_n and C_n subgenomes. The peak number of H3K4me3 modification increased, while the peak number of H3K27ac and H3K27me3 decreased in A_n and C_n subgenomes, and more genes that maintained parental histone modifications were identified in C_n subgenome. The differential multiples of DEGs in B. napus were positively correlated with DNA methylation levels in promoters and the gene body, and the differential multiples of these DEGs were also affected by the degree of variation in DNA methylation levels. Further analysis revealed that about 99% of DEGs were of DNA methylation, and about 68% of DEGs were modified by at least two types of DNA methylation and H3K4me3, H3K27ac, and H3K27me3 histone modifications. These results demonstrate that DNA methylation is crucial for gene expression regulation, and different epigenetic modifications have an essential function in regulating the differential expression of genes in B. napus.

Epigenetic Landscape Is Largely Shaped by Diversiform Transposons in Aegilops tauschii

Transposons (TEs) account for more than 80% of the wheat genome, the highest among all known crop species. They play an important role in shaping the elaborate genomic landscape, which is the key to the speciation of wheat. In this study, we analyzed the association between TEs, chromatin states, and chromatin accessibility in Aegilops tauschii, the D genome donor of bread wheat. We found that TEs contributed to the complex but orderly epigenetic landscape as chromatin states showed diverse distributions on TEs of different orders or superfamilies. TEs also contributed to the chromatin state and openness of potential regulatory elements, affecting the expression of TE-related genes. Some TE superfamilies, such as hAT-Ac, carry active/open chromatin regions. In addition, the histone mark H3K9ac was found to be associated with the accessibility shaped by TEs. These results suggest the role of diversiform TEs in shaping the epigenetic landscape and in gene expression regulation in Aegilops tauschii. This has positive implications for understanding the transposon roles in Aegilops tauschii or the wheat D genome.

Low impact of polyploidization on the transcriptome of synthetic allohexaploid wheat

BACKGROUND

Bread wheat is a recent allohexaploid (genomic constitution AABBDD) that emerged through a hybridization between tetraploid Triticum turgidum (AABB) and diploid Aegilops tauschii (DD) less than 10,000 years ago. The hexaploidization can be re-created artificially, producing synthetic wheat that has been used to study immediate genomic responses to polyploidization. The scale of the consequences of polyploidization, and their mechanism of establishment, remain uncertain.

RESULTS

Here we sampled several synthetic wheats from alternative parental genotypes and reciprocal crosses, and examined transcriptomes from two different tissues and successive generations. We did not detect any massive reprogramming in gene expression, with only around 1% of expressed genes showing significant differences compared to their lower-ploidy parents. Most of this differential expression is located on the D subgenome, without consistency in the direction of the expression change. Homoeolog expression bias in synthetic wheat is similar to the pattern observed in the parents. Both differential expression and homoeolog bias are tissue-specific. While up to three families of transposable elements became upregulated in wheat synthetics, their position and distance are not significantly associated with expression changes in proximal genes.

DISCUSSION

While only a few genes change their expression pattern after polyploidization, they can be involved in agronomically important pathways. Alternative parental combinations can lead to opposite changes on the same subset of D-located genes, which is relevant for harnessing new diversity in wheat breeding. Tissue specificity of the polyploidization-triggered expression changes indicates the remodelling of transcriptomes in synthetic wheat is plastic and likely caused by regulome interactions rather than permanent changes. We discuss the pitfalls of transcriptomic comparisons across ploidy levels that can inflate the de-regulation signal.

CONCLUSIONS

Transcriptomic response to polyploidization in synthetic AABBDD wheat is modest and much lower than some previous estimates. Homoeolog expression bias in wheat allohexaploids is mostly attributed to parental legacy, with polyploidy having a mild balancing effect.

Global transcriptome analysis of allopolyploidization reveals large-scale repression of the D-subgenome in synthetic hexaploid wheat

Synthetic hexaploid wheat (SHW) lines are created as pre-breeding germplasm to diversify the D subgenome of hexaploid wheat and capitalize upon the untapped genetic diversity of the Aegilops tauschii gene pool. However, the phenotypes observed in the Ae. tauschii parents are not always recovered in the SHW lines, possibly due to inter-subgenome interactions. To elucidate this post-polyploidization genome reprogramming phenomenon, we performed RNA-seq of four SHW lines and their corresponding tetraploid and diploid parents, across ten tissues and three biological replicates. Homoeologue expression bias (HEB) analysis using more than 18,000 triads suggests massive suppression of homoeoalleles of the D subgenome in SHWs. Comparative transcriptome analysis of the whole-genome gene set further corroborated this finding. Alternative splicing analysis of the high-confidence genes indicates an additional layer of complexity where all five splice events are identified, and retained intron is predominant. Homoeologue expression upon resynthesis of hexaploid wheat has implications to the usage and handling of this germplasm in breeding as it relates to capturing the effects of epistatic interaction across subgenomes upon polyploidization. Special considerations must be given to this germplasm in pre-breeding activities to consider the extent of the inter-subgenome interactions on gene expression and their impact on traits for crop improvement.

Does one subgenome become dominant in the formation and evolution of a polyploid?

BACKGROUND

Polyploids are common in flowering plants and they tend to have more expanded ranges of distributions than their diploid progenitors. Possible mechanisms underlying polyploid success have been intensively investigated. Previous studies showed that polyploidy generates novel changes and that subgenomes in allopolyploid species often differ in gene number, gene expression levels and levels of epigenetic alteration. It is widely believed that such differences are the results of conflicts among the subgenomes. These differences have been treated by some as subgenome dominance, and it is claimed that the magnitude of subgenome dominance increases in polyploid evolution.

SCOPE

In addition to changes which occurred during evolution, differences between subgenomes of a polyploid species may also be affected by differences between the diploid donors and changes which occurred during polyploidization. The variable genome components in many plant species are extensive, which would result in exaggerated differences between a subgenome and its progenitor when a single genotype or a small number of genotypes are used to represent a polyploid or its donors. When artificially resynthesized polyploids are used as surrogates for newly formed genotypes which have not been exposed to evolutionary selection, differences between diploid genotypes available today and those involved in the formation of the natural polyploid genotypes must also be considered.

CONCLUSIONS

Contrary to the now widely held views that subgenome biases in polyploids are the results of conflicts among the subgenomes and that one of the parental subgenomes generally retains more genes which are more highly expressed, available results show that subgenome biases mainly reflect legacy from the progenitors and that they can be detected before the completion of polyploidization events. Further, there is no convincing evidence that the magnitudes of subgenome biases have significantly changed during evolution for any of the allopolyploid species assessed.

Molecular prospective on the wheat grain development

Wheat grain development is an important biological process to determine grain yield and quality, which is controlled by the interplay of genetic, epigenetic, and environmental factors. Wheat grain development has been extensively characterized at the phenotypic and genetic levels. The advent of innovative molecular technologies allows us to characterize genes, proteins, and regulatory factors involved in wheat grain development, which have enhanced our understanding of the wheat seed development process. However, wheat is an allohexaploid with a large genome size, the molecular mechanisms underlying the wheat grain development have not been well understood as those in diploids. Understanding grain development, and how it is regulated, is of fundamental importance for improving grain yield and quality through conventional breeding or genetic engineering. Herein, we review the current discoveries on the molecular mechanisms underlying wheat grain development. Notably, only a handful of genes that control wheat grain development have, thus far, been well characterized, their interplay underlying the grain development remains elusive. The synergistic network-integrated genomics and epigenetics underlying wheat grain development and how the subgenome divergence dynamically and precisely regulates wheat grain development are unknown.

Recent allopolyploidy alters Spartina microRNA expression in response to xenobiotic-induced stress

Environmental contamination by xenobiotics represents a major threat for natural ecosystems and public health. In response, xenobiotic detoxification is a fundamental trait of organisms for developmental plasticity and stress tolerance, but the underlying molecular mechanisms remain poorly understood in plants. To decipher this process, we explored the consequences of allopolyploidy on xenobiotic tolerance in the genus Spartina Schreb. Specifically, we focused on microRNAs (miRNAs) owing to their central function in the regulation of gene expression patterns, including responses to stress. Small RNA-Seq was conducted on the parents S. alterniflora and S. maritima, their F1 hybrid S. x townsendii and the allopolyploid S. anglica under phenanthrene-induced stress (phe), a model Polycyclic Aromatic Hydrocarbon (PAH) compound. Differentially expressed miRNAs in response to phe were specifically identified within species. In complement, the respective impacts of hybridization and genome doubling were detected, through changes in miRNA expression patterns between S. x townsendii, S. anglica and the parents. The results support the impact of allopolyploidy in miRNA-guided regulation of plant response to phe. In total, we identified 17 phe-responsive miRNAs in Spartina among up-regulated MIR156 and down-regulated MIR159. We also describe novel phe-responsive miRNAs as putative Spartina-specific gene expression regulators in response to stress. Functional validation using Arabidopsis (L.) Heynh. T-DNA lines inserted in homologous MIR genes was performed, and the divergence of phe-responsive miRNA regulatory networks between Arabidopsis and Spartina was discussed.

Low-affinity SPL binding sites contribute to subgenome expression divergence in allohexaploid wheat

Expression divergence caused by genetic variation and crosstalks among subgenomes of the allohexaploid bread wheat (Triticum aestivum. L., BBAADD) is hypothesized to increase its adaptability and/or plasticity. However, the molecular basis of expression divergence remains unclear. Squamosa promoter-binding protein-like (SPL) transcription factors are critical for a wide array of biological processes. In this study, we constructed expression regulatory networks by combining DAP-seq for 40 SPLs, ATAC-seq, and RNA-seq. Our findings indicate that a group of low-affinity SPL binding regions (SBRs) were targeted by diverse SPLs and caused different sequence preferences around the core GTAC motif. The SBRs including the low-affinity ones are evolutionarily conserved, enriched GWAS signals related to important agricultural traits. However, those SBRs are highly diversified among the cis-regulatory regions (CREs) of syntenic genes, with less than 8% SBRs coexisting in triad genes, suggesting that CRE variations are critical for subgenome differentiations. Knocking out of TaSPL7A/B/D and TaSPL15A/B/D subfamily further proved that both high- and low-affinity SBRs played critical roles in the differential expression of genes regulating tiller number and spike sizes. Our results have provided baseline data for downstream networks of SPLs and wheat improvements and revealed that CRE variations are critical sources for subgenome divergence in the allohexaploid wheat.

Deleterious Mutations Accumulate Faster in Allopolyploid than Diploid Cotton (Gossypium) and Unequally between Subgenomes

Whole genome duplication (polyploidization) is among the most dramatic mutational processes in nature, so understanding how natural selection differs in polyploids relative to diploids is an important goal. Population genetics theory predicts that recessive deleterious mutations accumulate faster in allopolyploids than diploids due to the masking effect of redundant gene copies, but this prediction is hitherto unconfirmed. Here, we use the cotton genus (Gossypium), which contains seven allopolyploids derived from a single polyploidization event 1-2 million years ago, to investigate deleterious mutation accumulation. We use two methods of identifying deleterious mutations at the nucleotide and amino acid level, along with whole-genome resequencing of 43 individuals spanning six allopolyploid species and their two diploid progenitors, to demonstrate that deleterious mutations accumulate faster in allopolyploids than in their diploid progenitors. We find that, unlike what would be expected under models of demographic changes alone, strongly deleterious mutations show the biggest difference between ploidy levels, and this effect diminishes for moderately and mildly deleterious mutations. We further show that the proportion of nonsynonymous mutations that are deleterious differs between the two co-resident subgenomes in the allopolyploids, suggesting that homoeologous masking acts unequally between subgenomes. Our results provide a genome-wide perspective on classic notions of the significance of gene duplication that likely are broadly applicable to allopolyploids, with implications for our understanding of the evolutionary fate of deleterious mutations. Finally, we note that some measures of selection (e.g. dN/dS, πN/πS) may be biased when species of different ploidy levels are compared.

Asymmetric Divergence in Transmitted SNPs of DNA Replication/Transcription and Their Impact on Gene Expression in Polyploid Brassica napus

The marked increase in plant genomic data has provided valuable resources for investigating the dynamic evolution of duplicate genes in polyploidy. Brassica napus is an ideal model species for investigating polyploid genome evolution. The present study comprehensively analyzed DNA and RNA variation of two representative B. napus inbredlines, Zhongshuang11 and Zhongyou821, and we investigated gene expression levels of A_n and C_n subgenomes in multiple tissues of the two lines. The distribution of transmitted single nucleotide polymorphisms (SNPs) was significantly different in two subgenomes of B. napus. Gene expression levels were significantly negatively correlated with number of variations in replication and transcription of the corresponding genes, but were positively correlated with the ratios of transmitted SNPs from DNA to RNA. We found a higher density of SNP variation in A_n than that in C_n during DNA replication and more SNPs were transmitted to RNA during transcription, which may contribute to A_n expression dominance. These activities resulted in asymmetrical gene expression in polyploid B. napus. The SNPs transmitted from DNA to RNA could be an important complement feature in comparative genomics, and they may play important roles in asymmetrical genome evolution in polyploidy.

Histone acetyltransferase TaHAG1 acts as a crucial regulator to strengthen salt tolerance of hexaploid wheat

Polyploidy occurs prevalently and plays an important role during plant speciation and evolution. This phenomenon suggests polyploidy could develop novel features that enable them to adapt wider range of environmental conditions compared with diploid progenitors. Bread wheat (Triticum aestivum L., BBAADD) is a typical allohexaploid species and generally exhibits greater salt tolerance than its tetraploid wheat progenitor (BBAA). However, little is known about the underlying molecular basis and the regulatory pathway of this trait. Here, we show that the histone acetyltransferase TaHAG1 acts as a crucial regulator to strengthen salt tolerance of hexaploid wheat. Salinity-induced TaHAG1 expression was associated with tolerance variation in polyploidy wheat. Overexpression, silencing, and CRISPR-mediated knockout of TaHAG1 validated the role of TaHAG1 in salinity tolerance of wheat. TaHAG1 contributed to salt tolerance by modulating reactive oxygen species (ROS) production and signal specificity. Moreover, TaHAG1 directly targeted a subset of genes that are responsible for hydrogen peroxide production, and enrichment of TaHAG1 triggered increased H3 acetylation and transcriptional upregulation of these loci under salt stress. In addition, we found the salinity-induced TaHAG1-mediated ROS production pathway is involved in salt tolerance difference of wheat accessions with varying ploidy. Our findings provide insight into the molecular mechanism of how an epigenetic regulatory factor facilitates adaptability of polyploidy wheat and highlights this epigenetic modulator as a strategy for salt tolerance breeding in bread wheat.

An atlas of wheat epigenetic regulatory elements reveals subgenome divergence in the regulation of development and stress responses

Wheat (Triticum aestivum) has a large allohexaploid genome. Subgenome-divergent regulation contributed to genome plasticity and the domestication of polyploid wheat. However, the specificity encoded in the wheat genome determining subgenome-divergent spatio-temporal regulation has been largely unexplored. The considerable size and complexity of the genome are major obstacles to dissecting the regulatory specificity. Here, we compared the epigenomes and transcriptomes from a large set of samples under diverse developmental and environmental conditions. Thousands of distal epigenetic regulatory elements (distal-epiREs) were specifically linked to their target promoters with coordinated epigenomic changes. We revealed that subgenome-divergent activity of homologous regulatory elements is affected by specific epigenetic signatures. Subgenome-divergent epiRE regulation of tissue specificity is associated with dynamic modulation of H3K27me3 mediated by Polycomb complex and demethylases. Furthermore, quantitative epigenomic approaches detected key stress responsive cis- and trans-acting factors validated by DNA Affinity Purification and sequencing, and demonstrated the coordinated interplay between epiRE sequence contexts, epigenetic factors, and transcription factors in regulating subgenome divergent transcriptional responses to external changes. Together, this study provides a wealth of resources for elucidating the epiRE regulomics and subgenome-divergent regulation in hexaploid wheat, and gives new clues for interpreting genetic and epigenetic interplay in regulating the benefits of polyploid wheat.

Polyploidy: an evolutionary and ecological force in stressful times

Polyploidy has been hypothesized to be both an evolutionary dead-end and a source for evolutionary innovation and species diversification. Although polyploid organisms, especially plants, abound, the apparent nonrandom long-term establishment of genome duplications suggests a link with environmental conditions. Whole-genome duplications seem to correlate with periods of extinction or global change, while polyploids often thrive in harsh or disturbed environments. Evidence is also accumulating that biotic interactions, for instance, with pathogens or mutualists, affect polyploids differently than nonpolyploids. Here, we review recent findings and insights on the effect of both abiotic and biotic stress on polyploids versus nonpolyploids and propose that stress response in general is an important and even determining factor in the establishment and success of polyploidy.

New insights into homoeologous copy number variations in the hexaploid wheat genome

Bread wheat is an allohexaploid species originating from two successive and recent rounds of hybridization between three diploid species that were very similar in terms of chromosome number, genome size, TE content, gene content and synteny. As a result, it has long been considered that most of the genes were in three pairs of homoeologous copies. However, these so-called triads represent only one half of wheat genes, while the remaining half belong to homoeologous groups with various number of copies across subgenomes. In this study, we examined and compared the distribution, conservation, function, expression and epigenetic profiles of triads with homoeologous groups having undergone a deletion (dyads) or a duplication (tetrads) in one subgenome. We show that dyads and tetrads are mostly located in distal regions and have lower expression level and breadth than triads. Moreover, they are enriched in functions related to adaptation and more associated with the repressive H3K27me3 modification. Altogether, these results suggest that triads mainly correspond to housekeeping genes and are part of the core genome, while dyads and tetrads belong to the Triticeae dispensable genome. In addition, by comparing the different categories of dyads and tetrads, we hypothesize that, unlike most of the allopolyploid species, subgenome dominance and biased fractionation are absent in hexaploid wheat. Differences observed between the three subgenomes are more likely related to two successive and ongoing waves of post-polyploid diploidization, that had impacted A and B more significantly than D, as a result of the evolutionary history of hexaploid wheat.

Dynamic and reversible DNA methylation changes induced by genome separation and merger of polyploid wheat

Background

Wheat is a powerful genetic model for studying polyploid evolution and crop domestication. Hexaploid bread wheat was formed by two rounds of interspecific hybridization and polyploidization, processes which are often accompanied by genetic and epigenetic changes, including DNA methylation. However, the extent and effect of such changes during wheat evolution, particularly from tetraploid-to-hexaploid wheat, are currently elusive.

Results

Here we report genome-wide DNA methylation landscapes in extracted tetraploid wheat (ETW, AABB), natural hexaploid wheat (NHW, AABBDD), resynthesized hexaploid wheat (RHW, AABBDD), natural tetraploid wheat (NTW, AABB), and diploid (DD). In the endosperm, levels of DNA methylation, especially in CHG (H=A, T, or C) context, were dramatically decreased in the ETW relative to natural hexaploid wheat; hypo-differentially methylated regions (DMRs) (850,832) were 24-fold more than hyper-DMRs (35,111). Interestingly, those demethylated regions in ETW were remethylated in the resynthesized hexaploid wheat after the addition of the D genome. In ETW, hypo-DMRs correlated with gene expression, and TEs were demethylated and activated, which could be silenced in the hexaploid wheat. In NHW, groups of TEs were dispersed in genic regions of three subgenomes, which may regulate the expression of TE-associated genes. Further, hypo-DMRs in ETW were associated with reduced H3K9me2 levels and increased expression of histone variant genes, suggesting concerted epigenetic changes after separation from the hexaploid.

Conclusion

Genome merger and separation provoke dynamic and reversible changes in chromatin and DNA methylation. These changes correlate with altered gene expression and TE activity, which may provide insights into polyploid genome and wheat evolution.

Using asexual vertebrates to study genome evolution and animal physiology: Banded (Fundulus diaphanus) x Common Killifish (F. heteroclitus) hybrid lineages as a model system

Wild, asexual, vertebrate hybrids have many characteristics that make them good model systems for studying how genomes evolve and epigenetic modifications influence animal physiology. In particular, the formation of asexual hybrid lineages is a form of reproductive incompatibility, but we know little about the genetic and genomic mechanisms by which this mode of reproductive isolation proceeds in animals. Asexual lineages also provide researchers with the ability to produce genetically identical individuals, enabling the study of autonomous epigenetic modifications without the confounds of genetic variation. Here, we briefly review the cellular and molecular mechanisms leading to asexual reproduction in vertebrates and the known genetic and epigenetic consequences of the loss of sex. We then specifically discuss what is known about asexual lineages of Fundulus diaphanus x F. heteroclitus to highlight gaps in our knowledge of the biology of these clones. Our preliminary studies of F. diaphanus and F. heteroclitus karyotypes from Porter's Lake (Nova Scotia, Canada) agree with data from other populations, suggesting a conserved interspecific chromosomal arrangement. In addition, genetic analyses suggest that: (a) the same major clonal lineage (Clone A) of F. diaphanus x F. heteroclitus has remained dominant over the past decade, (b) some minor clones have also persisted, (c) new clones may have recently formed, and iv) wild clones still mainly descend from F. diaphanus ♀ x F. heteroclitus ♂ crosses (96% in 2017-2018). These data suggest that clone formation may be a relatively rare, but continuous process, and there are persistent environmental or genetic factors causing a bias in cross direction. We end by describing our current research on the genomic causes and consequences of a transition to asexuality and the potential physiological consequences of epigenetic variation.

Reduced chromatin accessibility underlies gene expression differences in homologous chromosome arms of diploid Aegilops tauschii and hexaploid wheat

Background

Polyploidy is centrally important in the evolution and domestication of plants because it leads to major genomic changes, such as altered patterns of gene expression, which are thought to underlie the emergence of new traits. Despite the common occurrence of these globally altered patterns of gene expression in polyploids, the mechanisms involved are not well understood.

Results

Using a precisely defined framework of highly conserved syntenic genes on hexaploid wheat chromosome 3DL and its progenitor 3 L chromosome arm of diploid Aegilops tauschii, we show that 70% of these gene pairs exhibited proportionately reduced gene expression, in which expression in the hexaploid context of the 3DL genes was ∼40% of the levels observed in diploid Ae tauschii. Several genes showed elevated expression during the later stages of grain development in wheat compared with Ae tauschii. Gene sequence and methylation differences probably accounted for only a few cases of differences in gene expression. In contrast, chromosome-wide patterns of reduced chromatin accessibility of genes in the hexaploid chromosome arm compared with its diploid progenitor were correlated with both reduced gene expression and the imposition of new patterns of gene expression.

Conclusions

Our pilot-scale analyses show that chromatin compaction may orchestrate reduced gene expression levels in the hexaploid chromosome arm of wheat compared to its diploid progenitor chromosome arm.

One subgenome to rule them all: underlying mechanisms of subgenome dominance

Allopolyploids, which are formed from the hybridization of two or more diploid progenitor species, often experience subgenome dominance, where one of the parental genomes (subgenomes) has higher levels of gene expression and ultimately greater gene retention compared to the other subgenomes. Low transposable element (TE) abundance near genes has been associated with the dominant subgenome in several allopolyploids, but TEs are unlikely to be the only causal factor responsible for subgenome expression dominance. In this review, we will examine the role of TEs in subgenome dominance as well as discuss how genetic incompatibilities among subgenomes likely contributes to the rapid emergence of subgenome dominance. Lastly, we highlight several burning questions about subgenome dominance that remain largely unanswered.

Evolution of small RNA expression following hybridization and allopolyploidization: insights from Spartina species (Poaceae, Chloridoideae)

Differential expression of mi-RNAs targeting developmental processes and progressive downregulation of repeat-associated siRNAs following genome merger and genome duplication in the context of allopolyploid speciation in Spartina. The role of small RNAs on gene expression regulation and genome stability is arousing increased interest and is being explored in various plant systems. In spite of prominence of reticulate evolution and polyploidy that affects the evolutionary history of all plant lineages, very few studies analysed RNAi mechanisms with this respect. Here, we explored small RNAs diversity and expression in the context of recent allopolyploid speciation, using the Spartina system, which offers a unique opportunity to explore the immediate changes following hybridization and genome duplication. Small RNA-Seq analyses were conducted on hexaploid parental species (S. alterniflora and S. maritima), their F1 hybrid S. x townsendii, and the neoallododecaploid S. anglica. We identified 594 miRNAs, 2197 miRNA-target genes, and 3730 repeat-associated siRNAs (mostly targeting Class I/Copia-Ivana- Copia-SIRE and LINEs elements). For both mi- and ra-siRNAs, we detected differential expression patterns following genome merger and genome duplication. These misregulations include non-additive expression of miRNAs in the F1 hybrid and additional changes in the allopolyploid targeting developmental processes. Expression of repeat-associated siRNAs indicates a strengthen of transposable element repression during the allopolyploidization process. Altogether, these results confirm the central role small RNAs play in shaping regulatory changes in naturally formed recent allopolyploids.

The allotetraploid origin and asymmetrical genome evolution of the common carp Cyprinus carpio

Common carp (Cyprinus carpio) is an allotetraploid species derived from recent whole genome duplication and provides a model to study polyploid genome evolution in vertebrates. Here, we generate three chromosome-level reference genomes of C. carpio and compare to related diploid Cyprinid genomes. We identify a Barbinae lineage as potential diploid progenitor of C. carpio and then divide the allotetraploid genome into two subgenomes marked by a distinct genome similarity to the diploid progenitor. We estimate that the two diploid progenitors diverged around 23 Mya and merged around 12.4 Mya based on the divergence rates of homoeologous genes and transposable elements in two subgenomes. No extensive gene losses are observed in either subgenome. Instead, we find gene expression bias across surveyed tissues such that subgenome B is more dominant in homoeologous expression. CG methylation in promoter regions may play an important role in altering gene expression in allotetraploid C. carpio.

The common carp is derived from recent whole genome duplication and represents a model for polyploid genome evolution, rare in vertebrates. Here, the authors generate and analyse chromosome-level reference genomes for common carp, and describe subgenome gene expression changes.

Evolutionary and ecological functional genomics, from lab to the wild

Plant phenotypes are the result of both genetic and environmental forces that act to modulate trait expression. Over the last few years, numerous approaches in functional genomics and systems biology have led to a greater understanding of plant phenotypic variation and plant responses to the environment. These approaches, and the questions that they can address, have been loosely termed evolutionary and ecological functional genomics (EEFG), and have been providing key insights on how plants adapt and evolve. In particular, by bringing these studies from the laboratory to the field, EEFG studies allow us to gain greater knowledge of how plants function in their natural contexts.

Integrating Networks, Phylogenomics, and Population Genomics for the Study of Polyploidy

Duplication events are regarded as sources of evolutionary novelty, but our understanding of general trends for the long-term trajectory of additional genomic material is still lacking. Organisms with a history of whole genome duplication (WGD) offer a unique opportunity to study potential trends in the context of gene retention and/or loss, gene and network dosage, and changes in gene expression. In this review, we discuss the prevalence of polyploidy across the tree of life, followed by an overview of studies investigating genome evolution and gene expression. We then provide an overview of methods in network biology, phylogenomics, and population genomics that are critical for advancing our understanding of evolution post-WGD, highlighting the need for models that can accommodate polyploids. Finally, we close with a brief note on the importance of random processes in the evolution of polyploids with respect to neutral versus selective forces, ancestral polymorphisms, and the formation of autopolyploids versus allopolyploids.

Genome-wide analysis of DNA polymorphisms, the methylome and transcriptome revealed that multiple factors are associated with low pollen fertility in autotetraploid rice

Autotetraploid rice is a useful germplasm with high biomass production; however, low fertility is the main barrier in commercial utilization. In our previous study, differential expression of meiosis-related miRNAs was found to be involved in the pollen sterility of autotetraploid rice. However, genome-wide DNA variations and methylomes associated with low fertility of autotetraploid rice are still poorly understood. Here, we measured both global DNA variations and the methylome and compared them with the transcriptome during pollen development in autotetraploid rice by high-throughput sequencing. A total of 34416 SNPs, 6993 InDels, 1003 SVs and 25 CNVs were detected, and 11367 and 41117 differentially methylated regions showed hypermethylation and hypomethylation in 02428-4x. In total, 1110 genes displayed differentially expression in 02428-4x during meiosis, of these six harbored CNVs, including four upregulated genes with gain CNVs, such as LOC_Os11g38620. We identified 122 genes by comparing with the previous data that might be associated with low fertility during pollen development in 02428-4x. Of the 122 gens, 98 were displayed methylation and differential expression, including OsMADS98, CYP703A3 and OsABCG26. The downregulation of these three genes were confirmed by qPCR during meiosis of 02428-4x, which played pivotal roles in pollen fertility. These results indicate that the low fertility of autotetraploid rice is not only caused by the differential expression of genes involved in pollen development, but also by sequence variation and differential methylation, suggesting that the reason for pollen sterility in autotetraploid rice is complex and might be affected by multiple factors.