Repeated reunions and splits feature the highly dynamic evolution of 5S and 35S ribosomal RNA genes (rDNA) in the Asteraceae family
© Garcia et al; licensee BioMed Central Ltd. 2010
Received: 24 May 2010
Accepted: 16 August 2010
Published: 16 August 2010
In flowering plants and animals the most common ribosomal RNA genes (rDNA) organisation is that in which 35S (encoding 18S-5.8S-26S rRNA) and 5S genes are physically separated occupying different chromosomal loci. However, recent observations established that both genes have been unified to a single 35S-5S unit in the genus Artemisia (Asteraceae), a genomic arrangement typical of primitive eukaryotes such as yeast, among others. Here we aim to reveal the origin, distribution and mechanisms leading to the linked organisation of rDNA in the Asteraceae by analysing unit structure (PCR, Southern blot, sequencing), gene copy number (quantitative PCR) and chromosomal position (FISH) of 5S and 35S rRNA genes in ~200 species representing the family diversity and other closely related groups.
Dominant linked rDNA genotype was found within three large groups in subfamily Asteroideae: tribe Anthemideae (93% of the studied cases), tribe Gnaphalieae (100%) and in the "Heliantheae alliance" (23%). The remaining five tribes of the Asteroideae displayed canonical non linked arrangement of rDNA, as did the other groups in the Asteraceae. Nevertheless, low copy linked genes were identified among several species that amplified unlinked units. The conserved position of functional 5S insertions downstream from the 26S gene suggests a unique, perhaps retrotransposon-mediated integration event at the base of subfamily Asteroideae. Further evolution likely involved divergence of 26S-5S intergenic spacers, amplification and homogenisation of units across the chromosomes and concomitant elimination of unlinked arrays. However, the opposite trend, from linked towards unlinked arrangement was also surmised in few species indicating possible reversibility of these processes.
Our results indicate that nearly 25% of Asteraceae species may have evolved unusual linked arrangement of rRNA genes. Thus, in plants, fundamental changes in intrinsic structure of rDNA units, their copy number and chromosomal organisation may occur within relatively short evolutionary time. We hypothesize that the 5S gene integration within the 35S unit might have repeatedly occurred during plant evolution, and probably once in Asteraceae.
It is generally considered that in prokaryotes and in some species of early diverging eukaryote groups e.g. yeast, Saccharomyces cerevisiae, rRNA genes are organised in a single operon, clustered in tandem and transcribed by the same RNA polymerase. In eukaryotes, the 35S (encoding 18S-5.8S-26S rRNA) and 5S genes are transcribed by different polymerases, RNA polymerase I and III, respectively. Independent control of transcription probably enabled physical separation of both loci in chromosomes, an arrangement that is typical for most eukaryotic organisms. Nevertheless, it seems that there are several exceptions to this rule. For example, 5S genes linkage to other repetitive sequences including 35S, histone genes or the trans-spliced leader has been demonstrated . These linked arrangements are found among diverse biological taxa including nematodes , fungi , crustaceans , slime moulds  or mosses [6, 7], and they are believed to represent transition states between linked (prokaryotic) and unlinked (eukaryotic) arrangements. However, our recent observations in a group of angiosperms (genus Artemisia, Asteraceae)  clearly point to the possibility that linked arrangements might not be restricted to prokaryotes and primitive eukaryotes but may occur throughout the tree of life.
The Asteraceae, also named Compositae, is the largest family of the angiosperms in terms of numbers of species, with 1,620 genera and 23,600 species , constituting approximately 8-10% of all flowering plants. Although many trees and shrubs exist, a majority of species of Asteraceae are herbaceous, and they are easily recognizable by a suite of characters including fused anthers, a fruit with a single ovule, and a specialized inflorescence termed the capitulum. The family occurs in all continents of the world except Antarctica  and includes many edible, medicinal, noxious, invasive or endangered species. Estimates indicate that the Asteraceae originated in the mid Eocene (45-49 Mya) and that most tribal splits occurred during the Oligocene (28-36 Mya). The family has been the subject of intensive phylogenetic analyses using all kinds of molecular and morphological data [9, 11] and it is currently considered to hold 12 major lineages , four of which (Asteroideae, Carduoideae, Cichorioideae and Mutisioideae) comprise 99% of its species diversity. Considering the size and importance of this family, however, relatively little has been published about rDNA position, organisation and structure in this group as a whole. In the past, the 5S and 35S genes have been mapped on chromosomes in some Asteraceae genera, e.g. Tragopogon , Centaurea , Helianthus , and Hypochaeris [16, 17] among others, showing one or several separate loci of each, mostly in different chromosomes. In contrast, initial FISH mapping on Artemisia showed colocalised 35S and 5S signals [18–20]. Employing molecular methods we have recently demonstrated that several Artemisia species have evolved unusual ribosomal units resembling the arrangement in yeast: the 5S genes locate within the 26S-18S intergenic spacer and are transcribed from the opposite strand of the 35S rDNA operon . In species displaying linked arrangement the homogenisation of linked rRNA genes most likely went to completion, since no unlinked genes were detected by FISH or Southern blot hybridisation. Taking into account the large size of the family it is not known how frequently such arrangement occurs, whether it is conserved in related species, its chromosomal dynamics and evolutionary success. With these issues in mind, we examined ribosomal DNA structure and organisation in selected representatives of the Asteraceae using molecular and cytogenetic methods. In addition, we searched for any nodes in extant Asteraceae phylogeny [21, 22] that mark a switch in the organisation of rDNA.
Leaf or seed material for 199 populations of the different species studied was obtained either from wild populations, botanical gardens/herbaria worldwide, or through plants grown at the greenhouses of the Institute of Biophysics (Brno) and of the Botanical Institute of Barcelona, from seeds obtained by Index Seminum or purchased. Additional file 1 summarizes the provenance of all the taxa investigated.
Southern blot hybridisation
Purified genomic DNAs and PCR products were digested with restriction enzymes and separated by gel electrophoresis on a 0.9% (w/v) agarose gel. After electrophoresis, the gels were alkali blotted onto Hybond-XL membranes (GE Healthcare, Little Chalfont, UK) and hybridised with 32P-labelled DNA probes (DekaLabel kit, MBI, Fermentas, Vilnius, Lithuania) as described in . After washing under high-stringency conditions, the hybridisation bands were visualized with a PhosphorImager (Storm, Molecular Dynamics, Sunnyvale, CA), and the data were processed by ImageQuant software (Molecular Dynamics). The probes were a 220-bp PCR product derived from the 3' end of the 26S rRNA gene of tobacco  and a cloned ~120 bp of the 5S genic region, also from tobacco .
Real time PCR
Quantification of the rRNA gene copies was done using the SYBR Green PCR Master Mix (SAbiosciences, Fredrick, MD, USA) by the 7300 Real Time PCR System (Applied Biosystems/Ambion, Austin, TX, USA). The 26S-5S unit was amplified with the Pr-1 (5'-AGACGACTTAAATACGCGAC-3') and 5SLf primers; the control P2P gene was amplified with the forward primer 5'-CTGGATTTGCTGGTGATGAT-3' and the reverse primer 5'-CYCTCTTGGATTGAGCTT-3'  in the same PCR cycle to get an estimate of DNA concentration (initial denaturation at 94°C for 7 min followed by 40 cycles of 20 s at 94°C, 20 s at 56°C and 30 s at 72°C; and a final extension step at 72°C, 10 min). The SYBR Green I fluorescence was monitored consecutively after the extension step. A standard curve was made of series dilutions of the pAA_IGS-12 plasmid containing a cloned IGS from A. absinthium [GenBank accession number: EU649669]. This curve was then used to estimate amounts of rDNA in each of the samples based on the cross point values. The number of genes were calculated according to the formula: N = m/M × C × 1/10 where m is the amount of rDNA in ng; M, amount of genomic DNA input in ng; C, genome size in kb; 1/10 is a constant related to average size of rDNA unit (10 kb). The genome size data were taken from [26–28] and Plant DNA C-values Database (http://data.kew.org/cvalues/).
To isolate intergenic spacer sequences we used Pr-1 and Pr-2 (5'-GGCTTAATCTTTGAGACAA-3') primers . The DyNAzyme™II DNA polymerase (Finnzymes, Espoo, Finland), was used in these PCRs; each reaction contained ~100 ng of genomic DNA. Following PCR, the products were gel-purified and the isolated fragments (QIAquick PCR purification kit, Hilden, Germany) were cloned into the pDrive vector (QIAGEN PCR Cloning Kit (Qiagen). Several positive clones were recovered from each transformation reaction and analysed for the lengths of insert by gel electrophoresis. Only clones containing inserts > 3 kb in size were further considered. The full length IGS constructs were obtained from Tagetes patula (pTPa_IGS, two clones), Matricaria matricarioides (pMMa_IGS, three clones), Gnaphalium luteoalbum (pGLu_IGS, three clones), Helichrysum cymosum (pHCy_IGS, two clones), Madia sativa (pMSa_IGS, two clones), Elachanthemum intricatum (pEIn_IGS, four clones) and Helianthus annuus (pHAn_IGS, two clones). The clones were sequenced both from the side of 26S and 18S genes (Eurofins, MWG Operon, Eberberg, Germany). The newly identified sequences were submitted to the EMBL/GenBank under the accession numbers [HM160523-29].
Sequences were assembled with the help of Bioedit software , manually edited and aligned. Phylogenetic relationships between clones were inferred from phylograms constructed based on the neighbor-joining (NJ) algorithm with the setting of the Kimura two-parameter substitution model. The statistical support for branching was obtained with a bootstrap analysis: 500 repetitions were carried out using the PHYLOWIN program . The distances between individual clones were calculated with the assistance of DnaSP 4.0 software .
Fluorescent in situ hybridisation
The fresh root tips were pretreated with an aqueous solution of colchicine 0.05% at room temperature, for 2.5 - 4 h and fixed in 3:1 (v/v) ethanol: acetic acid. Protoplasts were obtained using cellulolytic enzymes, dropped onto microscope slides, frozen and desiccated using liquid nitrogen and 70% ethanol. Before FISH, the slides were pre-treated with 50 μg mL-1 RNaseA for 1 h at 37°C in humid chamber. After washing three times in 2× SSC (2× standard saline citrate + 0.1% (w/v) sodium dodecyl sulfate), slides were dehydrated in an ethanol series and air-dried. Remnants of cytoplasm were removed with pepsin treatment (10 μg mL-1 in 10 mM HCl, 4 min room temperature). The slides were then washed, dehydrated in ethanol and fixed for 10 min in 3.7% formaldehyde in 1× PBS, washed three times in 2× SSC, dehydrated again and air dried. The hybridisation mix contained 50 ng μL-1 of Cy3-labelled (GE Healthcare, Chalfont, St Giles, UK) 5S probe, a 116 bp-long insert of the cloned tobacco 5S rRNA gene , and 20 ng μL-1 of 35S rDNA probe (a 2.5 kb fragment of 26S rRNA gene from tomato labelled with Spectrum Green, Abbott Molecular, IL, USA). The FISH hybridisation mixture (20 μL per slide) consisted of labelled DNA probes, 4 μL of a 50% solution of dextran sulphate, 10 μL formamide, 0.5 μL TE buffer and 2 μL 20 × SSC, and it was denatured at 75°C for 15 min and immediately cooled on ice. Slides were denatured in a thermocycler using a flat plate: 5 min at 75°C, 2 min at 65°C, 2 min at 55°C, 2 min at 45°C, and transferred into a pre-warmed moist chamber and put into an incubator. After overnight hybridisation at 37°C, the slides were washed with 2× SSC, 0.1× SSC (high stringency), at 42°C for 10 min each followed by washes with 2× SSC, 4× SSC + 0.1% Tween 20 at room temperature. Slides were rinsed in PBS and mounted in Vectashield (Vector Laboratories, Burlinghame, CA, USA) containing DAPI (1 μg/ml-1). FISH signals were observed using an Olympus AX 70 fluorescent microscope equipped with a digital camera. Images were analysed and processed using ISIS software (MetaSystems, Altlussheim, Germany).
We aimed to sample representative members of all subfamilies and most tribes within the Asteraceae (according to the recent phylogenetic revisions, see [22, 34]) In decreasing order, the largest groups within this family are: subfamily Asteroideae, 17,000 species (73% of the family) > subfamily Cichorioideae (12.4%) > subfamily Carduoideae (11%). Out of the selected 199 accessions, 72% belong to subfamily Asteroideae, 9% to Cichorioideae and 7% to Carduoideae. The remaining 12% of the sample is composed by the other 9 subfamilies that complete the Asteraceae and by a representation of the sister families Calyceraceae and Goodeniaceae, the closely related Campanulaceae and Menyanthaceae, and the Solanaceae as a control as done in . We consider that our sample selection is roughly reflecting distribution and relative abundance of species among the subfamilies and tribes. Multiple individuals per species were analysed and where possible different populations were scored. The complete list of taxa included in this study and details of specimens are given in Additional file 1.
Genomic organisation of rDNA analysed by PCR
DNAs producing amplification products with the 5SLf-5SLr primers but not with the 26Sf-5S primers (both combinations). These species were categorized as Group I species. Multiple regularly spaced bands in the 5S-5S profiles correspond to amplification of oligomeric units in tandems (mono-, di-, tri-, etc.). Thus Group I profile is consistent with a canonical separate unlinked tandem arrangement of 5S rRNA genes.
The DNA samples yielding strong 26Sf-5SLf PCR products and no 5SLf-5SLr products were categorized as Group II species. The patterns were interpreted as dominant linked 35S-5S arrangement in which the 5S genes are in an inverse orientation with respect to the direction of the 26S gene transcription.
There were cases of amplification products obtained with both 5SLf-5SLr and 26Sf-5SLf primer sets. This situation (Group III) might correspond to the presence of both linked and unlinked genes in the same genome, albeit in different proportions, as shown further.
Finally, few species such as Madia sativa (or Mutisia speciosa and also Tithonia rotundifolia, not shown) yielded amplification products with 26Sf-5SLr primer sets instead of the 26Sf-5SLf combination, that did not produce a product in these species), possibly indicating a direct orientation of the 26S-5S linkage at least in some units. These species (Group IV) also yielded positive PCR signals with the 5SLf-5SLr primer combination.
Genomic organisation of rDNA analysed by Southern blot hybridisation
The 5S probe hybridisation bands having no overlap with the 26S hybridisation signals. These ladder profiles are consistent with long tandems of unlinked 5S genes. The 26S probe usually hybridised to several high molecular weight bands indicating the presence multiple IGS families lacking 5S insertions. Contrast to the 5S profile the 26S hybridisation did not generate a regular ladder of fragments which is explained by the presence of multiple BamHI sites within the 35S unit (Figure 1A). These 5S and 26S hybridisation patterns were observed with all Group I, III and IV DNAs after the digestion with BamHI. Patterns do not support the presence of long arrays composed of linked units in Group III and IV species despite positive PCR signals with the 26S-5S primer sets.
The Group II plants showed extensive cohybridisation of 26S and 5S probes to the same restriction fragments, consistent with dominant linked arrangement of both genes in the genome. The sizes of major fragments (2.5-3.0 kb) differed slightly between the species probably reflecting sequence polymorphism in the spacer between the 26S and 5S genes (confirmed by sequencing - see further below). In Antennaria, Tagetes and Tanacetum the major fraction (2.5-3.0 kb) migrated as doublets or triplets, indicating the presence of multiple gene families. On the other hand Coreopsis, Evax, Matricaria and Pseudognaphalium, had single bands in this region showing a single dominant family. The minor higher molecular weight fragments (> 5 kb) resulted most likely from inhibition of enzyme digestion due to methylation or mutation. Few 5S-hybridising fragments (e.g. in Coreopsis, Matricaria, Pseudognaphalium and Tanacetum) did not hybridise with the 26S probe while they did hybridise to the IGS probe from Artemisia (not shown).
Structure of the 26S-18S intergenic spacers
Since 5S-homologous sequences have been found among the SINE3 retrotransposons in eukaryotic organisms , including plants , we searched the 5S flanking sequences for distinctive motifs typical of retroelements such as short repeats, long terminal repeats (LTR), transcription priming sites and duplications of 18S and 26S genes. All IGS analysed contained motifs resembling consensus sequences of terminal inverted repeats (TIR) of Cassandra retroelements (Figure 4B), which carry independently transcribed 5S rRNA genes . Their position was variable within <200 bp upstream or downstream from the 5S genic region (Figure 4A). In some cases homology extended to the LTRs. Duplications of a putative integration site as well as domains typical from autonomous retrotransposons (e.g. reverse transcriptase) were, however, not found.
The aligned 5S sequences comprising both genic and flanking regions are shown in Figure 5. In species that homogenised the 35S-5S units (Group II plants) the linked 5S copies contained all regulatory elements to execute transcription by polymerase III: the A and C boxes were found at conserved positions. While box A was invariant there was slight sequence variation in box C. The transcription starts (TSS) at G is preceded by a highly conserved C at -1. The TATA sequence appeared to be fairly conserved located about 24-30 bp upstream from the transcription start site. Most species had two oligo dT terminators in the 3' flanking region while Helichrysum cymosum sequence bears a single 16-17 bp-long T track immediately downstream from the genic region. In short, the species with homogenised linked arrangements bear apparently functional 5S genes in the 26S-18S intergenic spacer(s).
The sequence divergences along individual regions of units are given in Additional file 3. As expected, both 26S and 5S genic parts were highly conserved while the inter- and intra genic spacers were less conserved. There was a 10 fold difference in divergence levels between genic and non-genic regions. Thus, the 5S insertions evolve in concert with linked 18S-5.8S-26S genes.
Phylogenetic relationships between rRNA gene families
The phyogenetic tree constructed from the IGS1 sequences is shown in Figure 6B. The topology of major clades is reminiscent to that previously published [21, 37]. Well supported branches at the genus levels are of note.
Fluorescent in situ hybridisation
Figure 8 shows result of FISH analysis in representative Group I species with separate rDNA arrangement. Dahlia pinnata (Coreopsidae, from "Heliantheae alliance"), Helianthus annus (Heliantheae, also from "Heliantheae alliance"), Chrysanthemum zawadskii (Anthemideae), Aster alpinus (Astereae), Calendula officinalis (Calenduleae), and Tragopogon mirus (tribe Cichorieae, from subfamily Cichorioideae) present typical unlinked rDNA arrangement with 5S located interstitially and the 35S loci being terminal. One rDNA locus in Dahlia pinnata showed partial overlap of both 26S and 5S probes.
Summary of cytogenetic analysis by rDNA-FISH.
Linked2 35S-5S loci
Unlinked2 5S loci
Unlinked2 35S loci
Two contrasting modes of rDNA organisation in Asteraceae
The assorted situation here resembles that of primitive eukaryotes (yeast or other fungi) in which 5S genes display either linked arrangement with the other rDNAs, unlinked tandems, or even a dispersed organisation . The results obtained in this study are significant since they show variation in rDNA arrangements even between species assumed to diverge relatively recently, such as the mentioned case of the very closely related genera Elachanthemum and Artemisia. Based on molecular clock dating, it appears that the family is not older than 50 Myrs  and species within most tribes probably diverged recently . Together, linkage cannot be considered as an exclusive characteristic of the family, subfamily or tribe since significant rearrangements of rDNA probably occur at lower taxonomic levels such as the genus.
Mechanisms of rDNA rearrangements
The initial phase probably involved a retrotransposition event that inserted a 5S gene(s) into the IGS (Figure 10A). This assumption is supported by several retroelement features found in the IGS of many of the studied species: (i) imperfect short 4-7 bp terminal inverted repeats, resembling those of non-autonomous TRIM retroelements called Cassandra (carrying 5S-like sequences, ), were found in close proximity to the 5S insertions (Figure 4B). In some species (e.g. G. luteoalbum), the homology apparently included several additional base pairs of Cassandra LTR; (ii) motifs resembling potential primer binding sites (PBS) for reverse transcription (minus strand) occur not only in the IGS of some species  but, interestingly, also within the genic region (close to the 5'end) pointing to the retrotransposon origin of 5S genes, in general.
It also has to be mentioned that there are features that distinguish the 5S IGS insertions from Cassandra retroelements. First, the 5S insertions in IGS contain all structural features of functional 5S genes at least in species that homogenised linked units to completion. On the other hand the 5S-LTRs of the sequenced Cassandra clones are mutated and unlikely to encode functional RNA molecules. This is supported by the phylogenetic analysis that failed to distinguish between linked and unlinked cellular copies of 5S but clearly separated these from Cassandra sequences (Figure 6A). The second distinction is that while termination signals follow immediately after the 5S coding region in the genic 5S, absence of terminators in Cassandra sequences extends polymerization to LTR by >400 nt . It is unlikely that molecules of this size would be accurately processed and packed into ribosomes. Clearly, Cassandra elements, on their own, could not serve as 5S vehicles that carried functional 5S genes into the IGS. However, while there are thousands of similar copies in plant genomes we cannot exclude the possibility that there may be a yet undiscovered Cassandra locus bearing a functional 5S gene. Alternatively, a recombination event could have occurred between a retroelement and 5S tandems in a putative Asteroideae progenitor giving rise to a novel element. The presence of 5S-5S spacer sequences and duplication of 5S insertions in some IGS clones (Figure 4) could support the hypothesis.
Concerted evolution. FISH and quantitative PCR experiments showed that in many genomes concerted evolution has homogenised thousands of linked units indicating that insertion of a putative 5S-carrying retroelement into IGS does not prevent unit amplification. This may be consistent with situation in tomato , Drosophila  and nodding onion  in which units bearing retroelement insertions were observed in high copy. Classical theoretical models explained homogenisation of tandem repeats by multiple non-homologous recombination steps [ and references herein]. Computer simulations imply that homogenisation time of a tandemly arranged sequence is a quadratic function of copy numbers . It follows that 104 copies of a typical rDNA array homogenise within 108 years. However, the Asteroideae subfamily is probably not older than 26-39 million yrs [9, 10]. Many of its members with contrasting rDNA arrangements (e.g. Tagetes v. Helianthus or the previously mentioned Elachanthemum v. Artemisia) probably diverged more recently. In addition, some species, e.g. Tagetes patula, carry arrays of homogenised 35S-5S genes distributed across many chromosomes suggesting extensive interlocus recombination and/or recurrent polyploidy events. Thus, while the model based solely on non-homologous recombination between chromatids seems to work with yeast (mainly due to the lower number of rDNA copies, less than 200) and may elucidate the initial phase of 35S-5S repeat expansion (Figure 10B), an alternative model is needed to explain such rapid turnover of thousands of rRNA genes and their spreading across the chromosomes. Several mechanisms including rolling cycle replication and reintegration of covalently linked circles [46, 47] have been recently proposed to mediate homogenisation of satellite repeats in allopolyploids  and may be possibly applied to the rDNA switches as well (Figure 10C).
It is intriguing that we were unable to identify a species bearing equivalent proportion of linked and unlinked tandems, since FISH studies confirmed homogenisation across all rDNA loci (Figure 7 and Table 1). Either mixed genotypes are counterselected (intermediate stages, Figure 10B) may not be evolutionarily successful or the turnover (birth and extinction) of rDNA arrays in genomes is extremely fast so intermediate genotypes are short-lived and escape detection. The best evidence for the latter hypothesis is the organisation of rDNA in Elachanthemum intricatum, an Old World Anthemideae species. In this species, the functional rRNA genes are unlinked . However, a non-functional 5S remnant was found in the 26S-18S IGS. Its strand orientation, proximity to the 26S gene and flanking sequences were similar to functional 5S-IGS copies in the closely related Artemisia species (Figure 4A). These data suggest that a putative Elachanthemum progenitor did have linked arrangement of units that have become eventually lost. Hypothetically, a large deletion of locus or a chromosome loss could have promoted reamplification of unlinked 5S arrays and reversion of the unlinked rDNA genotype (Figure 10C). It will be interesting to analyse Asteraceae genomes (which are now largely unexplored) by deep sequencing in order to reveal intermediate stages of rDNA rearrangements, to map potentially active retroelement loci that mobilised 5S genes, and to determine the degree of tandem arrays homogeneity, in general.
The rDNA arrangement as a potential phylogenic tool
The intra- and interspecies divergence rates (Additional file 3) for 26S and 5S genes were comparable (low) confirming concerted evolution of the 35S-5S units and corroborating studies in yeasts . In contrast, the sequence between the 26S and 5S genes (termed IGS1) is not conserved and displays significant length and sequence polymorphisms (Figure 4). The cladogram based on these sequences showed highly supported branches (Figure 6B) and the tree topology presented good congruence with other phylogenetic reconstructions using different markers, as mentioned previously. The divergence of IGS1 was significantly higher than that of ITS1 (Additional file 3): between closely related species e.g. Artemisia tridentata and A. absinthium, the divergence of ITS1 was ~4% (10 substitutions per 250 bp) whereas for IGS1 the divergence was ~7% (35 substitutions per 488 bp). Thus IGS1 could be used as a suitable, alternative marker to infer phylogenetic relationships in taxonomically difficult taxa, e.g. Asteroideae. Besides, on the basis of the finding of linked arrangement in Anthemideae, Gnaphalieae and "Heliantheae alliance", maybe phylogenetic relationships between these groups should be reassessed, as these groups could be more closely related than currently considered [9, 21]. However, present nuclear and chloroplastic molecular phylogenies [ and references therein] confirm that, though related, these groups do not form a clade, but that Heliantheae alliance's sister group are tribes Athroismeae and Inuleae, and that Anthemideae are sister to Astereae (Figure 9). So in fact, our findings could reflect the lability of rDNA arrangements in these species.
The rearrangement of intrinsic locus structure in which thousands of copies are lost and replaced with new copies certainly represents a drastic genomic change with potential phylogenetic signal. Though a dual 5S and 35S FISH appears to be a reliable method to determine organisation of rDNA loci on chromosomes, it is somewhat impractical for large scale sampling as ours. In addition, the signals overlap may not necessarily mean physical linkage but rather juxtaposition of both arrays as demonstrated in the case of Dahlia pinnata (Figure 8). The simple PCR methodology described in Figure 1, based on primers amplifying 26S-5S and 5S-5S sequences, may represent a feasible approach.
Family Asteraceae nowadays occupy all kinds of habitats, from hot deserts to high mountain areas. Just in genus Artemisia there may have been 17-22 direct and 2-4 reverse migrations from Asia to America . We have shown that the unusual linked organisation of 35S-5S rDNA units occurs in three large groups of Asteraceae family comprising species-rich genera including abundant and economically important plants. The units were homogenised within the loci and across the chromosomes to several thousands of copies, with concomitant loss of the independent 5S and 35S tandems. Given that these tribes represent a large proportion of the species richness of the family we estimate that nearly a quarter of Asteraceae species could have evolved a linked arrangement of rRNA genes. We can therefore envisage that the linked rDNA arrangement in plants may be far more common than previously thought. Uncovering the adaptive advantages (if any) of a linked rDNA arrangement is the next step in the understanding of the dynamic evolution of rRNA genes in plants.
The authors wish to thank all people, Botanical Gardens, Herbaria and other institutions listed in Additional file 1 who generously contributed leaf or seed material to the present study. Prof. Joan Vallès (Universitat de Barcelona) and Dr. Teresa Garnatje (Institut Botànic de Barcelona, CSIC) are acknowledged for helpful suggestions that improved the work as well as for supplying many of the plant materials used. Also, personnel from the Departments of Molecular Epigenetics and Plant Developmental Genetics of the Institute of Biophysics, Academy of Sciences of the Czech Republic, are acknowledged for their help in different technical aspects. This research was funded by the Grant Agency of the Czech Republic (P501-10-0208), the Academy of Sciences of the Czech Republic (AVOZ50040507 and AVOZ50040702, MSMT LC0604), and the Ministry of Innovation and Science of Spain, MICINN (CGL2007-64839-C02/BOS). S. G. was granted a JAE-DOC contract from the CSIC (Superior Council of Scientific Investigations) and a José Castillejo Grant from the MICINN of the Spanish Government.
- Drouin G, de Sa MM: The concerted evolution of 5S ribosomal genes linked to the repeat units of other multigene families. Mol Biol Evol. 1995, 12 (3): 481-493.PubMedGoogle Scholar
- Vahidi H, Purac A, LeBlanc JM, Honda BM: Characterization of potentially functional 5S rRNA-encoding genes within ribosomal DNA repeats of the nematode Meloidogyne arenaria. Gene. 1991, 108 (2): 281-284. 10.1016/0378-1119(91)90446-I.PubMedView ArticleGoogle Scholar
- Bergeron J, Drouin G: The evolution of 5S ribosomal RNA genes linked to the rDNA units of fungal species. Curr Genet. 2008, 54 (3): 123-131. 10.1007/s00294-008-0201-2.PubMedView ArticleGoogle Scholar
- Drouin G, Hofman JD, Doolittle WF: Unusual ribosomal RNA gene organization in copepods of the genus Calanus. J Mol Biol. 1987, 196 (4): 943-946. 10.1016/0022-2836(87)90417-7.PubMedView ArticleGoogle Scholar
- Maizels N: Dictyostelium 17S, 25S, and 5S rDNAs lie within a 38,000 base pair repeated unit. Cell. 1976, 9 (3): 431-438. 10.1016/0092-8674(76)90088-X.PubMedView ArticleGoogle Scholar
- Capesius I: Analysis of the ribosomal RNA gene repeat from the moss Funaria hygrometrica. Plant Mol Biol. 1997, 33 (3): 559-564. 10.1023/A:1005740031313.PubMedView ArticleGoogle Scholar
- Sone T, Fujisawa M, Takenaka M, Nakagawa S, Yamaoka S, Sakaida M, Nishiyama R, Yamato KT, Ohmido N, Fukui K, et al: Bryophyte 5S rDNA was inserted into 45S rDNA repeat units after the divergence from higher land plants. Plant Mol Biol. 1999, 41 (5): 679-685. 10.1023/A:1006398419556.PubMedView ArticleGoogle Scholar
- Garcia S, Lim KY, Chester M, Garnatje T, Pellicer J, Valles J, Leitch AR, Kovarik A: Linkage of 35S and 5S rRNA genes in Artemisia (family Asteraceae): first evidence from angiosperms. Chromosoma. 2009, 118 (1): 85-97. 10.1007/s00412-008-0179-z.PubMedView ArticleGoogle Scholar
- Funk VA, Bayer RJ, Keeley S, Chan R, Watson L, Gemeinholzer B, Schilling E, Panrelo JL, Baldwin BG, Garcia-Jacas N, et al: Everywhere but Antartica: using a supertree to understand the diversity and distribution of the Compositae. Biol Skr. Edited by: Friis I, Balslev H. 2005, 55: 343-373.Google Scholar
- Kim KJ, Choi KS, Jansen RK: Two chloroplast DNA inversions originated simultaneously during the early evolution of the sunflower family (Asteraceae). Mol Biol Evol. 2005, 22 (9): 1783-1792. 10.1093/molbev/msi174.PubMedView ArticleGoogle Scholar
- Funk VA, Robinson H: Sampera, a new genus of Liabeae (Compositae or Asteraceae) from the northern Andes. P Biol Soc Wash. 2009, 122 (2): 155-161. 10.2988/08-31.1.View ArticleGoogle Scholar
- Fransz PF, Alonso-Blanco C, Liharska TB, Peeters AJ, Zabel P, de Jong JH: High-resolution physical mapping in Arabidopsis thaliana and tomato by fluorescence in situ hybridization to extended DNA fibres. Plant J. 1996, 9 (3): 421-430. 10.1046/j.1365-313X.1996.09030421.x.PubMedView ArticleGoogle Scholar
- Pires JC, Lim KY, Kovarík A, Matyasek R, Boyd A, Leitch AR, Leitch IJ, Bennett MD, Soltis PS, Soltis DE: Molecular cytogenetic analysis of recently evolved Tragopogon (Asteraceae) allopolyploids reveal a karyotype that is additive of the diploid progenitors. Am J Bot. 2004, 91 (7): 1022-1035. 10.3732/ajb.91.7.1022.PubMedView ArticleGoogle Scholar
- Dydak M, Kolano B, Nowak T, Siwinska D, Maluszynska J: Cytogenetic studies of three European species of Centaurea L. (Asteraceae). Hereditas. 2009, 146 (4): 152-161. 10.1111/j.1601-5223.2009.02113.x.PubMedView ArticleGoogle Scholar
- Cuellar T, Belhassen E, FernandezCalvin B, Orellana J, Bella JL: Chromosomal differentiation in Helianthus annuus var macrocarpus: Heterochromatin characterization and rDNA location. Heredity. 1996, 76: 586-591. 10.1038/hdy.1996.84.View ArticleGoogle Scholar
- Ruas CD, Vanzela ALL, Santos MO, Fregonezi JN, Ruas PM, Matzenbacher NI, de Aguiar-Perecin MLR: Chromosomal organization and phylogenetic relationships in Hypochaeris species (Asteraceae) from Brazil. Genet Mol Biol. 2005, 28 (1): 129-139. 10.1590/S1415-47572005000100023.View ArticleGoogle Scholar
- Weiss-Schneeweiss H, Stuessy TF, Siljak-Yakovlev S, Baeza CM, Parker J: Karyotype evolution in South American species of Hypochaeris (Asteraceae, Lactuceae). Plant Syst Evol. 2003, 241 (3-4): 171-184. 10.1007/s00606-003-0026-x.View ArticleGoogle Scholar
- Torrell M, Cerbah M, Siljak-Yakovlev S, Valles J: Molecular cytogenetics of the genus Artemisia (Asteraceae, Anthemideae): fluorochrome banding and fluorescence in situ hybridization. I. Subgenus Seriphidium and related taxa. Plant Syst Evol. 2003, 239 (1-2): 141-153. 10.1007/s00606-002-0259-0.View ArticleGoogle Scholar
- Garcia S, Garnatje T, Hidalgo O, McArthur ED, Siljak-Yakovlev S, Valles J: Extensive ribosomal DNA (18S-5.8S-26S and 5S) colocalization in the North American endemic sagebrushes (subgenus Tridentatae, Artemisia, Asteraceae) revealed by FISH. Plant Syst Evol. 2007, 267 (1-4): 79-92. 10.1007/s00606-007-0558-6.View ArticleGoogle Scholar
- Yoshikazu H, Matoba H, Kondo K: Physical mapping of ribosomal RNA genes in the genus Artemisia L. (Asteraceae). Caryologia. 2006, 59 (4): 312-318.View ArticleGoogle Scholar
- Panero JL, Funk VA: The value of sampling anomalous taxa in phylogenetic studies: Major clades of the Asteraceae revealed. Mol Phylogenet Evol. 2008, 47 (2): 757-782. 10.1016/j.ympev.2008.02.011.PubMedView ArticleGoogle Scholar
- Funk VA, Susanna A, Stuessy TF, Robinson H: Classification of Compositae. Systematics, evolution, and biogeography of Compositae. Edited by: Funk VA, Susanna A, Stuessy TF, Bayer RJ. 2009, Michigan, USA.: Sheridan Books Inc, 171-192.Google Scholar
- Lim KY, Kovarik A, Matyasek R, Bezdek M, Lichtenstein CP, Leitch AR: Gene conversion of ribosomal DNA in Nicotiana tabacum is associated with undermethylated, decondensed and probably active gene units. Chromosoma. 2000, 109 (3): 161-172. 10.1007/s004120050424.PubMedView ArticleGoogle Scholar
- Fulnecek J, Matyasek R, Kovarik A, Bezdek M: Mapping of 5-methylcytosine residues in Nicotiana tabacum 5S rRNA genes by genomic sequencing. Mol Gen Genet. 1998, 259 (2): 133-141. 10.1007/s004380050798.PubMedView ArticleGoogle Scholar
- Tate JA, Ni Z, Scheen A-C, Koh J, Gilbert CA, Lefkowitz D, Chen ZJ, Soltis PA, Soltis DE: Evolution and expression of homeologous loci in Tragopogon miscellus (Asteraceae), a recent and reciprocally formed allopolyploid. Genetics. 2006, 3: 1599-1611. 10.1534/genetics.106.057646.View ArticleGoogle Scholar
- Bancheva S, Greilhuber J: Genome size in Bulgarian Centaurea s.l. (Asteraceae). Plant Syst Evol. 2006, 257 (1-2): 95-117. 10.1007/s00606-005-0384-7.View ArticleGoogle Scholar
- Torrell M, Valles J: Genome size in 21 Artemisia L. species (Asteraceae, Anthemideae): Systematic, evolutionary, and ecological implications. Genome. 2001, 44 (2): 231-238. 10.1139/gen-44-2-231.PubMedView ArticleGoogle Scholar
- Garcia S, Inceer H, Garnatje T, Valles J: Genome size variation in some representatives of the genus Tripleurospermum. Biol Plantarum. 2005, 49 (3): 381-387. 10.1007/s10535-005-0011-z.View ArticleGoogle Scholar
- Zonneveld BJM, Leitch IJ, Bennett MD: First nuclear DNA amounts in more than 300 angiosperms. Ann Bot-London. 2005, 96 (2): 229-244. 10.1093/aob/mci170.View ArticleGoogle Scholar
- Komarova NY, Grabe T, Huigen DJ, Hemleben V, Volkov RA: Organization, differential expression and methylation of rDNA in artificial Solanum allopolyploids. Plant Mol Biol. 2004, 56 (3): 439-463. 10.1007/s11103-004-4678-x.PubMedView ArticleGoogle Scholar
- Hall TA: BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucl Acids Symp Ser. 1999, 41: 95-98.Google Scholar
- Galtier N, Gouy M, Gautier C: SEAVIEW and PHYLO_WIN: two graphic tools for sequence alignment and molecular phylogeny. Comput Appl Biosci. 1996, 12 (6): 543-548.PubMedGoogle Scholar
- Rozas J, Sanchez-DelBarrio JC, Messeguer X, Rozas R: DnaSP, DNA polymorphism analyses by the coalescent and other methods. Bioinformatics. 2003, 19 (18): 2496-2497. 10.1093/bioinformatics/btg359.PubMedView ArticleGoogle Scholar
- Panero JL, Funk VA: The value of sampling anomalous taxa in phylogenetic studies: major clades of the Asteraceae revealed. Mol Phylogenet Evol. 2008, 47 (2): 757-782. 10.1016/j.ympev.2008.02.011.PubMedView ArticleGoogle Scholar
- Kapitonov VV, Jurka J: A novel class of SINE elements derived from 5S rRNA. Mol Biol Evol. 2003, 20 (5): 694-702. 10.1093/molbev/msg075.PubMedView ArticleGoogle Scholar
- Kalendar R, Tanskanen J, Chang W, Antonius K, Sela H, Peleg O, Schulman AH: Cassandra retrotransposons carry independently transcribed 5S RNA. Proc Natl Acad Sci USA. 2008, 105 (15): 5833-5838. 10.1073/pnas.0709698105.PubMedPubMed CentralView ArticleGoogle Scholar
- Sanz M, Vilatersana R, Hidalgo O, Garcia-Jacas N, Susanna A, Schneeweiss GM, Valles J: Molecular phylogeny and evolution of floral characters of Artemisia and allies (Anthemideae, Asteraceae): evidence from nrDNA ETS and ITS sequences. Taxon. 2008, 57 (1): 66-78.Google Scholar
- Mazzella C, Rodriguez M, Vaio M, Gaiero P, Lopez-Carro B, Santinaque FF, Folle GA, Guerra M: Karyological Features of Achyrocline (Asteraceae, Gnaphalieae): Stable Karyotypes, Low DNA Content Variation and Linkage of rRNA Genes. Cytogenet Genome Res. 2010, 128 (1-3): 169-176. 10.1159/000290689.PubMedView ArticleGoogle Scholar
- Magallon S, Crane PR, Herendeen PS: Phylogenetic pattern, diversity, and diversification of eudicots. Ann Mo Bot Gard. 1999, 86 (2): 297-372. 10.2307/2666180.View ArticleGoogle Scholar
- Oberprieler C: Temporal and spatial diversification of Circum-Mediterranean Compositae-Anthemideae. Taxon. 2005, 54 (4): 951-966. 10.2307/25065480.View ArticleGoogle Scholar
- Jo SH, Koo DH, Kim JF, Hur CG, Lee S, Yang TJ, Kwon SY, Choi D: Evolution of ribosomal DNA-derived satellite repeat in tomato genome. BMC Plant Biol. 2009, 9: 42-10.1186/1471-2229-9-42.PubMedPubMed CentralView ArticleGoogle Scholar
- Lim KY, Skalicka K, Koukalova B, Volkov RA, Matyasek R, Hemleben V, Leitch AR, Kovarik A: Dynamic changes in the distribution of a satellite homologous to intergenic 26-18S rDNA spacer in the evolution of Nicotiana. Genetics. 2004, 166 (4): 1935-1946. 10.1534/genetics.166.4.1935.PubMedPubMed CentralView ArticleGoogle Scholar
- Eickbush TH, Eickbush DG: Finely orchestrated movements: evolution of the ribosomal RNA genes. Genetics. 2007, 175 (2): 477-485. 10.1534/genetics.107.071399.PubMedPubMed CentralView ArticleGoogle Scholar
- Chester M, Sykorova E, Fajkus J, Leitch AR: Single integration and spread of a copia-like sequence nested in rDNA Intergenic Spacers of Allium cernuum (Alliaceae). Cytogenet Genome Res. 2010, 129 (1-3): 35-46. 10.1159/000312959.PubMedView ArticleGoogle Scholar
- Smith GP: Unequal crossover and the evolution of multigene families. Cold Spring Harb Symp Quant Biol. 1974, 38: 507-513.PubMedView ArticleGoogle Scholar
- Navratilova A, Koblizkova A, Macas J: Survey of extrachromosomal circular DNA derived from plant satellite repeats. Bmc Plant Biol. 2008, 8: 90-10.1186/1471-2229-8-90.PubMedPubMed CentralView ArticleGoogle Scholar
- Cohen S, Agmon N, Sobol O, Segal D: Extrachromosomal circles of satellite repeats and 5S ribosomal DNA in human cells. Mob DNA. 2010, 1 (1): 11-10.1186/1759-8753-1-11.PubMedPubMed CentralView ArticleGoogle Scholar
- Koukalova B, Moraes AP, Renny-Byfield S, Matyasek R, Leitch AR, Kovarik A: Fall and rise of satellite repeats in allopolyploids of Nicotiana over c. 5 million years. New Phytol. 2010, 186 (1): 148-160. 10.1111/j.1469-8137.2009.03101.x.PubMedView ArticleGoogle Scholar
- Ganley AR, Kobayashi T: Highly efficient concerted evolution in the ribosomal DNA repeats: total rDNA repeat variation revealed by whole-genome shotgun sequence data. Genome Res. 2007, 17 (2): 184-191. 10.1101/gr.5457707.PubMedPubMed CentralView ArticleGoogle Scholar
- Tkach NV, Hoffmann MH, Roser M, Korobkov AA, von Hagen KB: Parallel evolutionary patterns in multiple lineages of arctic Artemisia L. (Asteraceae). Evolution. 2008, 62 (1): 184-198.PubMedGoogle Scholar
- Cloix C, Yukawa Y, Tutois S, Sugiura M, Tourmente S: In vitro analysis of the sequences required for transcription of the Arabidopsis thaliana 5S rRNA genes. Plant J. 2003, 35 (2): 251-261. 10.1046/j.1365-313X.2003.01793.x.PubMedView ArticleGoogle Scholar