Evolutionary trajectory of organelle-derived nuclear DNAs in the Triticum/Aegilops complex species.

Organelle-derived nuclear DNAs, nuclear plastid DNAs (NUPTs), and nuclear mitochondrial DNAs (NUMTs) have been identified in plants. Most, if not all, genes residing in NUPTs/NUMTs (NUPGs/NUMGs) are known to be inactivated and pseudogenized. However, the role of epigenetic control in silencing NUPGs/NUMGs and the dynamic evolution of NUPTs/NUMTs with respect to organismal phylogeny remain barely explored. Based on the available nuclear and organellar genomic resources of wheat (genus Triticum) and goat grass (genus Aegilops) within Triticum/Aegilops complex species, we investigated the evolutionary fates of NUPTs/NUMTs in terms of their epigenetic silencing and their dynamic occurrence rates in the nuclear diploid genomes and allopolyploid subgenomes. NUPTs and NUMTs possessed similar genomic atlas, including (i) predominantly located in intergenic regions and preferential integration to gene regulation regions and (ii) generating sequence variations in the nuclear genome. Unlike nuclear indigenous genes, the alien NUPGs/NUMGs were associated with repressive epigenetic signals, namely high levels of DNA methylation and low levels of active histone modifications. Phylogenomic analyses suggested that the species-specific and gradual accumulation of NUPTs/NUMTs accompanied the speciation processes. Moreover, based on further pan-genomic analyses, we found significant subgenomic asymmetry in the NUPT/NUMT occurrence, which accumulated during allopolyploid wheat evolution. Our findings provide insight into the dynamic evolutionary fates of organelle-derived nuclear DNA in plants.

Cell type-specific cytonuclear coevolution in three allopolyploid plant species.

Cytonuclear disruption may accompany allopolyploid evolution as a consequence of the merger of different nuclear genomes in a cellular environment having only one set of progenitor organellar genomes. One path to reconcile potential cytonuclear mismatch is biased expression for maternal gene duplicates (homoeologs) encoding proteins that target to plastids and/or mitochondria. Assessment of this transcriptional form of cytonuclear coevolution at the level of individual cells or cell types remains unexplored. Using single-cell (sc-) and single-nucleus (sn-) RNAseq data from eight tissues in three allopolyploid species, we characterized cell type-specific variations of cytonuclear coevolutionary homoeologous expression and demonstrated the temporal dynamics of expression patterns across development stages during cotton fiber development. Our results provide unique insights into transcriptional cytonuclear coevolution in plant allopolyploids at the single-cell level.

Meiotic pairing irregularity and homoeologous chromosome compensation cause rapid karyotype variation in synthetic allotetraploid wheat.

Allopolyploidization may initiate rapid evolution due to heritable karyotypic changes. The types and extents of these changes, the underlying causes, and their effects on phenotype remain to be fully understood. Here, we designed experimental populations suitable to address these issues using a synthetic allotetraploid wheat. We show that extensive variation in both chromosome number (NCV) and structure (SCV) accumulated in a selfed population of a synthetic allotetraploid wheat (genome Sb Sb DD). The combination of NCVs and SCVs generated massive organismal karyotypic heterogeneity. NCVs and SCVs were intrinsically correlated and highly variable across the seven sets of homoeologous chromosomes. Both NCVs and SCVs stemmed from meiotic pairing irregularity (presumably homoeologous pairing) but were also constrained by homoeologous chromosome compensation. We further show that homoeologous meiotic pairing was positively correlated with sequence synteny at the subtelomeric regions of both chromosome arms, but not with genic nucleotide similarity per se. Both NCVs and SCVs impacted phenotypic traits but only NCVs caused significant reduction in reproductive fitness. Our results implicate factors influencing meiotic homoeologous chromosome pairing and reveal the type and extent of karyotypic variation and its immediate phenotypic manifestation in synthetic allotetraploid wheat. This has relevance for our understanding of allopolyploid evolution.

Compensatory Genetic and Transcriptional Cytonuclear Coordination in Allopolyploid Lager Yeast (Saccharomyces pastorianus).

Cytonuclear coordination between biparental-nuclear genomes and uniparental-cytoplasmic organellar genomes in plants is often resolved by genetic and transcriptional cytonuclear responses. Whether this mechanism also acts in allopolyploid members of other kingdoms is not clear. Additionally, cytonuclear coordination of interleaved allopolyploid cells/individuals within the same population is underexplored. The yeast Saccharomyces pastorianus provides the opportunity to explore cytonuclear coevolution during different growth stages and from novel dimensions. Using S. pastorianus cells from multiple growth stages in the same environment, we show that nuclear mitochondria-targeted genes have undergone both asymmetric gene conversion and growth stage-specific biased expression favoring genes from the mitochondrial genome donor (Saccharomyces eubayanus). Our results suggest that cytonuclear coordination in allopolyploid lager yeast species entails an orchestrated and compensatory genetic and transcriptional evolutionary regulatory shift. The common as well as unique properties of cytonuclear coordination underlying allopolyploidy between unicellular yeasts and higher plants offers novel insights into mechanisms of cytonuclear evolution associated with allopolyploid speciation.

