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

Transcriptomics and Metabolomics Reveal Biosynthetic Pathways and Regulatory Mechanisms of Phenylpropanes in Different Ploidy of Capsicum frutescens

Pepper is a significant cash crop, and Capsicum frutescens is an exemplary variety of pepper cultivated for its distinctive flavor and substantial nutritional value. Polyploidization of plants often leads to an increase in their biomass and improved stress tolerance, and thus has important applications in plant breeding and improvement. In this study, germplasm innovation was carried out by polyploidy induction of C. frutescens by colchicine. To investigate the effects of polyploidization on C. frutescens, we conducted transcriptomic and metabolomic analyses of diploids and homotetraploids of C. frutescens to gain insights into the mechanisms of metabolite composition and molecular regulation of C. frutescens by polyploidization. Based on the analysis of metabolomics and transcriptomics data, a total of 551 differential metabolites were identified in the leaves of C. frutescens of different ploidy and 634 genes were significantly differentially expressed. In comparison, 241 differential metabolites and 454 genes were significantly differentially expressed in the mature fruits of C. frutescens of different ploidy. Analysis of KEGG enrichment of differentially expressed genes and differential metabolites revealed that both differential metabolites and differentially expressed genes were highly enriched in the phenylalanine metabolic pathway. It is worth noting that phenylpropanoids are highly correlated with capsaicin synthesis and also have an effect on fruit development. Therefore, we comprehensively analyzed the phenylalanine metabolic pathway and found that chromosome doubling significantly down-regulated the expression of genes upstream of phenylalanine (PAL, 4CL), which promoted lignin accumulation, and we suggested that this might have led to the enlargement of polyploid C. frutescens fruits. This study provides some references for further research on the phenotypic traits of different ploidy of C. frutescens, cloning of key regulatory genes, and using genetic engineering techniques in C. frutescens breeding for germplasm improvement.

Transcriptional Silencing of 35S rDNA in Tragopogon porrifolius Correlates with Cytosine Methylation in Sequence-Specific Manner

Despite the widely accepted involvement of DNA methylation in the regulation of rDNA transcription, the relative participation of different cytosine methylation pathways is currently described only for a few model plants. Using PacBio, Bisulfite, and RNA sequencing; PCR; Southern hybridizations; and FISH, the epigenetic consequences of rDNA copy number variation were estimated in two T. porrifolius lineages, por1 and por2, the latter with more than twice the rDNA copy numbers distributed approximately equally between NORs on chromosomes A and D. The lower rDNA content in por1 correlated with significantly reduced (>90%) sizes of both D-NORs. Moreover, two (L and S) prominent rDNA variants, differing in the repetitive organization of intergenic spacers, were detected in por2, while only the S-rDNA variant was detected in por1. Transcriptional activity of S-rDNA in por1 was associated with secondary constriction of both A-NORs. In contrast, silencing of S-rDNA in por2 was accompanied by condensation of A-NORs, secondary constriction on D-NORs, and L-rDNA transcriptional activity, suggesting (i) bidirectional nucleolar dominance and (ii) association of S-rDNAs with A-NORs and L-rDNAs with D-NORs in T. porrifolius. Each S- and L-rDNA array was formed of several sub-variants differentiating both genetically (specific SNPs) and epigenetically (transcriptional efficiency and cytosine methylation). The most significant correlations between rDNA silencing and methylation were detected for symmetric CWG motifs followed by CG motifs. No correlations were detected for external cytosine in CCGs or asymmetric CHHs, where methylation was rather position-dependent, particularly for AT-rich variants. We conclude that variations in rDNA copy numbers in plant diploids can be accompanied by prompt epigenetic responses to maintain an appropriate number of active rDNAs. The methylation dynamics of CWGs are likely to be the most responsible for regulating silent and active rDNA states.

Uniparental expression of ribosomal RNA in ×Festulolium grasses: a link between the genome and nucleolar dominance

