Early Evolution of Histone Genes: Prevalence of an ‘Orphon’ H1 Lineage in Protostomes and Birth-and-Death Process in the H2A Family

Evolution of methyltransferase like (METTL) proteins in Metazoa: A complex gene family involved in epitranscriptomic regulation and other epigenetic processes

The methyltransferase-like (METTL) proteins constitute a family of seven-beta-strand methyltransferases with S-adenosyl methionine-binding domains that modify DNA, RNA, and proteins. Methylation by METTL proteins contributes to the epigenetic, and in the case of RNA modifications, epitranscriptomic regulation of a variety of biological processes. Despite their functional importance, most investigations of the substrates and functions of METTLs within metazoans have been restricted to model vertebrate taxa. In the present work, we explore the evolutionary mechanisms driving the diversification and functional differentiation of 33 individual METTL proteins across Metazoa. Our results show that METTLs are nearly ubiquitous across the animal kingdom, with most having arisen early in metazoan evolution (i.e., occur in basal metazoan phyla). Individual METTL lineages each originated from single independent ancestors, constituting monophyletic clades, which suggests that each METTL was subject to strong selective constraints driving its structural and/or functional specialization. Interestingly, a similar process did not extend to the differentiation of nucleoside-modifying and protein-modifying METTLs (i.e., each METTL type did not form a unique monophyletic clade). The members of these two types of METTLs also exhibited differences in their rates of evolution. Overall, we provide evidence that the long-term evolution of METTL family members was driven by strong purifying selection, which in combination with adaptive selection episodes, led to the functional specialization of individual METTL lineages. This work contributes useful information regarding the evolution of a gene family that fulfills a variety of epigenetic functions, and can have profound influences on molecular processes and phenotypic traits.

A preliminary integrated genetic map distinguishes every chromosome pair and locates essential genes related to abiotic adaptation of Crassostrea angulata/gigas

BACKGROUND

The re-sequencing of C. angulata has revealed many polymorphisms in candidate genes related to adaptation to abiotic stress that are not present in C. gigas; these genes, therefore, are probably related to the ability of this oyster to retain high concentrations of toxic heavy metals. There is, in addition, an unresolved controversy as to whether or not C. angulata and C. gigas are the same species or subspecies. Both oysters have 20 metacentric chromosomes of similar size that are morphologically indistinguishable. From a genomic perspective, as a result of the great variation and selection for heterozygotes in C. gigas, the assembly of its draft genome was difficult: it is fragmented in more than seven thousand scaffolds.

RESULTS

In this work sixty BAC sequences of C. gigas downloaded from NCBI were assembled in BAC-contigs and assigned to BACs that were used as probes for mFISH in C. angulata and C. gigas. In addition, probes of H3, H4 histone, 18S and 5S rDNA genes were also used. Hence we obtained markers identifying 8 out the 10 chromosomes constituting the karyotype. Chromosomes 1 and 9 can be distinguished morphologically. The bioinformatic analysis carried out with the BAC-contigs annotated 88 genes. As a result, genes associated with abiotic adaptation, such as metallothioneins, have been positioned in the genome. The gene ontology analysis has also shown many molecular functions related to metal ion binding, a phenomenon associated with detoxification processes that are characteristic in oysters. Hence the provisional integrated map obtained in this study is a useful complementary tool for the study of oyster genomes.

CONCLUSIONS

In this study 8 out of 10 chromosome pairs of Crassostrea angulata/gigas were identified using BAC clones as probes. As a result all chromosomes can now be distinguished. Moreover, FISH showed that H3 and H4 co-localized in two pairs of chromosomes different that those previously escribed. 88 genes were annotated in the BAC-contigs most of them related with Molecular Functions of protein binding, related to the resistance of the species to abiotic stress. An integrated genetic map anchored to the genome has been obtained in which the BAC-contigs structure were not concordant with the gene structure of the C. gigas scaffolds displayed in the Genomicus database.

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.

Analysis of the histone cluster in Senegalese sole (Solea senegalensis): evidence for a divergent evolution of two canonical histone clusters