Chromatin architectural alterations due to null mutation of a major CG methylase in rice.

Associations between 3D chromatin architectures and epigenetic modifications have been characterized in animals. However, any impact of DNA methylation on chromatin architecture in plants is understudied, which is confined to Arabidopsis thaliana. Because plant species differ in genome size, composition, and overall chromatin packing, it is unclear to what extent findings from A. thaliana hold in other species. Moreover, the incomplete chromatin architectural profiles and the low-resolution high-throughput chromosome conformation capture (Hi-C) data from A. thaliana have hampered characterizing its subtle chromatin structures and their associations with DNA methylation. We constructed a high-resolution Hi-C interaction map for the null OsMET1-2 (the major CG methyltransferase in rice) mutant (osmet1-2) and isogenic wild-type rice (WT). Chromatin structural changes occurred in osmet1-2, including intra-/inter-chromosomal interactions, compartment transition, and topologically associated domains (TAD) variations. Our findings provide novel insights into the potential function of DNA methylation in TAD formation in rice and confirmed DNA methylation plays similar essential roles in chromatin packing in A. thaliana and rice.

Effects of homoeologous exchange on gene expression and alternative splicing in a newly formed allotetraploid wheat.

Homoeologous exchange (HE) is a major mechanism generating post-polyploidization genetic variation with important evolutionary consequences. However, the direct impacts of HE on gene expression and transcript diversity in allopolyploids without the intertwined evolutionary processes remain to be fully understood. Here, we analyzed high-throughput RNA-seq data of young leaves from plant groups of a synthetic allotetraploid wheat (AADD), which contained variable numbers of HEs. We aimed to investigate if and to which extent HE directly impacts gene expression and alternative splicing (AS). We found that HE impacts expression of genes located within HE regions primarily via a cis-acting dosage effect, which led to significant changes in the total expression level of homoeologous gene pairs, especially for homoeologs whose original expression was biased. In parallel, HE also influences expression of a large number of genes residing in non-HE regions by trans-regulation leading to convergent expression of homoeologs. Intriguingly, when taking the original relative homoeolog expression states into account, homoeolog pairs under trans-effect are more prone to manifesting a convergent response to the HEs whereas those under cis-regulation tended to show further exacerbated subgenome-biased expression. Moreover, HE-induced quantitative, largely individual-specific, changes of AS events were detected. Similar to homoeologous expression, homoeo-AS events under trans-effect were more responsive to HE. HE therefore exerts multifaceted immediate effects on gene expression and, to a less extent, on individualized transcript diversity in nascent allopolyploidy.

Mutation of a major CG methylase alters genome-wide lncRNA expression in rice.

Plant long non-coding RNAs (lncRNAs) function in diverse biological processes, and lncRNA expression is under epigenetic regulation, including by cytosine DNA methylation. However, it remains unclear whether 5-methylcytosine (5mC) plays a similar role in different sequence contexts (CG, CHG, and CHH). In this study, we characterized and compared the profiles of genome-wide lncRNA profiles (including long intergenic non-coding RNAs [lincRNAs] and long noncoding natural antisense transcripts [lncNATs]) of a null mutant of the rice DNA methyltransferase 1, OsMET1-2 (designated OsMET1-2-/-) and its isogenic wild type (OsMET1-2+/+). The En/Spm transposable element (TE) family, which was heavily methylated in OsMET1-2+/+, was transcriptionally de-repressed in OsMET1-2-/- due to genome-wide erasure of CG methylation, and this led to abundant production of specific lncRNAs. In addition, RdDM-mediated CHH hypermethylation was increased in the 5'-upstream genomic regions of lncRNAs in OsMET1-2-/-. The positive correlation between the expression of lincRNAs and that of their proximal protein-coding genes was also analyzed. Our study shows that CG methylation negatively regulates the TE-related expression of lncRNA and demonstrates that CHH methylation is also involved in the regulation of lncRNA expression.

Homoeologous exchanges occur through intragenic recombination generating novel transcripts and proteins in wheat and other polyploids.