Genome or genomic dominance (GD) is a phenomenon observed in hybrids when one parental genome becomes dominant over the other. It is manifested by the replacement of chromatin of the submissive genome by that of the dominant genome and by biased gene expression. Nucleolar dominance (ND) - the functional expression of only one parental set of ribosomal genes in hybrids - is another example of an intragenomic competitive process which, however, concerns ribosomal DNA only. Although GD and ND are relatively well understood, the nature and extent of their potential interdependence is mostly unknown. Here, we ask whether hybrids showing GD also exhibit ND and, if so, whether the dominant genome is the same. To test this, we used hybrids between Festuca and Lolium grasses (Festulolium), and between two Festuca species in which GD has been observed (with Lolium as the dominant genome in Festulolium and F. pratensis in interspecific Festuca hybrids). Using amplicon sequencing of ITS1 and ITS2 of the 45S ribosomal DNA (rDNA) cluster and molecular cytogenetics, we studied the organization and expression of rDNA in leaf tissue in five hybrid combinations, four generations and 31 genotypes [F. pratensis × L. multiflorum (F1, F2, F3, BC1), L. multiflorum × F. pratensis (F1), L. multiflorum × F. glaucescens (F2), L. perenne × F. pratensis (F1), F. glaucescens × F. pratensis (F1)]. We have found that instant ND occurs in Festulolium, where expression of Lolium-type rDNA reached nearly 100% in all F1 hybrids and was maintained through subsequent generations. Therefore, ND and GD in Festulolium are manifested by the same dominant genome (Lolium). We also confirmed the concordance between GD and ND in an interspecific cross between two Festuca species.

Switch them off or not: selective rRNA gene repression in grasses

Nucleolar dominance (ND) is selective epigenetic silencing of 35-48S rDNA loci. In allopolyploids, it is frequently manifested at the cytogenetic level by the inactivation of nucleolar organiser region(s) (NORs) inherited from one or several evolutionary ancestors. Grasses are ecologically and economically one of the most important land plant groups, which have frequently evolved through hybridisation and polyploidisation events. Here we review common and unique features of ND phenomena in this monocot family from cytogenetic, molecular, and genomic perspectives. We highlight recent advances achieved by using an allotetraploid model grass, Brachypodium hybridum, where ND commonly occurs at a population level, and we cover modern genomic approaches that decipher structural features of core arrays of NORs.

Structure and Methylation of 35S rDNA in Allopolyploids Anemone multifida (2n = 4x = 32, BBDD) and Anemone baldensis (2n = 6x = 48, AABBDD) and Their Parental Species Show Evidence of Nucleolar Dominance

Transcriptional silencing of 35S rDNA loci inherited from one parental species is occurring relatively frequently in allopolyploids. However, molecular mechanisms by which it is selected for transcriptional silencing remain unclear. We applied NGS, silver staining and bisulfite sequencing to study the structure, expression and methylation landscape of 35S rDNA in two allopolyploids of common origin, allotetraploid Anemone multifida (2n = 4x = 32, genome composition BBDD) and allohexaploid A. baldensis (2n = 6x = 48, AABBDD), and their genome donors, A. sylvestris (2n = 16, AA), A. cylindrica (2n = 16, BB) and A. parviflora (2n = 16, DD). The size of the recovered 35S rDNA units varied from 10,489 bp in A. cylindrica to 12,084 bp in A. sylvestris. Anemone showed an organization typical of most ribosomal 35S rDNA composed of NTS, ETS, rRNA genes, TTS and TIS with structural features of plant IGS sequences and all functional elements needed for rRNA gene activity. The NTS was more variable than the ETS and consisted of SRs which are highly variable among Anemone. Five to six CpG-rich islands were found within the ETS. CpG island located adjacent to the transcription initiation site (TIS) was highly variable regarding the sequence size and methylation level and exhibited in most of the species lower levels of methylation than CpG islands located adjacent to the 18S rRNA gene. Our results uncover hypomethylation of A. sylvestris- and A. parviflora-derived 35S rDNA units in allopolyploids A. multifida and A. baldensis. Hypomethylation of A. parviflora-derived 35S rDNA was more prominent in A. baldensis than in A. multifida. We showed that A. baldensis underwent coupled A. sylvestris-derived 35S rDNA array expansion and A. parviflora-derived 35S rDNA copy number decrease that was accompanied by lower methylation level of A. sylvestris-derived 35S rDNA units in comparison to A. parviflora-derived 35S rDNA units. These observations suggest that in A. baldensis nucleolar dominance is directed toward A. sylvestris-derived chromosomes. This work broadens our current knowledge of the 35S rDNA organization in Anemone and provides evidence of the progenitor-specific 35S rDNA methylation in nucleolar dominance.