The Senegalese sole (Solea senegalensis) is commercially very important and a priority species for aquaculture product diversification. The main histone cluster was identified within two BAC clones. However, two replacement histones (H1.0 and H3.3) were found in another BAC clone. Different types of canonical histones H2A and H2B were found within the same species for the first time. Phylogenetic analysis demonstrated that the different types of H1, H2A, and H2B histones were all more similar to each other than to canonical histones from other species. The canonical histone H3 of S. senegalensis differs from subtypes H3.1 and H3.2 in humans at the site of residue 96, where a serine is found instead of an alanine. This same polymorphism has been found only in Danio rerio. The karyotype of S. senegalensis comprises 21 pairs of chromosomes, distributed in 3 metacentric pairs, 2 submetacentric pairs, 4 subtelocentric pairs, and 12 acrocentric pairs. The two BAC clones that contain the clusters of canonical histones were both mapped on the largest metacentric pair, and mFISH analysis confirmed the co-location with the dmrt1 gene in that pair. Three chromosome markers have been identified which, in addition to those previously described, account for 18 chromosome pairs in S. senegalensis.

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.

Novel O-GlcNAcylation on Ser(40) of canonical H2A isoforms specific to viviparity

We report here newly discovered O-linked-N-acetylglucosamine (O-GlcNAc) modification of histone H2A at Ser⁴⁰ (H2AS40Gc). The mouse genome contains 18 H2A isoforms, of which 13 have Ser⁴⁰ and the other five have Ala⁴⁰. The combination of production of monoclonal antibody and mass spectrometric analyses with reverse-phase (RP)-high performance liquid chromatography (HPLC) fractionation indicated that the O-GlcNAcylation is specific to the Ser⁴⁰ isoforms. The H2AS40Gc site is in the L1 loop structure where two H2A molecules interact in the nucleosome. Targets of H2AS40Gc are distributed genome-wide and are dramatically changed during the process of differentiation in mouse trophoblast stem cells. In addition to the mouse, H2AS40Gc was also detected in humans, macaques and cows, whereas non-mammalian species possessing only the Ala⁴⁰ isoforms, such as silkworms, zebrafish and Xenopus showed no signal. Genome database surveys revealed that Ser⁴⁰ isoforms of H2A emerged in Marsupialia and persisted thereafter in mammals. We propose that the emergence of H2A Ser⁴⁰ and its O-GlcNAcylation linked a genetic event to genome-wide epigenetic events that correlate with the evolution of placental animals.

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.

The Evolutionary Dynamics of Ribosomal Genes, Histone H3, and Transposable Rex Elements in the Genome of Atlantic Snappers

Lutjanidae is a family of primarily marine and carnivorous fishes distributed in the Atlantic, Indian, and Pacific oceans, with enormous economic and ecological importance. In order to better clarify the conservative chromosomal evolution of Lutjanidae, we analyzed the evolutionary dynamics of 5 repetitive DNA classes in 5 Lutjanus and in 1 Ocyurus species from the Western Atlantic. The ribosomal 18S sites were generally located in a single chromosome pair, except for L. jocu and L. alexandrei where they are found in 2 pairs. In turn, the 5S rDNA sites are unique, terminal and nonsyntenic with the 18S rDNA sites. In 3 species analyzed, H3 hisDNA genes were found in 1 chromosomal pair. However, while L. jocu presented 2 H3 sites, O. chrysurus showed a noteworthy dispersion of this gene in almost all chromosomes of the karyotype. Retrotransposons Rex1 and Rex3 do not exhibit any association with the explosive distribution of H3 sequences in O. chrysurus. The low compartmentalization of Rex elements, in addition to the general nondynamic distribution of ribosomal and H3 genes, corroborate the karyotype conservatism in Lutjanidae species, also at the microstructural level. However, some "disturbing evolutionary waves" can break down this conservative scenario, as evidenced by the massive random dispersion of H3 hisDNA in the genome of O. chrysurus. The implication of the genomic expansion of H3 histone genes and their functionality remain unknown, although suggesting that they have higher evolutionary dynamics than previously thought.

