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

Assay Optimization Can Equalize the Sensitivity of Real-Time PCR with ddPCR for Detection of Helicoverpa armigera (Lepidoptera: Noctuidae) in Bulk Samples

Simple Summary

Invasive species are a constant threat to agriculture throughout the world against which early detection is one of the primary defenses. The Old World bollworm is one of the most important invasive agricultural pests in the world. While historically absent from the Americas, this species was first found in South America in 2013 and poses an ongoing threat of spreading into North America. Surveys are conducted each year, which result in hundreds or thousands of traps that must be screened for this species. Unfortunately, the most common non-target is the native corn earworm, which is nearly identical morphologically to the Old World bollworm and cannot be easily separated. Molecular methods have been developed to screen these trap samples, but the required equipment is expensive and not commonly available. This study details improvements to current molecular methods that will allow for screening of bulk trap samples using standard laboratory instruments and protocols. The ability to perform these methods in nearly any molecular biology lab will greatly enhance our ability to detect and exclude this important pest.

Abstract

Helicoverpa armigera (Hübner) is one of the most important agricultural pests in the world. This historically Old World species was first reported in Brazil in 2013 and has since spread throughout much of South America and into the Caribbean. Throughout North America, H. armigera surveys are ongoing to detect any incursions. Each trap is capable of capturing hundreds of native Helicoverpa zea (Boddie). The two species cannot be separated without genitalic dissection or molecular methods. A ddPCR assay is currently used to screen large trap samples, but this equipment is relatively uncommon and expensive. Here, we optimized a newly designed assay for accurate and repeatable detection of H. armigera in bulk samples across both ddPCR and less costly, and more common, real-time PCR methods. Improvements over previously designed assays were sought through multiple means. Our results suggest bulk real-time PCR assays can be improved through changes in DNA extraction and purification, so that real-time PCR can be substituted for ddPCR in screening projects. While ddPCR remains a more sensitive method for detection of H. armigera in bulk samples, the improvements in assay design, DNA extraction, and purification presented here also enhance assay performance over previous protocols.

Cytogenomic analysis unveils mixed molecular evolution and recurrent chromosomal rearrangements shaping the multigene families on Schistocerca grasshopper genomes

Multigene families are essential components of eukaryotic genomes and play key roles either structurally and functionally. Their modes of evolution remain elusive even in the era of genomics, because multiple multigene family sequences coexist in genomes, particularly in large repetitive genomes. Here, we investigate how the multigene families 18S rDNA, U2 snDNA, and H3 histone evolved in 10 species of Schistocerca grasshoppers with very large and repeat-enriched genomes. Using sequenced genomes and fluorescence in situ hybridization mapping, we find substantial differences between species, including the number of chromosomal clusters, changes in sequence abundance and nucleotide composition, pseudogenization, and association with transposable elements (TEs). The intragenomic analysis of Schistocerca gregaria using long-read sequencing and genome assembly unveils conservation for H3 histone and recurrent pseudogenization for 18S rDNA and U2 snDNA, likely promoted by association with TEs and sequence truncation. Remarkably, TEs were frequently associated with truncated copies, were also among the most abundant in the genome, and revealed signatures of recent activity. Our findings suggest a combined effect of concerted and birth-and-death models driving the evolution of multigene families in Schistocerca over the last 8 million years, and the occurrence of intra- and interchromosomal rearrangements shaping their chromosomal distribution. Despite the conserved karyotype in Schistocerca, our analysis highlights the extensive reorganization of repetitive DNAs in Schistocerca, contributing to the advance of comparative genomics for this important grasshopper genus.

Cytogenetic characterization and mapping of the repetitive DNAs in Cycloramphus bolitoglossus (Werner, 1897): More clues for the chromosome evolution in the genus Cycloramphus (Anura, Cycloramphidae)