Fine structure and transcription dynamics of bread wheat ribosomal DNA loci deciphered by a multi-omics approach

Three out of four RNA components of ribosomes are encoded by 45S ribosomal DNA (rDNA) loci, which are organized as long head-to-tail tandem arrays of nearly identical units, spanning several megabases of sequence. Due to this structure, the rDNA loci are the major sources of gaps in genome assemblies, and gene copy number, sequence composition, and expression status of particular arrays remain elusive, especially in complex genomes harboring multiple loci. Here we conducted a multi-omics study to decipher the 45S rDNA loci in hexaploid bread wheat. Coupling chromosomal genomics with optical mapping, we reconstructed individual rDNA arrays, enabling locus-specific analyses of transcription activity and methylation status from RNA- and bisulfite-sequencing data. We estimated a total of 6,650 rDNA units in the bread wheat genome, with approximately 2,321, 3,910, 253, and 50 gene copies located in short arms of chromosomes 1B, 6B, 5D, and 1A, respectively. Only 1B and 6B loci contributed substantially to rRNA transcription at a roughly 2:1 ratio. The ratio varied among five tissues analyzed (embryo, coleoptile, root tip, primary leaf, mature leaf), being the highest (2.64:1) in mature leaf and lowest (1.72:1) in coleoptile. Cytosine methylation was considerably higher in CHG context in the silenced 5D locus as compared with the active 1B and 6B loci. In conclusion, a fine genomic organization and tissue-specific expression of rDNA loci were deciphered, for the first time, in a complex polyploid species. The results are discussed in the context of wheat evolution and transcription regulation.

The Nuclear 35S rDNA World in Plant Systematics and Evolution: A Primer of Cautions and Common Misconceptions in Cytogenetic Studies

The ubiquitous presence of rRNA genes in nuclear, plastid, and mitochondrial genomes has provided an opportunity to use genomic markers to infer patterns of molecular and organismic evolution as well as to assess systematic issues throughout the tree of life. The number, size, location, and activity of the 35S rDNA cistrons in plant karyotypes have been used as conventional cytogenetic landmarks. Their scrutiny has been useful to infer patterns of chromosomal evolution and the data have been used as a proxy for assessing species discrimination, population differentiation and evolutionary relationships. The correct interpretation of rDNA markers in plant taxonomy and evolution is not free of drawbacks given the complexities derived from the lability of the genetic architecture, the diverse patterns of molecular change, and the fate and evolutionary dynamics of the rDNA units in hybrids and polyploid species. In addition, the terminology used by independent authors is somewhat vague, which often complicates comparisons. To date, no efforts have been reported addressing the potential problems and limitations involved in generating, utilizing, and interpreting the data from the 35S rDNA in cytogenetics. This review discusses the main technical and conceptual limitations of these rDNA markers obtained by cytological and karyological experimental work, in order to clarify biological and evolutionary inferences postulated in a systematic and phylogenetic context. Also, we provide clarification for some ambiguity and misconceptions in terminology usually found in published work that may help to improve the usage of the 35S ribosomal world in plant evolution.

The fate of 35S rRNA genes in the allotetraploid grass Brachypodium hybridum

