Protamines and the sperm nuclear basic proteins Pandora's Box of insects.

Insects are the largest group of animals when it comes to the number and diversity of species. Yet, with the exception of Drosophila, no information is currently available on the primary structure of their sperm nuclear basic proteins (SNBPs). This paper represents the first attempt in this regard and provides information about six species of Neoptera: Poecillimon thessalicus, Graptosaltria nigrofuscata, Apis mellifera, Nasonia vitripennis, Parachauliodes continentalis, and Tribolium castaneum. The SNBPs of these species were characterized by acetic acid urea gel electrophoresis (AU-PAGE) and high-performance liquid chromatography fractionated. Protein sequencing was obtained using a combination of mass spectrometry sequencing, Edman N-terminal degradation sequencing and genome mining. While the SNBPs of several of these species exhibit a canonical arginine-rich protamine nature, a few of them exhibit a protamine-like composition. They appear to be the products of extensive cleavage processing from a precursor protein which are sometimes further processed by other post-translational modifications that are likely involved in the chromatin transitions observed during spermiogenesis in these organisms.

Evolution of high mobility group nucleosome-binding proteins and its implications for vertebrate chromatin specialization.

High mobility group (HMG)-N proteins are a family of small nonhistone proteins that bind to nucleosomes (N). Despite the amount of information available on their structure and function, there is an almost complete lack of information on the molecular evolutionary mechanisms leading to their exclusive differentiation. In the present work, we provide evidence suggesting that HMGN lineages constitute independent monophyletic groups derived from a common ancestor prior to the diversification of vertebrates. Based on observations of the functional diversification across vertebrate HMGN proteins and on the extensive silent nucleotide divergence, our results suggest that the long-term evolution of HMGNs occurs under strong purifying selection, resulting from the lineage-specific functional constraints of their different protein domains. Selection analyses on independent lineages suggest that their functional specialization was mediated by bursts of adaptive selection at specific evolutionary times, in a small subset of codons with functional relevance-most notably in HMGN1, and in the rapidly evolving HMGN5. This work provides useful information to our understanding of the specialization imparted on chromatin metabolism by HMGNs, especially on the evolutionary mechanisms underlying their functional differentiation in vertebrates.

Characterization of mussel H2A.Z.2: a new H2A.Z variant preferentially expressed in germinal tissues from Mytilus.

Histones are the fundamental constituents of the eukaryotic chromatin, facilitating the physical organization of DNA in chromosomes and participating in the regulation of its metabolism. The H2A family displays the largest number of variants among core histones, including the renowned H2A.X, macroH2A, H2A.B (Bbd), and H2A.Z. This latter variant is especially interesting because of its regulatory role and its differentiation into 2 functionally divergent variants (H2A.Z.1 and H2A.Z.2), further specializing the structure and function of vertebrate chromatin. In the present work we describe, for the first time, the presence of a second H2A.Z variant (H2A.Z.2) in the genome of a non-vertebrate animal, the mussel Mytilus. The molecular and evolutionary characterization of mussel H2A.Z.1 and H2A.Z.2 histones is consistent with their functional specialization, supported on sequence divergence at promoter and coding regions as well as on varying gene expression patterns. More precisely, the expression of H2A.Z.2 transcripts in gonadal tissue and its potential upregulation in response to genotoxic stress might be mirroring the specialization of this variant in DNA repair. Overall, the findings presented in this work complement recent reports describing the widespread presence of other histone variants across eukaryotes, supporting an ancestral origin and conserved role for histone variants in chromatin.

Comparative structure of vertebrate sperm chromatin.

A consistent feature of sperm nuclei is its exceptionally compact state in comparison with somatic nuclei. Here, we have examined the structural organization of sperm chromatin from representatives of three vertebrate lineages, bony fish (Danio rerio), birds (Gallus gallus domesticus) and mammals (Mus musculus) using light and transmission electron microscopy (TEM). Although the three sperm nuclei are all highly compact, they differ in morphology and in the complement of compaction-inducing proteins. Whereas zebrafish sperm retain somatic histones and a nucleosomal organization, in the rooster and mouse, histones are largely replaced by small, arginine-rich protamines. In contrast to the mouse, the rooster protamine contains no cysteine residues and lacks the potential stabilizing effects of S-S bonds. Protamine driven chromatin compaction results in a stable, highly condensed chromatin, markedly different from the somatic nucleosome-based beads-on-a-string architecture, but its structure remains poorly understood. When prepared gently for whole mount TEM, the rooster and mouse sperm chromatin reveal striking rod-like units 40-50 nm in width. Also present in the mouse, which has very flattened sperm nuclei, but not rooster, where nuclei take the form of elongated cylinders, are toroidal shaped structures, with an external diameter of about 90 nm. In contrast, similarly prepared zebrafish sperm exhibit nucleosomal chromatin. We also examined the early stages in the binding of salmine (the salmon protamine) to defined sequence DNA. These images suggest an initial side-by-side binding of linear DNA-protamine complexes leading to the nucleation of thin, flexible rods with the potential to bend, allowing the ends to come into contact and fuse to form toroidal structures. We discuss the relationship between these in vitro observations and the rods and toroids seen in nuclei, and suggest an explanation for the apparent absence of these structures in TEM images of fully condensed sperm nuclei.

