- Open Access
Epigenetic regulation of the ribosomal cistron seasonally modulates enrichment of H2A.Z and H2A.Zub in response to different environmental inputs in carp (Cyprinus carpio)
Epigenetics & Chromatin volume 6, Article number: 22 (2013)
The specific deposition of histone variants into chromatin is an important epigenetic mechanism that contributes to gene regulation through chromatin architectural changes. The histone variant H2A.Z is essential in higher eukaryotes, and its incorporation within chromatin is a relevant process for gene expression and genome stability. However, the dual positive and negative roles of H2A.Z in gene regulation still remain unclear. We previously reported that acclimatization in common carp fish (Cyprinus carpio) involves cyclical seasonal gene reprogramming as an adaptation response to its natural environment, when rRNA synthesis and processing are profoundly affected. Epigenetic mechanisms primarily contribute to the transcriptional modulation of ribosomal genes concomitant with the acclimatization process, thus significantly regulating this process. The aim of this study was to describe the presence of several H2A.Z subtypes in carp, and assess the role of H2A.Z on the ribosomal cistron in summer- and winter-acclimatized carp.
This paper reports for the first time about the transcriptional expression of four different H2A.Z subtypes belonging to the same organism. Remarkably, a novel H2A.Z.7 was found, which corresponds to a tissue-specific histone subtype that contains seven amino acid residues longer than the canonical H2A.Z. Moreover, H2A.Z enrichment through the ribosomal cistron was significantly higher during summer, when rRNA transcription and processing are highly active, than it was in winter. Similar patterns of H2A.Z enrichment are found in two seasonally active promoters for genes transcribed by RNA polymerase II, the L41 and Δ9-desaturase genes. Interestingly, ubiquitylated-H2A.Z (H2A.Zub) was strongly enriched on regulatory regions of the ribosomal cistron in summer-acclimatized carp. Additionally, H2A.Z was present in both heterochromatin and euchromatin states on ribosomal cistron and RNA polymerase II promoters.
Our study revealed seasonally-dependent H2A.Z enrichment for active ribosomal cistron and RNA polymerase II promoters during the carp environmental adaptation. Moreover, seasonal H2A.Zub enrichment appears as a specific mechanism contributing to the regulation of chromatin architecture under natural conditions. The existence of several H2A.Z subtypes in carp suggests that the epigenetic regulation in this species constitutes a complex and finely tuned mechanism developed to cope with seasonal environmental changes that occur in its habitat.
Living organisms are constantly exposed to environmental fluctuations. In order to survive, they must adapt by inducing compensatory responses . In particular, teleost fish are an interesting model for studying environmental plasticity due to the fact that they are constantly exposed to several types of stressors in the environment. Thus, fish are living under large fluctuations of physical parameters, such as photoperiod, water temperature, nutrition availability, oxygen concentration, and salinity, among others . In this context, the common carp provides a powerful, natural, non-manipulated model for studying several molecular and cellular responses as part of the acclimatization process [2, 3].
It is well known that eukaryotic gene regulation and expression are connected to the structural organization of chromatin. Chromatin, through its capability to transiently pack and unpack the genome, plays a key role in controlling access to genetic information . The basic unit of the chromatin is the nucleosome, which is formed by two copies of each histone protein (H2A, H2B, H3, and H4) assembled into an octamer that has 145 to 147 base pairs (bp) of DNA wrapped around it to form a nucleosome core . In addition to packaging DNA, chromatin is a dynamic structure that adjusts to different cellular processes in response to biotic and abiotic environmental parameters . This flexible and efficient intercommunication between genomic DNA and exogenous influences can be explained by several epigenetic mechanisms, such as DNA methylation , ATP-dependent chromatin remodeling , non-coding RNAs , post-translational modifications (PTMs) of histones, and the replacement of canonical histones with non-allelic histone variants [10, 11]. Among these, the non-random deposition of histone variants into chromatin that occurs independently of DNA replication, substantially contributes to gene regulation and architectural changes in chromosomes .
H2A is the canonical histone that displays a great number of variants, such as macroH2A, H2A.Bbd, H2A.X, and H2A.Z . Among these, H2A.Z is the most evolutionarily ancient one between species, since it is highly conserved during evolution and is widely distributed in all eukaryotic lineages . The identity of H2A.Z compared with the canonical histone H2A is around 60%, suggesting a differential and particular role of this histone variant [15, 16]. Indeed, H2A.Z genes are essential for the development and viability of many organisms, such as Drosophila melanogaster, Tetrahymena thermophila, Xenopus laevis, and mice [20, 21]. Additionally, absence of H2A.Z in Saccharomyces cerevisiae leads to slow growth and a deficiency in transcriptional induction by RNA polymerase II . In vertebrates, two non-allelic genes exist (with a common origin in early chordate evolution) that code for two highly similar H2A.Z proteins, termed H2A.Z.1 and H2A.Z.2 . Even more recently in humans, two alternatively spliced variants for H2A.Z.2, named Z.2.1 and Z.2.2, have been described .
Several PTMs have been reported for H2A.Z, such as acetylation , sumoylation , and ubiquitylation . In particular, the ubiquitylated form of H2A.Z (H2A.Zub) is predominantly found in the inactive X chromosome of female cells, which suggests that this PTM is involved in the formation of silent chromatin. Therefore, the factors that control ubiquitylation of H2A.Z, such as RING1b E3 ligase  and Bmi1 , as well as those that regulate deubiquitylation, such as USP10 , have arisen as important regulators for the epigenetic control of H2A.Z.
rRNA transcription is a demanding and crucial mechanism for all eukaryotic cells because growing cells require a continuous ribosome supply . Particularly in carp, this process undergoes dramatic changes in the context of adaptation to summer and winter . Previously, we demonstrated that epigenetic mechanisms contribute to the transcriptional regulation of ribosomal genes concomitant with the acclimatization process. In this context, we have reported the hypermethylated state of ribosomal gene promoter in winter-acclimatized carp, as well the histone H2A replacement by variant subtypes of macroH2A, help the seasonal transcriptional activity of ribosomal genes [32, 33]. In addition to trans-acting factors, we also reported on the genomic organization of the carp ribosomal cistron , when cis-acting elements located in the intergenic spacer (IGS), named T0´ and T0, and core promoter (CP), can contribute to the modulation of rRNA synthesis.
Genome-wide location analyses have reported a non-random distribution of histone H2A.Z. However, despite intense efforts of many groups, H2A.Z function remains enigmatic. Previously, it was shown that H2A.Z is globally localized in promoter regions . In human and Drosophila cells, the enrichment of H2A.Z at gene promoters is positively correlated with transcription [36, 37], whereas in yeast cells H2A.Z content of the same regions is generally inversely correlated with transcription levels [38, 39].
Taking the above into account, and the limited evidence concerning the epigenetic response as a whole with a model under natural conditions, we evaluated the role of H2A.Z on the ribosomal cistron in summer- and winter-acclimatized carp. We found significant changes of H2A.Z deposition on ribosomal cistron depending on the seasonal transcriptional activity of rRNA genes. Furthermore, because PTMs of H2A.Z are associated with different chromatin states [27, 40, 41], we also found a differential enrichment of H2A.Zub in the IGS and CP regions. In addition, we report for the first time a tissue-specific subtype of H2A.Z, called H2A.Z.7, which is characterized by encoding seven additional amino acid residues, more than canonical H2A.Z.
Seasonal temperature and photoperiod variation in the natural habitat of the common carp in Santiago, Chile
In order to determine periods when abiotic conditions were extreme, we analyzed the temperature and photoperiod as parameters that continuously fluctuate during the year 2011. We found that the months between November and February (summer period) corresponded to the months with longer daylight (13 to 14 hours), and higher mean temperatures (17 to 22°C). Conversely, the period between May and August (winter period) corresponded with shorter daylight (10 to 11 hours), and lower mean temperatures (5 to 7°C) (Additional file 1: Figure S1).
Isolation of different H2A.Z subtypes, and evaluation of tissue-specific expression of H2A.Z and H2A.Z.7 subtypes
We isolated four different sequences that encode for different subtypes of H2A.Z in carp (Additional file 2: Figure S2). The deduced amino acid sequences were obtained and then compared with canonical H2A.Z sequences previously described in the GenBank (Figure 1). The first carp sequence corresponded to H2A.Z.2.1 (100% identity with human H2A.Z.2.1). The second and third sequences found corresponded to two new subtypes, which we named H2A.Z.3.1 and H2A.Z.3.2. Both subtypes differed in only one residue in position 11 (G and S, respectively). The fourth sequence was longer than H2A.Z and contained seven additional amino acids (positions 66 to 72). We named this subtype H2A.Z.7 according to new histone variant nomenclature .
With the purpose of evaluating tissue-specific gene expression of different H2A.Z subtypes in seasonally adapted carp, RT-PCR assays were performed on different tissues. As shown in Figure 2A,B, the presence of the H2A.Z.7 transcript was detected only in the brain and pituitary tissues. Particularly, in the brain, no significant differences in the transcriptional expression of H2A.Z.7 were observed between summer and winter seasons. In contrast, the pituitary H2A.Z.7 transcript was only detected during the winter. On the other hand, we performed RT-PCR assays with primers that did not discriminate between H2A.Z subtypes (H2A.Z.2.1, H2A.Z.3.1, and H2A.Z.3.2). This amplification was denoted as H2A.Z, and it was detected in all examined tissues (Figure 2B). Higher levels of the H2A.Z transcript were detected in testis of summer-acclimatized fish, while the same tendency was observed in the liver during winter-acclimatized fish. No significant variation in the content of H2A.Z was observed between seasons in the brain, pituitary, and gills.
The phylogenetic relationship of carp H2A.Z subtypes
In order to analyze the phylogenetic relationship among the new H2A.Z subtypes and the previously described H2A.Z, we constructed a phylogenetic tree using 29 coding DNA sequences (CDS) from the GenBank, including H2A.Z from three different kingdoms, such as Animalia, Fungi, and Protist (Figure 3). According to the phylogenetic reconstruction, subtypes of the histone variant H2A.Z isolated from carp, together with the group comprising H2A.Z variants belonging to the Actinopterygii class of fish, share a monophyletic origin. Furthermore, these sequences cluster together with the H2A.Z from other teleost fish, such as Danio rerio, Salmo salar, Pagrus major, Oncorhynchus mykiss, and Esox lucius.
The analysis revealed that H2A.Z subtypes from carp, and especially H2A.Z.7, are subject to a marked diversification process within this vertebrate group of animals.
Nuclear content of the histone H2A.Z varies seasonally in carp hepatocytes
In order to validate and confirm that commercial antibodies developed against mouse H2A.Z cross-react with samples of carp hepatocytes, comparative western blot with rat and human cells nuclear extracts was performed. We found a specific signal of approximately 14 kDa corresponding to carp H2A.Z, indicating that antibodies detect their orthologs in carp (Figure 4A).
Subsequently, we analyzed the protein content of H2A.Z in hepatocytes from seasonally adapted carp. We found that in winter, H2A.Z is 1.7 times higher compared to summer-adapted carp (Figure 4B,C), a result that is consistent with the liver tissue-specific RT-PCR assay shown in Figure 2.
Since the extracellular matrix contributes to the autofluorescence emission in hepatocyte tissue, we performed an H2A.Z immunodetection assay in pituitary tissue and found similar results as those previously observed in liver. Unexpectedly, nuclear protein content of H2A.Z from pituitary cells was much higher in winter-acclimatized fish (Figure 4D).
Seasonally-dependent distribution of H2A.Z throughout the ribosomal cistron in the carp
In order to evaluate the possible differences of H2A.Z enrichment along the carp ribosomal cistron during acclimatization, we used chromatin immunoprecipitation (ChIP) assays and subsequent quantification by quantitative PCR (qPCR) (Figure 5A). The results showed a sevenfold increase (± SD) of H2A.Z enrichment in IGS, T0´, T0, and CP regions, and a threefold increase (± SD) in the 5.8S region during summer (Figure 5B). A H3 ChIP control to determine nucleosome occupancy during seasons at investigated loci is shown in Figure 5D.
Additionally, with the purpose of correlating the differential seasonal enrichments of H2A.Z with general gene transcriptional activity, we evaluated the association of H2A.Z with the promoter regions of three RNA polymerase II-related genes. Two of them, L41 and Δ9-desaturase genes, exhibit seasonally-dependent expression in carp, whereas JAK1 gene has a constitutive expression between seasons (see Additional file 3: Figure S3) (Figure 5C). Previously, it has been described that L41 and Δ9-desaturase genes display increased transcriptional activity in hepatocytes from summer- and winter-acclimatized carp, respectively [43, 44]. In agreement with this observation, H2A.Z was systematically enriched (by approximately three times) on promoters when genes are seasonally activated. In contrast, on the promoter of the JAK1-constitutive gene, there is no significant difference of H2A.Z enrichment between seasons.
Deposition of H2A.Zub in summer-acclimatized carp throughout the ribosomal cistron
The monoubiquitylation of H2A.Z distinguishes its association with euchromatin or facultative heterochromatin . Consequently, we decided to study whether or not there is a relationship between the H2A.Zub deposition on ribosomal cistron with the transcriptional activity of carp rRNA. The results showed that during summer, there was an increase of fourfold, threefold and sixfold (± SD) of H2A.Zub enrichment in IGS, T0 and CP regions, respectively (Figure 6A); whereas, H2A.Zub enrichment in 5.8S was not significant. To determine whether these findings are exclusive in ribosomal cistron or not, we evaluated the association of H2A.Zub with the promoter regions of the ribosomal protein L41 and Δ9-desaturase genes. In both promoter regions, H2A.Zub enrichment was not significant (Figure 6B).
Moreover, we performed RT-qPCR assays against Bmi1, Ring1B, and USP10 ubiquitylation/deubiquitylation factors to reveal any possible variations in the ubiquitylation machinery concurrent to seasonal changes (Figure 6C). The mRNA levels of Bmi1 and RING1b were fourfold and sixfold higher, respectively, in the summer compared to winter. Conversely, there were no changes in the mRNA content of USP10 between seasons.
Correlation between the enrichment of H2A.Z and epigenetics marks
We define colocalization as the presence in the same chromatin region of H2A.Z, and either euchromatin (acetylated histone H4 at lysine 12 (H4K12ac)) or heterochromatin (trimethylated histone H3 at lysine 9 (H3K9me3)) markers. Using this colocalization, we can correlate H2A.Z occupancy with chromatin transcriptional states of the rDNA cistron and two RNA polymerase II promoters.
We found a similar proportion of colocalization between H2A.Z and H4K12ac between seasons along the rDNA cistron (Figure 7A). In contrast, when we analyzed the colocalization between H2A.Z and H3K9me3, this was significantly greater during winter on the T0´, T0, and CP regions (Figure 7B).
H2A.Z enrichment was higher on the promoter of the Δ9-desaturase gene during winter, co-localizing with the euchromatin marker H4K12ac. The same colocalization between H2A.Z and H4K12ac was found in the summer-active L41 promoter (Figure 7C). In contrast, H2A.Z enrichment co-localized with the heterochromatin marker (H3K9me3) when the transcriptional activity of gene expression is seasonally repressed (Figure 7D), as in the case of L41 during winter and Δ9-desaturase during summer.
In carp, seasonal acclimatization involves a reprogramming of molecular and cellular functions, which begin by sensing environmental clues (temperature and photoperiod) that occur during an annual cycle. This process entails a fine modulation of a vast array of genes in order to provide allostasis. At a cellular level, one of the most dramatic differences corresponds to the reorganization of the nucleolar components observed during seasonal adaptation, when rRNA transcription and processing are clearly affected [45, 46]. In this context, we previously reported that epigenetic mechanisms are clearly involved in the phenotypical plasticity that is achieved by the carp. In particular, the hypermethylated state of ribosomal gene promoter in winter-acclimatized carp, as well as histone H2A replacement by variant subtypes of macroH2A, is dependent on the seasonal transcriptional activity of ribosomal genes, corroborating that changes in the architecture of the chromatin are reversible and dynamically respond to environmental changes [32, 33].
In the present study, the expression of H2A.Z was detected in all tissues studied. However, we observed that the total content of H2A.Z transcripts can change in the liver and testis in a seasonally-dependent manner. The same difference was found in western blot experiments, when liver H2A.Z content tended to be higher during the winter than summer. On the other hand, in immunodetection assays from pituitary cells, protein content of H2A.Z was significantly higher in winter-acclimatized fish. Nevertheless, this result is not correlated with mRNA levels. It is important to notice that in the present study, polyclonal antibodies were used for immunodetection and these antibodies cannot differentiate among different H2A.Z subtypes because they are directed to the extremely conserved C-terminal region of the protein; thus, limiting us to relate mRNA levels with protein levels of different H2A.Z subtypes.
We report for the first time four distinct subtypes of histone H2A.Z in the same organism. One of the deduced amino acid sequences from carp shares a 100% identity with the H2A.Z.2.1 of humans , whereas the other two (H2A.Z.3.1 and H2A.Z.3.2) do not correspond to any H2A.Z described until now. Between them, they share a 99.2% identity, differing only in Gly11 to Ser. In carp, we could not identify the existence of histone H2A.Z.1, but we do not discard its existence in other tissues different to those studied in this work.
Unexpectedly, we identified a new subtype of the histone variant H2A.Z (which was designated as H2A.Z.7) that is similar to H2A.Z.3.2, but possesses seven additional residues, more than other H2A.Z variants. The extra codons in H2A.Z.7 mRNA are in the same location as an intron in human H2A.Z, most likely indicating that H2A.Z.7 is a splicing (or mis-splicing) variant of H2A.Z.3.2. These seven additional residues are located in the region that corresponds to α-helix 2 , which could deviate H2A.Z C-terminal docking domain, most likely affecting the interaction surface with (H3-H4)2 tetramer. Clearly, structural studies should be carried out to determine the functional implications of these additional residues. However, considering we have no indication that the H2A.Z.7 protein is translated, it could also be possible that this corresponds to a pseudogene.
H2A.Z.7 presents a tissue-specific expression, which was detected only in the brain (summer and winter) and pituitary tissues (winter). This is not the first time that an unusual length of H2A.Z has been described. In human tissues and cell lines, a shorter version named H2A.Z.2.2 was found, and whose incorporation leads to severely unstable nucleosomes due to it containing a characteristic-docking domain .
Compared to mammals, the presence of this particular and species-specific histone H2A.Z.7 could be due to the three whole genome duplications (WGDs) that teleost have experienced. In particular, carp are believed to have had another round of genome duplications (4R) and became an evolutionarily recent tetraploid fish . This suggests that the presence of H2A.Z.7, and the other forms of carp H2A.Z, could be a consequence of WGD as a selective and important mechanism for eukaryote genome evolution in response to biotic and abiotic stress in carp.
The clustering pattern of carp H2A.Z subtypes in the phylogenetic tree is consistent with the differentiations and diversifications of teleost fish. So far, H2A.Z.7 from carp has an intense evolutionary diversification between other H2A.Z from mammals and, even more, within teleost fish. The presence of more than one subtype of H2A.Z is not a whim of random genetic drift, but different forms may have acquired new or complementary functions . Future experiments must show if the new subtypes described in this work perform similar or different roles in chromatin.
The incorporation of histone H2A.Z within chromatin is crucial for proper gene regulation [17, 22]. Genome-wide studies in several organisms show that H2A.Z is found throughout the genome; however, H2A.Z deposition is not random on chromatin, but rather specific and highly regulated [35, 49]. In yeast, H2A.Z is distributed in active and inactive RNA polymerase II promoters [35, 38, 39], whereas in higher eukaryotes it is mostly in active genes . Interestingly, we report for the first time the seasonally-dependent enrichment of H2A.Z throughout the ribosomal cistron, concomitant with differential modulation of rDNA activity using a natural model of environmental plasticity, such as the carp fish. We observed that the levels of H2A.Z enrichment were higher in the ribosomal cistron during summer, a season in which rRNA transcription is highly active . This particular evidence confirms that the correlation between H2A.Z enrichment with the ribosomal transcriptional activity is not only associated to RNA polymerase I-related genes, since we detected that the enrichment of H2A.Z agrees with seasonally active gene promoters that are transcribed by RNA polymerase II.
In higher eukaryotes, p400 and SRCAP (SWI2/SNF2-related CBP activator protein) complexes are ATP-dependent chromatin remodelers specifically exchanging canonical H2A–H2B with H2A.Z–H2B dimers within the nucleosomes [51, 52]. In carp, both glycogen and lipid content in hepatocytes from acclimatized carp differ significantly between seasons . Taking into account that ribosome production is a major biosynthetic and energy-consuming activity of eukaryotic cells, which adapts rapidly to changes in intracellular energy status , it is plausible to speculate that variations of cellular AMP/ATP ratios could modulate H2A.Z–H2B incorporation by p400 and SRCAP complexes in an ATP-availability manner. Certainly, this speculation must be demonstrated experimentally.
Additionally, in order to corroborate if H2A.Z enrichment is consistent with a chromatin state, we evaluated H2A.Z colocalization with specific epigenetic markers of active (H4K12ac)  or inactive (H3K9me3)  gene expression. Our results showed that H2A.Z and the euchromatin marker (H4K12ac) co-localize along the ribosomal cistron but not in a seasonally-dependent manner. In contrast, we showed a significant colocalization between H3K9me3 and H2A.Z on rDNA cistron during winter, a season in which rRNA transcription and processing are highly inactive. Then, when we compared the colocalization of H2A.Z with these epigenetic markers on RNA polymerase II genes, they were consistent with the active or inactive gene states. These results agree with the controversial role of H2A.Z in the transcription of RNA polymerase II genes .
Histone tails are subjected to PTMs that alter the higher-order folding of chromatin. In particular, H2A.Z has been reported to be acetylated, sumoylated, and ubiquitylated. The H2A.Zub marker is involved in the maintenance or formation of heterochromatin because it was found predominantly on the inactive X chromosome . In this study, we investigated if H2A.Z ubiquitylation in the ribosomal cistron can contribute to regulate carp seasonal rRNA transcription. First, we evaluated transcriptional expression of the most relevant ubiquitylation factors. We found a greater gene expression of RING1b and Bmi1 in summer than in winter, whereas USP10 showed no change between seasons. These findings are consistent with the relevant enrichment of H2A.Zub on the ribosomal cistron of summer-acclimatized carp, when the activity of rRNA transcription is significantly higher. Unlike the uniform pattern of H2A.Z on the ribosomal cistron during summer, H2A.Zub enrichment in the same season was found mainly on the IGS and CP regions. These two locations of the rDNA gene might play a crucial role in rRNA transcription since H2A ubiquitylation (and consequently H2A.Zub) can serve as a docking site for recruiting other factors, or preventing chromatin access to transcriptional regulators . On the other hand, at least in hepatocytes isolated from carp, the two promoters of RNA polymerase II studied did not show significant enrichment of H2A.Zub in any season. Nevertheless, similar findings have been reported in Saccharomyces cerevisiae, where H2BK123ub has negative consequences on the assembly of RNA polymerase II at promoters and stimulates transcription elongation . Finally, it would certainly be desirable to obtain antibodies specific for carp H2A.Zub to unambiguously corroborate the ubiquitylation state of carp H2A.Z. Additionally, wide-genome screening of H2A.Zub should be performed in order to better understand the contribution of this PTM in transcriptional regulation during carp seasonal acclimatization.
Our findings reveal seasonal H2A.Z enrichment in active genes of the ribosomal cistron and RNA polymerase II promoters. In the same context, it seems that H2A.Zub enrichment is a specific marker of post-translational histone modification on active ribosomal cistron during the carp acclimatization process. Furthermore, we identified for the first time four distinct subtypes for the histone H2A.Z in the same organism. Further research is necessary to distinguish the particular contribution of each of the H2A.Z subtype within chromatin, and how the predicted monoubiquitylation of H2A.Z modulates gene transcription.
Animal and tissue preparation
All experimental work on fish described in this paper was approved by the Bioethics Committee of the Universidad Andres Bello, Santiago, Chile. Male carp were maintained in seasonal and environmental conditions as previously described . The fish were sacrificed, and tissue was collected from the liver, brain, pituitary, testis, and gills, washed with PBS, pH 7.4, and stored at −80°C until required.
Isolation of carp H2A.Z sequences
From Danio rerio, H2A.Z mRNA [NCBI Reference Sequence:NM_001201563.1] were deduced as the following oligonucleotides: Z-380pb, F 5′-GCAGGTGGAAAAGCAGGT AAAG-3′ and Z-380pb, R 5′-GCGGTTTTCTGCTGGCCCTTCTT-3′. A 380 bp fragment corresponding to a partial mRNA region of H2A. Z was amplified by RT-PCR using carp liver tissue as a template. The product was cloned in the pGEM-T vector (Promega, Madison, WI, USA) and sequenced. Full mRNA of different carp H2A.Z subtypes were isolated using FirstChoice RLM-RACE (Ambion, Foster City, CA, USA) with specific oligonucleotides deduced from the fragments described earlier. The 5′UTR region was completed with Z-380pb, R and 5′Z inner 5′-CCAGCTGTAGATGA CGAG-3′, whereas the 3′UTR region was completed with Z-380 pb, F and 3′Z inner 5′-GTCTACAGTGCAGCCATC-3′. All PCR reactions were performed with Phusion High-Fidelity DNA Polymerase (New England Biolabs, Ipswich, MA, USA) and DNA sequencing services were performed by Macrogen (Tokyo, Japan).
A total of 100 mg of pituitary tissue from summer- and winter-acclimatized carp were fixed using 4% paraformaldehyde (w/v). Sections were embedded in Paraplast and cut at 5 μm in a cryostat Leica CM1850 (Leica Microsystems, Wetzlar, Germany). The primary antibody used was anti-H2A.Z (1:500, ab4174; Abcam, Cambridge, MA, USA) incubated at 4°C for 16 hours and then washed three times with PBS. The secondary antibody was anti-rabbit IgG (1:400) conjugated to Alexa Fluor 647 (Cell Signaling Technology, Danvers, MA, USA) and 4',6-diamidino-2-phenylindole (DAPI) (1:1000), which was incubated for 1 hour at 25°C. Confocal microscopy was carried out in Axiovert 100 M (Zeiss, Oberkochen, Germany) and analyzed with ImageJ software (National Institutes of Health (NIH), Bethesda, MD, USA).
RT-PCR and RT-qPCR
Total RNA was purified from tissues, such as the liver, brain, pituitary, testis, and gills using TRIzol reagent (Invitrogen, Carlsbad, CA, USA). The quality of RNA was evaluated on 1.2% agarose gel electrophoresis and quantified spectrophotometrically at A260nm. The cDNA was synthesized from 2 μg of total RNA with M-MLV Reverse Transcriptase (Promega) and random primers according to the manufacturer’s instructions. The resulting cDNAs were PCR-amplified in GoTaq Green Master Mix (Promega) using oligonucleotides, showed in Additional file 4: Table S1.
RT-qPCR was performed with the appropriate oligonucleotides (Additional file 4: Table S1) and Brilliant SYBR Master Mix (Stratagene, La Jolla, CA, USA) in a final volume of 20 μl on a Mx3000P Real-Time PCR System (Stratagene). Seasonal average ratio of mRNA levels was calculated using the Pfaffl method  and carp 16S mitochondrial rRNA as a seasonally constitutive gene expression.
Nuclear extracts from HeLa cells, and rat and carp hepatocytes were prepared according to Araya et al. . For the immune detections, we used polyclonal anti-H2A.Z (1:1500; Abcam) and anti-H2B (1:5000; Abcam) antibodies. The peroxidase-conjugated anti-rabbit IgG (H+L) secondary antibodies were obtained from KPL (Gaithersburg, MD, USA). Protein-antibody complexes were visualized by an enhanced chemiluminescence (ECL) detection system (GE Healthcare, Piscataway, NJ, USA).
Chromatin immunoprecipitation combined with qPCR
ChIP and sequential chromatin immunoprecipitation (ChIP-reChIP) assays were undertaken as previously described . For each assay, we used 3 × 107 nuclei isolated from the liver of summer- and winter-acclimatized carp, respectively. ChIPs were performed with polyclonal anti-H2A.Z (ab4174; Abcam) and anti-H3 (ab1791; Abcam). The ChIP-reChIP assays were performed immunoprecipitating first with anti-H2A.Z antibodies, and then an additional immunoprecipitation with polyclonal histone H4 (acetyl K12), anti-histone H3 (Tri-methyl K9; Abcam), and anti-ubiquitin (Cell Signaling Technology). In the last case, although we did not have antibodies for carp H2A.Zub, our results were in agreement with H2A.Zub form because all data values obtained were higher than the unspecific IgG IP (normalizer). The immunoprecipitated DNA was amplified by qPCR on an Mx3000P Real-Time PCR System (Stratagene) using oligonucleotides corresponding to IGS, T0`, T0, CP and 5.8S regions of carp rDNA (Figure 5A), and the promoter regions of L41 [GenBank:AY115478], Δ9-desaturase [GenBank:BX571725.1] and JAK1 [GenBank:AH004872.1] genes (Additional file 4: Table S1).
ChIP-qPCR was performed using SYBR Master Mix (Stratagene). The specificity of the PCR reaction was assessed by melting-temperature profiles. The fold enrichment was calculated using the ΔΔCt method :
Where Ctx corresponds to the Ct value for each specific antibody used. All assays were performed in triplicate and repeated with three different animal samples for each season. Results represent the mean (± SD) of three experiments.
Phylogenetic analysis of H2A.Z variants
The CDS from different eukaryotes corresponding to H2A.Z were retrieved from the GenBank database (Additional file 5: Table S2). A total of 29 sequences were aligned using ClustalW alignment  and analyzed by phylogeny reconstruction using the minimum evolution (ME) method and the Close-Neighbor-Interchange (CNI) algorithm. The tree is drawn to scale, with branch lengths in the same units as those of the evolutionary distances used to infer the phylogenetic tree. The evolutionary distances were computed using the p-distance method. The bootstrap consensus tree inferred from 5,000 replicates was taken to represent the evolutionary history of the taxa analyzed. Evolutionary analyses were conducted in MEGA 5.05 software (Tokyo, Japan) .
Coding DNA sequences
Chromatin immunoprecipitation, ChIP-reChIP, sequential chromatin immunoprecipitation
Trimethylated histone H3 at lysine 9
acetylated histone H4 at lysine 12
Phosphate buffered saline
Polymerase chain reaction
Reverse transcription polymerase chain reaction
University College Dublin
Whole genome duplication.
López-Maury L, Marguerat S, Bähler J: Tuning gene expression to changing environments: from rapid responses to evolutionary adaptation. Nat Rev Genet. 2008, 9: 583-593.
Cossins AR, Crawford DL: Fish as models for environmental genomics. Nat Rev Genet. 2005, 6: 324-333. 10.1038/nrg1590.
Alvarez M, Molina A, Quezada C, Pinto R, Krauskopf M, Vera M: Eurythermal fish acclimatization and nucleolar function: a review. J Therm Biol. 2004, 29: 663-667. 10.1016/j.jtherbio.2004.08.036.
Vaquero A, Loyola A, Reinberg D: The constantly changing face of chromatin. Sci Aging Knowledge Environ. 2003, 2003: RE4.
Luger K, Mäder AW, Richmond RK, Sargent DF, Richmond TJ: Crystal structure of the nucleosome core particle at 2.8 Å resolution. Nature. 1997, 389: 251-260. 10.1038/38444.
Bollati V, Baccarelli A: Environmental epigenetics. Heredity. 2010, 105: 105-112. 10.1038/hdy.2010.2.
Rottach A, Leonhardt H, Spada F: DNA methylation-mediated epigenetic control. J Cell Biochem. 2009, 108: 43-51. 10.1002/jcb.22253.
Peterson CL: ATP-dependent chromatin remodeling: going mobile. FEBS Lett. 2000, 476: 68-72. 10.1016/S0014-5793(00)01673-2.
Esteller M: Non-coding RNAs in human disease. Nat Rev Genet. 2011, 12: 861-874.
Henikoff S, Ahmad K: Assembly of variant histones into chromatin. Annu Rev Cell Dev Biol. 2005, 21: 133-153. 10.1146/annurev.cellbio.21.012704.133518.
Kouzarides T: Chromatin modifications and their function. Cell. 2007, 128: 693-705. 10.1016/j.cell.2007.02.005.
Jin J, Cai Y, Li B, Conaway RC, Workman JL, Conaway JW, Kusch T: In and out: histone variant exchange in chromatin. Trends Biochem Sci. 2005, 30: 680-687. 10.1016/j.tibs.2005.10.003.
Kamakaka RT: Histone variants: deviants?. Genes Dev. 2005, 19: 295-316. 10.1101/gad.1272805.
Iouzalen N, Moreau J, Méchali M: H2A.ZI, a new variant histone expressed during Xenopus early development exhibits several distinct features from the core histone H2A. Nucleic Acids Res. 1996, 24: 3947-3952. 10.1093/nar/24.20.3947.
Thatcher TH, Gorovsky MA: Phylogenetic analysis of the core histones H2A, H2B, H3, and H4. Nucleic Acids Res. 1994, 22: 174-179. 10.1093/nar/22.2.174.
Jackson J, Gorovsky MA: Histone H2A.Z has a conserved function that is distinct from that of the major H2A sequence variants. Nucleic Acids Res. 2000, 28: 3811-3816. 10.1093/nar/28.19.3811.
Clarkson MJ, Wells JR, Gibson F, Saint R, Tremethick DJ: Regions of variant histone His2AvD required for Drosophila development. Nature. 1999, 399: 694-697. 10.1038/21436.
Liu X, Bowen J, Gorovsky MA: Either of the major H2A genes but not an evolutionarily conserved H2A.F/Z variant of Tetrahymenathermophila can function as the sole H2A gene in the yeast Saccharomyces cerevisiae. Mol Cell Biol. 1996, 16: 2878-2887.
Ridgway P, Brown KD, Rangasamy D, Svensson U, Tremethick DJ: Unique residues on the H2A.Z containing nucleosome surface are important for Xenopuslaevis development. J Biol Chem. 2004, 279: 43815-43820. 10.1074/jbc.M408409200.
Faast R, Thonglairoam V, Schulz TC, Beall J, Wells JR, Taylor H, Matthaei K, Rathjen PD, Tremethick DJ, Lyons I: Histone variant H2A.Z is required for early mammalian development. Curr Biol. 2001, 11: 1183-1187. 10.1016/S0960-9822(01)00329-3.
Nashun B, Yukawa M, Liu H, Akiyama T, Aoki F: Changes in the nuclear deposition of histone H2A variants during pre-implantation development in mice. Development. 2010, 137: 3785-3794. 10.1242/dev.051805.
Adam M, Robert F, Larochelle M, Gaudreau L: H2A.Z is required for global chromatin integrity and for recruitment of RNA polymerase II under specific conditions. Mol Cell Biol. 2001, 21: 6270-6279. 10.1128/MCB.21.18.6270-6279.2001.
Eirín-López JM, González-Romero R, Dryhurst D, Ishibashi T, Ausió J: 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. BMC Evol Biol. 2009, 9: 31-10.1186/1471-2148-9-31.
Bonisch C, Schneider K, Punzeler S, Wiedemann SM, Bielmeier C, Bocola M, Eberl HC, Kuegel W, Neumann J, Kremmer E, Leonhardt H, Mann M, Michaelis J, Schermelleh L, Hake SB: H2A.Z.2.2 is an alternatively spliced histone H2A.Z variant that causes severe nucleosome destabilization. Nucleic Acids Res. 2012, 40: 5951-5964. 10.1093/nar/gks267.
Millar CB, Xu F, Zhang K, Grunstein M: Acetylation of H2AZ Lys 14 is associated with genome-wide gene activity in yeast. Genes Dev. 2006, 20: 711-722. 10.1101/gad.1395506.
Kalocsay M, Hiller NJ, Jentsch S: Chromosome-wide Rad51 spreading and SUMO-H2A.Z-dependent chromosome fixation in response to a persistent DNA double-strand break. Mol Cell. 2009, 33: 335-343. 10.1016/j.molcel.2009.01.016.
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.
Ben-Saadon R, Zaaroor D, Ziv T, Ciechanover A: The polycomb protein Ring1B generates self atypical mixed ubiquitin chains required for its in vitro histone H2A ligase activity. Mol Cell. 2006, 24: 701-711. 10.1016/j.molcel.2006.10.022.
Cao R, Tsukada Y, Zhang Y: Role of Bmi-1 and Ring1A in H2A ubiquitylation and Hox gene silencing. Mol Cell. 2005, 20: 845-854. 10.1016/j.molcel.2005.12.002.
Draker R, Sarcinella E, Cheung P: USP10 deubiquitylates the histone variant H2A.Z and both are required for androgen receptor-mediated gene activation. Nucleic Acids Res. 2011, 39: 3529-3542. 10.1093/nar/gkq1352.
Grummt I, Pikaard CS: Epigenetic silencing of RNA polymerase I transcription. Nat Rev Mol Cell Biol. 2003, 4: 641-649. 10.1038/nrm1171.
Pinto R, Ivaldi C, Reyes M, Doyen C, Mietton F, Mongelard F, Alvarez M, Molina A, Dimitrov S, Krauskopf M, Vera MI, Bouvet P: Seasonal environmental changes regulate the expression of the histone variant macroH2A in an eurythermal fish. FEBS Lett. 2005, 579: 5553-5558. 10.1016/j.febslet.2005.09.019.
Araya I, Nardocci G, Morales J, Vera M, Molina A, Alvarez M: MacroH2A subtypes contribute antagonistically to the transcriptional regulation of the ribosomal cistron during seasonal acclimatization of the carp fish. Epigenetics Chromatin. 2010, 3: 14-10.1186/1756-8935-3-14.
Vera MI, Molina A, Pinto R, Reyes M, Alvarez M, Krauskopf E, Quezada C, Torres J, Krauskopf M: Genomic organization of the rDNAcistron of the teleost fish Cyprinuscarpio. Biol Res. 2003, 36: 241-251.
Zhang H, Roberts DN, Cairns BR: Genome-wide dynamics of Htz1, a histone H2A variant that poises repressed/basal promoters for activation through histone loss. Cell. 2005, 123: 219-231. 10.1016/j.cell.2005.08.036.
Barski A, Cuddapah S, Cui K, Roh T, 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.
Mavrich TN, Jiang C, Ioshikhes IP, Li X, Venters BJ, Zanton SJ, Tomsho LP, Qi J, Glaser RL, Schuster SC, Gilmour DS, Albert I, Pugh BF: Nucleosome organization in the Drosophila genome. Nature. 2008, 453: 358-362. 10.1038/nature06929.
Guillemette B, Bataille AR, Gévry N, Adam M, Blanchette M, Robert FO, Gaudreau L: Variant histone H2A.Z is globally localized to the promoters of inactive yeast genes and regulates nucleosome positioning. PLoS Biol. 2005, 3: e384-10.1371/journal.pbio.0030384.
Li B, Pattenden SG, Lee D, Gutiérrez J, Chen J, Seidel C, Gerton J, Workman JL: Preferential occupancy of histone variant H2AZ at inactive promoters influences local histone modifications and chromatin remodeling. Proc Natl Acad Sci U S A. 2005, 102: 18385-18390. 10.1073/pnas.0507975102.
Hardy S, Jacques P, Gévry N, Forest A, Fortin M, Laflamme L, Gaudreau L, Robert F: The Euchromatic and heterochromatic landscapes are shaped by antagonizing effects of transcription on H2A.Z Deposition. PLoS Genet. 2009, 5: e1000687-10.1371/journal.pgen.1000687.
Valdes-Mora F, Song JZ, Statham AL, Strbenac D, Robinson MD, Nair SS, Patterson KI, Tremethick DJ, Stirzaker C, Clark SJ: Acetylation of H2A.Z is a key epigenetic modification associated with gene deregulation and epigenetic remodeling in cancer. Genome Res. 2012, 22: 307-321. 10.1101/gr.118919.110.
Talbert PB, Ahmad K, Almouzni G, Ausió J, Berger F, Bhalla PL, Bonner WM, Cande WZ, Chadwick BP, Chan SWL, Cross GAM, Cui L, Dimitrov SI, Doenecke D, Eirín-López JM, Gorovsky MA, Hake SB, Hamkalo BA, Holec S, Jacobsen SE, Kamieniarz K, Khochbin S, Ladurner AG, Landsman D, Latham JA, Loppin B, Malik HS, Marzluff WF, Pehrson JR, Postberg J, et al: A unified phylogeny-based nomenclature for histone variants. Epigenetics Chromatin. 2012, 5: 7-10.1186/1756-8935-5-7.
Gracey A, Macartney A, Beynon R, Cossins A: Cold-induced expression of Δ9-desaturase in carp by transcriptional and posttranslational mechanisms. Science. 1996, 271: 815-818. 10.1126/science.271.5250.815.
Molina A, Corta A, Martin R, Alvarez M, Burzio LO, Krauskopf M, Vera MI: Gene structure of the carp fish ribosomal protein L41: seasonally regulated expression. Biochem Biophys Res Commun. 2002, 295: 582-586. 10.1016/S0006-291X(02)00715-5.
Vera MI, Norambuena L, Alvarez M, Figueroa J, Molina A, Leon G, Krauskopf M: Reprogramming of nucleolar gene expression during the acclimatization of the carp. Cell Mol Biol Res. 1993, 39: 665-674.
Alvarez M, Quezada C, Molina A, Krauskopf M, Vera MI, Thiry M: Ultrastructural changes of the carp (Cyprinus carpio) hepatocyte nucleolus during seasonal acclimatization. Biol Cell. 2006, 98: 457-463. 10.1042/BC20060006.
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.
Danzmann RG, Davidson EA, Ferguson MM, Gharbi K, Koop BF, Hoyheim B, Lien S, Lubieniecki KP, Moghadam HK, Park J, Phillips RB, Davidson WS: Distribution of ancestral proto-Actinopterygian chromosome arms within the genomes of 4R-derivative salmonid fishes (Rainbow trout and Atlantic salmon). BMC Genomics. 2008, 9: 557-10.1186/1471-2164-9-557.
Billon P, Côté J: Precise deposition of histone H2A.Z in chromatin for genome expression and maintenance. Biochim Biophys Acta. 2012, 1819: 290-302. 10.1016/j.bbagrm.2011.10.004.
Weber CM, Henikoff JG, Henikoff S: H2A.Z nucleosomes enriched over active genes are homotypic. Nat Struct Mol Biol. 2010, 17: 1500-1507. 10.1038/nsmb.1926.
Ruhl DD, Jin J, Cai Y, Swanson S, Florens L, Washburn MP, Conaway RC, Conaway JW, Chrivia JC: Purification of a human SRCAP complex that remodels chromatin by incorporating the Histone variant H2A.Z into nucleosomes. Biochemistry. 2006, 45: 5671-5677. 10.1021/bi060043d.
Cai Y: The mammalian YL1 protein is a shared subunit of the TRRAP/TIP60 histone acetyltransferase and SRCAP complexes. J Biol Chem. 2005, 280: 13665-13670. 10.1074/jbc.M500001200.
Sáez L, Zuvić T, Amthauer R, Rodríguez E, Krauskopf M: Fish liver protein synthesis during cold acclimatization: seasonal changes of the ultrastructure of the carp hepatocyte. J Exp Zool. 1984, 230: 175-186. 10.1002/jez.1402300203.
Murayama A, Ohmori K, Fujimura A, Minami H, Yasuzawa-Tanaka K, Kuroda T, Oie S, Daitoku H, Okuwaki M, Nagata K, Fukamizu A, Kimura K, Shimizu T, Yanagisawa J: Epigenetic control of rDNA loci in response to intracellular energy status. Cell. 2008, 133: 627-639. 10.1016/j.cell.2008.03.030.
Bernstein BE, Kamal M, Lindblad-Toh K, Bekiranov S, Bailey DK, Huebert DJ, McMahon S, Karlsson EK, Kulbokas EJ, Gingeras TR, Schreiber SL, Lander ES: Genomic maps and comparative analysis of histone modifications in human and mouse. Cell. 2005, 120: 169-181. 10.1016/j.cell.2005.01.001.
Berletch JB, Yang F, Disteche CM: Escape from X inactivation in mice and humans. Genome Biol. 2010, 11: 213-10.1186/gb-2010-11-6-213.
Marques M, Laflamme L, Gervais AL, Gaudreau L: Reconciling the positive and negative roles of histone H2A.Z in gene transcription. Epigenetics. 2010, 5: 267-272. 10.4161/epi.5.4.11520.
Braun S, Madhani HD: Shaping the landscape: mechanistic consequences of ubiquitin modification of chromatin. EMBO Rep. 2012, 13: 619-630. 10.1038/embor.2012.78.
Fleming AB, Kao CF, Hillyer C, Pikaart M, Osley MA: H2B ubiquitylation plays a role in nucleosome dynamics during transcription elongation. Mol Cell. 2008, 31: 57-66. 10.1016/j.molcel.2008.04.025.
Pfaffl MW: A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001, 29: 45e-45. 10.1093/nar/29.9.e45.
Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta C(T)) Method. Methods. 2001, 25: 402-408. 10.1006/meth.2001.1262.
Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22: 4673-4680. 10.1093/nar/22.22.4673.
Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011, 28: 2731-2739. 10.1093/molbev/msr121.
Financial support was received from Fondo Nacional de Desarrollo Científico y Tecnológico (FONDECYT): 1120873; Universidad Andrés Bello (UNAB): DI-4511/R. MA also acknowledges to the Comisión Nacional de Investigación Científica y Tecnológica (CONICYT)/Fondo de Financiamiento de Centros de Investigación en Areas Prioritarias (FONDAP): 15110027. We are grateful to Eduardo Fuentes for constructive comments on the manuscript, and to P Rudolffi and P Carvajal for technical support.
The authors declare that they have no competing interests.
NS carried out the experiments and wrote the manuscript. GN carried out some ChIP assays. MR carried out the immunofluorescence experiment. * MR is currently working at Departamento Ciencias Fundantes, Universidad Central de Chile. MA directed the research and co-wrote the manuscript. AM assisted with project design and development, and drafted the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Figure S1: Carp fish exposed to a wide range of temperature and photoperiod conditions. Graph of temperature and photoperiod changes during a seasonal annual cycle. (TIFF 1 MB)
Additional file 2: Figure S2: Identification and analysis of H2A.Z subtypes isolated from carp. Multiple nucleotide and amino acid sequence alignment of H2A.Z subtypes isolated from carp are shown. Black boxes represent the precise sites of nucleotide differences between CDS from H2A.Z.2.1, H2A.Z.3.1, H2A.Z.3.2, and H2A.Z.7. The deduced amino acid sequences of each H2A.Z subtype are shown in colors. The additional seven amino acids found in carp H2A.Z.7 are denoted with a red rectangle. CDS, coding DNA sequences. (PDF 2 MB)
Additional file 3: Figure S3: Constitutive gene expression of JAK1 during seasonal acclimatization of carp. Transcriptional expression analysis of JAK1 in liver tissues from summer- and winter-acclimatized carp. The mRNA levels of JAK1 were calculated using the Pfaffl method . Three different fish for each season were used. (TIFF 2 MB)
Additional file 4: Table S1: List of oligonucleotides used for RT-PCR, RT-qPCR, and ChIP assays. ChIP, chromatin immunoprecipitation; RT-PCR, reverse transcription polymerase chain reaction; RT-qPCR, reverse transcription quantitative PCR. (DOCX 21 KB)
Additional file 5: Table S2: GenBank Accession Numbers of CDS used for the phylogenetic analysis. The list depicts histone H2A.Z in eukaryotic organisms used in the phylogenetic analysis. CDS, coding DNA sequences. (DOCX 86 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Simonet, N.G., Reyes, M., Nardocci, G. et al. Epigenetic regulation of the ribosomal cistron seasonally modulates enrichment of H2A.Z and H2A.Zub in response to different environmental inputs in carp (Cyprinus carpio). Epigenetics & Chromatin 6, 22 (2013). https://doi.org/10.1186/1756-8935-6-22
- Carp fish
- Epigenetic mechanisms
- H2A.Z subtypes