Divergent evolutionary behavior of H3 histone gene and rDNA clusters in venerid clams

BACKGROUND

Histone H3 gene clusters have been described as highly conserved chromosomal markers in invertebrates. Surprisingly, in bivalves remarkable interspecific differences were found among the eight mussels and between the two clams in which histone H3 gene clusters have already been located. Although the family Veneridae comprises 10 % of the species of marine bivalves, their chromosomes are poorly studied. The clams belonging to this family present 2n = 38 chromosomes and similar karyotypes showing chromosome pairs gradually decreasing in length. In order to assess the evolutionary behavior of histone and rRNA multigene families in bivalves, we mapped histone H3 and ribosomal RNA probes to chromosomes of ten species of venerid clams.

RESULTS

In contrast with the reported conservation of histone H3 gene clusters and their intercalary location in invertebrates, these loci varied in number and were mostly subterminal in venerid clams. On the other hand, while a single 45S rDNA cluster, highly variable in location, was found in these organisms, 5S rDNA clusters showed interspecific differences in both number and location. The distribution patterns of these sequences were species-specific and mapped to different chromosomal positions in all clams but Ruditapes decussatus, in which one of the minor rDNA clusters and the major rDNA cluster co-located.

CONCLUSION

The diversity in the distribution patterns of histone H3 gene, 5S rDNA and 28S rDNA clusters found in venerid clams, together with their different evolutionary behaviors in other invertebrate taxa, strongly suggest that the control of the spreading of these multigene families in a group of organisms relies upon a combination of evolutionary forces that operate differently depending not only on the specific multigene family but also on the particular taxa. Our data also showed that H3 histone gene and rDNA clusters are useful landmarks to integrate nex-generation sequencing (NGS) and evolutionary genomic data in non-model species.

Bivalve omics: state of the art and potential applications for the biomonitoring of harmful marine compounds

The extraordinary progress experienced by sequencing technologies and bioinformatics has made the development of omic studies virtually ubiquitous in all fields of life sciences nowadays. However, scientific attention has been quite unevenly distributed throughout the different branches of the tree of life, leaving molluscs, one of the most diverse animal groups, relatively unexplored and without representation within the narrow collection of well established model organisms. Within this Phylum, bivalve molluscs play a fundamental role in the functioning of the marine ecosystem, constitute very valuable commercial resources in aquaculture, and have been widely used as sentinel organisms in the biomonitoring of marine pollution. Yet, it has only been very recently that this complex group of organisms became a preferential subject for omic studies, posing new challenges for their integrative characterization. The present contribution aims to give a detailed insight into the state of the art of the omic studies and functional information analysis of bivalve molluscs, providing a timely perspective on the available data resources and on the current and prospective applications for the biomonitoring of harmful marine compounds.

A unified phylogeny-based nomenclature for histone variants

Histone variants are non-allelic protein isoforms that play key roles in diversifying chromatin structure. The known number of such variants has greatly increased in recent years, but the lack of naming conventions for them has led to a variety of naming styles, multiple synonyms and misleading homographs that obscure variant relationships and complicate database searches. We propose here a unified nomenclature for variants of all five classes of histones that uses consistent but flexible naming conventions to produce names that are informative and readily searchable. The nomenclature builds on historical usage and incorporates phylogenetic relationships, which are strong predictors of structure and function. A key feature is the consistent use of punctuation to represent phylogenetic divergence, making explicit the relationships among variant subtypes that have previously been implicit or unclear. We recommend that by default new histone variants be named with organism-specific paralog-number suffixes that lack phylogenetic implication, while letter suffixes be reserved for structurally distinct clades of variants. For clarity and searchability, we encourage the use of descriptors that are separate from the phylogeny-based variant name to indicate developmental and other properties of variants that may be independent of structure.

Symbiodinium transcriptomes: genome insights into the dinoflagellate symbionts of reef-building corals