The CHROMEVALOA database: a resource for the evaluation of Okadaic Acid contamination in the marine environment based on the chromatin-associated transcriptome of the mussel Mytilus galloprovincialis.

Okadaic Acid (OA) constitutes the main active principle in Diarrhetic Shellfish Poisoning (DSP) toxins produced during Harmful Algal Blooms (HABs), representing a serious threat for human consumers of edible shellfish. Furthermore, OA conveys critical deleterious effects for marine organisms due to its genotoxic potential. Many efforts have been dedicated to OA biomonitoring during the last three decades. However, it is only now with the current availability of detailed molecular information on DNA organization and the mechanisms involved in the maintenance of genome integrity, that a new arena starts opening up for the study of OA contamination. In the present work we address the links between OA genotoxicity and chromatin by combining Next Generation Sequencing (NGS) technologies and bioinformatics. To this end, we introduce CHROMEVALOAdb, a public database containing the chromatin-associated transcriptome of the mussel Mytilus galloprovincialis (a sentinel model organism) in response to OA exposure. This resource constitutes a leap forward for the development of chromatin-based biomarkers, paving the road towards the generation of powerful and sensitive tests for the detection and evaluation of the genotoxic effects of OA in coastal areas.

Birth-and-death long-term evolution promotes histone H2B variant diversification in the male germinal cell line.

The rich diversity within each of the five histone families (H1, H2A, H2B, H3, and H4) can hardly be reconciled with the notion of homogenizing evolution. The prevalence of birth-and-death long-term evolution over concerted evolution has already been demonstrated in the linker histone H1 family as well as for the H2A, H3, and H4 core histone families. However, information about histone H2B is lacking. In the present work, we have analyzed the diversity of the members of this histone family across different eukaryotic genomes and have characterized the mechanisms involved in their long-term evolution. Our results reveal that, quite in contrast with other histones, H2B variants are subject to a very rapid process of diversification that primarily affects the male germinal cell lineage and involves their functional specialization probably as a consequence of neofunctionalization and subfunctionalization events after gene duplication. The overall parallelism observed between the molecular phylogenies and the relationships among the electrostatic potentials of the different variants suggests that the latter may have played a major structural selective constraint during H2B evolution. It thus seems that the reorganization of chromatin structure during spermiogenesis might have affected the evolutionary constraints driving histone H2B evolution, leading to an increase in diversity.

The evolutionary differentiation of two histone H2A.Z variants in chordates (H2A.Z-1 and H2A.Z-2) is mediated by a stepwise mutation process that affects three amino acid residues.

BACKGROUND: The histone H2A family encompasses the greatest number of core histone variants of which the replacement variant H2A.Z is currently one of the most heavily studied. No clear mechanism for the functional variability that H2A.Z imparts to chromatin has yet been proposed. While most of the past studies have referred to H2A.Z generically as a single protein, in vertebrates it is a mixture of two protein forms H2A.Z-1 (previously H2A.Z) and H2A.Z-2 (previously H2A.F/Z or H2A.V) that differ by three amino acids. RESULTS: We have performed an extensive study on the long-term evolution of H2A.Z across metazoans with special emphasis on the possible selective mechanisms responsible for the differentiation between H2A.Z-1 and H2A.Z-2. Our results reveal a common origin of both forms early in chordate evolution. The evolutionary process responsible for the differentiation involves refined stepwise mutation change within the codons of the three differential residues. This eventually led to differences in the intensity of the selective constraints acting upon the different H2A.Z forms in vertebrates. CONCLUSION: The results presented in this work definitively reveal that the existence of H2A.Z-1 and H2A.Z-2 is not a whim of random genetic drift. Our analyses demonstrate that H2A.Z-2 is not only subject to a strong purifying selection but it is significantly more evolutionarily constrained than H2A.Z-1. Whether or not the evolutionary drift between H2A.Z-1 and H2A.Z-2 has resulted in a functional diversification of these proteins awaits further research. Nevertheless, the present work suggests that in the process of their differently constrained evolutionary pathways, these two forms may have acquired new or complementary functions.

