Two classes of 5S rDNA unit arrays of the silver fir, Abies alba Mill.: structure, localization and evolution

Contrasting patterns of 5S rDNA repeats in European and Asian ecotypes of greater duckweed, Spirodela polyrhiza (Lemnaceae)

Ribosomal DNA (rDNA) contains highly conserved, specifically organized sequences encoding ribosomal RNAs (rRNAs) separated by variable non-transcribed intergenic spacers (NTSs) and is abundant in eukaryotic genomes. These characteristics make the rDNA an informative molecular target to study genome organization, molecular evolution, and phylogenetics. In this study, we characterized the 5S rDNA repeats in the greater duckweed Spiroldela polyrhiza, a species known for its small size, rapid growth, highly conserved genome organization, and low mutation rate. Sequence analysis of at least 12 individually cloned PCR fragments containing the 5S rDNA units for each of six ecotypes that originated from Europe (Ukraine) and Asia (China) revealed two distinct types of 5S rDNA repeats containing NTSs of different lengths and nucleotide compositions. The shorter 5S rDNA repeat units had a highly homogeneous 400-bp NTS, with few ecotype- or region-specific single-nucleotide polymorphisms (SNPs). The longer 5S rDNA units had NTSs of 1056-1084 bp with characteristic intra- and inter-genomic variants due to specific SNPs and insertions/deletions of 4-15-bp DNA elements. We also detected significant variability in the ratio of short/long 5S rDNA variants between ecotypes of S. polyrhiza. The contrasting dynamics of the two types of 5S rDNA units, combined with the unusually low repeat copy number (for plants) in S. polyrhiza (46-220 copies per genome), shows that this species could serve as an excellent model for examining the mechanisms of concerted evolution and functional significance of rDNA variability.

The Use of Ribosomal DNA for Comparative Cytogenetics

Fluorescence in situ hybridization (FISH) with ribosomal DNA (rDNA) sequences provides excellent chromosome markers for comparative cytogenetic analyses, especially in non-model plant species. The tandem repeat nature of a sequence and the presence of a highly conserved genic region make rDNA sequences relatively easy to isolate and clone. In this chapter, we describe the use of rDNA as markers for comparative cytogenetics studies. Traditionally, cloned probes labeled with Nick-translation have been used to detect rDNA loci. Recently, pre-labeled oligonucleotides are also employed quite frequently to detect both 35S and 5S rDNA loci. Ribosomal DNA sequences, together with other DNA probes in FISH/GISH or with fluorochromes such as CMA₃ banding or silver staining, are very useful tools in comparative analyses of plant karyotypes.

Molecular characterization of the 5S rDNA non-transcribed spacer and reconstruction of phylogenetic relationships in Capsicum

Abstract Capsicum includes ca. 41 species of chili peppers. In this original report we PCR amplified, cloned, sequenced and characterized the 5S rDNA non-transcribed spacer -NTS- in 23 taxa of nine clades of Capsicum, divergent at geographical origin and fruit and chromosome traits, and compared the NTS features throughout Solanaceae. According to GC content, inner variability and regulatory elements, the NTS organizes into three distinct structural regions; genetic variability at the NTS in Capsicum and related genus clusters into defined taxa hierarchies. Based on the reconstruction of a maximum-likelihood phylogenetic tree and phylogenetic networks, NTS sequences of Capsicum and related taxa grouped into well recognized categories -genus, section, clade, species, variety-. An evolutionary scenario arose from combined genetic and phylogenetic NTS data, in which monophyly and lineage diversification over time of Capsicum are addressed. Our analysis is original to include all domesticated species of Capsicum prevailing in germplasm collections and breeding programs, together with a large group of wild taxa that demanded further genetic characterization. The NTS set up as a double purpose marker in Capsicum, to directly evaluate genetic variability and reconstruct phylogenetic relationships to a broad extent, and constitutes a valuable tool for germplasm characterization and evolutionary studies within Solanaceae.

Molecular Evolution and Organization of Ribosomal DNA in the Hawkweed Tribe Hieraciinae (Cichorieae, Asteraceae)