Nucleolar dominance (ND) consists of the reversible silencing of 35S/45S rDNA loci inherited from one of the ancestors of an allopolyploid. The molecular mechanisms by which one ancestral rDNA set is selected for silencing remain unclear. We applied a combination of molecular (Southern blot hybridization and reverse-transcription cleaved amplified polymorphic sequence analysis), genomic (analysis of variants) and cytogenetic (fluorescence in situ hybridization) approaches to study the structure, expression and epigenetic landscape of 35S rDNA in an allotetraploid grass that exhibits ND, Brachypodium hybridum (genome composition DDSS), and its putative progenitors, Brachypodium distachyon (DD) and Brachypodium stacei (SS). In progenitor genomes, B. stacei showed a higher intragenomic heterogeneity of rDNA compared with B. distachyon. In all studied accessions of B. hybridum, there was a reduction in the copy number of S homoeologues, which was accompanied by their inactive transcriptional status. The involvement of DNA methylation in CG and CHG contexts in the silencing of the S-genome rDNA loci was revealed. In the B. hybridum allotetraploid, ND is stabilized towards the D-genome units, irrespective of the polyphyletic origin of the species, and does not seem to be influenced by homoeologous 35S rDNA ratios and developmental stage.

Competition of Parental Genomes in Plant Hybrids

Interspecific hybridization represents one of the main mechanisms of plant speciation. Merging of two genomes from different subspecies, species, or even genera is frequently accompanied by whole-genome duplication (WGD). Besides its evolutionary role, interspecific hybridization has also been successfully implemented in multiple breeding programs. Interspecific hybrids combine agronomic traits of two crop species or can be used to introgress specific loci of interests, such as those for resistance against abiotic or biotic stresses. The genomes of newly established interspecific hybrids (both allopolyploids and homoploids) undergo dramatic changes, including chromosome rearrangements, amplifications of tandem repeats, activation of mobile repetitive elements, and gene expression modifications. To ensure genome stability and proper transmission of chromosomes from both parental genomes into subsequent generations, allopolyploids often evolve mechanisms regulating chromosome pairing. Such regulatory systems allow only pairing of homologous chromosomes and hamper pairing of homoeologs. Despite such regulatory systems, several hybrid examples with frequent homoeologous chromosome pairing have been reported. These reports open a way for the replacement of one parental genome by the other. In this review, we provide an overview of the current knowledge of genomic changes in interspecific homoploid and allopolyploid hybrids, with strictly homologous pairing and with relaxed pairing of homoeologs.

The Hordeum bulbosum 25S-18S rDNA region: comparison with Hordeum vulgare and other Triticeae

