Repeated reunions and splits feature the highly dynamic evolution of 5S and 35S ribosomal RNA genes (rDNA) in the Asteraceae family

Background  

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.     

Asexuality Associated with Marked Genomic Expansion of Tandemly Repeated rRNA and Histone Genes

How does asexual reproduction influence genome evolution? Although is it clear that genomic structural variation is common and important in natural populations, we know very little about how one of the most fundamental of eukaryotic traits—mode of genomic inheritance—influences genome structure. We address this question with the New Zealand freshwater snail Potamopyrgus antipodarum, which features multiple separately derived obligately asexual lineages that coexist and compete with otherwise similar sexual lineages. We used whole-genome sequencing reads from a diverse set of sexual and asexual individuals to analyze genomic abundance of a critically important gene family, rDNA (the genes encoding rRNAs), that is notable for dynamic and variable copy number. Our genomic survey of rDNA in P. antipodarum revealed two striking results. First, the core histone and 5S rRNA genes occur between tandem copies of the 18S–5.8S–28S gene cluster, a unique architecture for these crucial gene families. Second, asexual P. antipodarum harbor dramatically more rDNA–histone copies than sexuals, which we validated through molecular and cytogenetic analysis. The repeated expansion of this genomic region in asexual P. antipodarum lineages following distinct transitions to asexuality represents a dramatic genome structural change associated with asexual reproduction—with potential functional consequences related to the loss of sexual reproduction.     

Linkage of 35S and 5S rRNA genes in <Emphasis Type="Italic">Artemisia</Emphasis> (family Asteraceae): first evidence from angiosperms

Typically in plants, the 5S and 35S ribosomal DNA (rDNA) encoding two major ribosomal RNA species occur at separate loci. However, in some algae, bryophytes and ferns, they are at the same locus (linked arranged). Southern blot hybridisation, polymerase chain reactions (PCR), fluorescent in situ hybridisation, cloning and sequencing were used to reveal 5S and 35S rDNA genomic organisation in Artemisia. We observed thousands of rDNA units at two–three loci containing 5S rDNA in an inverted orientation within the inter-genic spacer (IGS) of 35S rDNA. The sequenced clones of 26–18S IGS from Artemisia absinthium appeared to contain a conserved 5S gene insertion proximal to the 26S gene terminus (5S rDNA-1) and a second less conserved 5S insertion (5S rDNA-2) further downstream. Whilst the 5S rDNA-1 showed all the structural features of a functional gene, the 5S-rDNA-2 had a deletion in the internal promoter and probably represents a pseudogene. The linked arrangement probably evolved before the divergence of Artemisia from the rest of Asteraceae (>10 Myrs). This arrangement may have involved retrotransposons and once formed spread via mechanisms of concerted evolution. Heterogeneity in unit structure may reflect ongoing homogenisation of variant unit types without fixation for any particular variant. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Unique arrangement of coding sequences for 5 S, 5.8 S, 18 S and 25 S ribosomal RNA in Saccharomyces cerevisiae as determined by R-loop and hybridization analysis.

The arrangement of the coding sequences for the 5 S, 5.8 S, 18 S and 25 S ribosomal RNA from Saccharomyces cerevisiae was analyzed in λ-yeast hybrids containing repeating units of the ribosomal DNA. After mapping of restriction sites, the positions of the coding sequences were determined by hybridization of purified rRNAs to restriction fragments, by R-loop analysis in the electron microscope, and by electrophoresis of S₁ nuclease-treated rRNA/rDNA hybrids in alkaline agarose gels. The R-loop method was improved with respect to the length calibration of RNA/DNA duplexes and to the spreading conditions resulting in fully extended 18 S and 25 S rRNA R-loops. The qualitative results are: (1) the 5 S rRNA genes, unlike those in higher eukaryotes, alternate with the genes of the precursor for the 5.8 S, 18 S and 25 S rRNA; (2) the coding sequence for 5.8 S rRNA maps, as in higher eukaryotes, between the 18 S and 25 S rRNA coding sequences. The quantitative results are: (1) the tandemly repeating rDNA units have a constant length of 9060 ± 100 nucleotide pairs with one SstI, two HindIII and, dependent on the strain, six or seven EcoRI sites; (2) the 18 S and 25 S rRNA coding regions consist of 1710 ± 80 and 3360 ± 80 nucleotide pairs, respectively; (3) an 18 S rRNA coding region is separated by a 780 ± 70 nucleotide pairs transcribed spacer from a 25 S rRNA coding region. This is then followed by a 3210 ± 100 nucleotide pairs mainly non-transcribed spacer which contains a 5 S rRNA gene.     

Secondary structure prediction for complete rDNA sequences (18S, 5.8S,and 28S rDNA) of Demodex folliculorum, and comparison of divergent domains structures across Acari