Cycloramphus bolitoglossus (Werner, 1897) is a rare species with a low population density in the Serra do Mar region of Paraná and Santa Catarina, in southern Brazil. Currently, it has been assigned to the Near Threatened (NT) category in the Brazilian List of Endangered Animal Species. Here, we described the karyotype of this species for the first time and investigated the patterns of some repetitive DNA classes in the chromosomes using molecular cytogenetic approaches. We isolated, sequenced and mapped the 5S rDNA and the satellite DNA PcP190 of C. bolitoglossus, as well as mapped the telomeric sequences and seven microsatellites motifies [(GA)15, (CA)15, (GACA)4, (GATA)8, (CAG)10, (CGC)10, and (GAA)]10. Cycloramphus bolitoglossus has 2n = 26 chromosomes and a fundamental number (FN) equal to 52, with a highly conserved karyotype compared to other genus members. Comparative cytogenetic under the phylogenetic context of genus allowed evolutionary interpretations of the morphological changes in the homologs of pairs 1, 3, and 6 along with the evolutionary history of Cycloramphus. Two subtypes of 5S rDNA type II were isolated in C. bolitoglossus genome, and several comparative analysis suggests mixed effects of concerted and birth-and-death evolution acting in this repetitive DNA. The 5S rDNA II subtype "a" and "b" was mapped on chromosome 1. However, their different position along chromosome 1 provide an excellent chromosome marker for future studies. PcP190 satellite DNA, already reported for species of the families Hylidae, Hylodidae, Leptodactylidae, and Odontophrynidae, is scattered throughout the C. bolitoglossus genome, and even non-heterochromatic regions showed hybridization signals using the PcP190 probe. Molecular analysis suggests that PcP190 satellite DNA exhibit a high-level of homogenization of this sequence in the genome of C. bolitoglossus. The PcP190 satDNA from C. bolitoglossus represents a novel sequence group, compared to other anurans, based on its hypervariable region. Overall, the present data on repetitive DNA sequences showed pseudogenization evidence and corroborated the hypothesis of the emergence of satDNA from rDNA 5S clusters. These two arguments that reinforced the importance of the birth-and-death evolutionary model to explain 5S rDNA patterns found in anuran genomes.

Evolution and diversity of alpha-carbonic anhydrases in the mantle of the Mediterranean mussel (Mytilus galloprovincialis)

The α-carbonic anhydrases (α-CAs) are a large and ancient group of metazoan-specific enzymes. They generate bicarbonate from metabolic carbon dioxide and through calcium carbonate crystal formation play a key role in the regulation of mineralized structures. To better understand how α-CAs contribute to shell mineralization in the marine Mediterranean mussel (Mytilus galloprovincialis) we characterized them in the mantle. Phylogenetic analysis revealed that mollusc α-CA evolution was affected by lineage and species-specific events. Ten α-CAs were found in the Mediterranean mussel mantle and the most abundant form was named, MgNACR, as it grouped with oyster nacreins (NACR). Exposure of the Mediterranean mussel to reduced water salinity (18 vs 37 ppt), caused a significant reduction (p < 0.05) in mantle esterase activity and MgNACR transcript abundance (p < 0.05). Protonograms revealed multiple proteins in the mantle with α-CA hydratase activity and mapped to a protein with a similar size to that deduced for monomeric MgNACR. Our data indicate that MgNACR is a major α-CA enzyme in mantle and that by homology with oyster nacreins likely regulates mussel shell production. We propose that species-dependent α-CA evolution may contribute to explain the diversity of bivalve shell structures and their vulnerability to environmental changes.

Chromosomal and molecular characterization of 5S rRNA genes in the North American abalones Haliotis rufescens Swainson (red abalone) and H. fulgens Philippi (blue abalone)

Abalone is an extremely valuable food source derived from cultured and wild animals, the later from populations under intense fishing exploitation and of high conservation value. As part of a long-term study to characterize genes from abalone that can be used as markers for hybrids certification, we characterised 5S ribosomal DNA (5S rDNA) in red abalone (Haliotis rufescens) and blue abalone (H. fulgens). The 5S rDNA arrays occur to a single pair of metacentric chromosomes at interstitial positions in both species. Two types of 5S genes were found, named types I and II, each associated with different non-transcribed spacer (NTS) sequences. The structure of the 5S rRNA genes and the NTS indicate incomplete homogenisation of the 5S rDNA arrays. The divergence of the 5S genes between species provide polymorphisms which can be used to distinguish red from blue abalone in forensic analysis of commercial production.

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.