Histone genes of the razor clam Solen marginatus unveil new aspects of linker histone evolution in protostomes.

The association of DNA with histones results in a nucleoprotein complex called chromatin that consists of repetitive nucleosomal subunits. Nucleosomes are joined together in the chromatin fiber by short stretches of linker DNA that interact with a wide diversity of linker H1 histones involved in chromatin compaction and dynamics. Although the long-term evolution of the H1 family has been the subject of different studies during the last 5 years, the lack of molecular data on replication-independent (RI) H1 variants from protostomes has been hampering attempts to complete the evolutionary picture of this histone family in eukaryotes, especially as it pertains to the functional specialization they impart to the chromatin structure in members of this bilaterian lineage. In an attempt to fill this gap, the present work characterizes the histone gene complement from the razor clam Solen marginatus. Molecular evolutionary analyses reveal that the H1 gene from this organism represents one of the few protostome RI H1 genes known to date, a notion which is further supported by its location within the monophyletic group encompassing the RI H1 variants in the overall phylogeny of eukaryotic H1 proteins. Although the detailed characterization of the nucleotide substitution patterns in RI H1 variants agrees with the model of birth-and-death evolution under strong purifying selection, maximum-likelihood approaches unveil the presence of adaptive selection during at least part of the evolutionary differentiation between protostomes and deuterostomes. The presence of increased levels of specialization in RI H1 proteins from deuterostomes as well as the significant differences observed in electrostatic properties between protostome and deuterostome RI H1s represent novel and important preliminary results for future studies of the functional differentiation of this histone H1 lineage across bilaterians.

Quickly evolving histones, nucleosome stability and chromatin folding: all about histone H2A.Bbd.

Histone H2A.Bbd (Barr body-deficient) is a novel histone variant which is largely excluded from the inactive X chromosome of mammals. Discovered only 6 years ago, H2A.Bbd displays very unusual structural and functional properties, for instance, it is relatively shorter and only 48% identical compared to H2A, lacking both the typical C-terminal tail of the H2A family and the very last sequence of the docking domain, making it the most specialized among all histone variants known to date. Indeed, molecular evolutionary analyses have shown that H2A.Bbd is a highly hypervariable and quickly evolving protein exclusive to mammalian lineages, in striking contrast to all other histones. Different studies have described a deposition pattern of H2A.Bbd in the chromatin that overlaps with regions of histone H4 acetylation suggesting its association with transcriptionally active euchromatic regions of the genome. In this regard, it is believed that this histone variant plays an important role in determining such regions by destabilizing the nucleosome and locally unfolding the chromatin fiber. This review provides a concise, comprehensive and timely summary of the work published on H2A.Bbd structure and function. Special emphasis is placed on its chromatin deposition patterns in relation to gene expression profiles and its evolutionary history, as well as on the dynamics of H2A.Bbd-containing nucleosomes.

Molecular and Biochemical Methods Useful for the Epigenetic Characterization of Chromatin-Associated Proteins in Bivalve Molluscs.

Bivalve molluscs constitute a ubiquitous taxonomic group playing key functions in virtually all ecosystems, and encompassing critical commercial relevance. Along with a sessile and filter-feeding lifestyle in most cases, these characteristics make bivalves model sentinel organisms routinely used for environmental monitoring studies in aquatic habitats. The study of epigenetic mechanisms linking environmental exposure and specific physiological responses (i.e., environmental epigenetics) stands out as a very innovative monitoring strategy, given the role of epigenetic modifications in acclimatization and adaptation. Furthermore, the heritable nature of many of those modifications constitutes a very promising avenue to explore the applicability of epigenetic conditioning and selection in management and restoration strategies. Chromatin provides a framework for the study of environmental epigenetic responses. Unfortunately, chromatin and epigenetic information are very limited in most non-traditional model organisms and even completely lacking in most environmentally and ecologically relevant organisms. The present work aims to provide a comprehensive and reproducible experimental workflow for the study of bivalve chromatin. First, a series of guidelines for the molecular isolation of genes encoding chromatin-associated proteins is provided, including information on primers suitable for conventional PCR, Rapid Amplification of cDNA Ends (RACE), genome walking and quantitative PCR (qPCR) experiments. This section is followed by the description of methods specifically developed for the analysis of histone and SNBP proteins in different bivalve tissues, including protein extraction, purification, separation and immunodetection. Lastly, information about available antibodies, their specificity and performance is also provided. The tools and protocols described here complement current epigenetic analyses (usually limited to DNA methylation) by incorporating the study of structural elements modulating chromatin dynamics.

