The present work represents the first characterization of a clustered histone repetitive unit containing an H1 gene in a bivalve mollusk. To complete the knowledge on the evolutionary history of the histone multigene family in invertebrates, we undertake its characterization in five mussel Mytilus species, as an extension of our previous work on the H1 gene family. We report the quintet H4-H2B-H2A-H3-H1 as the major organization unit in the genome of Mytilus galloprovincialis with two 5S rRNA genes with interspersed nontranscribed spacer segments linked to the unit, which is not justified by their cotranscription with histone genes. Surprisingly, 3' UTR regions of histone genes show two different mRNA termination signals, a stem-loop and a polyadenylation signal, both related to the evolution of histone gene expression patterns throughout the cell cycle. The clustered H1 histones characterized share essential features with "orphon" H1 genes, suggesting a common evolutionary origin for both histone subtypes which is supported by the reconstructed phylogeny for H1 genes. The characterization of histone genes in four additional Mytilus species revealed the presence of strong purifying selection acting among the members of the family. The chromosomal location of most of the core histone genes studied was identified by FISH close to telomeric regions in M. galloprovincialis. Further analysis on nucleotide variation would be necessary to assess if H1 proteins evolve according to the birth-and-death model of evolution and if the effect of the strong purifying selection maintaining protein homogeneity could account for the homologies detected between clustered and "orphon" variants.

Histones are small basic nuclear proteins with critical structural and functional roles in eukaryotic genomes. The H1 multigene family constitutes a very interesting histone class gathering the greatest number of isoforms, with many different arrangements in the genome, including clustered and solitary genes, and showing replication-dependent (RD) or replication-independent (RI) expression patterns. The evolution of H1 histones has been classically explained by concerted evolution through a rapid process of interlocus recombination or gene conversion. Given such intriguing features, we have analyzed the long-term evolutionary pattern of the H1 multigene family through the evaluation of the relative importance of gene conversion, point mutation, and selection in generating and maintaining the different H1 subtypes. We have found the presence of an extensive silent nucleotide divergence, both within and between species, which is always significantly greater than the nonsilent variation, indicating that purifying selection is the major factor maintaining H1 protein homogeneity. The results obtained from phylogenetic analysis reveal that different H1 subtypes are no more closely related within than between species, as they cluster by type in the topologies, and that both RD and RI H1 variants follow the same evolutionary pattern. These findings suggest that H1 histones have not been subject to any significant effect of interlocus recombination or concerted evolution. However, the diversification of the H1 isoforms seems to be enhanced primarily by mutation and selection, where genes are subject to birth-and-death evolution with strong purifying selection at the protein level. This model is able to explain not only the generation and diversification of RD H1 isoforms but also the origin and long-term persistence of orphon RI H1 subtypes in the genome, something that is still unclear, assuming concerted evolution.

Linker histones are a divergent group of histone proteins with an independent evolutionary history in which, besides somatic subtypes, tissue- and differentiation-specific subtypes are included. In the present work H1 histone coding and noncoding segments from five Mytilus mussel species (Mollusca: Bivalvia) widely distributed throughout the world have been determined and characterized. Analysis of promoter regions shows clear homologies among Mytilus H1 genes, sea urchin H1 genes, and vertebrate differentiation-specific H1 subtypes (H5 and H1(o)), all having an H4 box motif in common. The amino acid sequence of the H1 protein central conserved domain is also closely related to that previously defined for the vertebrate divergent subtypes. A phylogenetic tree reconstructed from different H1 genes from several species strengthens the hypothesis of an "orphon" origin for the Mytilus H1 genes, as well as for the H1(o)/H5 genes from vertebrates and the H1D gene from the sea urchin Strongylocentrotus purpuratus, is suggested. As additional data, the average copy number of the H1 genes in the species analyzed was estimated as being 100 to 110 copies per haploid genome, where FISH revealed telomeric chromosomal location for several H1 copies in M. galloprovincialis. The contribution of such proximity to heterochromatic regions over the amount of codon bias detected for H1 genes is discussed.