According to base pairing, the rRNA folds into corresponding secondary structures, which contain additional phylogenetic information. On the basis of sequencing for complete rDNA sequences (18S, ITS1, 5.8S, ITS2 and 28S rDNA) of Demodex, we predicted the secondary structure of the complete rDNA sequence (18S, 5.8S, and 28S rDNA) of Demodex folliculorum, which was in concordance with that of the main arthropod lineages in past studies. And together with the sequence data from GenBank, we also predicted the secondary structures of divergent domains in SSU rRNA of 51 species and in LSU rRNA of 43 species from four superfamilies in Acari (Cheyletoidea, Tetranychoidea, Analgoidea and Ixodoidea). The multiple alignment among the four superfamilies in Acari showed that, insertions from Tetranychoidea SSU rRNA formed two newly proposed helixes, and helix c3-2b of LSU rRNA was absent in Demodex (Cheyletoidea) taxa. Generally speaking, LSU rRNA presented more remarkable differences than SSU rRNA did, mainly in D2, D3, D5, D7a, D7b, D8 and D10.     

Physical mapping of 5S and 18S-26S ribosomal RNA gene families inAllium victorialis var.platyphyllum

A physical map of the 5S and 18S–26S rRNA genes was determined using bi-color fluorescencein situ hybridization technique inA. victorialis var.platyphyllum. 5S rRNA genes were positioned in the intercalary regions of the short arms in homologous chromosomes 6. Two major loci of the 18S-26S rRNA genes were detected in the secondary constrictions flanking with a pair of satellite and terminal region of short arm in chromosome 4. And two additional minor loci were heterotype, representing one signal on the terminal region of the short arm in one homolog of chromsome 2, and other on one homolog of chromosome 6 with linked 5S rRNA loci. In addition chromomycin A₃ (CMA,) fluorescent banding method was used to identify the relation between Nucleolus Organizer Region (NOR) sites and CMA, positive heterochromatin sites. In homologous chromosome 4 showing 18S–26S rDNA hybridization signals revealed also distinct CMA, positive band.     

Organization of 5S ribosomal DNA of Melitaea trivia

Two length variants of 5S rDNA repeated units were detected in the genome of East European butterfly Melitaea trivia. Both repeat variants contain the 5S rRNA coding region of the same length of 120 bp, but possess the intergenic spacer region (IGS) of different size, 78 and 125 bp, respectively. The level of sequence similarity between the two 5S rDNA variants amounts to 43.9-45.5% in the IGS, whereas the coding region appears to be more conservative. In the IGS, microsatellite sequence motives were found; amplification of these motives could be involved in the evolution of the 5S rDNA.     

Structural motifs of syringyl peroxidases predate not only the gymnosperm-angiosperm divergence but also the radiation of tracheophytes

* The most distinctive variation in the monomer composition of lignins in vascular land plants is that found between the two main groups of seed plants. Thus, while gymnosperm lignins are typically composed of guaiacyl (G) units, angiosperm lignins are largely composed of similar levels of G and syringyl (S) units. * However, and contrary to what might be expected, peroxidases isolated from basal (Cycadales and Ginkgoales) and differentially evolved (Coniferales and Gnetales) gymnosperms are also able to oxidize S moieties, and this ability is independent of the presence or absence of S-type units in their lignins. * The results obtained led us to look at the protein database to search for homologies between gymnosperm peroxidases and true eudicot S-peroxidases, such as the Zinnia elegans peroxidase. * The findings showed that certain structural motifs characteristic of eudicot S-peroxidases (certain amino acid sequences and beta-sheet secondary structures) predate the gymnosperm-angiosperm divergence and the radiation of tracheophytes, since they are found not only in peroxidases from basal gymnosperms, ferns and lycopods, but also in peroxidases from the moss Physcomitrella patens (Bryopsida) and the liverwort Marchantia polymorpha (Marchantiopsida), which, as typical of bryophytes, do not have xylem tissue nor lignins.     

Restriction site and length polymorphism of the rDNA unit in the cultivated basidiomycetePleurotus cornucopiae

In the ribosomal DNA unit ofPleurotus cornucopiae, the rDNA coding regions are in the order 5, 5S-18S-5.8S-25S, 3, with the 5 location of the 5S gene differing from its 3 location found in other basidiomycetes. The most discriminating probe used to study the rDNA polymorphism consisted of a fragment that included the 5S, 18S and part of the 5.8S and 25S genes flanking three intergenic sequences. A high degree of rDNA polymorphism was observed in the sevenP. cornucopiae dikaryons studied. For the first time within a basidiomycete species, the restrictions maps distinguished two types of rDNA units (I and II). In each rDNA type, length variations in the external intergenic sequence IGS 1 located between the 25S and 5S genes allowed characterization of two different rDNA units in type I and four rDNA units in type II. This suggested that theP. cornucopiae rDNA units were derived from two kinds of ancestors (type I and II) by insertion or deletion events (100–700 bp) in the IGS 1. In four dikaryotic strains, two rDNA units of the same type (I or II) differing only by the IGS 1 length, were found in a similar number of copies, and presented a meiotic segregation in homokaryotic progeny. In one progeny, some homokaryotic strains possessed two different rDNA units: one with a high copy number and another with a lower one, showing that two different rDNA units could coexist in a single nucleus.     

