Uniparental silencing of 5S rRNA genes in plant allopolyploids - insights from Cardamine (Brassicaceae).

While the phenomenon of uniparental silencing of 35S rDNA in interspecific hybrids and allopolyploids is well documented, there is a notable absence of information regarding whether such silencing extends to the 5S RNA component of ribosomes. To address this gap in knowledge, we analyzed the 5S and 35S rDNA expression in Cardamine (Brassicaceae) allopolyploids, namely C. × insueta (2n = 3x = 24, genome composition RRA), C. flexuosa (2n = 4x = 32, AAHH), and C. scutata (2n = 4x = 32, PPAA) which share a common diploid ancestor (AA). We employed high-throughput sequencing of transcriptomes and genomes and phylogenetic analyses of 5S rRNA variants. The genomic organization of rDNA was further scrutinized through clustering and fluorescence in situ hybridization. In the C. × insueta allotriploid, we observed uniparental dominant expression of 5S and 35S rDNA loci. In the C. flexuosa and C. scutata allotetraploids, the expression pattern differed, with the 35S rDNA being expressed from the A subgenome, whereas the 5S rDNA was expressed from the partner subgenome. Both C. flexuosa and C. scutata but not C. × insueta showed copy and locus number changes. We conclude that in stabilized allopolyploids, transcription of ribosomal RNA components occurs from different subgenomes. This phenomenon appears to result in the formation of chimeric ribosomes comprising rRNA molecules derived from distinct parental origins. We speculate that the interplay of epigenetic silencing and rDNA rearrangements introduces an additional layer of variation in multimolecule ribosomal complexes, potentially contributing to the evolutionary success of allopolyploids.

The Dynamic Interplay Between Ribosomal DNA and Transposable Elements: A Perspective From Genomics and Cytogenetics.

Although both are salient features of genomes, at first glance ribosomal DNAs and transposable elements are genetic elements with not much in common: whereas ribosomal DNAs are mainly viewed as housekeeping genes that uphold all prime genome functions, transposable elements are generally portrayed as selfish and disruptive. These opposing characteristics are also mirrored in other attributes: organization in tandem (ribosomal DNAs) versus organization in a dispersed manner (transposable elements); evolution in a concerted manner (ribosomal DNAs) versus evolution by diversification (transposable elements); and activity that prolongs genomic stability (ribosomal DNAs) versus activity that shortens it (transposable elements). Re-visiting relevant instances in which ribosomal DNA-transposable element interactions have been reported, we note that both repeat types share at least four structural and functional hallmarks: (1) they are repetitive DNAs that shape genomes in evolutionary timescales, (2) they exchange structural motifs and can enter co-evolution processes, (3) they are tightly controlled genomic stress sensors playing key roles in senescence/aging, and (4) they share common epigenetic marks such as DNA methylation and histone modification. Here, we give an overview of the structural, functional, and evolutionary characteristics of both ribosomal DNAs and transposable elements, discuss their roles and interactions, and highlight trends and future directions as we move forward in understanding ribosomal DNA-transposable element associations.

Epigenetic histone H3 phosphorylation marks discriminate between univalent- and bivalent-forming chromosomes during canina asymmetrical meiosis.

BACKGROUND AND AIMS: Dogroses (Rosa sect. Caninae) are mostly pentaploid, bearing 2n = 5x = 35 chromosomes in somatic cells. They evolved a unique form of asymmetrical meiosis characterized by two types of chromosomes: (1) chromosomes forming bivalents and distributed in the normal sexual way; and (2) chromosomes occurring as univalents and transferred by a female gamete only. In the mature pollen of pentaploid species, seven bivalent-derived chromosomes are transmitted to offspring, and 21 unpaired univalent chromosomes are eliminated during microsporogenesis. To discriminate between bivalent- and univalent-forming chromosomes, we studied histone H3 phosphorylation patterns regulating meiotic chromosome condensation and segregation. METHODS: We analysed histone modification patterns during male canina meiosis in two representative dogrose species, 5x Rosa canina and 5x Rosa rubiginosa, by immunohistochemical and molecular cytogenetics approaches. Immunostaining of meiotic cells included α-tubulin, histone H3 phosphorylation (H3S10p, H3S28p and H3T3p) and methylation (H3K4me3 and H3K27me3) marks. In addition, fluorescent in situ hybridization was carried out with an 18S rDNA probe. KEY RESULTS: In the first meiotic division, univalent chromosomes underwent equational division into chromatids, while homologues in bivalents were segregated as regular dyads. In diakinesis, bivalent chromosomes displayed strong H3 phosphorylation signals in proximal regions, spreading to the rest of the chromosome. In contrast, in univalents, the H3 phosphorylation signals were weaker, occurring mostly outside proximal regions largely overlapping with the H3K4me3 signals. Reduced phosphorylation was associated with relative under-condensation of the univalent chromosomes, particularly at early diakinesis. CONCLUSIONS: We hypothesize that the absence of pairing and/or recombination in univalent chromosomes negatively affects the histone H3 phosphorylation of their chromatin and perhaps the loading of meiotic-specific cohesins. This apparently destabilizes cohesion of sister chromatids, leading to their premature split in the first meiotic division.

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.