Evolutionary Dynamics of Copy Number and Meiotic Recombination in Murine 5S rDNA: Possible Involvement of Natural Selection

We investigated evolutionary trends of the 5S ribosomal RNA gene in the house mouse, Mus musculus. First, we assessed the 5S cluster and copy numbers in eight laboratory strains by pulsed-field gel electrophoresis. The copy numbers in seven lines were estimated to be around 130-170 copies per cluster, with 63 copies in the remaining strain, implying that the copy number can change drastically and has been maintained under certain evolutionary constraints at ~ 140 copies. Second, we addressed the frequency of meiotic recombination mediated by the 5S cluster by performing a mating experiment with laboratory strains, and found that the 5S cluster did not accelerate recombination events. Third, we surveyed recombination events of the 5S-containing chromosome region in wild mice from the Japanese Islands, where the two subspecies lineages, M. m. castaneus and M. m. musculus, are historically mingled, and found that the influence of the 5S cluster on meiotic recombination was limited. Finally, we examined the nucleotide diversity of six genes in the neighboring regions of the 5S cluster and found reduced genetic diversity in the regions on both sides of the cluster, suggesting the involvement of either positive or background selection in the population-level sequence similarity of the 5S clusters. Therefore, the mouse 5S genes are considered to be evolving toward sequence similarity within a given cluster by certain intrachromosomal mechanisms and toward sharing of a specific 5S cluster within a population by certain selective processes.

Evolution of blue-flowered species of genus Linum based on high-throughput sequencing of ribosomal RNA genes

BACKGROUND

The species relationships within the genus Linum have already been studied several times by means of different molecular and phylogenetic approaches. Nevertheless, a number of ambiguities in phylogeny of Linum still remain unresolved. In particular, the species relationships within the sections Stellerolinum and Dasylinum need further clarification. Also, the question of independence of the species of the section Adenolinum still remains unanswered. Moreover, the relationships of L. narbonense and other species of the section Linum require further clarification. Additionally, the origin of tetraploid species of the section Linum (2n = 30) including the cultivated species L. usitatissimum has not been explored. The present study examines the phylogeny of blue-flowered species of Linum by comparisons of 5S rRNA gene sequences as well as ITS1 and ITS2 sequences of 35S rRNA genes.

RESULTS

High-throughput sequencing has been used for analysis of multicopy rRNA gene families. In addition to the molecular phylogenetic analysis, the number and chromosomal localization of 5S and 35S rDNA sites has been determined by FISH. Our findings confirm that L. stelleroides forms a basal branch from the clade of blue-flowered flaxes which is independent of the branch formed by species of the sect. Dasylinum. The current molecular phylogenetic approaches, the cytogenetic analysis as well as different genomic DNA fingerprinting methods applied previously did not discriminate certain species within the sect. Adenolinum. The allotetraploid cultivated species L. usitatissimum and its wild ancestor L. angustifolium (2n = 30) could originate either as the result of hybridization of two diploid species (2n = 16) related to the modern L. gandiflorum and L. decumbens, or hybridization of a diploid species (2n = 16) and a diploid ancestor of modern L. narbonense (2n = 14).

CONCLUSIONS

High-throughput sequencing of multicopy rRNA gene families allowed us to make several adjustments to the phylogeny of blue-flowered flax species and also reveal intra- and interspecific divergence of the rRNA gene sequences.

Variations in 5S rDNAs in diploid and tetraploid offspring of red crucian carp × common carp

BACKGROUND