Effects of Florida Red Tides on histone variant expression and DNA methylation in the Eastern oyster Crassostrea virginica.

Massive algal proliferations known as Harmful Algal Blooms (HABs) represent one of the most important threats to coastal areas. Among them, the so-called Florida Red Tides (FRTs, caused by blooms of the dinoflagellate Karenia brevis and associated brevetoxins) are particularly detrimental in the southeastern U.S., causing high mortality rates and annual losses in excess of $40 million. The ability of marine organisms to cope with environmental stressors (including those produced during HABs) is influenced by genetic and epigenetic mechanisms, the latter resulting in phenotypic changes caused by heritable modifications in gene expression, without involving changes in the genetic (DNA) sequence. Yet, studies examining cause-effect relationships between environmental stressors, specific epigenetic mechanisms and subsequent responses are still lacking. The present work contributes to increase this knowledge by investigating the effects of Florida Red Tides on two types of mechanisms participating in the epigenetic memory of Eastern oysters: histone variants and DNA methylation. For that purpose, a HAB simulation was conducted in laboratory conditions, exposing oysters to increasing concentrations of K. brevis. The obtained results revealed, for the first time, the existence of H2A.X, H2A.Z and macroH2A genes in this organism, encoding histone variants potentially involved in the maintenance of genome integrity during responses to the genotoxic effect of brevetoxins. Additionally, an increase in H2A.X phosphorylation (γH2A.X, a marker of DNA damage) and a decrease in global DNA methylation were observed as the HAB simulation progressed. Overall, the present work provides a basis to better understand how epigenetic mechanisms participate in responses to environmental stress in marine invertebrates, opening new avenues to incorporate environmental epigenetics approaches into management and conservation programs.

Basic surface features of nuclear FKBPs facilitate chromatin binding.

The nucleoplasmin family of histone chaperones is identified by a pentamer-forming domain and multiple acidic tracts that mediate histone binding and chaperone activity. Within this family, a novel domain organization was recently discovered that consists of an N-terminal nucleoplasmin-like (NPL) domain and a C-terminal FKBP peptidyl-proline isomerase domain. Saccharomyces cerevisiae Fpr4 is one such protein. Here we report that in addition to its known histone prolyl isomerase activities, the Fpr4 FKBP domain binds to nucleosomes and nucleosome arrays in vitro. This ability is mediated by a collection of basic patches that enable the enzyme to stably associate with linker DNA. The interaction of the Fpr4 FKBP with recombinant chromatin complexes condenses nucleosome arrays independently of its catalytic activity. Based on phylogenetic comparisons we propose that the chromatin binding ability of 'basic' FKBPs is shared amongst related orthologues present in fungi, plants, and insects. Thus, a subclass of FKBP prolyl isomerase enzymes is recruited to linker regions of chromatin.

The characterization of macroH2A beyond vertebrates supports an ancestral origin and conserved role for histone variants in chromatin.

Histone variants play a critical role in chromatin structure and epigenetic regulation. These "deviant" proteins have been historically considered as the evolutionary descendants of ancestral canonical histones, helping specialize the nucleosome structure during eukaryotic evolution. Such view is now challenged by 2 major observations: first, canonical histones present extremely unique features not shared with any other genes; second, histone variants are widespread across many eukaryotic groups. The present work further supports the ancestral nature of histone variants by providing the first in vivo characterization of a functional macroH2A histone (a variant long defined as a specific refinement of vertebrate chromatin) in a non-vertebrate organism (the mussel Mytilus) revealing its recruitment into heterochromatic fractions of actively proliferating tissues. Combined with in silico analyses of genomic data, these results provide evidence for the widespread presence of macroH2A in metazoan animals, as well as in the holozoan Capsaspora, supporting an evolutionary origin for this histone variant lineage before the radiation of Filozoans (including Filasterea, Choanoflagellata and Metazoa). Overall, the results presented in this work help configure a new evolutionary scenario in which histone variants, rather than modern "deviants" of canonical histones, would constitute ancient components of eukaryotic chromatin.

