5S rDNA and U2 snRNA are linked in the genome of Crassostrea angulata and Crassostrea gigas oysters: does the (CT)n.(GA)n microsatellite stabilize this novel linkage of large tandem arrays?

Cytogenetic markers using single-sequence probes reveal chromosomal locations of tandemly repetitive genes in scleractinian coral Acropora pruinosa

The short and similar sized chromosomes of Acropora pose a challenge for karyotyping. Conventional methods, such as staining of heterochromatic regions, provide unclear banding patterns that hamper identification of such chromosomes. In this study, we used short single-sequence probes from tandemly repetitive 5S ribosomal RNA (rRNA) and core histone coding sequences to identify specific chromosomes of Acropora pruinosa. Both the probes produced intense signals in fluorescence in situ hybridization, which distinguished chromosome pairs. The locus of the 5S rDNA probe was on chromosome 5, whereas that of core histone probe was on chromosome 8. The sequence of the 5S rDNA probe was composed largely of U1 and U2 spliceosomal small nuclear RNA (snRNA) genes and their interspacers, flanked by short sequences of the 5S rDNA. This is the first report of a tandemly repetitive linkage of snRNA and 5S rDNA sequences in Cnidaria. Based on the constructed tentative karyogram and whole genome hybridization, the longest chromosome pair (chromosome 1) was heteromorphic. The probes also hybridized effectively with chromosomes of other Acropora species and population, revealing an additional core histone gene locus. We demonstrated the applicability of short-sequence probes as chromosomal markers with potential for use across populations and species of Acropora.

Sequence characterization of the 5S ribosomal DNA and the internal transcribed spacer (ITS) region in four European Donax species (Bivalvia: Donacidae)

Background

The whole repeat unit of 5S rDNA and the internal transcribed spacer (ITS) of four European Donax species were analysed. After amplifying, cloning and sequencing several 5S and ITS units, their basic features and their variation were described. The phylogenetic usefulness of 5S and ITS sequences in the inference of evolutionary relationships among these wedge clams was also investigated.

Results

The length of the 5S repeat presented little variation among species, except D. trunculus that differed from the rest of the Donax species in 170–210 bp. The deduced coding region covered 120 bp, and showed recognizable internal control regions (ICRs) involved in the transcription. The length of non-transcribed spacer region (NTS) ranged from 157 bp to 165 bp in Donax trunculus and from 335 bp to 367 bp in the other three species. The conservation degree of transcriptional regulatory regions was analysed revealing a conserved TATA-like box in the upstream region. Regarding ITS sequences, the four Donax species showed slight size differences among clones due to the variation occurring in the ITS1 and ITS2, except Donax variegatus did not display size differences in the ITS2. The total length of the ITS sequence ranged between 814 and 1014 bp. Resulting phylogenetic trees display that the two ribosomal DNA regions provide well-resolved phylogenies where the four European Donax species form a single clade receiving high support in nodes. The topology obtained with 5S sequences was in agreement with Donax evolutionary relationships inferred from several sequences of different nature in previous studies.

Conclusions

This is not only a basic research work, where new data and new knowledge is provided about Donax species, but also have allowed the authentication of these wedge clams and offers future applications to provide other genetic resources.

Electronic supplementary material

The online version of this article (10.1186/s12863-018-0684-x) contains supplementary material, which is available to authorized users.

Dispersion of transposable elements and multigene families: Microstructural variation in Characidium (Characiformes: Crenuchidae) genomes

Eukaryotic genomes consist of several repetitive DNAs, including dispersed DNA sequences that move between chromosome sites, tandem repeats of DNA sequences, and multigene families. In this study, repeated sequences isolated from the genome of Characidium gomesi were analyzed and mapped to chromosomes in Characidium zebra and specimens from two populations of C. gomesi. The sequences were transposable elements (TEs) named retroelement of Xiphophorus (Rex); multigene families of U2 small nuclear RNA (U2 snRNA); and histones H1, H3, and H4. Sequence analyses revealed that U2 snRNA contains a major portion corresponding to the Tx1-type non-LTR retrotransposon Keno, the preferential insertion sites of which are U2 snRNA sequences. All histone sequences were found to be associated with TEs. In situ localization revealed that these DNA sequences are dispersed throughout the autosomes of the species, but they are not involved in differentiation of the specific region of the W sex chromosome in C. gomesi. We discuss mechanisms of TE invasion into multigene families that lead to microstructural variation in Characidium genomes.