Molecular and chromosomal organization of DNA sequences coding for the ribosomal RNAs in cereals

The chromosomal locations of ribosomal DNA in wheat, rye and barley have been determined by in situ hybridization using high specific activity ¹²⁵I-rRNA. The 18S-5.8S-26S rRNA gene repeat units in hexaploid wheat (cv. Chinese Spring) are on chromosomes 1B, 6B and 5D. In rye (cv. Imperial) the repeat units occur at a single site on chromosome 1R(E), while in barley (cv. Clipper) they are on both the chromosomes (6 and 7) which show secondary constrictions. In wheat and rye the major 5S RNA gene sites are close to the cytological secondary constrictions where the 18S-5.8S-26S repeating units are found, but in barley the site is on a chromosome not carrying the other rDNA sequences. — Restriction enzyme and R-loop analyses showed the 18S-5.8S-26S repeating units to be approximately 9.5 kb long in wheat, 9.0 kb in rye and barley to have two repeat lengths of 9.5 kb and 10 kb. Electron microscopic and restriction enzyme data suggest that the two barley forms may not be interpersed. Digestion with EcoR1 gave similar patterns in the three species, with a single site in the 26S gene. Bam H1 digestion detected heterogeneity in the spacer regions of the two different repeats in barley, while in rye and wheat heterogeneity was shown within the 26S coding sequence by an absence of an effective Bam H1 site in some repeat units. EcoR1 and Bam H1 restriction sites have been mapped in each species. — The repeat unit of the 5S RNA genes was approximately 0.5 kb in wheat and rye and heterogeneity was evident. The analysis of the 5S RNA genes emphasizes the homoeology between chromosomes 1B of wheat and 1R of rye since both have these genes in the same position relative to the secondary constriction. In barley we did not find a dominant monomer repeat unit for the 5S genes.     

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

To study the relationship between uniparental rDNA (encoding 18S, 5.8S and 26S ribosomal RNA) silencing (nucleolar dominance) and rRNA gene dosage, we studied a recently emerged (within the last 80 years) allotetraploid Tragopogon mirus (2n=24), formed from the diploid progenitors T. dubius (2n=12, D-genome donor) and T. porrifolius (2n=12, P-genome donor). Here, we used molecular, cytogenetic and genomic approaches to analyse rRNA gene activity in two sibling T. mirus plants (33A and 33B) with widely different rRNA gene dosages. Plant 33B had ~400 rRNA genes at the D-genome locus, which is typical for T. mirus, accounting for ~25% of total rDNA. We observed characteristic expression dominance of T. dubius-origin genes in all organs. Its sister plant 33A harboured a homozygous macrodeletion that reduced the number of T. dubius-origin genes to about 70 copies (~4% of total rDNA). It showed biparental rDNA expression in root, flower and callus, but not in leaf where D-genome rDNA dominance was maintained. There was upregulation of minor rDNA variants in some tissues. The RNA polymerase I promoters of reactivated T. porrifolius-origin rRNA genes showed reduced DNA methylation, mainly at symmetrical CG and CHG nucleotide motifs. We hypothesise that active, decondensed rDNA units are most likely to be deleted via recombination. The silenced homeologs could be used as a ‘first reserve'' to ameliorate mutational damage and contribute to evolutionary success of polyploids. Deletion and reactivation cycles may lead to bidirectional homogenisation of rRNA arrays in the long term.     

A novel clade of Prochlorococcus found in high nutrient low chlorophyll waters in the South and Equatorial Pacific Ocean

A novel high-light (HL)-adapted Prochlorococcus clade was discovered in high nutrient and low chlorophyll (HNLC) waters in the South Pacific Ocean by phylogenetic analyses of 16S ribosomal RNA (rRNA) and 16S–23S internal transcribed spacer (ITS) sequences. This clade, named HNLC fell within the HL-adapted Prochlorococcus clade with sequences above 99% similarity to one another, and was divided into two subclades, HNLC1 and HNLC2. The distribution of the whole HNLC clade in a northwest to southeast transect in the South Pacific (HNLC-to-gyre) and two 8°N to 8°S transects in the Equatorial Pacific was determined by quantitative PCR using specific primers targeting ITS regions. HNLC was the dominant HL Prochlorococcus clade (2–9% of bacterial 16S rRNA genes) at the three westernmost stations in the South Pacific but decreased to less than 0.1% at the other stations being replaced by the eMIT9312 ecotype in the hyperoligotrophic gyre. The highest contributions of HNLC Prochlorococcus in both Equatorial Pacific transects along the latitudinal lines of 170°W and 155°W were observed at the southernmost stations, reaching 16 and 6% of bacterial 16S rRNA genes, respectively, whereas eMIT9312 dominated near the Equator. Spearman Rank Order correlation analysis indicated that although both the HNLC clade and eMIT9312 were correlated with temperature, they showed different correlations with regard to nutrients. HNLC only showed significant correlations to ammonium uptake and regeneration rates, whereas eMIT9312 was negatively correlated with inorganic nutrients.