ALL Metrics
-
Views
-
Downloads
Get PDF
Get XML
Cite
Export
Track
Research Note

Analysis of distribution of chromatin marks across "divergence islands" in three-spined stickleback (Gasterosteus aculeatus)

[version 1; peer review: 2 approved with reservations]
PUBLISHED 20 Dec 2016
Author details Author details
OPEN PEER REVIEW
REVIEWER STATUS

Abstract

The three-spined stickleback (Gasterosteus aculeatus) is a well-known model organism for studying adaptations to water salinity. In this work, we investigate the dynamics of an epigenetic landscape of water salinity adaptation using three chromatin marks: H3K27ac, H3K4me1 and H3K4me3. The choice of marks was determined by the fact that some adaptive genomic loci are situated in gene-free regions, suggesting their regulatory role as enhancers. Histone modifications seem to be a promising mechanism that could regulate such regions. Difference between histone modifications in sea and freshwater - both in genes and intergenic enhancers - may contribute to epigenetic plasticity of stickleback adaptation. As a result of this study, we found differential chromatin peaks in "divergence islands" at enhancer elements and promoters of genes, which are responsible for stress adaptation and homeostasis. However, a full genome study analysis is required to fully understand mechanism of adaptation to water salinity.

Keywords

divergence islands, chromatin marks, Gasterosteus aculeatus

Introduction

The three-spined stickleback (Gasterosteus aculeatus) is a model organism that can be used to elucidate molecular mechanisms of adaption to various salinities. In this work, we analyze epigenetic signatures that can be specific for various salinities. We study three chromatin marks H3K27ac, H3K4me1 and H3K4me3, which were chosen as these marks label active enhancers and promoters1, in marine and fresh water sticklebacks in their natural habitats, as well as in foreign (change of salinity) environments. These marks were studied in 19 “divergence islands”2, short genomic regions, which are highly diverged between marine and fresh water species. In the present study, 17 divergence islands overlapped in protein-coding genes, which are relevant to fresh water adaptation. In addition, we report the changes of histone modifications between marine and freshwater fish in promoters of protein-coding genes and enhancers.

Methods

Sample preparation and ChIP experiment

In this study, we used six marine and six fresh water sticklebacks, which were collected in august of 2015 from the White Sea (near the Pertsov White Sea Biological Station of Lomonosov Moscow State University, Murmanskaya Oblast, Russia) and Mashinnoe Lake (located near the village of Chkalvosky, Repuyblic of Karelia, Russia), respectively. The fish were placed for four days in various water tanks: half of the sticklebacks from each group was kept in the water with their natural salinity (FF and MM for fresh and marine water, respectively) and the other half was kept in a modified environment (i.e. three fresh water fish was placed in salty water (FM) and three marine fish were placed in fresh water (MF)). The live fish were transferred to the laboratory to Moscow.