The allotetraploid hybrid fish (4nAT) that was created in a previous study through an intergeneric cross between red crucian carp (Carassius auratus red var., ♀) and common carp (Cyprinus carpio L., ♂) provided an excellent platform to investigate the effect of hybridization and polyploidization on the evolution of 5S rDNA. The 5S rDNAs of paternal common carp were made up of a coding sequence (CDS) and a non-transcribed spacer (NTS) unit, and while the 5S rDNAs of maternal red crucian carp contained a CDS and a NTS unit, they also contained a variable number of interposed regions (IPRs). The CDSs of the 5S rDNAs in both parental fishes were conserved, while their NTS units seemed to have been subjected to rapid evolution.

RESULTS

The diploid hybrid 2nF₁ inherited all the types of 5S rDNAs in both progenitors and there were no signs of homeologous recombination in the 5S rDNAs of 2nF₁ by sequencing of PCR products. We obtained two segments of 5S rDNA with a total length of 16,457 bp from allotetraploid offspring 4nAT through bacterial artificial chromosome (BAC) sequencing. Using this sequence together with the 5S rDNA sequences amplified from the genomic DNA of 4nAT, we deduced that the 5S rDNAs of 4nAT might be inherited from the maternal progenitor red crucian carp. Additionally, the IPRs in the 5S rDNAs of 4nAT contained A-repeats and TA-repeats, which was not the case for the IPRs in the 5S rDNAs of 2nF₁. We also detected two signals of a 200-bp fragment of 5S rDNA in the chromosomes of parental progenitors and hybrid progenies by fluorescence in situ hybridization (FISH).

CONCLUSIONS

We deduced that during the evolution of 5S rDNAs in different ploidy hybrid fishes, interlocus gene conversion events and tandem repeat insertion events might occurred in the process of polyploidization. This study provided new insights into the relationship among the evolution of 5S rDNAs, hybridization and polyploidization, which were significant in clarifying the genome evolution of polyploid fish.

U1 snDNA clusters in grasshoppers: chromosomal dynamics and genomic organization

The spliceosome, constituted by a protein set associated with small nuclear RNA (snRNA), is responsible for mRNA maturation through intron removal. Among snRNA genes, U1 is generally a conserved repetitive sequence. To unveil the chromosomal/genomic dynamics of this multigene family in grasshoppers, we mapped U1 genes by fluorescence in situ hybridization in 70 species belonging to the families Proscopiidae, Pyrgomorphidae, Ommexechidae, Romaleidae and Acrididae. Evident clusters were observed in all species, indicating that, at least, some U1 repeats are tandemly arrayed. High conservation was observed in the first four families, with most species carrying a single U1 cluster, frequently located in the third or fourth longest autosome. By contrast, extensive variation was observed among Acrididae, from a single chromosome pair carrying U1 to all chromosome pairs carrying it, with occasional occurrence of two or more clusters in the same chromosome. DNA sequence analysis in Eyprepocnemis plorans (species carrying U1 clusters on seven different chromosome pairs) and Locusta migratoria (carrying U1 in a single chromosome pair) supported the coexistence of functional and pseudogenic lineages. One of these pseudogenic lineages was truncated in the same nucleotide position in both species, suggesting that it was present in a common ancestor to both species. At least in E. plorans, this U1 snDNA pseudogenic lineage was associated with 5S rDNA and short interspersed elements (SINE)-like mobile elements. Given that we conclude in grasshoppers that the U1 snDNA had evolved under the birth-and-death model and that its intragenomic spread might be related with mobile elements.

Karyotypic diversification in Mytilus mussels (Bivalvia: Mytilidae) inferred from chromosomal mapping of rRNA and histone gene clusters

Background

Mussels of the genus Mytilus present morphologically similar karyotypes that are presumably conserved. The absence of chromosome painting probes in bivalves makes difficult verifying this hypothesis. In this context, we comparatively mapped ribosomal RNA and histone gene families on the chromosomes of Mytilus edulis, M. galloprovincialis, M. trossulus and M. californianus by fluorescent in situ hybridization (FISH).

Results