Molecular evolution of ribosomal DNA can be highly dynamic. Hundreds to thousands of copies in the genome are subject to concerted evolution, which homogenizes sequence variants to different degrees. If well homogenized, sequences are suitable for phylogeny reconstruction; if not, sequence polymorphism has to be handled appropriately. Here we investigate non-coding rDNA sequences (ITS/ETS, 5S-NTS) along with the chromosomal organization of their respective loci (45S and 5S rDNA) in diploids of the Hieraciinae. The subtribe consists of genera Hieracium, Pilosella, Andryala, and Hispidella and has a complex evolutionary history characterized by ancient intergeneric hybridization, allele sharing among species, and incomplete lineage sorting. Direct or cloned Sanger sequences and phased alleles derived from Illumina genome sequencing were subjected to phylogenetic analyses. Patterns of homogenization and tree topologies based on the three regions were compared. In contrast to most other plant groups, 5S-NTS sequences were generally better homogenized than ITS and ETS sequences. A novel case of ancient intergeneric hybridization between Hispidella and Hieracium was inferred, and some further incongruences between the trees were found, suggesting independent evolution of these regions. In some species, homogenization of ITS/ETS and 5S-NTS sequences proceeded in different directions although the 5S rDNA locus always occurred on the same chromosome with one 45S rDNA locus. The ancestral rDNA organization in the Hieraciinae comprised 4 loci of 45S rDNA in terminal positions and 2 loci of 5S rDNA in interstitial positions per diploid genome. In Hieracium, some deviations from this general pattern were found (3, 6, or 7 loci of 45S rDNA; three loci of 5S rDNA). Some of these deviations concerned intraspecific variation, and most of them occurred at the tips of the tree or independently in different lineages. This indicates that the organization of rDNA loci is more dynamic than the evolution of sequences contained in them and that locus number is therefore largely unsuitable to inform about species relationships in Hieracium. No consistent differences in the degree of sequence homogenization and the number of 45S rDNA loci were found, suggesting interlocus concerted evolution.

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.

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.

Non-Random Distribution of 5S rDNA Sites and Its Association with 45S rDNA in Plant Chromosomes

5S and 45S rDNA sites are the best mapped chromosome regions in eukaryotic chromosomes. In this work, a database was built gathering information about the position and number of 5S rDNA sites in 784 plant species, aiming to identify patterns of distribution along the chromosomes and its correlation with the position of 45S rDNA sites. Data revealed that in most karyotypes (54.5%, including polyploids) two 5S rDNA sites (a single pair) are present, with 58.7% of all sites occurring in the short arm, mainly in the proximal region. In karyotypes of angiosperms with only 1 pair of sites (single sites) they are mostly found in the proximal region (52.0%), whereas in karyotypes with multiple sites the location varies according to the average chromosome size. Karyotypes with multiple sites and small chromosomes (<3 µm) often display proximal sites, while medium-sized (between 3 and 6 µm) and large chromosomes (>6 µm) more commonly show terminal or interstitial sites. In species with holokinetic chromosomes, the modal value of sites per karyotype was also 2, but they were found mainly in a terminal position. Adjacent 5S and 45S rDNA sites were often found in the short arm, reflecting the preferential distribution of both sites in this arm. The high frequency of genera with at least 1 species with adjacent 5S and 45S sites reveals that this association appeared several times during angiosperm evolution, but it has been maintained only rarely as the dominant array in plant genera.

Triparental origin of triploid onion, Allium × cornutum (Clementi ex Visiani, 1842), as evidenced by molecular, phylogenetic and cytogenetic analyses

BACKGROUND

Reconstruction of the parental origins of cultivated plants from wild relatives, especially after long periods of domestication, is not a trivial task. However, recent advances in molecular phylogenetics, among other approaches, have proved to be very informative in analyses of the origin and evolution of polyploid genomes. An established minor garden crop, triploid onion Allium × cornutum (Clementi ex Visiani, 1842) (2n = 3x = 24), is widespread in southeastern Asia and Europe. Our previous cytogenetic analyses confirmed its highly heterozygous karyotype and indicated its possible complex triparental genome origin. Allium cepa L. and Allium roylei Stearn were suggested as two putative parental species of A. × cornutum, whereas the third parental species remained hitherto unknown.

RESULTS

