Loss of 5S rDNA units in the evolution of Agropyron, Pseudoroegneria, and Douglasdeweya

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.

Identification of 5S and 45S rDNA sites in Chrysanthemum species by using oligonucleotide fluorescence in situ hybridization (Oligo-FISH)

Fluorescence in situ hybridization (FISH) is a conventional method used to visualize the distribution of DNA elements within a genome. To examine the relationships within the Chrysanthemum genus, ribosomal DNA (rDNA), a popular cytogenetic marker, was utilized as a probe for FISH within this genus. Based on the genome data of Chrysanthemum nankingense, C. seticuspe and its allied genera in the Compositae(Asteraceae), we explored rDNA sequences to design oligonucleotide probes and perform oligonucleotide fluorescence in situ hybridization (Oligo-FISH) in eight Chrysanthemum accessions. The results showed that the majority of 5S rDNA signals were located in subterminal chromosome regions and that the number of 5S rDNA sites might be tightly associated with ploidy. For 45S rDNA sites, the number and intensity of signals differed from those of previously investigated Chrysanthemum resources. These findings may provide an optimally reliable method of examining the chromosome composition and structural variation of Chrysanthemum and its related species and allow researchers to understand the evolutionary history and phylogenetic relationships of Chrysanthemum.

The Utility of Graph Clustering of 5S Ribosomal DNA Homoeologs in Plant Allopolyploids, Homoploid Hybrids, and Cryptic Introgressants

INTRODUCTION

Ribosomal DNA (rDNA) loci have been widely used for identification of allopolyploids and hybrids, although few of these studies employed high-throughput sequencing data. Here we use graph clustering implemented in the RepeatExplorer (RE) pipeline to analyze homoeologous 5S rDNA arrays at the genomic level searching for hybridogenic origin of species. Data were obtained from more than 80 plant species, including several well-defined allopolyploids and homoploid hybrids of different evolutionary ages and from widely dispersed taxonomic groups.

RESULTS

(i) Diploids show simple circular-shaped graphs of their 5S rDNA clusters. In contrast, most allopolyploids and other interspecific hybrids exhibit more complex graphs composed of two or more interconnected loops representing intergenic spacers (IGS). (ii) There was a relationship between graph complexity and locus numbers. (iii) The sequences and lengths of the 5S rDNA units reconstituted in silico from k-mers were congruent with those experimentally determined. (iv) Three-genomic comparative cluster analysis of reads from allopolyploids and progenitor diploids allowed identification of homoeologous 5S rRNA gene families even in relatively ancient (c. 1 Myr) Gossypium and Brachypodium allopolyploids which already exhibit uniparental partial loss of rDNA repeats. (v) Finally, species harboring introgressed genomes exhibit exceptionally complex graph structures.

CONCLUSION

We found that the cluster graph shapes and graph parameters (k-mer coverage scores and connected component index) well-reflect the organization and intragenomic homogeneity of 5S rDNA repeats. We propose that the analysis of 5S rDNA cluster graphs computed by the RE pipeline together with the cytogenetic analysis might be a reliable approach for the determination of the hybrid or allopolyploid plant species parentage and may also be useful for detecting historical introgression events.

Diversity within the genus Elymus (Poaceae: Triticeae) as investigated by the analysis of the nr5S rDNA variation in species with St and H haplomes

The genus Elymus ("Ryegrass") is a repository for a range of species with a variety of haplome contents; hence the pejorative name "dustbin" genus. We have analyzed 1,059 sequences from 128 accessions representing 24 species to investigate the relationships among the StH haplomes-containing species described by Yen and Yang (Genus Elymus Beijing 5:58-362, 2013). Sequences were assigned to "unit classes" of orthologous sequences and subjected to a suite of analyses including BLAST (Basic Local Alignment Search Tool) searches, phylogenetic analysis and population genetic analysis to estimate species diversity. Our results support the genome analyses in Yen and Yang (Genus Elymus Beijing 5:58-362, 2013), i.e., genomic constitution StStHH including variants restricted to Elymus. Population genetic analysis of the 5S nrDNA sequence data revealed that the within-species variance component is roughly ±89 %; thus, we were unable to identify molecular markers capable to separate the 24 species analyzed. Separate phylogenetic analyses of the two unit classes and of all the data exhibit a trend only of the species to cluster on the phylograms. Finally, the analysis provides evidence for the multiple origins of American and Eurasian species.

What does the nr5S DNA multigene family tell us about the genomic relationship between Dasypyrum breviaristatum and D. villosum (Triticeae: Poaceae)?

The genus Dasypyrum contains two species: the annual and widespread D. villosum (2x = 2n = 14) and the perennial and generally rare D. breviaristatum (2x = 2n = 14 and 4x = 2n = 28). The origin of the latter and its genome constitution have been subject of several studies. There is agreement that the genome of the diploid D. villosum (VV) is different from the diploid cytotype of D. breviaristatum (VbVb), but there is no agreement of the constitution of the tetraploid cytotype, specifically whether is it an autotetraploid or an allotetraploid. This is a long-standing disagreement that this study aims to resolve using the 5S nrDNA as a genomic marker. Our studies suggest that the 4x D. breviaristatum is an allotetraploid (VVVbVb).