Major rRNA, core and linker histone gene clusters mapped to different chromosome pairs in the four taxa. In contrast, minor rRNA gene clusters showed a different behavior. In all Mytilus two of the 5S rDNA clusters mapped to the same chromosome pair and one of them showed overlapping signals with those corresponding to one of the histone H1 gene clusters. The overlapping signals on mitotic chromosomes became a pattern of alternate 5S rRNA and linker histone gene signals on extended chromatin fibers. Additionally, M. trossulus showed minor and major rDNA clusters on the same chromosome pair.

Conclusion

The results obtained suggest that at least some of the chromosomes bearing these sequences are orthologous and that chromosomal mapping of rRNA and histone gene clusters could be a good tool to help deciphering some of the many unsolved questions in the systematic classification of Mytilidae.

Disparate molecular evolution of two types of repetitive DNAs in the genome of the grasshopper Eyprepocnemis plorans

Wide arrays of repetitive DNA sequences form an important part of eukaryotic genomes. These repeats appear to evolve as coherent families, where repeats within a family are more similar to each other than to other orthologous representatives in related species. The continuous homogenization of repeats, through selective and non-selective processes, is termed concerted evolution. Ascertaining the level of variation between repeats is crucial to determining which evolutionary model best explains the homogenization observed for these sequences. Here, for the grasshopper Eyprepocnemis plorans, we present the analysis of intragenomic diversity for two repetitive DNA sequences (a satellite DNA (satDNA) and the 45S rDNA) resulting from the independent microdissection of several chromosomes. Our results show different homogenization patterns for these two kinds of paralogous DNA sequences, with a high between-chromosome structure for rDNA but no structure at all for the satDNA. This difference is puzzling, considering the adjacent localization of the two repetitive DNAs on paracentromeric regions in most chromosomes. The disparate homogenization patterns detected for these two repetitive DNA sequences suggest that several processes participate in the concerted evolution in E. plorans, and that these mechanisms might not work as genome-wide processes but rather as sequence-specific ones.

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.

Diversity and recombination of dispersed ribosomal DNA and protein coding genes in microsporidia

Microsporidian strains are usually classified on the basis of their ribosomal DNA (rDNA) sequences. Although rDNA occurs as multiple copies, in most non-microsporidian species copies within a genome occur as tandem arrays and are homogenised by concerted evolution. In contrast, microsporidian rDNA units are dispersed throughout the genome in some species, and on this basis are predicted to undergo reduced concerted evolution. Furthermore many microsporidian species appear to be asexual and should therefore exhibit reduced genetic diversity due to a lack of recombination. Here, DNA sequences are compared between microsporidia with different life cycles in order to determine the effects of concerted evolution and sexual reproduction upon the diversity of rDNA and protein coding genes. Comparisons of cloned rDNA sequences between microsporidia of the genus Nosema with different life cycles provide evidence of intragenomic variability coupled with strong purifying selection. This suggests a birth and death process of evolution. However, some concerted evolution is suggested by clustering of rDNA sequences within species. Variability of protein-coding sequences indicates that considerable intergenomic variation also occurs between microsporidian cells within a single host. Patterns of variation in microsporidian DNA sequences indicate that additional diversity is generated by intragenomic and/or intergenomic recombination between sequence variants. The discovery of intragenomic variability coupled with strong purifying selection in microsporidian rRNA sequences supports the hypothesis that concerted evolution is reduced when copies of a gene are dispersed rather than repeated tandemly. The presence of intragenomic variability also renders the use of rDNA sequences for barcoding microsporidia questionable. Evidence of recombination in the single-copy genes of putatively asexual microsporidia suggests that these species may undergo cryptic sexual reproduction, a possibility with profound implications for the evolution of virulence, host range and drug resistance in these species.

Evidence for 5S rDNA horizontal transfer in the toadfish Halobatrachus didactylus (Schneider, 1801) based on the analysis of three multigene families

BACKGROUND