Chromosomal Mapping of Repeat DNA in Bergiaria westermanni (Pimelodidae, Siluriformes): Localization of 45S rDNA in B Chromosomes

The occurrence of repetitive DNA in autosomes and B chromosomes of Bergiaria westermanni was examined using conventional and molecular cytogenetic techniques. This species exhibited 2n = 56 chromosomes, with intra- and interindividual variation in the number of heterochromatic B chromosomes (from 0 to 4). The 5S rDNA was localized in pairs 1 and 5, and histone probes (H1, H3, and H4) and U2 small nuclear RNA were syntenic with 5S rDNA in pair 5. Histone sequences were also located in chromosome pair 14. The (GATA)n sequence was dispersed throughout the autosomes and B chromosomes, with clusters (microsatellite accumulation) in some chromosome regions. The telomeric probe revealed no signs of chromosomal rearrangements in the genome of B. westermanni. The 45S rDNA sites were detected in the terminal region of pair 27; these sites corresponded to a GC-rich heterochromatin block. In addition, 3 of the 4 B chromosomes also contained 45S rDNA copies. Silver nitrate staining in interphase nuclei provided indirect evidence of the expression of these rRNA genes in B chromosomes, indicating the probable origin of these elements. This report shows plasticity in the chromosomal localization of repeat DNA in B. westermanni and features a discussion of genomic diversification.

The 5S rDNA in two Abracris grasshoppers (Ommatolampidinae: Acrididae): molecular and chromosomal organization

The 5S ribosomal DNA (rDNA) sequences are subject of dynamic evolution at chromosomal and molecular levels, evolving through concerted and/or birth-and-death fashion. Among grasshoppers, the chromosomal location for this sequence was established for some species, but little molecular information was obtained to infer evolutionary patterns. Here, we integrated data from chromosomal and nucleotide sequence analysis for 5S rDNA in two Abracris species aiming to identify evolutionary dynamics. For both species, two arrays were identified, a larger sequence (named type-I) that consisted of the entire 5S rDNA gene plus NTS (non-transcribed spacer) and a smaller (named type-II) with truncated 5S rDNA gene plus short NTS that was considered a pseudogene. For type-I sequences, the gene corresponding region contained the internal control region and poly-T motif and the NTS presented partial transposable elements. Between the species, nucleotide differences for type-I were noticed, while type-II was identical, suggesting pseudogenization in a common ancestor. At chromosomal point to view, the type-II was placed in one bivalent, while type-I occurred in multiple copies in distinct chromosomes. In Abracris, the evolution of 5S rDNA was apparently influenced by the chromosomal distribution of clusters (single or multiple location), resulting in a mixed mechanism integrating concerted and birth-and-death evolution depending on the unit.

Microsatellite Organization in the B Chromosome and A Chromosome Complement in Astyanax (Characiformes, Characidae) Species