Here we report the phylogenetic analyses of the internal transcribed spacers ITS1-5.8S-ITS2 of 35S rDNA and the non-transcribed spacer (NTS) region of 5S rDNA of A. × cornutum and its relatives of the section Cepa. Both ITS and NTS sequence data revealed intra-individual variation in triploid onion, and these data clustered into the three main clades, each with high sequence homology to one of three other species of section Cepa: A. cepa, A. roylei, and unexpectedly, the wild Asian species Allium pskemense B. Fedtsh. Allium pskemense is therefore inferred to be the third, so far unknown, putative parental species of triploid onion Allium × cornutum. The 35S and 5S rRNA genes were found to be localised on somatic chromosomes of A. × cornutum and its putative parental species by double fluorescent in situ hybridisation (FISH). The localisation of 35S and 5S rDNA in A. × cornutum chromosomes corresponded to their respective positions in the three putative parental species, A. cepa, A. pskemense, and A. roylei. GISH (genomic in situ hybridisation) using DNA of the three putative parental diploids corroborated the results of the phylogenetic study.

CONCLUSIONS

The combined molecular, phylogenetic and cytogenetic data obtained in this study provided evidence for a unique triparental origin of triploid onion A. × cornutum with three putative parental species, A. cepa, A. pskemense, and A. roylei.

Dancing together and separate again: gymnosperms exhibit frequent changes of fundamental 5S and 35S rRNA gene (rDNA) organisation

In higher eukaryotes, the 5S rRNA genes occur in tandem units and are arranged either separately (S-type arrangement) or linked to other repeated genes, in most cases to rDNA locus encoding 18S–5.8S–26S genes (L-type arrangement). Here we used Southern blot hybridisation, PCR and sequencing approaches to analyse genomic organisation of rRNA genes in all large gymnosperm groups, including Coniferales, Ginkgoales, Gnetales and Cycadales. The data are provided for 27 species (21 genera). The 5S units linked to the 35S rDNA units occur in some but not all Gnetales, Coniferales and in Ginkgo (∼30% of the species analysed), while the remaining exhibit separate organisation. The linked 5S rRNA genes may occur as single-copy insertions or as short tandems embedded in the 26S–18S rDNA intergenic spacer (IGS). The 5S transcript may be encoded by the same (Ginkgo, Ephedra) or opposite (Podocarpus) DNA strand as the 18S–5.8S–26S genes. In addition, pseudogenised 5S copies were also found in some IGS types. Both L- and S-type units have been largely homogenised across the genomes. Phylogenetic relationships based on the comparison of 5S coding sequences suggest that the 5S genes independently inserted IGS at least three times in the course of gymnosperm evolution. Frequent transpositions and rearrangements of basic units indicate relatively relaxed selection pressures imposed on genomic organisation of 5S genes in plants.

Evolution of the tetraploid Anemone multifida (2n = 32) and hexaploid A. baldensis (2n = 48) (Ranunculaceae) was accompanied by rDNA loci loss and intergenomic translocation: evidence for their common genome origin

BACKGROUND AND AIMS

In the genus Anemone two small groups of taxa occur with the highest ploidy levels 2n = 6x = 48, belonging to the closely related clades: the montane/alpine Baldensis clade and the more temperate Multifida clade. To understand the formation of polyploids within these groups, the evolution of allohexaploid A. baldensis (AABBDD, 2n = 6x = 48) from Europe and allotetraploid Anemone multifida (BBDD, 2n = 4x = 32) from America was analysed.

METHODS

Internal transcribed spacer and non-transcribed spacer sequences were used as molecular markers for phylogenetic analyses. Cytogenetic studies, including genomic in situ hybridization with genomic DNA of potential parental species as probe, fluorescence in situ hybridization with 5S and 18S rDNA as probes and 18S rDNA restriction analyses, were used to identify the parental origin of chromosomes and to study genomic changes following polyploidization.

KEY RESULTS