Mega-sized pericentromeric blocks of simple telomeric repeats and their variants reveal patterns of chromosome evolution in ancient Cycadales genomes.

Simple telomeric repeats composed of six to seven iterating nucleotide units are important sequences typically found at the ends of chromosomes. Here we analyzed their abundance and homogeneity in 42 gymnosperm (29 newly sequenced), 29 angiosperm (one newly sequenced), and eight bryophytes using bioinformatics, conventional cytogenetic and molecular biology approaches to explore their diversity across land plants. We found more than 10 000-fold variation in the amounts of telomeric repeats among the investigated taxa. Repeat abundance was positively correlated with increasing intragenomic sequence heterogeneity and occurrence at non-telomeric positions, but there was no correlation with genome size. The highest abundance/heterogeneity was found in the gymnosperm genus Cycas (Cycadaceae), in which megabase-sized blocks of telomeric repeats (i.e., billions of copies) were identified. Fluorescent in situ hybridization experiments using variant-specific probes revealed canonical Arabidopsis-type telomeric TTTAGGG repeats at chromosome ends, while pericentromeric blocks comprised at least four major telomeric variants with decreasing abundance: TTTAGGG>TTCAGGG >TTTAAGG>TTCAAGG. Such a diversity of repeats was not found in the sister cycad family Zamiaceae or in any other species analyzed. Using immunocytochemistry, we showed that the pericentromeric blocks of telomeric repeats overlapped with histone H3 serine 10 phosphorylation signals. We show that species of Cycas have amplified their telomeric repeats in centromeric and telomeric positions on telocentric chromosomes to extraordinary high levels. The ancestral chromosome number reconstruction suggests their occurrence is unlikely to be the product of ancient Robertsonian chromosome fusions. We speculate as to how the observed chromosome dynamics may be associated with the diversification of cycads.

Intragenomic rDNA variation - the product of concerted evolution, mutation, or something in between?

The classical model of concerted evolution states that hundreds to thousands of ribosomal DNA (rDNA) units undergo homogenization, making the multiple copies of the individual units more uniform across the genome than would be expected given mutation frequencies and gene redundancy. While the universality of this over 50-year-old model has been confirmed in a range of organisms, advanced high throughput sequencing techniques have also revealed that rDNA homogenization in many organisms is partial and, in rare cases, even apparently failing. The potential underpinning processes leading to unexpected intragenomic variation have been discussed in a number of studies, but a comprehensive understanding remains to be determined. In this work, we summarize information on variation or polymorphisms in rDNAs across a wide range of taxa amongst animals, fungi, plants, and protists. We discuss the definition and description of concerted evolution and describe whether incomplete concerted evolution of rDNAs predominantly affects coding or non-coding regions of rDNA units and if it leads to the formation of pseudogenes or not. We also discuss the factors contributing to rDNA variation, such as interspecific hybridization, meiotic cycles, rDNA expression status, genome size, and the activity of effector genes involved in genetic recombination, epigenetic modifications, and DNA editing. Finally, we argue that a combination of approaches is needed to target genetic and epigenetic phenomena influencing incomplete concerted evolution, to give a comprehensive understanding of the evolution and functional consequences of intragenomic variation in rDNA.

Intact ribosomal DNA arrays of Potentilla origin detected in Erythronium nucleus suggest recent eudicot-to-monocot horizontal transfer.

During our initial phylogenetic study of the monocot genus Erythronium (Liliaceae), we observed peculiar eudicot-type internal transcribed spacer (ITS) sequences in a dataset derived from genomic DNA of Erythronium dens-canis. This raised the possibility of horizontal transfer of a eudicot alien ribosomal DNA (rDNA) into the Erythronium genome. In this work we aimed to support this hypothesis by carrying out genomic, molecular, and cytogenetic analyses. Genome skimming coupled by PacBio HiFi sequencing of a bacterial artificial chromosome clone derived from flow-sorted nuclei was used to characterise the alien 45S rDNA. Integration of alien rDNA in the recipient genome was further proved by Southern blotting and fluorescence in situ hybridization using specific probes. Alien rDNA, nested among Potentilla species in phylogenetic analysis, likely entered the Erythronium lineage in the common ancestor of E. dens-canis and E. caucasicum. Transferred eudicot-type rDNA preserved its tandemly arrayed feature on a single chromosome and was found to be transcribed in the monocot host, albeit much less efficiently than the native counterpart. This study adds a new example to the rarely documented nuclear-to-nuclear jumps of DNA between eudicots and monocots while holding the scientific community continually in suspense about the mode of DNA transfer.

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.