The organization of microsatellites in B and sex chromosomes has been linked to chromosomal evolution in a number of animal groups. Here, the chromosomal organizations of (CA)15, (GA)15, (CG)15, (GACA)4, and (GATA)8 microsatellites were examined in several Astyanax species with different diploid numbers: Astyanax mexicanus (2n = 50 + 1 B chromosome), A. altiparanae (2n = 50), A. marionae (2n = 48), A. fasciatus (2n = 46), and A. schubarti (2n = 36). The (CA)15 and (GA)15 microsatellites were dispersed across the chromosomes of A. altiparanae and A. fasciatus but were also observed as clusters (CA and GA for A. altiparanae, and CA for A. fasciatus). In A. marionae and A. schubarti, the (CA)15 and (GA)15 microsatellites were dispersed but were also observed as clustered signals and coincident with heterochromatic regions. In all 4 of these species, the (CG)15 and (GACA)4 microsatellites were dispersed across chromosomes, and the (GATA)8 microsatellite was co-localized with 5S rDNA. In A. mexicanus, the (CA)15, (GA)15, (CG)15, (GATA)8, and (GACA)4 microsatellites were weakly detected and dispersed across the chromosomes of the A complement. On the B chromosome, signals for the different microsatellites were weak, strong, absent, weak, and absent, respectively. The distribution of microsatellites and the locational relationship between microsatellites and 5S rDNA are discussed, and a possible evolutionary pathway is proposed for microsatellites in Astyanax.

Chromosomal Mapping of Repetitive DNA Sequences in Five Species of Astyanax (Characiformes, Characidae) Reveals Independent Location of U1 and U2 snRNA Sites and Association of U1 snRNA and 5S rDNA

Astyanax is a genus of Characidae fishes currently composed of 155 valid species. Previous cytogenetic studies revealed high chromosomal diversification among them, and several studies have been performed using traditional cytogenetic techniques to investigate karyotypes and chromosomal locations of 18S and 5S rDNA genes. However, only a few studies are currently available about other repetitive sequences. Here, the chromosomal location of small nuclear RNA genes, identified as U1 and U2 snRNA clusters, was established and compared to the distribution of 5S rDNA and histone clusters in 5 Astyanax species (A. paranae, A. fasciatus, A. bockmanni, A. altiparanae, and A. jordani) using FISH. The cytogenetic mapping of U1 and U2 snRNA demonstrated a conserved pattern in the number of sites per genome independent of the location in Astyanax species. The location of the U1 snRNA gene was frequently associated with 5S rDNA sequences, indicating a possible interaction between the distinct repetitive DNA families. Finally, comparisons involving the location of U1 and U2 snRNA clusters in the chromosomes of Astyanax species revealed a very diverse pattern, suggesting that many rearrangements have occurred during the diversification process of this group.

Structure and Organization of the Engraulidae Family U2 snRNA: An Evolutionary Model Gene?

The U2 snRNA multigene family has been analyzed in four species of the Engraulidae family--Engraulis encrasicolus, Engraulis mordax, Engraulis ringens, and Engraulis japonicas--with the object of understanding more about the structure of this multigene family in these pelagic species and studying their phylogenetic relationships. The results showed that the cluster of this gene family in the Engraulis genus is formed by the U2-U5 snRNA with highly conserved sequences of mini- and micro-satellites, such as (CTGT)n, embedded downstream of the transcription unit; findings indicate that this gene family evolved following the concerted model. The phylogenetic analysis of the non-transcribed spacer of cluster U2-U5 snDNA in the 4 species showed that the sequences of the species E. encrasicolus and E. japonicus are closely related; these two are genetically close to E. mordax and slightly more distant from E. ringens. The data obtained by molecular analysis of U2-U5 snDNA and their secondary structure, with the presence of the micro-satellite (CTGT)n and mini-satellites, show clearly that the species E. encrasicolus and E. japonicus are closely related and would be older than E. mordax and E. ringens.

Adaptation to abiotic stress in the oyster Crassostrea angulata relays on genetic polymorphisms