The Batrachoididae family is a group of marine teleosts that includes several species with more complicated physiological characteristics, such as their excretory, reproductive, cardiovascular and respiratory systems. Previous studies of the 5S rDNA gene family carried out in four species from the Western Atlantic showed two types of this gene in two species but only one in the other two, under processes of concerted evolution and birth-and-death evolution with purifying selection. Here we present results of the 5S rDNA and another two gene families in Halobatrachus didactylus, an Eastern Atlantic species, and draw evolutionary inferences regarding the gene families. In addition we have also mapped the genes on the chromosomes by two-colour fluorescence in situ hybridization (FISH).

RESULTS

Two types of 5S rDNA were observed, named type α and type β. Molecular analysis of the 5S rDNA indicates that H. didactylus does not share the non-transcribed spacer (NTS) sequences with four other species of the family; therefore, it must have evolved in isolation. Amplification with the type β specific primers amplified a specific band in 9 specimens of H. didactylus and two of Sparus aurata. Both types showed regulatory regions and a secondary structure which mark them as functional genes. However, the U2 snRNA gene and the ITS-1 sequence showed one electrophoretic band and with one type of sequence. The U2 snRNA sequence was the most variable of the three multigene families studied. Results from two-colour FISH showed no co-localization of the gene coding from three multigene families and provided the first map of the chromosomes of the species.

CONCLUSIONS

A highly significant finding was observed in the analysis of the 5S rDNA, since two such distant species as H. didactylus and Sparus aurata share a 5S rDNA type. This 5S rDNA type has been detected in other species belonging to the Batrachoidiformes and Perciformes orders, but not in the Pleuronectiformes and Clupeiformes orders. Two hypotheses have been outlined: one is the possible vertical permanence of the shared type in some fish lineages, and the other is the possibility of a horizontal transference event between ancient species of the Perciformes and Batrachoidiformes orders. This finding opens a new perspective in fish evolution and in the knowledge of the dynamism of the 5S rDNA. Cytogenetic analysis allowed some evolutionary trends to be roughed out, such as the progressive change in the U2 snDNA and the organization of (GATA)n repeats, from dispersed to localized in one locus. The accumulation of (GATA)n repeats in one chromosome pair could be implicated in the evolution of a pair of proto-sex chromosomes. This possibility could situate H. didactylus as the most highly evolved of the Batrachoididae family in terms of sex chromosome biology.

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 and functional nuclear rDNA genes in the abalone Haliotis tuberculata: tissue differential expression

Analysis of the 18S rDNA sequences of Haliotis tuberculata tuberculata and H. t. coccinea subtaxa identified two different types of 18S rDNA genes and ITS1 regions. These two different genes were also detected in H. marmorata, H. rugosa and H. diversicolor that are separated from H. tuberculata by 5-65 mya. The mean divergence value between type I and type II sequences ranged from 7.25% for 18S to 80% for ITS1. ITS1 type II is homologous with the ITS1 consensus sequences published for many abalone species, whereas ITS1 type I presented only minor homology with a unique database entry for H. iris ITS1. A phylogenetic analysis makes a clear separation between type I and type II ITS1 sequences and supports grouping H. t. tuberculata, H. t. coccinea and H. marmorata together. The two subtaxa do not show any significant differences between the homologous 18S rDNA sequences. A general structure of the ITS1 transcript was proposed, with four major helices for the two types. The two genes were expressed and, for the first time, a putative differential expression of ITS1 type I was detected in the gills, digestive gland and gonads whereas ITS1 type II was expressed in all tissues.

Molecular organization and phylogenetic analysis of 5S rDNA in crustaceans of the genus Pollicipes reveal birth-and-death evolution and strong purifying selection

BACKGROUND

The 5S ribosomal DNA (5S rDNA) is organized in tandem arrays with repeat units that consist of a transcribing region (5S) and a variable nontranscribed spacer (NTS), in higher eukaryotes. Until recently the 5S rDNA was thought to be subject to concerted evolution, however, in several taxa, sequence divergence levels between the 5S and the NTS were found higher than expected under this model. So, many studies have shown that birth-and-death processes and selection can drive the evolution of 5S rDNA. In analyses of 5S rDNA evolution is found several 5S rDNA types in the genome, with low levels of nucleotide variation in the 5S and a spacer region highly divergent. Molecular organization and nucleotide sequence of the 5S ribosomal DNA multigene family (5S rDNA) were investigated in three Pollicipes species in an evolutionary context.

