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
© Eirín-López et al; licensee BioMed Central Ltd. 2009
Received: 30 October 2008
Accepted: 04 February 2009
Published: 04 February 2009
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.
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.
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.
In eurkaryotic organisms, DNA is found associated with histone proteins constituting a nucleoprotein complex called chromatin. Approximately 146 base pairs of DNA wrap around a core histone octamer to form a nucleosome which is the basic subunit of chromatin. This nucleoprotein complex allows for the high extent of compaction of genomic DNA within the cell nucleus and provides the support on which most DNA metabolic processes take place . There are five histone families which can be classified into core histones (H2A, H2B, H3, and H4) and linker histones (H1) according to structural and functional features. The histone H1 and H2A families show the most diversity of isoforms that have dedicated functions in many cellular processes including organization of chromatin structure in somatic and germinal cells, gene transcription, DNA replication, and DNA repair among others [2–9].
The histone H2A family contains the greatest number of variants among the core histones, some of which are essential for the maintenance of genome integrity and viability such as H2A.Z and H2A.X [2, 10, 11]. At present, H2A.Z is one of the most heavily studied histone variants and it has been ascribed multiple functions that may differ among species. In yeast, H2A.Z (Htz1) is present at active and inactive gene promoters in euchromatin, it is depleted at the silenced subtelomeric heterochromatin and it is enriched at the boundaries between euchromatin and heterochromatin . Although studies concerning the function of H2A.Z in mammalian cells have always yielded results that seem difficult to reconcile, a growing body of evidence suggests that H2A.Z is present at gene promoters and that in an acetylated form, its presence correlates with gene expression [13, 14]. However, it appears there are also at least two populations of H2A.Z present in heterochromatin. Greaves and colleagues show that H2A.Z is a feature of pericentric heterochromatin and contributes to the structure of the centromere . Another fraction of H2A.Z stains the length of the inactive X chromosome though intriguingly this fraction can be distinguished by monoubiquitination at K120, K121 or K125 .
Analysis of the H2A.Z-containing nucleosome has also yielded conflicting results [3, 17, 18]. The crystal structure of this nucleosome initially suggested a subtle destabilization between the H2A.Z-H2B dimer and the H3-H4 tetramer . However, FRET and analytical ultracentrifuge analysis using native H2A.Z have indicated that the H2A.Z nucleosome is in fact slightly more stable than the canonical H2A nucleosome [20, 21]. Also, when H2A.Z is present in nucleosome arrays it facilitates the formation of the 30 nm chromatin fiber . The recent study by Sarcinella and colleagues showed that it is a monoubiquitinated form of H2A.Z that is present on the inactive X chromosome . This study is a clear demonstration that a post-translational modification has the potential to define a subpopulation of H2A.Z. Indeed, a similar situation can be seen with H2A.Z N-terminal acetylation and active gene transcription . The difference in PTMs could reflect the different functional constraints of the H2A.Z-containing mono-nucleosomes in the in vivo setting. Conversely, the tendency to fold the chromatin fiber in arrays consisting of contiguous H2A.Z-containing nucleosomes may account for the presence of this variant in physiologically relevant situations such as that found in association with PcG proteins in the polycomb genes  or at the flanking sites of the insulator protein CTCF .
Mass spectrometry analysis showed that purified H2A.Z consists of an almost equimolar amount of two similar yet distinct proteins that differ by three amino acids  which are labeled here as H2A.Z-1 and H2A.Z-2. In the present work we have explored the long-term evolutionary pathway of H2A.Z-1 and H2A.Z-2 differentiation across metazoans and have analyzed the possible selective mechanisms and the constraints responsible for the potential functionalization of both variants. Our results have important implications for histone evolution and function as they show for the first time that H2A.Z-1 and H2A.Z-2 represent two very closely related variants that share a common evolutionary origin early in chordate evolution. Furthermore, our results show that the evolutionary constraints leading to the differentiation of both variants are primarily acting at the nucleotide level. This involves a refined process of stepwise mutation change within the codons of their three characteristic amino acid residues. Finally, we show that H2A.Z-2 is more tightly controlled (constrained) by selection than H2A.Z-1.
Results and Discussion
The phylogenetic context of H2A histone variants
Therefore, a phylogeny based on the complete nucleotide coding regions was reconstructed in order to further investigate the mechanisms of differentiation between H2A.Z-1 and H2A.Z-2. In contrast to the protein phylogeny, this defines a clear pattern between H2A.Z-1 and H2A.Z-2 in vertebrates (Fig. 2), with a monophyletic origin which is already independent from that of H2A.Z-e genes in early chordates, protostomes and fungi.
Two exceptions to this clustering pattern were detected. In the first instance, histones H2A.Z-1 and H2A.Z-2 from zebrafish are included in the same group. This could be due to a recent gene duplication within a short period of time (not enough time elapsed to allow for the accumulation of nucleotide substitutions). Under such circumstances, the presence of gene conversion would be unlikely because H2A.Z-1 and H2A.Z-2 genes are located on different chromosomes in zebrafish. Another exception was observed for the H2A.Z-1 genes from X. laevis which fall into an independent group outside both the H2A.Z-1 and H2A.Z-2 clusters. Given that the H2A.Z-1 gene from X. tropicalis and those of X. laevis and X. tropicalis H2A.Z-2 fall within their corresponding groups in the tree, the independent position of X. laevis H2A.Z-1 could be due to the extensive nucleotide variation exhibited among H2A genes.
The phylogeny from Fig. 2 represents a very useful tool in defining the identity of vertebrate H2A variants as either H2A.Z-1 or H2A.Z-2 on the basis of their nucleotide sequences. Indeed, the current topology allowed us to define 24 sequences previously annotated as H2A either as H2A.Z-1 or H2A.Z-2 following their position on the tree. Furthermore, 20 of those sequences were previously annotated as unknown in the databases, and thus the present analysis helped to reveal their true identity. In two instances (chicken and rhesus monkey) these sequences were wrongly defined as H2A.Z-2 in the databases, and our analyses unveiled their H2A.Z-1 identity (see Additional file 2). Such a clear differentiation between both histone H2A.Z forms using nucleotide phylogeny suggests not only a common phylogenetic origin early in metazoan evolution, but also the presence of a process of functional differentiation similar to that described for other histone multigene families [28–31]. However, the differentiation among other histone family members commonly encompasses extensive synonymous divergence under strong purifying selection at the protein level [28, 29, 32–34]. In contrast, the differentiation between H2A.Z-1 and H2A.Z-2 seems to primarily involve variation at the nucleotide level which in vertebrates resulted in a subtle protein differentiation encompassing three different amino acid residues. Such structural amino acid refinement represents a new mechanism among those previously known in histone diversification and is further investigated below.
Protein and nucleotide variation in H2A.Z-1 and H2A.Z-2 variants
The amino acid sequence alignments shown in Fig. 3A and Additional file 4 highlight the high extent of protein similarity between H2A.Z-1 and H2A.Z-2. This is especially evident when comparing the sequences in the same species, as in the case of human H2A.Z-1 and H2A.Z-2 (Fig. 3A). The alignment reveals the presence in vertebrates of only three residue differences (triresidue) defined as: Thr → Ala (pos. 15, N-terminal region), Ser → Thr (pos. 39, central globular region), and Val → Ala (pos.128, C-terminal region). The amino acid residues involved in the acidic patch of H2A  that are presumably involved in interchromosomal contacts within the chromatin fibre  are maintained in both H2A.Z-1 and H2A.Z-2 in contrast to what has been observed in other H2A variants such as H2ABbd [26, 37]. Although the second amino acid of the triresidue falls within the a-helix 1 domain of the histone fold , a computer analysis using molecular dynamic simulation (see materials and methods) indicated that such variation has no structural implications for secondary or tertiary structure of these histone variants (Additional file 4).
Given the highly charged nature of histones and their involvement in protein-DNA interactions that modulate chromatin dynamics, the electrostatic interaction properties of the different H2A variants were analyzed. Electrostatic potentials and the corresponding similarity indices were calculated for all H2A.Z-1 and H2A.Z-2 proteins listed in Additional file 2, allowing us to calculate the electrostatic distances between proteins (Additional file 5). As expected from Additional file 4, this analysis showed that the differentiation process between H2A.Z-1 and H2A.Z-2 does not involve significant differences in electrostatic potentials within different taxonomic groups
Nucleotide and protein variation in the H2A variants analyzed.
Nucleotide and protein variation between the H2A variants analyzed and between different taxonomic groups.
The nature of the selective constraints operating on H2A.Z-1 and H2A.Z-2 histones
At the protein level, the evolutionary mechanisms leading to the differential identity of H2A.Z-1 and H2A.Z-2 seem to be operating beyond a process of purifying selection. Therefore, we decided to investigate the levels of codon usage bias in both genes in order to analyze the effect of selection at the nucleotide level. Differences in codon bias are common between different histone multigene families, more or less independently of the particular organisms studied [28, 40, 41]. Indeed, our results reveal significant differences between the H2A variants analyzed, with H2A.Z-2 (47.122 ± 2.456) exhibiting a significantly higher bias than H2A.Z-1 (55.106 ± 1.098; t-test = -6.799, P = 0.000). Such results have important implications for the evolutionary constraints affecting both histone types at the nucleotide level. Interestingly, no such significant difference in codon bias has ever been previously reported between histone variants belonging to the same histone multigene family, due to the presence of a strong purifying selection acting at the protein level as a major evolutionary force that leads to an extensive and homogeneous silent variation at the nucleotide level [28, 29, 32–34]. Also, the codon bias analyses are in agreement with the previous results at the amino acid level indicating also that H2A.Z-1 is less constrained than H2A.Z-2 at the nucleotide level. Furthermore, histone H2A.Z-e from early chordates displays an intermediate overall codon bias magnitude (49.666 ± 2.546) which is not significantly different from either H2A.Z-1 or from H2A.Z-2 (P > 0.05 in Duncan multiple range-test). These results suggest that the differences in the intensity of the selective constraints operating on H2A.Z-1 and H2A.Z-2 probably arose during the differentiation of both variants from a common ancestor, represented here by H2A.Z-e from early chordates.
Correlations between GC content and the frequency of GC-rich and GC-poor amino acids
Spearman rank correlation coefficient, rS
Genomic GC vs. GAPW (GC-rich)
Genomic GC vs. FYMINK (GC-poor)
Genomic GC vs. Alanine
Genomic GC vs. Lysine
Genomic GC vs. GAPW (GC-rich)
Genomic GC vs. FYMINK (GC-poor)
Genomic GC vs. Alanine
Genomic GC vs. Lysine
An additional method to gauge the significance of mutation bias and selection at the nucleotide level involves the comparison of changes at first codon positions (all nonsynonymous) with changes at fourfold positions (all synonymous) in the most frequent residues in H2A.Z-1 and H2A.Z-2. Under the neutral model the nucleotide frequencies at both positions should not be significantly different . Codons for glycine and alanine (GC-rich) contain G at first codon positions, whereas codons for lysine and isoleucine (GC-poor) have A at first codon positions. Analysis of the mean G + A content at first codon positions in H2A.Z-1 (74.94 ± 0.67) and H2A.Z-2 (77.49 ± 0.84) showed that their values were significantly larger than the mean G + A content at fourfold degenerate positions in H2A.Z-1 (41.76 ± 4.73) and H2A.Z-2 (44.88 ± 2.86), (H2A.Z-1, t-test = 44.627, P = 0.000; H2A.Z-2, t-test = 37.901, P = 0.000). The values are significantly different in all species analyzed (see Additional files 6 and 7). While the neutral model of molecular evolution predicts that amino acid and nucleotide compositions are driven by the underlying GC content as a result of mutation bias, our results strongly suggest that selection has acted to maintain high levels of glycine, alanine, lysine and isoleucine in H2A.Z-1 and H2A.Z-2 variants, biasing their nucleotide composition. Few studies have shown that natural selection is more important than mutation bias in determining amino acid composition of proteins [28, 42, 47–49]. In this regard, our observations with H2A.Z stand in contrast to the neutral model.
The progressive differentiation of H2A variants is mediated by stepwise mutations
The analyses presented in this work indicate that selective constraints governing H2A.Z-1 and H2A.Z-2 evolution go far beyond the protein level, as shown by the significant differences detected in codon usage bias and the presence of selection for highly biased amino acid composition that influences the nucleotide composition. However, the specific mechanisms responsible for the subtle differentiation between both forms and the functional meaning of this process remain obscure. In order to define possible functional selective targets in these proteins we decided to look at the codon usage of the amino acids of the triresidue that defines the identity of vertebrate H2A.Z-1 (Thr/Ser/Val) and H2A.Z-2 (Ala/Thr/Ala). The codons involved, which are indicated near each of the H2A sequences analyzed in the phylogeny shown in Fig. 2, show a high degree of conservation of codon in each histone form and within each taxonomic group.
Maximum composite likelihood estimation of the probability of substitution in H2A.Z-1 and H2A.Z-2.
The function of histone H2A.Z in gene activation/silencing is still an important topic in chromatin research, as no clear mechanism for its structural and functional variability has yet been proposed. In this regard, the presence of two different H2A.Z forms is especially interesting . Although very little is known about H2A.Z-2, the results presented in this work definitively reveal that its existence is not a whim of random genetic drift. The functional significance of H2A.Z-2 is still obscure, however our group has been able to demonstrate the coexistence of both H2A.Z-1 and H2A.Z-2 in chicken and human tissues, and that significant differences in their mRNA expression levels exist and in this regard, it is very likely that the key to the existence of these two functional H2A.Z forms resides within their promoter regions (manuscript in preparation). Our analyses demonstrate not only that H2A.Z-2 is subject to a strong purifying selection (as most histones are) but that in fact it is significantly more evolutionarily constrained than H2A.Z-1.
Nevertheless, it appears that this selection does not proceed in conventional ways. While phylogenetic and evolutionary analyses reveal a typical process of birth-and-death evolution with strong purifying selection leading to the differentiation of H2A family members , an almost identical primary structure has been conserved between H2A.Z-1 and H2A.Z-2 except for three amino acid differences. This is surprising considering that the two forms occupy different chromosomal locations (as revealed by the in silico analyses performed in the present work) and that they have resulted from a progressive differentiation across vertebrates starting from a common ancestor early in chordate evolution. The main evolutionary constraints directing the limited amino acid variation between H2A.Z-1 and H2A.Z-2 are primarily acting at the nucleotide level. This defines a process of stepwise mutation change in the codons constituting the triresidue which mirrors H2A.Z-1 and H2A.Z-2 evolution.
According to Clapier et al., the amino acid sequence changes observed in the protein variants throughout the highly constrained evolution of histones, are of little structural but decisive functional consequences . Indeed, is worthwhile to mention that knocking out H2A.Z-1 in mice results in lethality and therefore (at least during early development), H2A.Z-2 cannot replace H2A.Z-1 (either in terms of abundance or function). In the instance of H2A.Z-1 and H2A.Z-2, the difference in amino acid sequence variability is minimal and it affects only three residues. Yet, we have observed a much closer proximity of the constraints imposed at the nucleotide level between the sequence of the genes encoding H2A.Z-2 and the histone H2A.Z ancestor (H2A.Z-e) when compared to H2AZ-1. This suggests that in the transition from chordates to vertebrates, the H2A.Z-1 has arisen to acquire a novel, or most likely complementary functions.
Mining of H2A nucleotide data
A total of 109 nonredundant H2A nucleotide coding sequences available from eukaryotes was collected from the histone database  and GenBank through BLAST searches, including 64 canonical H2A sequences, 30 H2A.X sequences, 8 macro H2A sequences, and 7 H2A.Bbd sequences (see Additional file 1). Additional data mining performed on complete and draft genome databases resulted in the identification of 69 nucleotide sequences encoding H2A.Z-1 and H2A.Z-2 histone variants from eukaryotes (see Additional file 2). In vertebrates, these variants were identified based on the differences shown in three residues (triresidue) characteristic either of H2A.Z-1 (21 sequences identified) or H2A.Z-2 (26 sequences identified), correcting sequence nomenclature when necessary (see Additional file 2). Given that sequences from early chordates contain mixed characteristics of both H2A.Z-1 and H2A.Z-2, these were called H2A.Z-e (4 sequences).
Variation in H2A.Z-1 and H2A.Z-2 histone variants
Nucleotide coding sequences were aligned on the basis of their translated amino acid sequences using the BioEdit and CLUSTAL_X programs with the default parameters [52, 53]. A bar chart representation was used in order show the frequency of each residue at every position of the alignment of vertebrate H2A.Z-1 and H2A.Z-2 forms using the LogoBar program . The 3D structures of H2A.Z-1 and H2A.Z-2 proteins from vertebrates as well as H2A.Z-e from early chordates were modeled using the coordinates determined for the crystal structure of a nucleosome particle containing the variant histone H2A.Z-1 from human (PDB accession code 1F66) as a reference . Evaluation of model qualities in homology modeling was performed by two approaches: 1) GROMOS empirical force energy to estimate the local quality of the predicted structure, with the y-axis representing the energy for each amino acid of the protein (negative and positive energy values represent favorable and unfavorable energy environments, respectively, for a given amino acid); 2) Verify3D to analyze the compatibility of an atomic model with its own amino acid sequence, in which the y-axis represents the average profile score for each residue in a 21-residue sliding window with scores ranging from -1 (bad score) to +1 (good score). All modeling and evaluation analyses were performed using the SWISS-MODEL workspace  and structures were rendered using the MacPyMOL program .
The comparisons between H2A protein variants with respect to their electrostatic properties were conducted in the webPIPSA pipeline , starting with a set of 3D structures modeled for all proteins listed in Additional file 2 and using human H2A.Z-1 as a reference. Electrostatic potentials were calculated using the University of Houston Brownian Dynamics (UHBD) program , and the absolute distances calculated from the similarity indices for the electrostatic potentials were represented in a colorized matrix and an epogram (tree representation of the relationships among potentials). The representation of the electrostatic potentials in the modeled 3D structures was implemented with the VMD program .
The extent of nucleotide and amino acid divergence between sequences was estimated using uncorrected differences (p-distance) as this distance is known to give better results than more complicated methods when the number of sequences is large and the number of positions used is relatively small, because of its smaller variance . The numbers of synonymous and nonsynonymous nucleotide differences per site were computed by means of the modified Nei-Gojobori method . Distances were estimated using the complete-deletion option in all cases and standard errors were calculated by the bootstrap method (1000 replicates).
All molecular and evolutionary analyses in this work were conducted using the program MEGA ver. 4.0 , as well as the calculation of amino acid and nucleotide frequencies, the relative synonymous codon usage (RSCU) and the maximum composite likelihood estimation of the nucleotide substitution patterns. The codon usage bias in H2A variants was referred to as the effective number of codons (ENC), which ranges from 61 (no bias) to 20 (maximum bias) and does not need any previous information on codon usage preferences in the genomes analyzed . The analysis of the nucleotide variation across the different protein domains of H2A variants was performed by estimating the proportion (p) of nucleotide sites at which the two sequences being compared are different and the numbers of synonymous substitutions per site (p S ), following a sliding-window approach (window length of 20 bp and step size of 5 bp for p, window length of 5 bp and step size of 1 bp for p S ) implemented in the program DnaSP ver. 4.0 .
Gauge of selection and selective constraints acting on H2A.Z-1 and H2A.Z-2 variants
The presence and nature of selection was studied following two strategies: 1) using the codon-based Z-test for selection  comparing the numbers of synonymous and nonsynonymous substitutions per site  in H2A genes, establishing the alternative hypothesis as H1: p N <p S and the null hypothesis as H0: p N = p S ; 2) analyzing deviations from neutrality following two different approaches. First, the influence of selection on overrepresented amino acids was revealed by determining the correlation between the genomic GC content and the proportion of GC-rich (GAPW) and GC-poor amino acids (FYMINK). In this case, while GC-rich amino acids will be positively correlated with genomic GC content and vice-versa under the neutral model [43–46], they will not show any correlation if they are influenced by selection. In both cases, the GC content at fourfold degenerate positions was assumed to represent the genomic GC content, given that the latter has already been shown to be a good approximation of the former  and was also used as a good approximation to the neutral expectation. Correlations were computed using Spearman rank correlation analyses and statistical significance was assessed using standard regression analyses. Second, the effect of mutation and selection bias at the nucleotide level was studied by comparing nucleotide frequencies at first codon positions (always nonsynonymous in the case of the residues studied here) and at fourfold positions (always synonymous). Under the neutral model, nucleotide frequencies should not be significantly different between both positions .
Inference of the phylogenetic relationships among H2A variants
Phylogenetic trees were reconstructed from the obtained p-distances using the neighbor-joining method . To assess that our results are not dependent on this choice, phylogenetic inference analyses were completed by the reconstruction of maximum-parsimony trees  using the close-neighbor-interchange (CNI) search method. The bootstrap  and the interior branch-test [67, 69] methods were combined in order to test the reliability of the trees, producing the bootstrap probability (BS) and the confidence probability (CP) values for each internal branch, assuming BP > 80% and CP ≥ 95% as statistically significant . Histone H2A sequences of the diplomonad protist Giardia were used as outgroups, given that this lineage is believed to be the first to diverge from all other eukaryotes .
This work was funded in part by a Canadian Institutes of Health Research (CIHR) Grant MOP-57718 to J.A. and by a Marie Curie Outgoing International Fellowship (MOIF-CT-2005-021900) within the 6th Framework Programme (European Union) and by a contract within the Isidro Parga Pondal Program (Xunta de Galicia) to J.M.E.-L. R.G.-R. Holds a fellowship from the Diputacion de A Coruña (Spain) and a Predoctoral fellowship from the Universidade da Coruña and T.I. holds a Michael Smith Foundation for Health Research (MSFHR) Postdoctoral Fellowship. J.M.E.-L. would like to thank all the faculty members and students in the EMBO Course 'Computational Biology: From genomes to cells and systems' held in Singapore during August 2008 for fruitful discussions and feedback on an early version of this work.
- van Holde KE: Chromatin. 1988, New York, NY: Springer-VerlagGoogle Scholar
- Eirín-López JM, Li A, Ausió J: H2AX: tailoring histone H2A for chromatin-dependent genomic integrity. Biochem Cell Biol. 2005, 83: 505-515. 10.1139/o05-114.View ArticlePubMedGoogle Scholar
- Dryhurst DD, Thambirajah AA, Ausió J: New twists on H2A.Z: a histone variant with a controversial structural and functional past. Biochem Cell Biol. 2004, 82: 490-497. 10.1139/o04-043.View ArticlePubMedGoogle Scholar
- Abbott DW, Chadwick BP, Thambirajah AA, Ausió J: Beyond the Xi: macroH2A chromatin distribution and post-translational modification in an avian system. J Biol Chem. 2005, 280: 16437-16435. 10.1074/jbc.M500170200.View ArticlePubMedGoogle Scholar
- González-Romero R, Méndez J, Ausió J, Eirín-López JM: Quickly evolving histones, nucleosome stability and chromatin folding: All about histone H2A.Bbd. Gene. 2008, 413: 1-7. 10.1016/j.gene.2008.02.003.View ArticlePubMedGoogle Scholar
- Martianov I, Brancorini S, Catena R, Gansmuller A, Kotaja N, Parvinen M, Sassone-Corsi P, Davidson I: Polar nuclear localization of H1T2, a histone H1 variant, required for spermatid elongation and DNA condensation during spermiogenesis. Proc Natl Acad Sci USA. 2005, 102: 2808-2813. 10.1073/pnas.0406060102.PubMed CentralView ArticlePubMedGoogle Scholar
- Konishi A, Shimizu S, Hirota J, Takao T, Fan y, Matsuoka Y, Zhang L, Yoneda Y, Fuji Y, Skoultchi AI, et al: Involvement of histone H1.2 in apoptosis induced by DNA double-strand breaks. Cell. 2003, 114: 673-688. 10.1016/S0092-8674(03)00719-0.View ArticlePubMedGoogle Scholar
- Tanaka M, Hennebold JD, MacFarlane J, Adashi EY: A mammalian oocyte-specific linker histone gene H1oo: homology with the genes for oocyte-specific cleavage stage histones (cs-H1) of sea urchin and the B4/H1M histone of the frog. Development. 2001, 128: 655-664.PubMedGoogle Scholar
- Parseghian MH, Hamkalo BA: A compendium of the H1 family of somatic subtypes: an elusive cast of characters and their characteristics. Biochem Cell Biol. 2001, 79: 289-304. 10.1139/bcb-79-3-289.View ArticlePubMedGoogle Scholar
- Eirín-López JM, Ausió J: H2A.Z-mediated genome-wide chromatin specialization. Curr Genomics. 2007, 8: 59-66. 10.2174/138920207780076965.PubMed CentralView ArticlePubMedGoogle Scholar
- Malik HS, Henikoff S: Phylogenomics of the nucleosome. Nat Struct Biol. 2003, 10: 882-891. 10.1038/nsb996.View ArticlePubMedGoogle Scholar
- Raisner RM, Madhani HD: Patterning chromatin: form and function for H2A.Z variant nucleosomes. Curr Opin Genet Dev. 2006, 16: 119-124. 10.1016/j.gde.2006.02.005.View ArticlePubMedGoogle Scholar
- Barski A, Cuddapah S, Cui K, Roth TY, Schones DE, Wang Z, Wei G, Chepelev I, Zhao K: High-resolution profiling of histone methylations in the human genome. Cell. 2007, 129: 823-837. 10.1016/j.cell.2007.05.009.View ArticlePubMedGoogle Scholar
- Bruce K, Myers FA, Mantouvalou E, Lefevre P, Greaves I, Bonifer C, Tremethick DJ, Thorne AW, Crane-Robinson C: The replacement histone H2A.Z in hyperacetylated form is a feature of active genes in chicken. Nucl Acids Res. 2005, 33: 5633-5639. 10.1093/nar/gki874.PubMed CentralView ArticlePubMedGoogle Scholar
- Greaves IK, Rangasamy D, Ridgway P, Tremethick DJ: H2A.Z contributes to the unique 3D structure of the centromere. Proc Natl Acad Sci USA. 2007, 104: 525-530. 10.1073/pnas.0607870104.PubMed CentralView ArticlePubMedGoogle Scholar
- Sarcinella E, Zuzarte PC, Lau PN, Draker R, Cheung P: Monoubiquitylation of H2A.Z distinguishes its association with euchromatin or facultative heterochromatin. Mol Cell Biol. 2007, 27: 6457-6468. 10.1128/MCB.00241-07.PubMed CentralView ArticlePubMedGoogle Scholar
- Guillemette B, Gaudreau L: Reuniting the contrasting functions of H2A.Z. Biochem Cell Biol. 2006, 84: 528-535. 10.1139/O06-077.View ArticlePubMedGoogle Scholar
- Zlatanova J, Thakar A: H2A.Z: view from the top. Structure. 2008, 16: 166-179. 10.1016/j.str.2007.12.008.View ArticlePubMedGoogle Scholar
- Suto RK, Clarkson MJ, Tremethick DJ, Luger K: Crystal structure of a nucleosome core particle containing the variant histone H2A.Z. Nat Struct Biol. 2000, 7: 1121-1124. 10.1038/81971.View ArticlePubMedGoogle Scholar
- Park Y-J, Dyer PN, Tremethick DJ, Luger K: A new fluorescence resonance energy transfer approach demonstrates that the histone variant H2AZ stabilizes the histone octamer within the nucleosome. J Biol Chem. 2004, 279: 24274-24282. 10.1074/jbc.M313152200.View ArticlePubMedGoogle Scholar
- Thambirajah AA, Dryhurst DD, Ishibashi T, Li A, Maffey AH, Ausió J: H2A.Z stabilizes chromatin in a way that is dependent on core histone acetylation. J Biol Chem. 2005, 281: 20036-20044. 10.1074/jbc.M601975200.View ArticleGoogle Scholar
- Fan JY, Gordon F, Luger K, Hansen JC, Tremethick DJ: The essential histone variant H2A.Z regulates the equilibrium between different chromatin conformational states. Nat Struct Biol. 2002, 9: 172-176. 10.1038/nsb0402-316b.View ArticlePubMedGoogle Scholar
- Creyghton MP, Markoulaki S, Levine SS, Hanna J, Lodato MA, Sha K, Young RA, Jaenisch R, Boyer LA: H2AZ is enriched at polycomb complex target genes in ES cells and is necessary for lineage commitment. Cell. 2008, 135: 649-661. 10.1016/j.cell.2008.09.056.PubMed CentralView ArticlePubMedGoogle Scholar
- Fu Y, Sinha M, Peterson CL, Weng Z: The insulator binding protein CCF positions 20 nucleosomes around its binding sites across the human genome. PLoS Genet. 2008, 4: e1000138-10.1371/journal.pgen.1000138.PubMed CentralView ArticlePubMedGoogle Scholar
- Coon JJ, Ueberheide B, Syka JE, Dryhurst DD, Ausió J, Shabanowitz J, Hunt DF: Protein identification using sequential ion/ion reactions and tandem mass spectrometry. Proc Natl Acad Sci USA. 2005, 102: 9463-9468. 10.1073/pnas.0503189102.PubMed CentralView ArticlePubMedGoogle Scholar
- Eirín-López JM, Ishibashi T, Ausió J: H2A.Bbd: a quickly evolving hypervariable mammalian histone that destabilizes nucleosomes in an acetylation-independent way. FASEB J. 2008, 22: 316-326. 10.1096/fj.07-9255com.View ArticlePubMedGoogle Scholar
- Chadwick BP, Willard HF: Histone H2A variants and the inactive X chromosome: identification of a second macroH2A variant. Hum Mol Genet. 2001, 10: 1101-1113. 10.1093/hmg/10.10.1101.View ArticlePubMedGoogle Scholar
- González-Romero R, Ausió J, Méndez J, Eirín-López JM: Early evolution of histone genes: prevalence of an 'orphon' H1 lineage in protostomes and birth-and-death process in the H2A family. J Mol Evol. 2008, 66: 505-518. 10.1007/s00239-008-9109-1.View ArticlePubMedGoogle Scholar
- Eirín-López JM, González-Tizón AM, Martínez A, Méndez J: Birth-and-death evolution with strong purifying selection in the histone H1 multigene family and the origin of orphon H1 genes. Mol Biol Evol. 2004, 21 (10): 1992-2003. 10.1093/molbev/msh213.View ArticlePubMedGoogle Scholar
- Eirín-López JM, Lewis JD, Howe L, Ausió J: Common phylogenetic origin of protamine-like (PL) proteins and histone H1: evidence from bivalve PL genes. Mol Biol Evol. 2006, 23: 1304-1317. 10.1093/molbev/msk021.View ArticlePubMedGoogle Scholar
- Eirín-López JM, Ruiz MF, González-Tizón AM, Martínez A, Ausió J, Sánchez L, Méndez J: Common evolutionary origin and birth-and-death process in the replication-independent histone H1 isoforms from vertebrate and invertebrate genomes. J Mol Evol. 2005, 61: 398-407. 10.1007/s00239-004-0328-9.View ArticlePubMedGoogle Scholar
- Nei M, Rooney AP: Concerted and birth-and-death evolution in multigene families. Annu Rev Genet. 2006, 39: 121-152. 10.1146/annurev.genet.39.073003.112240.View ArticleGoogle Scholar
- Piontkivska H, Rooney AP, Nei M: Purifying selection and birth-and-death evolution in the histone H4 gene family. Mol Biol Evol. 2002, 19: 689-697.View ArticlePubMedGoogle Scholar
- Rooney AP, Piontkivska H, Nei M: Molecular evolution of the nontandemly repeated genes of the histone 3 multigene family. Mol Biol Evol. 2002, 19: 68-75.View ArticlePubMedGoogle Scholar
- Luger K, Mäder AW, Richmond RK, Sargent DF, Richmond TJ: Crystal structure of the nucleosome core particle at 2.8 A resolution. Nature. 1997, 389: 251-260. 10.1038/38444.View ArticlePubMedGoogle Scholar
- Dorigo B, Schalch T, Bystricky K, Richmond TJ: Chromatin fiber folding: requirements for the histone H4 N-terminal tail. J Mol Biol. 2003, 327: 85-96. 10.1016/S0022-2836(03)00025-1.View ArticlePubMedGoogle Scholar
- Bao Y, Konesky K, Park Y-J, Rosu S, Dyer PN, Rangasamy D, Tremethick DJ, Laybourn PJ, Luger K: Nucleosomes containing the histone variant H2A.Bbd organize only 118 base pairs of DNA. EMBO J. 2004, 23: 3314-3324. 10.1038/sj.emboj.7600316.PubMed CentralView ArticlePubMedGoogle Scholar
- Arents G, Burlingame RW, Wang BC, Love WE, Moudrianakis E: The nucleosomal core histone octamer at 3.1 A resolution: a tripartite protein assembly and a left-handed superhelix. Proc Natl Acad Sci USA. 1991, 88: 10148-10152. 10.1073/pnas.88.22.10148.PubMed CentralView ArticlePubMedGoogle Scholar
- Eirín-López JM, Abbott DW, Boraston AB: Insight into Ligand Diversity and Novel Biological Roles for Family 32 Carbohydrate Binding Modules. Mol Biol Evol. 2008, 25: 155-157. 10.1093/molbev/msn121.View ArticlePubMedGoogle Scholar
- Eirín-López JM, González-Tizón AM, Martínez A, Méndez J: Molecular and evolutionary analysis of mussel histone genes (Mytilus spp.): possible evidence of an "orphon origin" for H1 histone genes. J Mol Evol. 2002, 55: 272-283. 10.1007/s00239-002-2325-1.View ArticlePubMedGoogle Scholar
- Eirín-López JM, Ruiz MF, González-Tizón AM, Martínez A, Sánchez L, Méndez J: Molecular evolutionary characterization of the mussel Mytilus histone multigene family: first record of a tandemly repeated unit of five histone genes containing an H1 subtype with "orphon" features. J Mol Evol. 2004, 58: 131-144. 10.1007/s00239-003-2531-5.View ArticlePubMedGoogle Scholar
- Eirín-López JM, Frehlick LJ, Ausió J: Long-term evolution and functional diversification in the members of the nucleophosmin/nucleoplasmin family of nuclear chaperones. Genetics. 2006, 173: 1835-1850. 10.1534/genetics.106.058990.PubMed CentralView ArticlePubMedGoogle Scholar
- Jukes TH, Bhushan V: Silent nucleotide substitutions and G+C content of some mitochondrial and bacterial genes. J Mol Evol. 1986, 24: 39-44. 10.1007/BF02099949.View ArticlePubMedGoogle Scholar
- Kimura M: The Neutral Theory of Molecular Evolution. 1983, Cambridge: Cambridge University PressView ArticleGoogle Scholar
- Sueoka N: Correlation between base composition of deoxyribonucleotic acid and composition of proteins. Proc Natl Acad Sci USA. 1961, 47: 1141-1149. 10.1073/pnas.47.8.1141.PubMed CentralView ArticlePubMedGoogle Scholar
- Sueoka N: Directional mutation pressure and neutral molecular evolution. Proc Natl Acad Sci USA. 1988, 85: 2653-2657. 10.1073/pnas.85.8.2653.PubMed CentralView ArticlePubMedGoogle Scholar
- Rooney AP, Zhang J, Nei M: An unusual form of purifying selection in a sperm protein. Mol Biol Evol. 2000, 17: 278-283.View ArticlePubMedGoogle Scholar
- Akashi H, Gojobori T: Metabolic efficiency and amino acid composition in the proteomes of Escherichia coli and Bacillus subtilis. Proc Natl Acad Sci USA. 2002, 99: 3695-36700. 10.1073/pnas.062526999.PubMed CentralView ArticlePubMedGoogle Scholar
- Rooney AP: Selection for highly biased amino acid frequency in the TolA cell envelope protein of proteobacteria. J Mol Evol. 2003, 57: 731-736. 10.1007/s00239-003-2530-6.View ArticlePubMedGoogle Scholar
- Clapier CR, Chakravarthy S, Petosa C, Fernandez-Tornero C, Luger K, Muller CW: Structure of the Drosophila nucleosome core particle highlights evolutionary constraints on the H2A-H2B histone dimer. Proteins. 2008, 71:Google Scholar
- Marino-Ramirez L, Hsu B, Baxevanis AD, Landsman D: The Histone Database: a comprehensive resource for histones and histone fold-containing proteins. Proteins. 2006, 62 (4): 838-842. 10.1002/prot.20814.PubMed CentralView ArticlePubMedGoogle Scholar
- Hall TA: BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucl Acids Symp Ser. 1999, 41: 95-98.Google Scholar
- Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG: The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 1997, 25: 4876-4882. 10.1093/nar/25.24.4876.PubMed CentralView ArticlePubMedGoogle Scholar
- Perez-Bercoff A, Koch J, Burglin TR: LogoBar: bar graph visualization of protein logos with gaps. Bioinformatics. 2006, 22: 112-114. 10.1093/bioinformatics/bti761.View ArticlePubMedGoogle Scholar
- Arnold K, Bordoli L, Kopp J, Schwede T: The SWISS-MODEL Workspace: A web-based environment for protein structure homology modelling. Bioinformatics. 2006, 22: 195-201. 10.1093/bioinformatics/bti770.View ArticlePubMedGoogle Scholar
- DeLano WL: MacPyMOL: A PyMOL-based Molecular Graphics Application for MacOS X. 2007, Palo Alto, CA: DeLano Scientific LLCGoogle Scholar
- Richter S, Wenzel A, Stein M, Gabdoulline RR, Wade RC: webPIPSA: a web server for the comparison of protein interaction properties. Nucleic Acids Res. 2008, 36: W276-280. 10.1093/nar/gkn181.PubMed CentralView ArticlePubMedGoogle Scholar
- Madura JD, Briggs JM, Wade RC, Davis ME, Luty BA, Ilin A, Antosiewicz J, Gilson MK, Bagheri B, Scott LR, et al: Electrostatics and difusion of molecules in solution: simulations with the University of Houston Brownian Dynamics Program. Comp Phys Commun. 1995, 91: 57-95. 10.1016/0010-4655(95)00043-F.View ArticleGoogle Scholar
- Humphrey W, Dalke A, Schulten K: VMD – Visual Molecular Dynamics. J Molec Graph. 1996, 14: 33-38. 10.1016/0263-7855(96)00018-5.View ArticleGoogle Scholar
- Nei M, Kumar S: Molecular Evolution and Phylogenetics. 2000, New York: Oxford University PressGoogle Scholar
- Zhang J, Rosenberg HF, Nei M: Positive Darwinian selection after gene duplication in primate ribonuclease genes. Proc Natl Acad Sci USA. 1998, 95: 3708-3713. 10.1073/pnas.95.7.3708.PubMed CentralView ArticlePubMedGoogle Scholar
- Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol. 2007, 24: 1596-1599. 10.1093/molbev/msm092.View ArticlePubMedGoogle Scholar
- Wright F: The 'effective number of codons' used in a gene. Gene. 1990, 87: 23-29. 10.1016/0378-1119(90)90491-9.View ArticlePubMedGoogle Scholar
- Rozas J, Sánchez-del Barrio JC, Messeguer X, Rozas P: DnaSP, DNA polymorphism analyses by the coalescent and other methods. Bioinformatics. 2003, 19: 2496-2497. 10.1093/bioinformatics/btg359.View ArticlePubMedGoogle Scholar
- Li WH: Molecular Evolution. 1997, Sunderland, MA: SinauerGoogle Scholar
- Saitou N, Nei M: The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987, 4: 406-425.PubMedGoogle Scholar
- Rzhetsky A, Nei M: A simple method for estimating and testing minimum-evolution trees. Mol Biol Evol. 1992, 9: 945-967.Google Scholar
- Felsenstein J: Confidence limits on phylogenies: an approach using the bootstrap. Evolution Int J Org Evolution. 1985, 39: 783-791.View ArticleGoogle Scholar
- Sitnikova T: Bootstrap method of interior-branch test for phylogenetic trees. Mol Biol Evol. 1996, 13: 605-611.View ArticlePubMedGoogle Scholar
- Sitnikova T, Rzhetsky A, Nei M: Interior-branch and bootstrap tests of phylogenetic trees. Mol Biol Evol. 1995, 12: 319-333.PubMedGoogle Scholar
- Roger AJ, Svard SG, Tovar J, Clark CG, Smith MW, Gillin FD, Sogin ML: A mitochondrial-like chaperonin 60 gene in Giardia lamblia: Evidence that diplomonads once harbored an endosymbiont related to the progenitor of mitochondria. Proc Natl Acad Sci USA. 1998, 95: 229-234. 10.1073/pnas.95.1.229.PubMed CentralView ArticlePubMedGoogle Scholar