Here we describe the whole genome re-sequencing of the Portuguese oyster Crassostrea angulata, an edible cupped oyster of major commercial importance with an important role as biosensor of coastal water pollution. We sequenced the genome of the C. angulata to 29.3-fold coverage using ABI SOLID system. Comparisons of the sequences with the reference assembly of the Pacific oyster (Crassostrea gigas), yielded 129 million SNPs, 151,620 from which were located in 20,908 genes from the C. gigas database. The analysis of Gene Ontology (GO) terms associated with gene regions containing SNPs, revealed that significant GO terms showing differences between the two oyster species, were related to activities of response to stress caused both by drying and by metal contamination. In the Biological Process domain, the GO terms ion transport, phosphorylation and proteolysis processes, among others, showed many polymorphic genes in C. angulata. These processes are related to combating genotoxic and hypo-osmotic stress in the oyster. It is noteworthy that more than 200 polymorphic genes were associated with DNA repair processes. These results reveal that most of the gene polymorphisms observed in C. angulata are associated with processes related to genome adaptation to abiotic stress in estuarine regions and support that genetic polymorphisms may be the base to the observed ability of C. angulata to retain the phenomenally high concentrations of toxic heavy metals. Our results also provide the framework for future investigations to establish the molecular basis of phenotypic variation of adaptive traits and should contribute to the management of the species' genetic resources.

High evolutionary dynamism in 5S rDNA of fish: state of the art

The 5S ribosomal DNA (rDNA) consists of one transcriptional unit of about 120 base pairs, which is separated from the next unit by a non-transcribed spacer (NTS). The coding sequence and the NTS together form a repeat unit which can be found in hundreds to thousands of copies tandemly repeated in the genomes. The NTS regions seem to be subject to rapid evolution. The first general model of evolution of these multigene families was referred to as divergent evolution, based on studies using hemoglobin and myoglobin as model systems. Later studies showed that nucleotide sequences of different multigene family members are more closely related within species than between species. This observation led to a new model of multigene family evolution, termed concerted evolution. Another model of evolution, named the birth-and-death model, has been found to be more suitable to explain the long-term evolution of these multigene families. According to this model, new genes originate by successive duplications, and these new genes are either maintained for a long time or are lost, or else degenerate into pseudogenes. In this review we describe different sources of variability in the 5S rDNA genes observed in several distinct fish species. This variability is mainly referred to NTSs and includes the presence of other multigene families (mainly LINEs, SINEs, non-LTR retrotransposons, and U snRNA families). Different types of microsatellites have also been found to contribute to the increase of variability in this region. Our recent results suggest that horizontal transfer contributes to the increase of diversity in the NTSs of some species. Variability in the 5S rDNA coding region affecting the stability of the structure, but without effects on the function of the 5S rRNA, is also described. Retrotransposons seem to be responsible for the high dynamism of 5S rDNA, while microsatellites acting as recombination hot spots could stabilize a wide variety of unusual DNA structures, affecting DNA replication and enhancing or decreasing promoter activity in gene expression. The relationship between the high variability found at molecular level and the low variability found at chromosomal level is also discussed.

Systematic analysis and evolution of 5S ribosomal DNA in metazoans

Several studies on 5S ribosomal DNA (5S rDNA) have been focused on a subset of the following features in mostly one organism: number of copies, pseudogenes, secondary structure, promoter and terminator characteristics, genomic arrangements, types of non-transcribed spacers and evolution. In this work, we systematically analyzed 5S rDNA sequence diversity in available metazoan genomes, and showed organism-specific and evolutionary-conserved features. Putatively functional sequences (12,766) from 97 organisms allowed us to identify general features of this multigene family in animals. Interestingly, we show that each mammal species has a highly conserved (housekeeping) 5S rRNA type and many variable ones. The genomic organization of 5S rDNA is still under debate. Here, we report the occurrence of several paralog 5S rRNA sequences in 58 of the examined species, and a flexible genome organization of 5S rDNA in animals. We found heterogeneous 5S rDNA clusters in several species, supporting the hypothesis of an exchange of 5S rDNA from one locus to another. A rather high degree of variation of upstream, internal and downstream putative regulatory regions appears to characterize metazoan 5S rDNA. We systematically studied the internal promoters and described three different types of termination signals, as well as variable distances between the coding region and the typical termination signal. Finally, we present a statistical method for detection of linkage among noncoding RNA (ncRNA) gene families. This method showed no evolutionary-conserved linkage among 5S rDNAs and any other ncRNA genes within Metazoa, even though we found 5S rDNA to be linked to various ncRNAs in several clades.