Acid-Catalyzed RNA-Oligomerization from 3',5'-cGMP.

The assembly of ancient informational polymers from nucleotide precursors is the central challenge of life's origin on our planet. Among the possible solutions, dry polymerization of 3',5'-cyclic guanosine monophosphate (3',5'-cGMP) has been proposed as a candidate to create oligonucleotides of 15-20 units in length. However, the reported sensitivity of the reaction to the presence of cations raised questions of whether this chemistry could be relevant in a geological context. The experiments in this study show that the presence of cations is not restrictive as long as the reaction is conducted in an acidic environment, in contrast to previous reports that suggested optimal conditions at pH 9.

Analyses of the Updated "Animal rDNA Loci Database" with an Emphasis on Its New Features.

We report on a major update to the animal rDNA loci database, which now contains cytogenetic information for 45S and 5S rDNA loci in more than 2600 and 1000 species, respectively.The data analyses show the following: (i) A high variability in 5S and 45S loci numbers, with both showing 50-fold or higher variability. However, karyotypes with an extremely high number of loci were rare, and medians generally converged to two 5S sites and two 45S rDNA sites per diploid genome. No relationship was observed between the number of 5S and 45S loci. (ii) The position of 45S rDNA on sex chromosomes was relatively frequent in some groups, particularly in arthropods (14% of karyotypes). Furthermore, 45S rDNA was almost exclusively located in microchromosomes when these were present (in birds and reptiles). (iii) The proportion of active NORs (positively stained with silver staining methods) progressively decreased with an increasing number of 45S rDNA loci, and karyotypes with more than 12 loci showed, on average, less than 40% of active loci. In conclusion, the updated version of the database provides some new insights into the organization of rRNA genes in chromosomes. We expect that its updated content will be useful for taxonomists, comparative cytogeneticists, and evolutionary biologists. .

Asymmetrical canina meiosis is accompanied by the expansion of a pericentromeric satellite in non-recombining univalent chromosomes in the genus Rosa.

BACKGROUND AND AIMS: Despite their abundant odd-ploidy (2n = 5x = 35), dogroses (Rosa sect. Caninae) are capable of sexual reproduction due to their unique meiosis. During canina meiosis, two sets of chromosomes form bivalents and are transmitted by male and female gametes, whereas the remaining chromosomes form univalents and are exclusively transmitted by the egg cells. Thus, the evolution of chromosomes is expected to be driven by their behaviour during meiosis. METHODS: To gain insight into differential chromosome evolution, fluorescence in situ hybridization was conducted for mitotic and meiotic chromosomes in four dogroses (two subsections) using satellite and ribosomal DNA probes. By exploiting high-throughput sequencing data, we determined the abundance and diversity of the satellite repeats in the genus Rosa by analysing 20 pentaploid, tetraploid and diploid species in total. KEY RESULTS: A pericentromeric satellite repeat, CANR4, was found in all members of the genus Rosa, including the basal subgenera Hulthemia and Hesperhodos. The satellite was distributed across multiple chromosomes (5-20 sites per mitotic cell), and its genomic abundance was higher in pentaploid dogroses (2.3 %) than in non-dogrose species (1.3 %). In dogrose meiosis, univalent chromosomes were markedly enriched in CANR4 repeats based on both the number and the intensity of the signals compared to bivalent-forming chromosomes. Single-nucleotide polymorphisms and cluster analysis revealed high intragenomic homogeneity of the satellite in dogrose genomes. CONCLUSIONS: The CANR4 satellite arose early in the evolution of the genus Rosa. Its high content and extraordinary homogeneity in dogrose genomes is explained by its recent amplification in non-recombining chromosomes. We hypothesize that satellite DNA expansion may contribute to the divergence of univalent chromosomes in Rosa species with non-symmetrical meiosis.

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.

Evolutionary trends in animal ribosomal DNA loci: introduction to a new online database.