Dinoflagellates are unicellular algae that are ubiquitously abundant in aquatic environments. Species of the genus Symbiodinium form symbiotic relationships with reef-building corals and other marine invertebrates. Despite their ecologic importance, little is known about the genetics of dinoflagellates in general and Symbiodinium in particular. Here, we used 454 sequencing to generate transcriptome data from two Symbiodinium species from different clades (clade A and clade B). With more than 56,000 assembled sequences per species, these data represent the largest transcriptomic resource for dinoflagellates to date. Our results corroborate previous observations that dinoflagellates possess the complete nucleosome machinery. We found a complete set of core histones as well as several H3 variants and H2A.Z in one species. Furthermore, transcriptome analysis points toward a low number of transcription factors in Symbiodinium spp. that also differ in the distribution of DNA-binding domains relative to other eukaryotes. In particular the cold shock domain was predominant among transcription factors. Additionally, we found a high number of antioxidative genes in comparison to non-symbiotic but evolutionary related organisms. These findings might be of relevance in the context of the role that Symbiodinium spp. play as coral symbionts.

Our data represent the most comprehensive dinoflagellate EST data set to date. This study provides a comprehensive resource to further analyze the genetic makeup, metabolic capacities, and gene repertoire of Symbiodinium and dinoflagellates. Overall, our findings indicate that Symbiodinium possesses some unique characteristics, in particular the transcriptional regulation in Symbiodinium may differ from the currently known mechanisms of eukaryotic gene regulation.

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.

Evolution in the block: common elements of 5S rDNA organization and evolutionary patterns in distant fish genera

The 5S rDNA is organized in the genome as tandemly repeated copies of a structural unit composed of a coding sequence plus a nontranscribed spacer (NTS). The coding region is highly conserved in the evolution, whereas the NTS vary in both length and sequence. It has been proposed that 5S rRNA genes are members of a gene family that have arisen through concerted evolution. In this study, we describe the molecular organization and evolution of the 5S rDNA in the genera Lepidorhombus and Scophthalmus (Scophthalmidae) and compared it with already known 5S rDNA of the very different genera Merluccius (Merluccidae) and Salmo (Salmoninae), to identify common structural elements or patterns for understanding 5S rDNA evolution in fish. High intra- and interspecific diversity within the 5S rDNA family in all the genera can be explained by a combination of duplications, deletions, and transposition events. Sequence blocks with high similarity in all the 5S rDNA members across species were identified for the four studied genera, with evidences of intense gene conversion within noncoding regions. We propose a model to explain the evolution of the 5S rDNA, in which the evolutionary units are blocks of nucleotides rather than the entire sequences or single nucleotides. This model implies a "two-speed" evolution: slow within blocks (homogenized by recombination) and fast within the gene family (diversified by duplications and deletions).

Chromosomal mapping of rRNA genes, core histone genes and telomeric sequences in Brachidontes puniceus and Brachidontes rodriguezi (Bivalvia, Mytilidae)

Background

Chromosome rearrangements are an important part of the speciation process in many taxa. The study of chromosome evolution in bivalves is hampered by the absence of clear chromosomal banding patterns and the similarity in both chromosome size and morphology. For this reason, obtaining good chromosome markers is essential for reliable karyotypic comparisons. To begin this task, the chromosomes of the mussels Brachidontes puniceus and B. rodriguezi were studied by means of fluorochrome staining and fluorescent in situ hybridization (FISH).

Results

Brachidontes puniceus and B. rodriguezi both have 2n = 32 chromosomes but differing karyotype composition. Vertebrate-type telomeric sequences appear at both ends of every single chromosome. B. puniceus presents a single terminal major rRNA gene cluster on a chromosome pair while B. rodriguezi shows two. Both mussels present two 5S rDNA and two core histone gene clusters intercalary located on the long arms of two chromosome pairs. Double and triple-FISH experiments demonstrated that one of the 5S rDNA and one of the major rDNA clusters appear on the same chromosome pair in B. rodriguezi but not in B. puniceus. On the other hand, the second 5S rDNA cluster is located in one of the chromosome pairs also bearing one of the core histone gene clusters in the two mussel species.

Conclusion

Knowledge of the chromosomal distribution of these sequences in the two species of Brachidontes is a first step in the understanding of the role of chromosome changes on bivalve evolution.

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.

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.

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.