Hordeum vulgare and Hordeum bulbosum are two closely related barley species, which share a common H genome. H. vulgare has two nucleolar organizer regions (NORs), while the NOR of H. bulbosum is only one. We sequenced the 2.5 kb 25S-18S region in the rDNA of H. bulbosum and compared it to the same region in H. vulgare as well as to the other Triticeae. The region includes an intergenic spacer (IGS) with a number of subrepeats, a promoter, and an external transcribed spacer (5'ETS). The IGS of H. bulbosum downstream of 25S rRNA contains two 143-bp repeats and six 128-bp repeats. In contrast, the IGS in H. vulgare contains an array of seven 79-bp repeats and a varying number of 135-bp repeats. The 135-bp repeats in H. vulgare and the 128-bp repeats in H. bulbosum show similarity. Compared to H. vulgare, the 5'ETS of H. bulbosum is shorter. Additionally, the 5'ETS regions in H. bulbosum and H. vulgare diverged faster than in other Triticeae genera. Alignment of the Triticeae promoter sequences suggests that in Hordeum, as in diploid Triticum, transcription starts with guanine and not with adenine as it is in many other plants.

Unique Epigenetic Features of Ribosomal RNA Genes (rDNA) in Early Diverging Plants (Bryophytes)

Introduction: In plants, the multicopy genes encoding ribosomal RNA (rDNA) typically exhibit heterochromatic features and high level of DNA methylation. Here, we explored rDNA methylation in early diverging land plants from Bryophyta (15 species, 14 families) and Marchantiophyta (4 species, 4 families). DNA methylation was investigated by methylation-sensitive Southern blot hybridization in all species. We also carried out whole genomic bisulfite sequencing in Polytrichum formosum (Polytrichaceae) and Dicranum scoparium (Dicranaceae) and used available model plant methyloms (Physcomitrella patents and Marchantia polymorpha) to determine rDNA unit-wide methylation patterns. Chromatin structure was analyzed using fluorescence in situ hybridization (FISH) and immunoprecipitation (CHIP) assays. Results: In contrast to seed plants, bryophyte rDNAs were efficiently digested with methylation-sensitive enzymes indicating no or low levels of CG and CHG methylation in these loci. The rDNA methylom analyses revealed variation between species ranging from negligible (<3%, P. formosum, P. patens) to moderate (7 and 17% in M. polymorpha and D. scoparium, respectively) methylation levels. There were no differences between coding and noncoding parts of rDNA units and between gametophyte and sporophyte tissues. However, major satellite repeat and transposable elements were heavily methylated in P. formosum and D. scoparium. In P. formosum rDNA, the euchromatic H3K4m3 and heterochromatic H3K9m2 histone marks were nearly balanced contrasting the angiosperms data where H3K9m2 typically dominates rDNA chromatin. In moss interphase nuclei, rDNA was localized at the nucleolar periphery and its condensation level was high. Conclusions: Unlike seed plants, the rRNA genes seem to escape global methylation machinery in bryophytes. Distinct epigenetic features may be related to rDNA expression and the physiology of these early diverging plants that exist in haploid state for most of their life cycles.

Intragenomic heterogeneity of intergenic ribosomal DNA spacers in Cucurbita moschata is determined by DNA minisatellites with variable potential to form non-canonical DNA conformations

The intergenic spacer (IGS) of rDNA is frequently built of long blocks of tandem repeats. To estimate the intragenomic variability of such knotty regions, we employed PacBio sequencing of the Cucurbita moschata genome, in which thousands of rDNA copies are distributed across a number of loci. The rRNA coding regions are highly conserved, indicating intensive interlocus homogenization and/or high selection pressure. However, the IGS exhibits high intragenomic structural diversity. Two repeated blocks, R1 (300-1250 bp) and R2 (290-643 bp), account for most of the IGS variation. They exhibit minisatellite-like features built of multiple periodically spaced short GC-rich sequence motifs with the potential to adopt non-canonical DNA conformations, G-quadruplex-folded and left-handed Z-DNA. The mutual arrangement of these motifs can be used to classify IGS variants into five structural families. Subtle polymorphisms exist within each family due to a variable number of repeats, suggesting the coexistence of an enormous number of IGS variants. The substantial length and structural heterogeneity of IGS minisatellites suggests that the tempo of their divergence exceeds the tempo of the homogenization of rDNA arrays. As frequently occurring among plants, we hypothesize that their instability may influence transcription regulation and/or destabilize rDNA units, possibly spreading them across the genome.

Ribosomal DNA loci derived from Brachypodium stacei are switched off for major parts of the life cycle of Brachypodium hybridum

Nucleolar dominance is an epigenetic phenomenon that occurs in some plant and animal allopolyploids and hybrids, whereby only one ancestral set of 35S rRNA genes retains the ability to form the nucleolus while the rDNA loci derived from the other progenitor are transcriptionally silenced. There is substantial evidence that nucleolar dominance is regulated developmentally. This study focuses upon the establishment and/or maintenance of nucleolar dominance during different stages of development in the model grass allotetraploid Brachypodium hybridum. Fluorescence in situ hybridization with a 25S rDNA probe to cells in three-dimensional cytogenetic preparations showed that nucleolar dominance is present not only in root meristematic and differentiated cells of this species, but also in male meiocytes at prophase I, tetrads of microspores, and different embryonic tissues. The inactive state of Brachypodium stacei-originated rDNA loci was confirmed by silver staining. Only B. distachyon-derived 35S rDNA loci formed nucleoli in the aforementioned tissues, whereas B. stacei-like loci remained highly condensed and thus transcriptionally suppressed. The establishment of nucleolar dominance during earlier stages of B. hybridum embryo development cannot be ruled out. However, we propose that gradual pseudogenization of B. stacei-like loci in the evolution of the allotetraploid seems to be more likely.

Origin and parental genome characterization of the allotetraploid Stylosanthes scabra Vogel (Papilionoideae, Leguminosae), an important legume pasture crop

BACKGROUNDS AND AIMS

The genus Stylosanthes includes nitrogen-fixing and drought-tolerant species of considerable economic importance for perennial pasture, green manure and land recovery. Stylosanthes scabra is adapted to variable soil conditions, being cultivated to improve pastures and soils worldwide. Previous studies have proposed S. scabra as an allotetraploid species (2n = 40) with a putative diploid A genome progenitor S. hamata or S. seabrana (2n = 20) and the B genome progenitor S. viscosa (2n = 20). We aimed to provide conclusive evidence for the origin of S. scabra.

METHODS

We performed fluorescence in situ hybridization (FISH) and genomic in situ hybridization (GISH) experiments and Illumina paired-end sequencing of S. scabra, S. hamata and S. viscosa genomic DNA, to assemble and compare complete ribosomal DNA (rDNA) units and chloroplast genomes. Plastome- and genome-wide single nucleotide variation detection was also performed.

KEY RESULTS

GISH and phylogenetic analyses of plastid DNA and rDNA sequences support that S. scabra is an allotetraploid formed from 0.63 to 0.52 million years ago (Mya), from progenitors with a similar genome structure to the maternal donor S. hamata and the paternal donor S. viscosa. FISH revealed a non-additive number of 35S rDNA sites in S. scabra compared with its progenitors, indicating the loss of one locus from A genome origin. In S. scabra, most 5S rDNA units were similar to S. viscosa, while one 5S rDNA site of reduced size most probably came from an A genome species as revealed by GISH and in silico analysis.

CONCLUSIONS

Our approach combined whole-plastome and rDNA assembly with additional cytogenetic analysis to shed light successfully on the allotetraploid origin of S. scabra. We propose a Middle Pleistocene origin for S. scabra involving species with maternal A and paternal B genomes. Our data also suggest that variation found in rDNA units in S. scabra and its progenitors reveals differences that can be explained by homogenization, deletion and amplification processes that have occurred since its origin.

Re-exploration of U's Triangle Brassica Species Based on Chloroplast Genomes and 45S nrDNA Sequences

The concept of U's triangle, which revealed the importance of polyploidization in plant genome evolution, described natural allopolyploidization events in Brassica using three diploids [B. rapa (A genome), B. nigra (B), and B. oleracea (C)] and derived allotetraploids [B. juncea (AB genome), B. napus (AC), and B. carinata (BC)]. However, comprehensive understanding of Brassica genome evolution has not been fully achieved. Here, we performed low-coverage (2-6×) whole-genome sequencing of 28 accessions of Brassica as well as of Raphanus sativus [R genome] to explore the evolution of six Brassica species based on chloroplast genome and ribosomal DNA variations. Our phylogenomic analyses led to two main conclusions. (1) Intra-species-level chloroplast genome variations are low in the three allotetraploids (2~7 SNPs), but rich and variable in each diploid species (7~193 SNPs). (2) Three allotetraploids maintain two 45SnrDNA types derived from both ancestral species with maternal dominance. Furthermore, this study sheds light on the maternal origin of the AC chloroplast genome. Overall, this study clarifies the genetic relationships of U's triangle species based on a comprehensive genomics approach and provides important genomic resources for correlative and evolutionary studies.

The fate of ribosomal RNA genes in spontaneous polyploid dogrose hybrids [Rosa L. sect. Caninae (DC.) Ser.] exhibiting non-symmetrical meiosis

Dogroses represent an exceptional system for studying the effects of genome doubling and hybridization: their asymmetrical meiosis enables recombination in bi-parentally inherited chromosomes but prevents it in maternally inherited ones. We employed fluorescent in situ hybridization, genome skimming, amplicon sequencing of genomic and cDNA as well as conventional cloning of nuclear ribosomal DNA in two phylogenetically distinct pentaploid (2n = 5x = 35) species, Rosa canina and Rosa inodora, and their naturally occurring reciprocal hybrids, Rosa dumalis (5x) and Rosa agrestis (5x, 6x). Both progenitor species differed in composition, meiotic behaviour and expression of rDNA loci: R. canina (five 18S and 5-8 5S loci) was dominated by the Canina ribotypes, but R. inodora (four 18S loci and 7-8 5S loci) by the Rubiginosa ribotype. The co-localized 5S/18S loci occurred on either bivalent-forming (R. canina) or univalent-forming (R. inodora) chromosomes. Ribosomal DNA loci were additively inherited; however, the Canina ribotypes were dominantly expressed, even in genotypes with relatively low copy number of these genes. Moreover, we observed rDNA homogenization towards the paternally transmitted Canina ribotype in 6x R. agrestis. The here-observed variation in arrangement and composition of rDNA types between R. canina and R. inodora suggests the involvement of different genomes in bivalent formation. This results supports the hypothesis that the asymmetrical meiosis arose at least twice by independent ancient hybridization events.

Unstable Inheritance of 45S rRNA Genes in Arabidopsis thaliana

The considerable genome size variation in Arabidopsis thaliana has been shown largely to be due to copy number variation (CNV) in 45S ribosomal RNA (rRNA) genes. Surprisingly, attempts to map this variation by means of genome-wide association studies (GWAS) failed to identify either of the two likely sources, namely the nucleolus organizer regions (NORs). Instead, GWAS implicated a trans-acting locus, as if rRNA gene CNV was a phenotype rather than a genotype. To explain these results, we investigated the inheritance and stability of rRNA gene copy number using the variety of genetic resources available in A. thaliana - F2 crosses, recombinant inbred lines, the multiparent advanced-generation inter-cross population, and mutation accumulation lines. Our results clearly show that rRNA gene CNV can be mapped to the NORs themselves, with both loci contributing equally to the variation. However, NOR size is unstably inherited, and dramatic copy number changes are visible already within tens of generations, which explains why it is not possible to map the NORs using GWAS. We did not find any evidence of trans-acting loci in crosses, which is also expected since changes due to such loci would take very many generations to manifest themselves. rRNA gene copy number is thus an interesting example of "missing heritability"-a trait that is heritable in pedigrees, but not in the general population.

Phenotypic reversion in fas mutants of Arabidopsis thaliana by reintroduction of FAS genes: variable recovery of telomeres with major spatial rearrangements and transcriptional reprogramming of 45S rDNA genes

Arabidopsis thaliana mutants dysfunctional in the evolutionarily conserved protein complex chromatin assembly factor-1 (CAF-1), which deposits the canonical histone H3 variant H3.1 during DNA synthesis-dependent chromatin assembly, display complex phenotypic changes including meristem and growth alterations, sensitivity to DNA-damaging agents, and reduced fertility. We reported previously that mutants in the FAS1 subunit of CAF-1 progressively lose telomere and 45S rDNA repeats. Here we show that multiple aspects of the fas phenotype are recovered immediately on expression of a reintroduced FAS1 allele, and are clearly independent of the recovery of rDNA copy-numbers and telomeres. In reverted lines, 45S rDNA genes are recovered to diverse levels with a strikingly different representation of their variants, and the typical association of nucleolar organizing region 4 with the nucleolus is perturbed. One of 45S rDNA variants (VAR1), which is silenced in wild-type (WT) plants without mutation history (Col-0 WT), dominates the expression pattern, whereas VAR2 is dominant in Col-0 WT plants. We propose an explanation for the variability of telomere and 45S rDNA repeats associated with CAF-1 function, suggesting that the differences in nuclear partitioning and expression of the rDNA variants in fas mutants and their revertants provide a useful experimental system to study genetic and epigenetic factors in gene dosage compensation.

Interpopulation hybridization generates meiotically stable rDNA epigenetic variants in allotetraploid Tragopogon mirus

Uniparental silencing of 35S rRNA genes (rDNA), known as nucleolar dominance (ND), is common in interspecific hybrids. Allotetraploid Tragopogon mirus composed of Tragopogon dubius (d) and Tragopogon porrifolius (p) genomes shows highly variable ND. To examine the molecular basis of such variation, we studied the genetic and epigenetic features of rDNA homeologs in several lines derived from recently and independently formed natural populations. Inbred lines derived from T. mirus with a dominant d-rDNA homeolog transmitted this expression pattern over generations, which may explain why it is prevalent among natural populations. In contrast, lines derived from the p-rDNA dominant progenitor were meiotically unstable, frequently switching to co-dominance. Interpopulation crosses between progenitors displaying reciprocal ND resulted in d-rDNA dominance, indicating immediate suppression of p-homeologs in F1 hybrids. Original p-rDNA dominance was not restored in later generations, even in those segregants that inherited the corresponding parental rDNA genotype, thus indicating the generation of additional p-rDNA and d-rDNA epigenetic variants. Despite preserved intergenic spacer (IGS) structure, they showed altered cytosine methylation and chromatin condensation patterns, and a correlation between expression, hypomethylation of RNA Pol I promoters and chromatin decondensation was apparent. Reversion of such epigenetic variants occurred rarely, resulting in co-dominance maintained in individuals with distinct genotypes. Generally, interpopulation crosses may generate epialleles that are not present in natural populations, underlying epigenetic dynamics in young allopolyploids. We hypothesize that highly expressed variants with distinct IGS features may induce heritable epigenetic reprogramming of the partner rDNA arrays, harmonizing the expression of thousands of genes in allopolyploids.

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.

Cytomolecular Analysis of Ribosomal DNA Evolution in a Natural Allotetraploid Brachypodium hybridum and Its Putative Ancestors-Dissecting Complex Repetitive Structure of Intergenic Spacers

Nucleolar dominance is an epigenetic phenomenon associated with nuclear 35S rRNA genes and consists in selective suppression of gene loci inherited from one of the progenitors in the allopolyploid. Our understanding of the exact mechanisms that determine this process is still fragmentary, especially in case of the grass species. This study aimed to shed some light on the molecular basis of this genome-specific inactivation of 35S rDNA loci in an allotetraploid Brachypodium hybridum (2n = 30), which arose from the interspecific hybridization between two diploid ancestors that were very similar to modern B. distachyon (2n = 10) and B. stacei (2n = 20). Using fluorescence in situ hybridization with 25S rDNA and chromosome-specific BAC clones as probes we revealed that the nucleolar dominance is present not only in meristematic root-tip cells but also in differentiated cell fraction of B. hybridum. Additionally, the intergenic spacers (IGSs) from both of the putative ancestors and the allotetraploid were sequenced and analyzed. The presumptive transcription initiation sites, spacer promoters and repeated elements were identified within the IGSs. Two different length variants, 2.3 and 3.5 kb, of IGSs were identified in B. distachyon and B. stacei, respectively, however only the IGS that had originated from B. distachyon-like ancestor was present in the allotetraploid. The amplification pattern of B. hybridum IGSs suggests that some genetic changes occurred in inactive B. stacei-like rDNA loci during the evolution of the allotetraploid. We hypothesize that their preferential silencing is an effect of structural changes in the sequence rather than just the result of the sole inactivation at the epigenetic level.