Evolutionary dynamics of rRNA gene clusters in cichlid fish

Background

Among multigene families, ribosomal RNA (rRNA) genes are the most frequently studied and have been explored as cytogenetic markers to study the evolutionary history of karyotypes among animals and plants. In this report, we applied cytogenetic and genomic methods to investigate the organization of rRNA genes among cichlid fishes. Cichlids are a group of fishes that are of increasing scientific interest due to their rapid and convergent adaptive radiation, which has led to extensive ecological diversity.

Results

The present paper reports the cytogenetic mapping of the 5S rRNA genes from 18 South American, 22 African and one Asian species and the 18S rRNA genes from 3 African species. The data obtained were comparatively analyzed with previously published information related to the mapping of rRNA genes in cichlids. The number of 5S rRNA clusters per diploid genome ranged from 2 to 15, with the most common pattern being the presence of 2 chromosomes bearing a 5S rDNA cluster. Regarding 18S rDNA mapping, the number of sites ranged from 2 to 6, with the most common pattern being the presence of 2 sites per diploid genome. Furthermore, searching the Oreochromis niloticus genome database led to the identification of a total of 59 copies of 5S rRNA and 38 copies of 18S rRNA genes that were distributed in several genomic scaffolds. The rRNA genes were frequently flanked by transposable elements (TEs) and spread throughout the genome, complementing the FISH analysis that detect only clustered copies of rRNA genes.

Conclusions

The organization of rRNA gene clusters seems to reflect their intense and particular evolutionary pathway and not the evolutionary history of the associated taxa. The possible role of TEs as one source of rRNA gene movement, that could generates the spreading of ribosomal clusters/copies, is discussed. The present paper reinforces the notion that the integration of cytogenetic data and genomic analysis provides a more complete picture for understanding the organization of repeated sequences in the genome.

Genetic characterization of Plectorhinchus mediterraneus yields important clues about genome organization and evolution of multigene families

Background

Molecular and cytogenetic markers are of great use for to fish characterization, identification, phylogenetics and evolution. Multigene families have proven to be good markers for a better understanding of the variability, organization and evolution of fish species. Three different tandemly-repeated gene families (45S rDNA, 5S rDNA and U2 snDNA) have been studied in Plectorhinchus mediterraneus (Teleostei: Haemulidae), at both molecular and cytogenetic level, to elucidate the taxonomy and evolution of these multigene families, as well as for comparative purposes with other species of the family.

Results

Four different types of 5S rDNA were obtained; two of them showed a high homology with that of Raja asterias, and the putative implication of a horizontal transfer event and its consequences for the organization and evolution of the 5S rDNA have been discussed. The other two types do not resemble any other species, but in one of them a putative tRNA-derived SINE was observed for the first time, which could have implications in the evolution of the 5S rDNA. The ITS-1 sequence was more related to a species of another different genus than to that of the same genus, therefore a revision of the Hamulidae family systematic has been proposed. In the analysis of the U2 snDNA, we were able to corroborate that U2 snDNA and U5 snDNA were linked in the same tandem array, and this has interest for tracing evolutionary lines. The karyotype of the species was composed of 2n = 48 acrocentric chromosomes, and each of the three multigene families were located in different chromosome pairs, thus providing three different chromosomal markers.

Conclusions

Novel data can be extracted from the results: a putative event of horizontal transfer, a possible tRNA-derived SINE linked to one of the four 5S rDNA types characterized, and a linkage between U2 and U5 snDNA. In addition, a revision of the taxonomy of the Haemulidae family has been suggested, and three cytogenetic markers have been obtained. Some of these results have not been described before in any other fish species. New clues about the genome organization and evolution of the multigene families are offered in this study.

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.

Analysis of three multigene families as useful tools in species characterization of two closely-related species, Dicentrarchus labrax, Dicentrarchus punctatus and their hybrids