This study shows that A. multifida (BBDD, 2n= 4x = 32) and A. baldensis (AABBDD, 2n = 6x = 48) are allopolyploids originating from the crosses of diploid members of the Multifida (donor of the A and B subgenomes) and Baldensis groups (donor of the D subgenome). The A and B subgenomes are closely related to the genomes of A. sylvestris, A. virginiana and A. cylindrica, indicating that these species or their progeny might be the ancestral donors of the B subgenome of A. multifida and A and B subgenomes of A. baldensis. Both polyploids have undergone genomic changes such as interchromosomal translocation affecting B and D subgenomes and changes at rDNA sites. Anemone multifida has lost the 35S rDNA loci characteristic of the maternal donor (B subgenome) and maintained only the rDNA loci of the paternal donor (D subgenome).

CONCLUSIONS

It is proposed that A. multifida and A. baldensis probably had a common ancestor and their evolution was facilitated by vegetation changes during the Quaternary, resulting in their present disjunctive distribution.

Cytogenetic and phylogenetic studies of diploid and polyploid members of tribe Anemoninae (Ranunculaceae)

The ancestry, phylogenetic differentiation and systematic classification of the worldwide-distributed genus Anemone have been debated for many years. In this paper 11 Anemone, three Pulsatilla species and Hepatica nobilis were subjected to detailed karyotype analysis with the aim of obtaining new cytogenetic data that will contribute to karyotype evolutionary studies of the tribe Anemoninae. The results are interpreted in a phylogenetic context, established from the intergenic nontranscribed spacer (NTS) of 5S rDNA and internal transcribed spacer (ITS) of 35S rDNA. One to three 35S and one to three 5S rDNA loci are present in diploid and polyploid taxa. The 35S rDNA loci are located terminally on the short arm of acrocentric chromosomes, while for 5S rDNA there is no preferential chromosomal position as it exhibits terminal, subterminal, interstitial or pericentromeric positions, and is located either on acrocentric or metacentric chromosomes. The karyotype of hexaploid A. baldensis (2n = 6x = 48) is presented for the first time, and A. sylvestris is proposed as one of its putative parental species. Chromosome fusion/translocation is proposed as the key mechanism involved in reduction of the basic chromosome number from 8 in the Anemone subgenus to 7 in the Anemonidium subgenus. The cytogenetic data obtained are mainly supported by ITS and NTS phylogeny. Diversification of the genus Anemone was accompanied by a large reduction of heterochromatin, from the Mediterranean anemones that have large amounts of heterochromatin to the New World anemones without any detectable heterochromatic blocks.

Two different size classes of 5S rDNA units coexisting in the same tandem array in the razor clam Ensis macha: is this region suitable for phylogeographic studies?

For a study of 5S ribosomal genes (rDNA) in the razor clam Ensis macha, the 5S rDNA region was amplified and sequenced. Two variants, so-called type I or short repeat (approximately 430 bp) and type II or long repeat (approximately 735 bp), appeared to be the main components of the 5S rDNA of this species. Their spacers differed markedly, both in length and nucleotide composition. The organization of the two variants was investigated by amplifying the genomic DNA with primers based on the sequence of the type I and type II spacers. PCR amplification products with primers EMLbF and EMSbR showed that the long and short repeats are associated within the same tandem array, suggesting an intermixed arrangement of both spacers. Nevertheless, amplifications carried out with inverse primers EMSinvF/R and EMLinvF/R revealed that some short and long repeats are contiguous in the same tandem array. This is the first report of the coexistence of two variable spacers in the same tandem array in bivalve mollusks.

Molecular organization of 5S rDNA in sharks of the genus Rhizoprionodon: insights into the evolutionary dynamics of 5S rDNA in vertebrate genomes

In this study, we attempted a molecular characterization of the 5S rDNA in two closely related species of carcharhiniform sharks, Rhizoprionodon lalandii and Rhizoprionodon porosus, as well as a further comparative analysis of available data on lampreys, several fish groups and other vertebrates. Our data show that Rhizoprionodon sharks carry two 5S rDNA classes in their genomes: a short repeat class (termed class I) composed of approximately 185 bp repeats, and a large repeat class (termed class II) arrayed in approximately 465 bp units. These classes were differentiated by several base substitutions in the 5S coding region and by completely distinct non-transcribed spacers (NTS). In class II, both species showed a similar composition for both the gene coding region and the NTS region. In contrast, class I varied extensively both within and between the two shark species. A comparative analysis of 5S rRNA gene sequences of elasmobranchs and other vertebrates showed that class I is closely related to the bony fishes, whereas the class II gene formed a separate cartilaginous clade. The presence of two variant classes of 5S rDNA in sharks likely maintains the tendency for dual ribosomal classes observed in other fish species. The present data regarding the 5S rDNA organization provide insights into the dynamics and evolution of this multigene family in the fish genome, and they may also be useful in clarifying aspects of vertebrate genome evolution.