Recombination between homeologous chromosomes, also known as homeologous exchange (HE), plays a significant role in shaping genome structure and gene expression in interspecific hybrids and allopolyploids of several plant species. However, the molecular mechanisms that govern HEs are not well understood. Here, we studied HE events in the progeny of a nascent allotetraploid (genome AADD) derived from two diploid progenitors of hexaploid bread wheat using cytological and whole-genome sequence analyses. In total, 37 HEs were identified and HE junctions were mapped precisely. HEs exhibit typical patterns of homologous recombination hotspots, being biased toward low-copy, subtelomeric regions of chromosome arms and showing association with known recombination hotspot motifs. But, strikingly, while homologous recombination preferentially takes place upstream and downstream of coding regions, HEs are highly enriched within gene bodies, giving rise to novel recombinant transcripts, which in turn are predicted to generate new protein fusion variants. To test whether this is a widespread phenomenon, a dataset of high-resolution HE junctions was analyzed for allopolyploid Brassica, rice, Arabidopsis suecica, banana, and peanut. Intragenic recombination and formation of chimeric genes was detected in HEs of all species and was prominent in most of them. HE thus provides a mechanism for evolutionary novelty in transcript and protein sequences in nascent allopolyploids.

Transcriptome Changes during Major Developmental Transitions Accompanied with Little Alteration of DNA Methylome in Two Pleurotus Species.

Pleurotus tuoliensis (Pt) and P. eryngii var. eryngii (Pe) are important edible mushrooms. The epigenetic and gene expression signatures characterizing major developmental transitions in these two mushrooms remain largely unknown. Here, we report global analyses of DNA methylation and gene expression in both mushrooms across three major developmental transitions, from mycelium to primordium and to fruit body, by whole-genome bisulfite sequencing (WGBS) and RNA-seq-based transcriptome profiling. Our results revealed that in both Pt and Pe the landscapes of methylome are largely stable irrespective of genomic features, e.g., in both protein-coding genes and transposable elements (TEs), across the developmental transitions. The repressive impact of DNA methylation on expression of a small subset of genes is likely due to TE-associated effects rather than their own developmental dynamics. Global expression of gene orthologs was also broadly conserved between Pt and Pe, but discernible interspecific differences exist especially at the fruit body formation stage, and which are primarily due to differences in trans-acting factors. The methylome and transcriptome repertories we established for the two mushroom species may facilitate further studies of the epigenetic and transcriptional regulatory mechanisms underpinning gene during development in Pleurotus and related genera.

The evolution of genomic and epigenomic features in two Pleurotus fungi.

Pleurotus tuoliensis (Bailinggu, designated Pt) and P. eryngii var. eryngii (Xingbaogu, designated Pe) are highly valued edible mushrooms. We report de novo assemblies of high-quality genomes for both mushrooms based on PacBio RS II sequencing and annotation of all identified genes. A comparative genomics analysis between Pt and Pe with P. ostreatus as an outgroup taxon revealed extensive genomic divergence between the two mushroom genomes primarily due to the rapid gain of taxon-specific genes and disruption of synteny in either taxon. The re-appraised phylogenetic relationship between Pt and Pe at the genome-wide level validates earlier proposals to designate Pt as an independent species. Variation of the identified wood-decay-related gene content can largely explain the variable adaptation and host specificity of the two mushrooms. On the basis of the two assembled genome sequences, methylomes and the regulatory roles of DNA methylation in gene expression were characterized and compared. The genome, methylome and transcriptome data of these two important mushrooms will provide valuable information for advancing our understanding of the evolution of Pleurotus and related genera and for facilitating genome- and epigenome-based strategies for mushroom breeding.

Transgenerationally Precipitated Meiotic Chromosome Instability Fuels Rapid Karyotypic Evolution and Phenotypic Diversity in an Artificially Constructed Allotetraploid Wheat (AADD).

Although a distinct karyotype with defined chromosome number and structure characterizes each biological species, it is intrinsically labile. Polyploidy or whole-genome duplication has played a pervasive and ongoing role in the evolution of all eukaryotes, and is the most dramatic force known to cause rapid karyotypic reconfiguration, especially at the initial stage. However, issues concerning transgenerational propagation of karyotypic heterogeneity and its translation to phenotypic diversity in nascent allopolyploidy, at the population level, have yet to be studied in detail. Here, we report a large-scale examination of transgenerationally propagated karyotypic heterogeneity and its phenotypic manifestation in an artificially constructed allotetraploid with a genome composition of AADD, that is, involving two of the three progenitor genomes of polyploid wheat. Specifically, we show that 1) massive organismal karyotypic heterogeneity is precipitated after 12 consecutive generations of selfing from a single euploid founder individual, 2) there exist dramatic differences in aptitudes between subgenomes and among chromosomes for whole-chromosome gain and/or loss and structural variations, 3) majority of the numerical and structural chromosomal variations are concurrent due to mutual contingency and possible functional constraint, 4) purposed and continuous selection and propagation for euploidy over generations did not result in enhanced karyotype stabilization, and 5) extent of karyotypic variation correlates with variability of phenotypic manifestation. Together, our results document that allopolyploidization catalyzes rampant and transgenerationally heritable organismal karyotypic heterogeneity that drives population-level phenotypic diversification, which lends fresh empirical support to the still contentious notion that whole-genome duplication enhances organismal evolvability.