By analyzing three multigene families, two closely related and commercially important species, Dicentrarchus labrax and Dicentrarchus punctatus, were characterized by cytogenetic and molecular methods. The interspecies hybrid Dicentrarchus labrax (♀) × Dicentrarchus punctatus (♂) was also analyzed. The multigene families studied were the 5S rDNA, 45S rDNA and the U2 snRNA. A microsatellite GTT motif was found within the non transcribed spacers (NTS) of the 5S rDNA from the two species. However, hexanucleotide duplication next to this microsatellite was observed in the D. labrax and hybrid clones, but not in D. punctatus. The U2 snRNA appeared to be linked to the U5 gene and showed two variant sequences, in both D. labrax and D. punctatus. They differed in one insertion/deletion of 7 nucleotides. The first internal transcribed spacer (ITS-1) region showed higher nucleotide variability in D. punctatus than in D. labrax. Nucleotide polymorphism within species and also nucleotide divergence between species were determined in the different gene regions. In a FISH analysis we obtained three chromosomal markers, because the 5S rDNA, 18S rDNA and U2 snRNA probes hybridized each in three different chromosome pairs. Hence none of them was co-localized. The 5S rDNA cluster and U2 snRNA were localized in acrocentric chromosome pairs, while the 18S rRNA gene probe hybridized in a subtelocentric pair. Finally, the usefulness of the results in developing tools for phylogenetic analysis and species identification are discussed in relation to other fish species.

The linked units of 5S rDNA and U1 snDNA of razor shells (Mollusca: Bivalvia: Pharidae)

The linkage between 5S ribosomal DNA and other multigene families has been detected in many eukaryote lineages, but whether it provides any selective advantage remains unclear. In this work, we report the occurrence of linked units of 5S ribosomal DNA (5S rDNA) and U1 small nuclear DNA (U1 snDNA) in 10 razor shell species (Mollusca: Bivalvia: Pharidae) from four different genera. We obtained several clones containing partial or complete repeats of both multigene families in which both types of genes displayed the same orientation. We provide a comprehensive collection of razor shell 5S rDNA clones, both with linked and nonlinked organisation, and the first bivalve U1 snDNA sequences. We predicted the secondary structures and characterised the upstream and downstream conserved elements, including a region at -25 nucleotides from both 5S rDNA and U1 snDNA transcription start sites. The analysis of 5S rDNA showed that some nontranscribed spacers (NTSs) are more closely related to NTSs from other species (and genera) than to NTSs from the species they were retrieved from, suggesting birth-and-death evolution and ancestral polymorphism. Nucleotide conservation within the functional regions suggests the involvement of purifying selection, unequal crossing-overs and gene conversions. Taking into account this and other studies, we discuss the possible mechanisms by which both multigene families could have become linked in the Pharidae lineage. The reason why 5S rDNA is often found linked to other multigene families seems to be the result of stochastic processes within genomes in which its high copy number is determinant.

Sequence characterization and phylogenetic analysis of the 5S ribosomal DNA in species of the family Batrachoididae

5S ribosomal DNA (rDNA) sequences were analyzed in four species belonging to different genera of the fish family Batrachoididae. Several 5S rDNA variants differing in their non-transcribed spacers (NTSs) were found and were grouped into two main types. Two species showed both types of 5S rDNA, whereas the other two species showed only one type. One type of NTS of Amphichthys cryptocentrus showed a high polymorphism due to several deletions and insertions, and phylogenetic analysis showed a between-species clustering of this type of NTS in Amphichthys cryptocentrus. These results suggest a clear differentiation in the model of 5S rDNA evolution of these four species of Batrachoididae, which appear to have been subject to processes of concerted evolution and birth-and-death evolution with purifying selection.

Evolutionary dynamics of the 5S rDNA gene family in the mussel Mytilus: mixed effects of birth-and-death and concerted evolution