Unique yeast histone sequences influence octamer and nucleosome stability.

Yeast nucleosomes are known to be intrinsically less stable than those from higher eukaryotes. This difference presents significant challenges for the production of yeast nucleosome core particles (NCPs) and chromatin for in vitro analyses. Using recombinant yeast, human, and chimeric histone proteins, we demonstrate that three divergent amino acids in histone H3 (Q120 K121 K125 ) are responsible for the poor reconstitution of yeast histones into octamers. This QKK motif is only found in Fungi, and is located at the nucleosome dyad axis. Yeast-to-human changes at these positions render yeast histones amenable to well-established octamer reconstitution and salt dialysis methods for generating nucleosomal and longer chromatin templates. By contrast, the most divergent yeast core histones, H2A and H2B, affect the biophysical properties of NCP but not their stability. An evolutionary analysis of H3 sequences shows that a gradual divergence in H3 sequences occurred in Fungi to yield QKK in budding yeast. This likely facilitates the highly euchromatic nature of yeast genomes. Our results provide an explanation for the long recognized difference in yeast nucleosome stability and they offer a simple method to generate yeast chromatin templates for in vitro studies.

Environmental epigenetics: A promising venue for developing next-generation pollution biomonitoring tools in marine invertebrates.

Environmental epigenetics investigates the cause-effect relationships between specific environmental factors and the subsequent epigenetic modifications triggering adaptive responses in the cell. Given the dynamic and potentially reversible nature of the different types of epigenetic marks, environmental epigenetics constitutes a promising venue for developing fast and sensible biomonitoring programs. Indeed, several epigenetic biomarkers have been successfully developed and applied in traditional model organisms (e.g., human and mouse). Nevertheless, the lack of epigenetic knowledge in other ecologically and environmentally relevant organisms has hampered the application of these tools in a broader range of ecosystems, most notably in the marine environment. Fortunately, that scenario is now changing thanks to the growing availability of complete reference genome sequences along with the development of high-throughput DNA sequencing and bioinformatic methods. Altogether, these resources make the epigenetic study of marine organisms (and more specifically marine invertebrates) a reality. By building on this knowledge, the present work provides a timely perspective highlighting the extraordinary potential of environmental epigenetic analyses as a promising source of rapid and sensible tools for pollution biomonitoring, using marine invertebrates as sentinel organisms. This strategy represents an innovative, groundbreaking approach, improving the conservation and management of natural resources in the oceans.

Interaction of chromatin with a histone H1 containing swapped N- and C-terminal domains.

Although the details of the structural involvement of histone H1 in the organization of the nucleosome are quite well understood, the sequential events involved in the recognition of its binding site are not as well known. We have used a recombinant human histone H1 (H1.1) in which the N- and C-terminal domains (NTD/CTD) have been swapped and we have reconstituted it on to a 208-bp nucleosome. We have shown that the swapped version of the protein is still able to bind to nucleosomes through its structurally folded wing helix domain (WHD); however, analytical ultracentrifuge analysis demonstrates its ability to properly fold the chromatin fibre is impaired. Furthermore, FRAP analysis shows that the highly dynamic binding association of histone H1 with the chromatin fibre is altered, with a severely decreased half time of residence. All of this suggests that proper binding of histone H1 to chromatin is determined by the simultaneous and synergistic binding of its WHD-CTD to the nucleosome.

dBigH1, a second histone H1 in Drosophila, and the consequences for histone fold nomenclature.

Recently, Pérez-Montero and colleagues (Developmental cell, 26: 578-590, 2013) described the occurrence of a new histone H1 variant (dBigH1) in Drosophila. The presence of unusual acidic amino acid patches at the N-terminal end of dBigH1 is in contrast to the arginine patches that exist at the N- and C-terminal domains of other histone H1-related proteins found in the sperm of some organisms. This departure from the strictly lysine-rich composition of the somatic histone H1 raises a question about the true definition of its protein members. Their minimal essential requirements appear to be the presence of a lysine- and alanine-rich, intrinsically disordered C-terminal domain, with a highly helicogenic potential upon binding to the linker DNA regions of chromatin. In metazoans, specific targeting of these regions is further achieved by a linker histone fold domain (LHFD), distinctively different from the characteristic core histone fold domain (CHFD) of the nucleosome core histones.

Chromatin specialization in bivalve molluscs: a leap forward for the evaluation of Okadaic Acid genotoxicity in the marine environment.