What does the 5S rRNA multigene family tell us about the origin of the annual Triticeae (Poaceae)?

We have investigated the complex relationships among the annual genera within the tribe Triticeae through phylogenetic analyses of the 5S rRNA multigene family. Cloned sequences were assigned to groups of orthologous sequences, called unit classes, that were subjected to several analyses including BLAST (Basic Local Alignment Search Tool) searches to assess possible ancestral relationships with perennial genera; phylogenetic analyses using parsimony (Pars), maximum likelihood (ML), and Bayesian methods; and minimum reticulation networks from the Pars, ML, and Bayesian trees. In this study, we included genera with both annual and perennial species, such as Dasypyrum, Hordeum, and Secale. BLAST pointed to Pseudoroegneria (carrier of the St genome) and possibly Thinopyrum (carrier of the J genome) as the potential next of kin. However, Thinopyrum and Pseudoroegneria have never fallen together on the individual trees with the former generally associated with Crithopsis, Aegilops, Triticum, and Dasypyrum, while the latter is usually associated with the rest of the genera within Triticeae. The "long" unit classes placed Dasypyrum breviaristatum together with Dasypyrum villosum, whereas the "short" unit classes put them far apart on the trees. None of the gene trees alone was able to summarize the complex relationships among the genera, in line with previous results in the Triticeae. However, the application of tools designed to display phylogenetic networks was able to depict the complex links among the genera based on the short and the long gene trees, including the close link between Thinopyrum and Pseudoroegneria suggested by the phylogenetic analyses. In addition, our analyses provide support for the hypothesis that at least some annual Triticeae taxa are derived from their perennial relatives.

Phylogenetic relationships among the polyploid and diploid Aegilops species inferred from the nuclear 5S rDNA sequences (Poaceae: Triticeae)

Phylogenetic inferences of the polyploid Aegilops taxa were drawn based upon the analysis of 909 nuclear 5S rDNA sequences obtained from 15 Aegilops polyploid taxa (531 sequences new to this paper) and 378 sequences from our previous study on the diploid taxa. The 531 sequences can be split into two orthologous groups (unit classes), the long AE1 and short AE1 previously identified in the diploid set. An examination of the relationships between unit classes and their associated haplomes suggests that U haplome sequences found in Ae. umbellulata are the closest to the T sequences found in Amblyopyrum muticum and that sequences of the polyploid species expected to be the M type found in Ae. comos are more similar to the T haplome sequences, except in the three hexaploids Ae. glumiaristata, Ae. juvenalis, and Ae. vavilovii and the tetraploid Ae. crassa where they are found to be similar to the M haplome sequences. These three hexaploid taxa likely originated from the tetraploid Ae. crassa (DM), while the closest taxon to the fourth hexaploid, Ae. recta, is the tetraploid Ae. neglecta (UM). Based upon the distribution of the unit classes, several reticulate phylogenies depicting evolutionary relationships among diploid, tetraploid, and hexaploid taxa were constructed; however, none of these widely used methods could depict the expected reticulate relationship as previously drawn from cytogenetic analyses in this group of allopolyploid species. These results suggest that evolutionary relationships derived from models based upon the assumption of bifurcating species require careful interpretation when these same models are applied to species with reticulate evolution.

Codependence of repetitive sequence classes in genomes: phylogenetic analysis of 5S rDNA families in Hordeum (Triticeae: Poaceae)

To complete our study of the genus Hordeum and to elaborate a phylogeny of species based upon 5S rDNA sequences, we have cloned and sequenced PCR amplicons from seven American polyploid species to generate 164 new 5S rRNA gene sequences. These sequences were analysed along with the more than 2000 5S rDNA sequences previously generated from the majority of species in Hordeum to provide a comprehensive picture of the distribution (presence or absence) of 5S rDNA unit classes (orthologous groups) in this genus as well as insights into the phylogeny of Hordeum. Testing of substitution models for each unit class based upon the consensus sequences of all the taxa as well as for each unit class within the genus found that the general best fit was TPM3uf+G, from which a maximum-likelihood tree was calculated. A novel application of cophylogenetic analysis, where relationships among unit classes were treated as host-parasite interactions, depicted some significant pair links under tests of randomness indicative of nonrandom codivergence among several unit classes within the same taxon. The previous classification of four genomic groups is reflected in combinations of unit classes, and it is proposed that current taxa developed from ancient diploidized paleopolyploids and that some were subjected to gene loss, i.e., unit class loss. Finally, separate phylogenetic analyses performed for the tetraploid and hexaploid species were used to derive a working model describing the phylogeny of the polyploid taxa from their putative diploid ancestry.