In higher eukaryotes, the gene family encoding the 5S ribosomal RNA (5S rRNA) has been used (together with histones) to showcase the archetypal example of a gene family subject to concerted evolution. However, recent studies have revealed conspicuous features challenging the predictions of this model, including heterogeneity of repeat units, the presence of functional 5S gene variants as well as the existence of 5S rDNA divergent pseudogenes lacking traces of homogenization. In the present work, we have broadened the scope in the evolutionary study of ribosomal gene families by studying the 5S rRNA family in mussels, a model organism which stands out among other animals due to the heterogeneity it displays regarding sequence and organization. To this end, 48 previously unknown 5S rDNA units (coding and spacer regions) were sequenced in five mussel species, leading to the characterization of two new types of units (referred to here as small-beta 5S rDNA and gamma-5S rDNA) coexisting in the genome with alpha and beta rDNA units. The intense genetic dynamics of this family is further supported by the first description of an association between gamma-5S rDNA units and tRNA genes. Molecular evolutionary and phylogenetic analyses revealed an extensive lack of homology among spacer sequences belonging to different rDNA types, suggesting the presence of independent evolutionary pathways leading to their differentiation. Overall, our results suggest that the long-term evolution of the 5S rRNA gene family in mussels is most likely mediated by a mixed mechanism involving the generation of genetic diversity through birth-and-death, followed by a process of local homogenization resulting from concerted evolution in order to maintain the genetic identities of the different 5S units, probably after their transposition to independent chromosomal locations.

Evolution of spliceosomal snRNA genes in metazoan animals

While studies of the evolutionary histories of protein families are commonplace, little is known on noncoding RNAs beyond microRNAs and some snoRNAs. Here we investigate in detail the evolutionary history of the nine spliceosomal snRNA families (U1, U2, U4, U5, U6, U11, U12, U4atac, and U6atac) across the completely or partially sequenced genomes of metazoan animals. Representatives of the five major spliceosomal snRNAs were found in all genomes. None of the minor splicesomal snRNAs were detected in nematodes or in the shotgun traces of Oikopleura dioica, while in all other animal genomes at most one of them is missing. Although snRNAs are present in multiple copies in most genomes, distinguishable paralogue groups are not stable over long evolutionary times, although they appear independently in several clades. In general, animal snRNA secondary structures are highly conserved, albeit, in particular, U11 and U12 in insects exhibit dramatic variations. An analysis of genomic context of snRNAs reveals that they behave like mobile elements, exhibiting very little syntenic conservation.

Patterns of evolution in the unique tRNA gene arrays of the genus Entamoeba

Genome sequencing of the protistan parasite Entamoeba histolytica HM-1:IMSS revealed that almost all the tRNA genes are organized into tandem arrays that make up over 10% of the genome. The 25 distinct array units contain up to 5 tRNA genes each and some also encode the 5S RNA. Between adjacent genes in array units are complex short tandem repeats (STRs) resembling microsatellites. To investigate the origins and evolution of this unique gene organization, we have undertaken a genome survey to determine the array unit organization in 4 other species of Entamoeba-Entamoeba dispar, Entamoeba moshkovskii, Entamoeba terrapinae, and Entamoeba invadens-and have explored the STR structure in other isolates of E. histolytica. The genome surveys revealed that E. dispar has the same array unit organization as E. histolytica, including the presence and numerical variation of STRs between adjacent genes. However, the individual repeat sequences are completely different to those in E. histolytica. All other species of Entamoeba studied also have tandem arrays of clustered tRNA genes, but the gene composition of the array units often differs from that in E. histolytica/E. dispar. None of the other species' arrays exhibit the complex STRs between adjacent genes although simple tandem duplications are occasionally seen. The degree of similarity in organization reflects the phylogenetic relationships among the species studied. Within individual isolates of E. histolytica most copies of the array unit are uniform in sequence with only minor variation in the number and organization of the STRs. Between isolates, however, substantial differences in STR number and organization can exist although the individual repeat sequences tend to be conserved. The origin of this unique gene organization in the genus Entamoeba clearly predates the common ancestor of the species investigated to date and their function remains unclear.