Marine biotoxins synthesized by Harmful Algal Blooms (HABs) represent one of the most important sources of contamination in marine environments as well as a serious threat to fisheries and aquaculture-based industries in coastal areas. Among these biotoxins Okadaic Acid (OA) is of critical interest as it represents the most predominant Diarrhetic Shellfish Poisoning biotoxin in the European coasts. Furthermore, OA is a potent tumor promoter with aneugenic and clastogenic effects on the hereditary material, most notably DNA breaks and alterations in DNA repair mechanisms. Therefore, a great effort has been devoted to the biomonitoring of OA in the marine environment during the last two decades, mainly based on physicochemical and physiological parameters using mussels as sentinel organisms. However, the molecular genotoxic effects of this biotoxin make chromatin structure a good candidate for an alternative strategy for toxicity assessment with faster and more sensitive evaluation. To date, the development of chromatin-based studies to this purpose has been hampered by the complete lack of information on chromatin of invertebrate marine organisms, especially in bivalve molluscs. Our preliminary results have revealed the presence of histone variants involved in DNA repair and chromatin specialization in mussels and clams. In this work we use this information to put forward a proposal focused on the development of chromatin-based tests for OA genotoxicity in the marine environment. The implementation of such tests in natural populations has the potential to provide an important leap in the biomonitoring of this biotoxin. The outcome of such monitoring may have critical implications for the evaluation of DNA damage in these marine organisms. They will provide as well important tools for the optimization of their harvesting and for the elaboration of additional tests designed to evaluate the safety of their consumption and potential implications for consumer's health.

Histone H2A (H2A.X and H2A.Z) variants in molluscs: molecular characterization and potential implications for chromatin dynamics.

Histone variants are used by the cell to build specialized nucleosomes, replacing canonical histones and generating functionally specialized chromatin domains. Among many other processes, the specialization imparted by histone H2A (H2A.X and H2A.Z) variants to the nucleosome core particle constitutes the earliest response to DNA damage in the cell. Consequently, chromatin-based genotoxicity tests have been developed in those cases where enough information pertaining chromatin structure and dynamics is available (i.e., human and mouse). However, detailed chromatin knowledge is almost absent in most organisms, specially protostome animals. Molluscs (which represent sentinel organisms for the study of pollution) are not an exception to this lack of knowledge. In the present work we first identified the existence of functionally differentiated histone H2A.X and H2A.Z variants in the mussel Mytilus galloprovincialis (MgH2A.X and MgH2A.Z), a marine organism widely used in biomonitoring programs. Our results support the functional specialization of these variants based on: a) their active expression in different tissues, as revealed by the isolation of native MgH2A.X and MgH2A.Z proteins in gonad and hepatopancreas; b) the evolutionary conservation of different residues encompassing functional relevance; and c) their ability to confer specialization to nucleosomes, as revealed by nucleosome reconstitution experiments using recombinant MgH2A.X and MgH2A.Z histones. Given the seminal role of these variants in maintaining genomic integrity and regulating gene expression, their preliminary characterization opens up new potential applications for the future development of chromatin-based genotoxicity tests in pollution biomonitoring programs.

Early evolution of histone genes: prevalence of an 'orphon' H1 lineage in protostomes and birth-and-death process in the H2A family.

The study of histone evolution has experienced a rebirth, for two main reasons: the identification of new essential histone variants responsible for regulating chromatin dynamics and the subsequent contradictions posed by this variability as it pertains to their long-term evolution process. Although different evolutionary models (e.g., birth-and-death evolution, concerted evolution) may account for the observed divergence of histone genes, conclusive evidence is lacking (e.g., histone H1) or totally nonexistent (e.g., histone H2A). While most of the published work has focused on deuterostomes, very little is known about the diversification and functional differentiation mechanisms followed by histone protein subtypes in protostomes, for which histone variants have only been recently described. In this study, we identify linker and core histone genes in three clam species. Our results demonstrate the prevalence of an 'orphon' H1 lineage in molluscs, a group in which the protostome H1 and sperm nuclear basic proteins are on the verge of diversification. They share an early monophyletic origin with vertebrate-specific variants prior to the differentiation between protostomes and deuterostomes. Given the intringuing evolutionary features of the histone H1 family, we have evaluated the relative importance of gene conversion, point mutation, and selection in maintaining the diversity found among H2A subtypes in eukaryotes. We show evidence for the first time that the long-term evolution of this family is not subject to concerted evolution but, rather, to a gradual evolution following a birth-and-death model under a strong purifying selection at the protein level.