Ribosomal DNA (rDNA) loci encoding 5S and 45S (18S-5.8S-28S) rRNAs are important components of eukaryotic chromosomes. Here, we set up the animal rDNA database containing cytogenetic information about these loci in 1343 animal species (264 families) collected from 542 publications. The data are based on in situ hybridisation studies (both radioactive and fluorescent) carried out in major groups of vertebrates (fish, reptiles, amphibians, birds, and mammals) and invertebrates (mostly insects and mollusks). The database is accessible online at www.animalrdnadatabase.com . The median number of 45S and 5S sites was close to two per diploid chromosome set for both rDNAs despite large variation (1-74 for 5S and 1-54 for 45S sites). No significant correlation between the number of 5S and 45S rDNA loci was observed, suggesting that their distribution and amplification across the chromosomes follow independent evolutionary trajectories. Each group, irrespective of taxonomic classification, contained rDNA sites at any chromosome location. However, the distal and pericentromeric positions were the most prevalent (> 75% karyotypes) for 45S loci, while the position of 5S loci was more variable. We also examined potential relationships between molecular attributes of rDNA (homogenisation and expression) and cytogenetic parameters such as rDNA positions, chromosome number, and morphology.

Parental transposable element loads influence their dynamics in young Nicotiana hybrids and allotetraploids.

The genomic shock hypothesis suggests that allopolyploidy is associated with genome changes driven by transposable elements, as a response to imbalances between parental insertion loads. To explore this hypothesis, we compared three allotetraploids, Nicotiana arentsii, N. rustica and N. tabacum, which arose over comparable time frames from hybridisation between increasingly divergent diploid species. We used sequence-specific amplification polymorphism (SSAP) to compare the dynamics of six transposable elements in these allopolyploids, their diploid progenitors and in corresponding synthetic hybrids. We show that element-specific dynamics in young Nicotiana allopolyploids reflect their dynamics in diploid progenitors. Transposable element mobilisation is not concomitant with immediate genome merger, but occurs within the first generations of allopolyploid formation. In natural allopolyploids, such mobilisations correlate with imbalances in the repeat profile of the parental species, which increases with their genetic divergence. Other restructuring leading to locus loss is immediate, nonrandom and targeted at specific subgenomes, independently of cross orientation. The correlation between transposable element mobilisation in allopolyploids and quantitative imbalances in parental transposable element loads supports the genome shock hypothesis proposed by McClintock.

Remarkable variation of ribosomal DNA organization and copy number in gnetophytes, a distinct lineage of gymnosperms.

INTRODUCTION: Gnetophytes, comprising the genera Ephedra, Gnetum and Welwitschia, are an understudied, enigmatic lineage of gymnosperms with a controversial phylogenetic relationship to other seed plants. Here we examined the organization of ribosomal DNA (rDNA) across representative species. METHODS: We applied high-throughput sequencing approaches to isolate and reconstruct rDNA units and to determine their intragenomic homogeneity. In addition, fluorescent in situ hybridization and Southern blot hybridization techniques were used to reveal the chromosome and genomic organization of rDNA. KEY RESULTS: The 5S and 35S rRNA genes were separate (S-type) in Gnetum montanum, Gnetum gnemon and Welwitschia mirabilis and linked (L-type) in Ephedra altissima. There was considerable variability in 5S rDNA abundance, ranging from as few as ~4000 (W. mirabilis) to >100 000 (G. montanum) copies. A similar large variation was also observed in 5S rDNA locus numbers (two to 16 sites per diploid cell). 5S rRNA pseudogenes were interspersed between functional genes forming a single unit in E. altissima and G. montanum. Their copy number was comparable or even higher than that of functional 5S rRNA genes. In E. altissima internal transcribed spacers of 35S rDNA were long and intrinsically repetitive while in G. montanum and W. mirabilis they were short without the subrepeats. CONCLUSIONS: Gnetophytes are distinct from other gymnosperms and angiosperms as they display surprisingly large variability in rDNA organization and rDNA copy and locus numbers between genera, with no relationship between copy numbers and genome sizes apparent. Concerted evolution of 5S rDNA units seems to have led to the amplification of 5S pseudogenes in both G. montanum and E. altissima. Evolutionary patterns of rDNA show both gymnosperm and angiosperm features underlining the diversity of the group.

Cytogenetic features of rRNA genes across land plants: analysis of the Plant rDNA database.