For the chromatin immunoprecipitation (ChIP)-seq experiments, gills were collected from all 12 sticklebacks. Chromatin was prepared from gills, as described by Cell Signalling (https://www.cellsignal.com/common/content/content.jsp?id=chip-agarose) and ChIP was performed as described by Filion et al.3.

Bioinformatics analysis

Reads were aligned to gasAcu1 reference genome from the UCSC Genome Browser Gateway (https://genome-euro.ucsc.edu/) by Bowtie2 (http://bowtie-bio.sourceforge.net/bowtie2/index.shtml) with a "very-sensitive-local" parameter4,5. For peaks calling we used MACS version 1.4.26. For the intersection of peaks with divergence islands and genes we used bedtools version 2.26.07.

Results

In our study, we determined between 7138 and 20828 histone modification peaks (Table 1). We selected histone modification peaks in “divergence islands” regions, which are highly divergent between marine and freshwater populations of sticklebacks2. We observed that the majority (17 out of 19) of the islands showed the same chromatin marks (H3K4me1, H3K4me3 and H3K27ac) in fresh water species in their natural salinity (FF) and fresh water species placed into a marine environment (FM). The same was true for marine water species in their natural salinity (MM) and marine water species placed into a fresh water environment (MF). In addition, the majority of the islands (14 out of 19) showed the same histone modifications between the whole set of marine and fresh water species (MM+MF vs. FF+FM; Table 2).

Table 1. Results of reads mapping and peaks calling.

MM - marine species in natural environment; MF - marine species in fresh water; FF - fresh water species in natural environment; FM - fresh water species in marine water.

SamFle NameAll readsMapped readsRatioPeaks
MMH3K27ac16,319,16213,242,76381.1485479511906
MMH3K4me115,805,73512,244,40177.4680899113289
MMH3K4me37,354,1666,872,44293.4496447312944
MFH3K27ac16,917,5839,700,33857.3387936113964
MFH3K4me112,773,84910,303,45580.6605354420828
MFH3K4me314,099,2897,836,36955.5798877518851
FFH3K27ac14,567,42511,367,24478.03193779130
FFH3K4me112,243,30510,737,65387.702242168825
FFH3K4me315,176,6979,476,32662.4399762314672
FMH3K27ac14,850,80011,887,59780.0468459612889
FMH3K4me113,213,51210,426,78678.910027867138
FMH3K4me314,512,63011,658,22580.3315801515723

Table 2. Results of intersection of chromatin peaks with divergence islands and differential chromatin peaks analysis between MM and MF, and between FF and FM.

“0” - no intersection; “1” – intersection. Regions with differential chromatin peaks are highlighted. MM - marine species in natural environment; MF - marine species in fresh water; FF - fresh water species in natural environment; FM - fresh water species in marine water. Differential peaks islands are highlighted in red.

Divergence islands coordinatesMMMFFFFM
ChromosomeStartEndH3K27acH3K4me1H3K4me3H3K27acH3K4me1H3K4me3H3K27acH3K4me1H3K4me3H3K27acH3K4me1H3K4me3
chrI21,487,99821,960,119111111111111
chrII14,874,36614,898,826000000000000
chrIV12,803,78012,881,296111111111111
chrIV13,930,00213,959,331010011011011
chrIV19,811,92219,914,666111111111111
chrIV23,954,63423,981,981101111101101
chrIV26,016,95526,166,536111111111111
chrV2,482,2092,501,295111111101101
chrVII17,982,35118,002,671111111111111
chrIX8,521,9358,537,559111111111101
chrIX8,901,8168,910,115111111001101
chrIX9,208,1589,227,809000000000000
chrIX1033410110353801111111001011
chrXI54457575855124111111111111
chrXII1433822914358336111111101101
chrXII1652202816538810000000000000
chrXIX24499032581858111111111111
chrXIX1478790414799088000000000000
chrXXI57598797486635111111111111

Nevertheless, we found 3 out of 19 “divergence islands” demonstrating differential chromatin marks in cases of short-term adaptation to water salinity for fish placed into foreign environments. Interestingly, one “divergence island” gained H3K4me1 in the promoter of the RPTOR gene in FM. The other two islands gained H3K27ac and lost H3K4me1 inf FM, suggesting their role as enhancers for genes outside the island. Also, two islands gained H3K4me3 and H3K4me1 at STC2 and PNPLA3 genes, respectively, in MF.

Finally, 5 out of 19 islands demonstrated differential histone modifications between fresh water species and marine species placed into fresh water (FF vs MF) (Table 3). In all these cases, MF gained a mark, which was absent in a FF. For example, H3K4me1 was gained in LRRC59 and BDH2 genes, which is involved in the adaption to stress and homeostasis8, suggesting that these genes might be activated after the placement of marine fish into a fresh water environment. In addition, 6 islands out of 19 demonstrated differential histone modifications between marine and fresh water species placed in marine water, with both gains and losses of marks (Table 3). Among these, we found that H3K4me3 was gained at STC2 in FM* compared to MM*, H3K4me1 at LRRC59 and BDH2 in MM* compared with FM*, and H3K27ac at RPTR in MM* compared with FM*.

Table 3. Results of intersection of chromatin peaks with divergence islands and differential chromatin peaks analysis between FM and MM, and between MF and FF.

“0” - no intersection; “1” – intersection. Regions with differential chromatin peaks are highlighted. MM - marine species in natural environment; MF - marine species in fresh water; FF - fresh water species in natural environment; FM - fresh water species in marine water. Differential peaks islands are highlighted in red.

Divergence islands coordinatesFMMMMFFF
ChromosomeStartEndH3K27acH3K4me1H3K4me3H3K27acH3K4me1H3K4me3H3K27acH3K4me1H3K4me3H3K27acH3K4me1H3K4me3
chrI21,487,99821,960,119111111111111
chrII14,874,36614,898,826000000000000
chrIV12,803,78012,881,296111111111111
chrIV13,930,00213,959,331011010011011
chrIV19,811,92219,914,666111111111111
chrIV23,954,63423,981,981101101111101
chrIV26,016,95526,166,536111111111111
chrV2,482,2092,501,295101111111101
chrVII17,982,35118,002,671111111111111
chrIX8,521,9358,537,559101111111111
chrIX8,901,8168,910,115101111111001
chrIX9,208,1589,227,809000000000000
chrIX1033410110353801011111111001
chrXI54457575855124111111111111
chrXII1433822914358336101111111101
chrXII1652202816538810000000000000
chrXIX24499032581858111111111111
chrXIX1478790414799088000000000000
chrXXI57598797486635111111111111

Conclusions

In this study, we analyzed the epigenetic profile of "divergence islands" with three chromatin marks, H3K4me1, H3K4me3 and H3K27ac. We report differential histone modifications that might be involved in the regulation of promoters and enhancers located in "divergent islands", and therefore contribute to adaptation to water salinity. Furthermore, we found differential chromatin peaks at promoters of genes that are responsible for stress adaptation and homeostasis. The results of this study contributes to our understanding of molecular mechanisms of adaptation to water salinity. However, a full genome histone modification analysis is required in order to further understand these mechanisms of adaptation.

Data availability

Fastq files can be found in the SRA archive under accession number SRX2403902 (https://www.ncbi.nlm.nih.gov/sra/?term=SRX2403902).

Comments on this article Comments (0)

Version 1
VERSION 1 PUBLISHED 20 Dec 2016
Comment
Author details Author details
Competing interests
Grant information
Copyright
Download
 
Export To
metrics
Views Downloads
F1000Research - -
PubMed Central
Data from PMC are received and updated monthly.
- -
Citations
CITE
how to cite this article
Sokolov A, Zhenilo S, Rastorguev S et al. Analysis of distribution of chromatin marks across "divergence islands" in three-spined stickleback (Gasterosteus aculeatus) [version 1; peer review: 2 approved with reservations]. F1000Research 2016, 5:2880 (https://doi.org/10.12688/f1000research.10428.1)
NOTE: If applicable, it is important to ensure the information in square brackets after the title is included in all citations of this article.
track
receive updates on this article
Track an article to receive email alerts on any updates to this article.

Open Peer Review

Current Reviewer Status: ?
Key to Reviewer Statuses VIEW
ApprovedThe paper is scientifically sound in its current form and only minor, if any, improvements are suggested
Approved with reservations A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit.
Not approvedFundamental flaws in the paper seriously undermine the findings and conclusions
Version 1
VERSION 1
PUBLISHED 20 Dec 2016
Views
11
Cite
Reviewer Report 20 Mar 2017
Ilkka Kronholm, Centre of Excellence in Biological Interactions, Department of Biological and Environmental Sciences, University of Jyväskylä, Jyväskylä, Finland 
Approved with Reservations
VIEWS 11
The study by Sokolov et al. investigates changes in chromatin marks among saline and freshwater populations of sticklebacks. While I like the idea of the study and I do appreciate that this manuscript is intended as a short note, the ... Continue reading
CITE
CITE
HOW TO CITE THIS REPORT
Kronholm I. Reviewer Report For: Analysis of distribution of chromatin marks across "divergence islands" in three-spined stickleback (Gasterosteus aculeatus) [version 1; peer review: 2 approved with reservations]. F1000Research 2016, 5:2880 (https://doi.org/10.5256/f1000research.11237.r20811)
NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article.
Views
20
Cite
Reviewer Report 03 Jan 2017
Alexey Ruzov, Division of Cancer and Stem Cells, Centre for Biomolecular Sciences (CBS), School of Medicine, University of Nottingham, Nottingham, UK 
Approved with Reservations
VIEWS 20
In my opinion, the results on differential distribution of histone modifications in sea and freshwater three-spined sticklebacks may potentially be interesting. Despite this, I strongly believe that text of the manuscript is currently beyond the minimal standards required for paper indexing.
... Continue reading
CITE
CITE
HOW TO CITE THIS REPORT
Ruzov A. Reviewer Report For: Analysis of distribution of chromatin marks across "divergence islands" in three-spined stickleback (Gasterosteus aculeatus) [version 1; peer review: 2 approved with reservations]. F1000Research 2016, 5:2880 (https://doi.org/10.5256/f1000research.11237.r18627)
NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article.

Comments on this article Comments (0)

Version 1
VERSION 1 PUBLISHED 20 Dec 2016
Comment
Alongside their report, reviewers assign a status to the article:
Approved - the paper is scientifically sound in its current form and only minor, if any, improvements are suggested
Approved with reservations - A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit.
Not approved - fundamental flaws in the paper seriously undermine the findings and conclusions
Sign In
If you've forgotten your password, please enter your email address below and we'll send you instructions on how to reset your password.

The email address should be the one you originally registered with F1000.

Email address not valid, please try again

You registered with F1000 via Google, so we cannot reset your password.

To sign in, please click here.

If you still need help with your Google account password, please click here.

You registered with F1000 via Facebook, so we cannot reset your password.

To sign in, please click here.

If you still need help with your Facebook account password, please click here.

Code not correct, please try again
Email us for further assistance.
Server error, please try again.