RESULTS

The nucleotide sequence variation revealed that several 5S rDNA variants occur in Pollicipes genomes. They are clustered in up to seven different types based on differences in their nontranscribed spacers (NTS). Five different units of 5S rDNA were characterized in P. pollicipes and two different units in P. elegans and P. polymerus. Analysis of these sequences showed that identical types were shared among species and that two pseudogenes were present. We predicted the secondary structure and characterized the upstream and downstream conserved elements. Phylogenetic analysis showed an among-species clustering pattern of 5S rDNA types.

CONCLUSIONS

These results suggest that the evolution of Pollicipes 5S rDNA is driven by birth-and-death processes with strong purifying selection.

Structure and molecular evolution of the ribosomal DNA external transcribed spacer in the cockroach genus Blattella

The ribosomal DNA (rDNA) cluster of insects contains several hundred repeating structural-functional units and, therefore, is a typical example of a multigene family. Eukaryotic ribosomal RNA (rRNA) genes (18S, 5.8S, and 28S like) are arranged in tandemly repeated clusters in the nucleolus organizers, separated by several spacers, namely the nontranscribed spacer, the external transcribed spacer (ETS), and the internal transcribed spacers. The nucleotide sequences of the ETS of the three closely related Blattella cockroach species, Blattella germanica (Linnaeus, 1767), Blattella asahinai (Mizukubo, 1981), and Blattella lituricollis (Walker, 1868), were determined and compared. The three species had relatively similar ETS lengths, and sequence differences among them could be explained by two types of rearrangements, namely deletions of subrepeats and nucleotide substitutions. Minor ETS variants in B. germanica differed from the major variant in the same way that the major ETS variants of the three Blattella species differed from each other. Concerted evolution and the birth-and-death models, which are often invoked to explain the diversity and evolution of the multigene families of rDNA clusters, are discussed in the light of our data. A new model is proposed to explain the evolutionary reorganization of the ETS region: evolution of rDNA by "magnification-and-fixation" is characterized by magnification of minor subrepeats, which become adaptive in a new rapidly changed environment, and subsequent fixation of this variant type as a major component of the multigene family of a new species.

The 5S rDNA family evolves through concerted and birth-and-death evolution in fish genomes: an example from freshwater stingrays

Background

Ribosomal 5S genes are well known for the critical role they play in ribosome folding and functionality. These genes are thought to evolve in a concerted fashion, with high rates of homogenization of gene copies. However, the majority of previous analyses regarding the evolutionary process of rDNA repeats were conducted in invertebrates and plants. Studies have also been conducted on vertebrates, but these analyses were usually restricted to the 18S, 5.8S and 28S rRNA genes. The recent identification of divergent 5S rRNA gene paralogs in the genomes of elasmobranches and teleost fishes indicate that the eukaryotic 5S rRNA gene family has a more complex genomic organization than previously thought. The availability of new sequence data from lower vertebrates such as teleosts and elasmobranches enables an enhanced evolutionary characterization of 5S rDNA among vertebrates.

Results

We identified two variant classes of 5S rDNA sequences in the genomes of Potamotrygonidae stingrays, similar to the genomes of other vertebrates. One class of 5S rRNA genes was shared only by elasmobranches. A broad comparative survey among 100 vertebrate species suggests that the 5S rRNA gene variants in fishes originated from rounds of genome duplication. These variants were then maintained or eliminated by birth-and-death mechanisms, under intense purifying selection. Clustered multiple copies of 5S rDNA variants could have arisen due to unequal crossing over mechanisms. Simultaneously, the distinct genome clusters were independently homogenized, resulting in the maintenance of clusters of highly similar repeats through concerted evolution.

Conclusions

We believe that 5S rDNA molecular evolution in fish genomes is driven by a mixed mechanism that integrates birth-and-death and concerted evolution.

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.