The online resource http://www.plantrdnadatabase.com/ stores information on the number, chromosomal locations and structure of the 5S and 18S-5.8S-26S (35S) ribosomal DNAs (rDNA) in plants. This resource was exploited to study relationships between rDNA locus number, distribution, the occurrence of linked (L-type) and separated (S-type) 5S and 35S rDNA units, chromosome number, genome size and ploidy level. The analyses presented summarise current knowledge on rDNA locus numbers and distribution in plants. We analysed 2949 karyotypes, from 1791 species and 86 plant families, and performed ancestral character state reconstructions. The ancestral karyotype (2n = 16) has two terminal 35S sites and two interstitial 5S sites, while the median (2n = 24) presents four terminal 35S sites and three interstitial 5S sites. Whilst 86.57% of karyotypes show S-type organisation (ancestral condition), the L-type arrangement has arisen independently several times during plant evolution. A non-terminal position of 35S rDNA was found in about 25% of single-locus karyotypes, suggesting that terminal locations are not essential for functionality and expression. Single-locus karyotypes are very common, even in polyploids. In this regard, polyploidy is followed by subsequent locus loss. This results in a decrease in locus number per monoploid genome, forming part of the diploidisation process returning polyploids to a diploid-like state over time.

Organization and evolution of two repetitive sequences, 18-24J and 12-13P, in the genome of Chenopodium (Amaranthaceae).

The abundance and chromosomal organization of two repetitive sequences named 12-13P and 18-24J were analyzed in 24 diploid and nine polyploid species of Chenopodium s.l., with special attention to Chenopodium s.s. Both sequences were predominantly present in species of Chenopodium s.s.; however, differences in the amplification levels were observed among the species. The 12-13P repeat was highly amplified in all of the analyzed Eurasian species, whereas the American diploids showed a marked variation in the amplification levels. The 12-13P repeat contains a tandemly arranged 40 bp minisatellite element forming a large proportion of the genome of Chenopodium (up to 3.5%). FISH revealed its localization to the pericentromeric regions of the chromosomes. The chromosomal distribution of 12-13P delivered additional chromosomal marker for B-genome diploids. The 18-24J repeat showed a dispersed organization in all of the chromosomes of the analyzed diploid species and the Eurasian tetraploids. In the American allotetraploids (C. quinoa, C. berlandieri) and Eurasian allohexaploids (e.g., C. album) very intense hybridization signals of 18-24J were observed only on 18 chromosomes that belong to the B subgenome of these polyploids. Combined cytogenetic and molecular analyses suggests that reorganization of these two repeats accompanied the diversification and speciation of diploid (especially A genome) and polyploid species of Chenopodium s.s.

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.

Higher-order organisation of extremely amplified, potentially functional and massively methylated 5S rDNA in European pikes (Esox sp.).

BACKGROUND: Pikes represent an important genus (Esox) harbouring a pre-duplication karyotype (2n = 2x = 50) of economically important salmonid pseudopolyploids. Here, we have characterized the 5S ribosomal RNA genes (rDNA) in Esox lucius and its closely related E. cisalpinus using cytogenetic, molecular and genomic approaches. Intragenomic homogeneity and copy number estimation was carried out using Illumina reads. The higher-order structure of rDNA arrays was investigated by the analysis of long PacBio reads. Position of loci on chromosomes was determined by FISH. DNA methylation was analysed by methylation-sensitive restriction enzymes. RESULTS: The 5S rDNA loci occupy exclusively (peri)centromeric regions on 30-38 acrocentric chromosomes in both E. lucius and E. cisalpinus. The large number of loci is accompanied by extreme amplification of genes (>20,000 copies), which is to the best of our knowledge one of the highest copy number of rRNA genes in animals ever reported. Conserved secondary structures of predicted 5S rRNAs indicate that most of the amplified genes are potentially functional. Only few SNPs were found in genic regions indicating their high homogeneity while intergenic spacers were more heterogeneous and several families were identified. Analysis of 10-30 kb-long molecules sequenced by the PacBio technology (containing about 40% of total 5S rDNA) revealed that the vast majority (96%) of genes are organised in large several kilobase-long blocks. Dispersed genes or short tandems were less common (4%). The adjacent 5S blocks were directly linked, separated by intervening DNA and even inverted. The 5S units differing in the intergenic spacers formed both homogeneous and heterogeneous (mixed) blocks indicating variable degree of homogenisation between the loci. Both E. lucius and E. cisalpinus 5S rDNA was heavily methylated at CG dinucleotides. CONCLUSIONS: Extreme amplification of 5S rRNA genes in the Esox genome occurred in the absence of significant pseudogenisation suggesting its recent origin and/or intensive homogenisation processes. The dense methylation of units indicates that powerful epigenetic mechanisms have evolved in this group of fish to silence amplified genes. We discuss how the higher-order repeat structures impact on homogenisation of 5S rDNA in the genome.