Molecular and cytological characterization of 5S rDNA inOryzaspecies: genomic organization and phylogenetic implications

We investigated the molecular characteristics and chromosomal organization of 5S rDNA in the genus Oryza , including diploid and tetraploid species. A phylogenetic tree of Oryza species was constructed based on the non-transcribed spacer sequences of 5S rDNA, and some novel relationships were discovered. Specifically, comparative sequence analysis of 5S rDNA in several wild rice species showed unique characteristics inconsistent with the model of concerted evolution: (1) multiple distinct 5S rDNA types were detected within a species, leading to intraspecific divergence of 5S rDNA; (2) multiple identical 5S rDNA types were shared among species, resulting in interspecies clustering of 5S rDNA types; and (3) intraspecific nucleotide diversity was detected within a 5S rDNA class. Our results obtained by fluorescence in situ hybridization revealed that each rice species studied contained only one 5S rDNA locus with two hybridization sites, which were located on either chromosome 7 or chromosome 11. These results suggest that different 5S rDNA classes within the rice genome were arranged together and that one pair of 5S rDNA loci from a diploid progenitor of the tetraploid species might have been lost during evolution. Taken together, our data show that 5S rDNA in rice species is more informative at the gene level than at the chromosome level.

Cytogenetic and molecular characterization of the Abies alba genome and its relationship with other members of the Pinaceae

Genome size, karyotype structure, heterochromatin distribution, position and number of ribosomal genes, as well as the ITS2 sequence of the internal transcribed spacer (ITS) were analysed in silver fir (Abies alba Mill.). The analysis also included characterization of the Arabidopsis-type of telomeric repeats in silver fir and in related species. The results were compared with results from other species of the Pinaceae, to evaluate phylogeny and chromosomal and molecular evolution in the Pinaceae. Integrated chromosomal data provided insights into chromosome and karyotype evolution in the Pinaceae. The evolutionary trend for GC-rich heterochromatic blocks seems to involve loss of blocks that are not associated with rDNA. Similarly, numerous large blocks of interstitial plant telomeric repeats that are typical for all analysed species of the genus Pinus were not observed in the evolutionarily younger genera, such as Abies, Picea and Larix. On the contrary, the majority of telomeric sequences in these three genera appeared confined to the chromosome ends. We confirmed the current position of Abies and Tsuga in subfamily Abietoideae and the position of Pinus in the subfamily Pinoideae based on ITS2 sequences. Pseudotsuga is placed together with Larix into the subfamily Laricoideae. We conclude that the current position of the genus Picea in the subfamily Abietoideae should be reconsidered and, possibly, the genus Picea should be reclassified as a separate subfamily, Piceoideae, as recently proposed.

Investigations of 5S rDNA of Vitis vinifera L.: sequence analysis and physical mapping

Here we report the first results of a study of 5S rDNA of Vitis vinifera . 5S rDNA sequences from seven genotypes were amplified by PCR, cloned, and sequenced. Three types of repeats were found. Two variants, denominated long repeat and short repeat, appeared to be the main components of the 5S rDNA of this species, since they were found in all genotypes analyzed. They differed markedly from each other in both the length and the nucleotide composition of the spacers. The third variant, classified as DEL short repeat, differs from the short repeat owing to a large deletion in the spacer region. It appears to be the most recent repeat type, since it was identified in only one genotype. The organization of the 5S rDNA repeat unit variants was investigated by amplifying the genomic DNA with primers designed on the sequence of the long and short spacers. The PCR-amplified fragments showed that the long repeat is associated with the other two repeats, indicating that in V. vinifera different repeat units coexist within the same tandem array. FISH analysis demonstrated that 5S rRNA genes are localized at a single locus. The variability of 5S rDNA repeats is discussed in relation to the putative allopolyploid origin of V. vinifera.