The hexamer repeat sequence (TTAGGG) n , found at the ends of all vertebrate chromosomes, was previously identified as the main building element of one member of a HindIII satellite DNA family characterized in the genome of the bivalve mollusc Donax trunculus. It was also found in 22 perfect tandem repeats in a cloned junction region juxtaposed to the proper satellite sequence, from which the DNA tract encompassing the clustered tandem copies was excised and subcloned. Here, the chromosomal distribution of (TTAGGG) n sequences in the Donax was studied by the sensitivity to Bal31 exonuclease digestion, fluorescence in situ hybridization (FISH) on metaphase chromosomes and rotating-field gel electrophoresis. To verify the occurrence of the hexamer repeat in the genomes of taxonomically related molluscs and other marine invertebrates, genomic DNA from the mussel Mytilus galloprovincialis and the echinoderm Holothuria tubulosa was also analyzed. The kinetics of Bal31 hydrolysis of high molecular mass DNA from the three marine invertebrates revealed a marked decrease over time of the hybridization with the cloned (TTAGGG) 22 sequence, concomitantly with a progressive shortening of the positively reacting DNA fragments. This revealed a marked susceptibility to exonuclease consistent with terminal positioning on the respective chromosomal DNAs. In full agreement, FISH results with the (TTAGGG) 22 probe showed that the repeat appears located in telomeric regions in all chromosomes of both bivalve molluscs. The presence of (TTAGGG) n repeat tracts in marine invertebrate telomeres points to its wider distribution among eukaryotic organisms and suggests an ancestry older than originally presumed from its vertebrate distinctiveness.

A study of nucleotide sequence variation of 5S ribosomal DNA from six Ensis species revealed that several 5S ribosomal DNA variants, based on differences in their nontranscribed spacers (NTS), occur in Ensis genomes. The 5S rRNA gene was not very polymorphic, compared with the NTS region. The phylogenetic analyses performed showed a between-species clustering of 5S ribosomal DNA variants. Sequence divergence levels between variants were very large, revealing a lack of sequence homogenization. These results strongly suggest that the long-term evolution of Ensis 5S ribosomal DNA is driven by birth-and-death processes and selection.

Satellite DNAs represent a large portion of all high eukaryotic genomes. They consist of numerous very similar repeated sequences, tandemly arranged in large clusters up to 100 million base pairs in length, usually located in the heterochromatic parts of chromosomes. The biological significance of satDNAs is still under discussion, but most of their proposed functions are related to heterochromatin and/or centromere formation and function. Because information about the structure of reptilian satDNA is far from exhaustive, we present a molecular and cytogenetic characterization of two satDNA families in four lacertid species. Two families of tandemly repeated DNAs, namely TaqI and HindIII satDNAs, have been cloned and sequenced from four species belonging to the genus Iberolacerta. These satDNAs are characterized by a monomer length of 171-188 and 170-172 bp, and by an AT content of 60.5% and 58.1%, respectively. FISH experiments with TaqI satDNA probe produced bright signals in pericentromeric regions of a subset of chromosomes whereas all the centromeres were marked by HindIII probe. The results obtained in this study suggest that chromosome location and abundance of satDNAs influence the evolution of these elements, with centromeric families evolving tenfold faster than interstitial/pericentromeric ones. Such different rates render different satellites useful for phylogenetic investigation at different taxonomic ranks.

Satellite DNAs compose a large portion of all higher eukaryotic genomes. The turnover of these highly repetitive sequences is an important element in genome organization and evolution. However, information about the structure and dynamics of reptilian satellite DNA is still scarce. Two satellite DNA families, HindIII and TaqI, have been previously characterized in four species of the genus Iberolacerta. These families showed different chromosomal locations, abundances, and evolutionary rates. Here, we extend the study of both satellite DNAs (satDNAs) to the remaining Iberolacerta species, with the aim to investigate the patterns of variability and factors influencing the evolution of these repetitive sequences. Our results revealed disparate patterns but also common traits in the evolutionary histories of these satellite families: (i) each satellite DNA is made up of a library of monomer variants or subfamilies shared by related species; (ii) species-specific profiles of satellite repeats are shaped by expansions and/or contractions of different variants from the library; (iii) different turnover rates, even among closely related species, result in great differences in overall sequence homogeneity and in concerted or non-concerted evolution patterns, which may not reflect the phylogenetic relationships among taxa. Contrasting turnover rates are possibly related to genomic constraints such as karyotype architecture and the interspersed organization of diverging repeat variants in satellite arrays. Moreover, rapid changes in copy number, especially in the centromeric HindIII satDNA, may have been associated with chromosomal rearrangements and even contributed to speciation within Iberolacerta.

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.