Rare genetic variation at transcription factor binding sites modulates local DNA methylation profiles
Authors:
Alejandro Martin-Trujillo aff001; Nihir Patel aff001; Felix Richter aff001; Bharati Jadhav aff001; Paras Garg aff001; Sarah U. Morton aff002; David M. McKean aff003; Steven R. DePalma aff004; Elizabeth Goldmuntz aff006; Dorota Gruber aff008; Richard Kim aff009; Jane W. Newburger aff010; George A. Porter, Jr. aff012; Alessandro Giardini aff013; Daniel Bernstein aff014; Martin Tristani-Firouzi aff015; Jonathan G. Seidman aff004; Christine E. Seidman aff003; Wendy K. Chung aff016; Bruce D. Gelb aff001; Andrew J. Sharp aff001
Authors place of work:
The Mindich Child Health and Development Institute and Department of Genetics & Genomic Sciences, Icahn School of Medicine at Mount Sinai, New York, New York, United States of America
aff001; Department of Newborn Medicine, Boston Children’s Hospital, Boston, Massachusetts, United States of America
aff002; Division of Cardiovascular Medicine, Brigham and Women’s Hospital, Boston, Massachusetts, United States of America
aff003; Department of Genetics, Harvard Medical School, Boston, Massachusetts, United States of America
aff004; Howard Hughes Medical Institute, Harvard University, Boston, Massachusetts, United States of America
aff005; Division of Cardiology, Children’s Hospital of Philadelphia, Philadelphia, PA, United States of America
aff006; Department of Pediatrics, University of Pennsylvania Perlman School of Medicine, Philadelphia, PA, United States of America
aff007; Department of Pediatrics, Cohen Children’s Medical Center, Northwell Health, New Hyde Park, NY, Unites States of America
aff008; Department of Pediatrics, Keck School of Medicine, University of Southern California, Los Angeles, California, United States of America
aff009; Department of Cardiology, Boston Children's Hospital, Boston, Massachusetts, United States of America
aff010; Department of Pediatrics, Harvard Medical School, Boston, Massachusetts, United States of America
aff011; Department of Pediatrics, University of Rochester Medical Center, Rochester, NY, United States of America
aff012; Cardiothoracic Unit, Great Ormond Street Hospital, London, England
aff013; Department of Pediatrics, Stanford University School of Medicine, Stanford, CA, United States of America
aff014; Department of Pediatrics, University of Utah School of Medicine, Salt Lake City, UT, United States of America
aff015; Departments of Pediatrics and Medicine, Columbia University, New York, NY, United States of America
aff016; Department of Pediatrics, Icahn School of Medicine at Mount Sinai, New York, United States of America
aff017; Graduate School of Biomedical Sciences, Icahn School of Medicine at Mount Sinai, New York, New York, United States of America
aff018
Published in the journal:
Rare genetic variation at transcription factor binding sites modulates local DNA methylation profiles. PLoS Genet 16(11): e1009189. doi:10.1371/journal.pgen.1009189
Category:
Research Article
doi:
https://doi.org/10.1371/journal.pgen.1009189
Summary
Although DNA methylation is the best characterized epigenetic mark, the mechanism by which it is targeted to specific regions in the genome remains unclear. Recent studies have revealed that local DNA methylation profiles might be dictated by cis-regulatory DNA sequences that mainly operate via DNA-binding factors. Consistent with this finding, we have recently shown that disruption of CTCF-binding sites by rare single nucleotide variants (SNVs) can underlie cis-linked DNA methylation changes in patients with congenital anomalies. These data raise the hypothesis that rare genetic variation at transcription factor binding sites (TFBSs) might contribute to local DNA methylation patterning.
In this work, by combining blood genome-wide DNA methylation profiles, whole genome sequencing-derived SNVs from 247 unrelated individuals along with 133 predicted TFBS motifs derived from ENCODE ChIP-Seq data, we observed an association between the disruption of binding sites for multiple TFs by rare SNVs and extreme DNA methylation values at both local and, to a lesser extent, distant CpGs. While the majority of these changes affected only single CpGs, 24% were associated with multiple outlier CpGs within ±1kb of the disrupted TFBS. Interestingly, disruption of functionally constrained sites within TF motifs lead to larger DNA methylation changes at nearby CpG sites. Altogether, these findings suggest that rare SNVs at TFBS negatively influence TF-DNA binding, which can lead to an altered local DNA methylation profile. Furthermore, subsequent integration of DNA methylation and RNA-Seq profiles from cardiac tissues enabled us to observe an association between rare SNV-directed DNA methylation and outlier expression of nearby genes.
In conclusion, our findings not only provide insights into the effect of rare genetic variation at TFBS on shaping local DNA methylation and its consequences on genome regulation, but also provide a rationale to incorporate DNA methylation data to interpret the functional role of rare variants.
Keywords:
DNA methylation – Gene disruption – Gene expression – Genetic polymorphism – Genetics – Genomics – Methylation – Sequence motif analysis
Introduction
Over the last decade, genomic DNA sequence variation has been associated with quantitative changes in multiple molecular phenotypes, including variation in gene expression and epigenetic marks such as chromatin accessibility and DNA methylation [1–5]. These variants are commonly referred to as quantitative trait loci (QTLs), the mapping of which has allowed us to unravel the primary mechanisms by which cis-regulatory variants can influence phenotypic variation. Single nucleotide variants (SNVs) present within cis-regulatory regions such as transcription factor binding sites (TFBSs) are enriched for different types of QTLs, such as those influencing chromatin structure/accessibility [1,2] and DNA methylation (meQTL) [3], suggesting these as functional regulatory variants. However, to date, most studies have focused on common variants, and the functional consequences of rare variants (minor allele frequency (MAF) ≤1%) have not been systematically interrogated.
Large population-scale projects using deep sequencing have provided extensive catalogues of human genomic variation, showing that rare genetic variants are abundant in the human genome [6–8]. It has been hypothesized that rare variants likely also contribute to the genetic architecture of complex diseases, which cannot be entirely explained by common genetic variation. Recent studies have demonstrated that some rare SNVs induce dramatic changes to the expression levels of nearby genes [9–11], which likely contribute to phenotypic variation. However, there is limited insight about how these variants can regulate these transcriptional changes. An important regulatory mechanism for controlling transcriptional activity is DNA methylation [12].
DNA methylation is an epigenetic modification crucial for mammalian development, playing a critical role in many cellular processes such as X-chromosome inactivation [13], genomic imprinting or maintenance of genome stability [12,14]. Although DNA methylation has been extensively characterized in mammals, the mechanisms by which it is targeted to specific genomic regions remains unclear. Recent data have shown that cis-acting methylation-determining regions (MDRs), short motifs that appear to encode epigenetic patterns in the local region, are sufficient for the de novo establishment of local DNA methylation patterns during embryonic development [15]. Interestingly, deleterious mutations at DNA binding motifs within these MDRs result in methylation changes at the adjacent CpGs. Furthermore, in line with this, we have recently shown that single base mutations at the canonical binding site for CCCTC-Binding Factor (CTCF) are enriched around cis-linked DNA methylation changes in both patients with congenital anomalies (CA) and in the normal population [16]. These findings suggest the potential for regulatory genetic variation to have profound effects on individual epigenetic profiles. However, it is still unknown whether this effect is limited to a few TFs or is a more general mechanism underlying the regulation of DNA methylation. Supporting the latter, experiments including targeted gene disruption of TFs in mouse stem cells have implicated DNA-binding factors in the regulation of local DNA methylation profiles [17]. Altogether, these data support an emerging role of DNA-binding TFs in regulating local DNA methylation profiles and suggest that this is likely influenced by genetic variation at their transcription factor binding sites (TFBSs).
Based on the potential role of cis-regulatory elements in modulating DNA methylation profiles locally [15,16], we hypothesized that the disruption of TFBSs by rare SNVs can modify local DNA methylation profiles and lead to altered transcription of nearby genes. To test this hypothesis, we combined blood DNA methylation profiles and SNVs derived from whole genome sequencing (WGS) of 247 unrelated individuals together with TFBS motifs for 133 different TFs predicted from chromatin immunoprecipitation sequencing (ChIP-Seq) data [18,19]. Subsequent integration of transcriptome and DNA methylation profiles from cardiac tissue enabled the study of the functional consequences of the TF-mediated DNA methylation at the transcriptional level.
Results
Rare variants at canonical TF motifs affect local DNA methylation profiles
The flow diagram presented in Fig 1A represents our approach to integrate, filter and analyze data to determine the influence of rare regulatory genetic variation on local DNA methylation profiles. Briefly, after identification of rare SNVs that lie within TFBSs, we extracted β-values corresponding to the Illumina Infinium MethylationEPIC BeadChip (EPIC) microarray probes located within the proximity of the SNV-disrupted TFBSs and ranked them from lowest (1) to highest (247) (see methods for further details). Following this approach, we expected to identify rare regulatory SNVs that result in extreme high or low DNA methylation values in individuals carrying the tested SNV as compared to controls, i.e, individuals who do not carry the tested variant (Fig 1B).
The intersection between our genotype and DNA methylation datasets yielded a total of 127,335 SNV:TFBS pairs including 91,356 rare SNVs (MAF≤1%) within 120,096 different TFBS motifs that were informative for DNA methylation, i.e., EPIC probes mapping within 1 kilobase (kb)- flanks of SNV-targeted TFBS (S1 Table). To detect the effects of rare genetic variation on local DNA methylation profiles, we focused on the 5% tails of the rank distribution that showed an absolute minimum β-value difference of 0.05 between SNV carriers and controls. We observed that carriers of rare SNVs that disrupt TFBSs (SNV-TFBSs) show an increased burden of extreme methylation values as compared to individuals without rare variants (Fig 1C), suggesting a potential effect of rare regulatory SNVs in modulating local DNA methylation profiles.
To determine the extent of the effect of rare genetic variation within TFBS on DNA methylation, we explored whether a single or multiple CpG sites show extreme methylation values in the proximity of the mutated TFBS. Here, we considered those TFBS with 3 or more CpG sites located within 1kb and, at least one of them showing extreme methylation (n = 16,708). Of these, we observed that 12,682 (76%) had only a single outlier CpG, while 4,026 (24%) had between 2 and 13 CpG sites within 1kb with outlier DNA methylation values. To formally assess whether loci with outlier DNA methylation could be identified using an analysis that searches for regions containing multiple outlier CpGs, termed differentially methylated regions (DMRs), we performed a prospective screen for DMRs in each sample using a sliding window algorithm, similar to that described in Barbosa et al. [16]. Using this approach, we identified a total of 9,689 different DMRs in our cohort (median of 202 DMRs per sample). Because our algorithm was only able to detect DMRs at loci containing ≥3 probes within a 1kb interval, we intersected our DMR list with the set of 16,708 SNV-TFBS described above, and observed 744 loci (4.5%) of SNV-TFBS that showed a DMR within 1kb of the mutated TFBS (S2 Table).
We used permutation testing to assess the association between the disruption of TFBSs by rare SNVs and extreme CpG methylation (see Methods), and observed that 46 of the 133 tested TFBSs were significantly associated with extreme methylation values levels (Fig 2) (S3 and S4 Tables), suggesting that the sequence-specific binding of TFs to these 46 TFBSs could be involved in shaping local DNA methylation patterns. Amongst these TFBSs, 36 correspond to the canonical motif for well-known human TFs such as Yin and Yang 1 (YY1), CTCF and Specificity Protein 1 (SP1), while the remaining 10 correspond to anonymous motifs for which the potential binding factor has not yet been characterized.
As the direction of DNA methylation change can have different potential functional consequences on genome regulation, we next assessed the direction of the effect of rare SNVs on local DNA methylation changes. We observed a 1.45-fold enrichment for methylation gains across all tested SNV-TFBS pairs. For example, both CTCF and SP1 showed an enrichment for gains of methylation, which is in accordance with the known protective role of these two factors against DNA methylation [20].
Minimal intersection between extreme DNA methylation values and DNA methylation quantitative trait loci
As variation in DNA methylation can also be driven by common genetic variation, we next assessed whether our results might be confounded by known meQTL or regions of haplotype-specific methylation (HSMs). We overlapped our results with CpGs where it is known that DNA methylation levels are associated with common genetic variation [21,22], and observed that only a very small fraction of CpG sites showing extreme DNA methylation and the disruption of TFBS by rare SNV (n = 189, or ~1% of the total) correspond to known meQTL or HSM regions. While these numbers do represent a significant enrichment for meQTL and HSM regions when compared to CpG sites that do not show outlier DNA methylation, these results indicate that the vast majority of our results are likely not driven by effects of common genetic variation.
Disruption of functionally constrained sites within TF binding motifs is associated with larger methylation changes
Since the binding affinity for many TFs relies on the recognition of specific DNA sequences [23], a possible direct consequence of DNA sequence variation at TFBSs could be alterations of the TF-DNA interaction. The binding affinity preferences of TFs for a specific DNA sequence can be represented by position weight matrices (PWM). These matrices denote the nucleotide frequencies at each position of a given DNA sequence motif, allowing an estimation of the effect of motif variation on TF binding affinity in a site-specific manner.
To infer the effect of rare SNV-TFBSs on TF-DNA binding, we computed separately for each SNV-TFBS the difference between the PWM score of the reference and alternate alleles (ΔPWM) [19] (http://hgdownload.soe.ucsc.edu/goldenPath/hg19/database/factorbookMotifPwm.txt.gz). Smaller values of ΔPWM correspond to SNVs at less constrained sites, which we hypothesized might be less likely to have an effect on TF binding efficiency (Fig 3A). Conversely, larger values of ΔPWM correspond to SNVs at more conserved positions within the motif that we hypothesized were more likely to alter the TF binding, and consequently impact DNA methylation (Fig 3B and 3C).
We analyzed 45 TFBSs, at which SNVs were associated with extreme DNA methylation changes and for which PWMs were available. Consistent with our hypothesis, we observed that rare SNVs disrupting those TFBSs associated with outlier methylation values showed higher ΔPWM scores compared to SNVs that are not associated with extreme methylation changes (Wilcoxon paired signed rank test, p = 0.0026) (Fig 3D). These results are consistent with rare variants altering local DNA methylation through altered binding of TFs to their motifs.
Rare SNVs within TFBS exert stronger effects on local DNA methylation than SNVs within flanking regions
To confirm the key role of TFs in modulating local DNA methylation profiles, we compared the effect on local DNA methylation between SNV-TFBSs and SNVs that lie adjacent to but outside of annotated TFBS (SNV-noTFBS). Following the same approach as above, we assessed the effect on local DNA methylation for 571,137 different rare SNV-noTFBSs.
Analogous to our previous analysis of rare SNV-TFBS, we observed a genome-wide enrichment for extreme methylation values for rare SNV-noTFBS, suggesting that local DNA methylation changes can also be driven by variants outside of annotated TFBSs. However, we observed that the effect of SNV-TFBSs on DNA methylation was 1.18 fold greater when compared to SNV-noTFBS (Paired t test, P-value<0.0001, Fig 4), indicating that rare genetic variants at TFBSs are enriched for effects on local epigenetics.
Rare genetic variation within TFBS can influence DNA methylation over larger distances in cis
So far, we have shown that the disruption of TFBSs by rare SNVs is associated with altered DNA methylation profiles at the local level. However, we hypothesized that the effect of genetic variation on DNA methylation might also occur across wider regions [4,5]. To address this, we followed a similar approach to that described above, but selected β-values corresponding to CpGs located between 1–100 kb around each disrupted TFBS. We first binned β-values according to their separation from the disrupted TFBS, ranked and categorized them to define those showing extreme methylation, and then performed permutation testing for each TFBS per bin. After applying a multiple testing correction, we observed many significant associations between extreme DNA methylation and different TFBSs. However, there was a clear effect of physical distance, with the number of TFBS associations reducing markedly once separation of TFBS and CpG was >30 kb (Fig 5), suggesting that the majority of TF binding events modulate DNA methylation over relatively short distances (S5 Table).
Disruption of DNA methylation by rare genetic variation at TFBS alters expression levels of nearby genes
Next, based on the key role of DNA methylation in the control of transcription [12], we explored the impact of DNA methylation changes associated with the disruption of TFBSs on the expression of nearby genes. As both DNA methylation and gene expression profiles are often tissue-specific, we used available cardiac DNA methylation and RNA sequencing (RNA-seq) expression profiles from 20 unrelated individuals for whom WGS data were also available.
After quality control (QC) procedures and filtering (see Methods), we identified 12,903 rare SNVs that overlapped with 17,620 TFBS motifs. From these, we selected 4,944 rare SNVs that disrupted TFBSs within gene promoter regions (TSS±2kb), intersected these with normalized gene expression data, allowing evaluation of the impact of 3,978 rare regulatory SNVs on 3,478 independent transcripts. In comparison to gene promoters with a normal DNA methylation profile, this analysis revealed a 2.58-fold enrichment for extreme expression levels (>2 Z-scores from the mean) associated with genes with altered promoter methylation and a rare SNV in a nearby TFBS (p = 0.0363, two-tailed Fisher’s exact test, Fig 6A) (S6 Table). An example is shown in Fig 6B and 6C, where an individual with a rare variant within the canonical TFBS for UA4 (predicted to be THAP1) shows increased methylation levels at the promoter of the nearby GUSBP11 gene (Fig 6B). Consistent with the known repressive effects of DNA methylation at gene promoters, this is associated with reduced expression of GUSBP11 (Fig 6C). This result suggests that a fraction of rare genetic variation at TFBS contribute to local transcriptional regulation by modulating DNA methylation profiles.
Discussion
Here, we evaluated the association of rare regulatory genetic variation on DNA methylation profiles by focusing on extreme methylation values in a cohort of 247 individuals for whom both WGS and methylation profiles were available. Our results indicate that a fraction of rare SNV-TFBSs are associated with outlier DNA methylation profiles, which we suggest is potentially due to alterations in TF binding, and these events are associated with altered expression of nearby genes. Furthermore, our results highlight the value of incorporating DNA methylation data to interpret the functional consequences of rare regulatory genetic variation.
Analogous to our results, it has been recently shown that rare genetic variation is sometimes associated with extreme expression levels of nearby genes [9]. Interestingly, these regions are strongly enriched for variants at CpG-rich promoter regions and TFBSs. Based on our results, we hypothesize that a certain fraction of this altered gene expression might be due to changes in DNA methylation profiles as a consequence of rare genetic variation, especially involving SNVs within TFBSs.
We observed that a total of 46 canonical TFBSs are significantly associated with local DNA methylation patterns. Consistent with previous reports showing that genetic variation at CTCF binding sites have been involved in modulating, creating and propagating DNA methylation profiles in mammals [16,17,21], we observed a significant association between local DNA methylation profiles and the presence of rare SNVs in CTCF-binding sites. In addition to CTCF, disruption of the binding sites of ten other TFs in our list of 46 has recently been shown to cause alterations of local allelic DNA methylation [24]. In line with our results, Onuchic et al. showed that the resulting DNA methylation patterns are more often linked to rare, rather than common, genetic variation, reinforcing our approach of focusing on rare genetic variation. Altogether, these studies not only validate our approach to detect rare SNV-TFBS-directed DNA methylation but also suggest that the modulation of local DNA methylation by TFs represents a general mechanism involved in the regulation of genome function, rather than being limited to a few TFs.
In addition to local effects on DNA methylation, we also observed that, in some instances, rare genetic variation at TFBSs is also associated with changes in DNA methylation up to distances of 100 kb. This is in agreement with the potential role of some TFs in DNA loop formation [25–27], which can bring regions located hundreds of kilobases away into contact. For instance, we observed that CTCF, a structural protein involved in the establishment and maintenance of the 3D conformation of the genome [27], is consistently associated with extreme methylation values up to 50 kb from the disrupted TFBS. Similarly to CTCF, we also observed that AP-1 has effects on DNA methylation over large distances; notably, AP-1 has also been implicated in the formation of DNA loops during development [28].
According to our results, TFs seem to play a key role in preferentially dictating local DNA methylation profiles. Nevertheless, we still observed an enrichment for extreme methylation values associated with rare genetic variation outside of TFBSs. An explanation for this enrichment could be that, despite being outside of these regulatory elements, SNVs can still have an impact on TF activity, as recently shown by both computational and experimental approaches [29–32]. For instance, mutations at positions immediately adjacent to GATA1 binding motifs result in reduced binding of this TF [31], indicating that TF binding is not only regulated by genetic variation within the TF motif. Another possible explanation could be that these variants might lie within unannotated binding motifs for TFs that were not included in this study. For example, it is estimated that approximately one third of human TFs lack DNA-binding affinity models [33]. Furthermore, we cannot rule out that local DNA methylation profiles could be modulated by additional mechanisms.
The mechanisms by which TF binding mediates local DNA methylation are still unclear. Since TFs lack enzymatic activity, it is unlikely that they directly catalyze the gain or loss of DNA methylation at CpG sites. One possible mechanism would involve the formation of complexes with DNA methylation machinery, such as DNA methyltransferases (DNMTs) and ten-eleven translocation (TET) proteins [34–36]. For example, it has recently been shown that nuclear receptor subfamily 6 group A member 1 (Nr6a1) can interact with Dnmt3a and Dnmt3b proteins and induce DNA methylation at the Oct4 gene promoter that contains Nr6a1 binding sites [36,37]. Conversely, local DNA demethylation could result from the interaction between TFs and TET proteins in mammalian cells [34,35]. However, mutations at CTCF binding sites result in abnormal methylation of the maternal unmethylated allele of the H19 Imprinted Control Region in mice [38], suggesting that the binding of the TF on the DNA itself can protect against DNA methylation, likely by preventing accessibility to DNMTs at these sites. Altogether these data provide a possible mechanistic link between the alteration of TFBS by rare SNVs and the resulting DNA methylation profile.
While previous studies have shown that epigenetic changes due to regulatory mutations often lead to regional changes affecting clusters of multiple CpGs, it should be noted that for the majority of SNV-TFBS in our study (74%), we identified only a single associated outlier CpG. There are several possible explanations for this observation: (i) the EPIC array that we used to interrogate DNA methylation only has relatively sparse genomic coverage, sampling only a small subset of CpGs at most genomic loci; (ii) due to the large number of SNV-TFBS we interrogated, some of the pairwise associations with outlier methylation occurred by chance, and (iii) differential methylation of single CpGs might sometimes result from inaccurate measurements due to poorly performing probes or other technical artefacts, and does not represent true epigenetic variation [39]. To minimize this latter possibility, we performed pre-processing of the methylation dataset including the removal of potentially confounded probes and normalization. Despite this, we observed a 1.7-fold enrichment (p<0.0001, two-tailed Fisher’s Exact test) for multi-mapping probes reported by Pidsley et al. [40] in our set of CpGs showing extreme DNA methylation. However, it is important to note that these CpGs only account for a very small fraction of the overall set of outlier CpGs we reported (approximately 500, or <3% of outlier CpGs). Therefore, although we cannot rule out that a proportion of our results might be driven by technical artefacts, we believe that this fraction is likely very small and does not affect the overall conclusions of this study.
Conclusions
Our data suggest that a meaningful fraction of rare genetic variation at TFBSs can play a role in shaping DNA methylation profiles in cis, which in turn can result in altered expression of the nearby genes. Furthermore, we provide a rationale for integrating DNA methylation data to identify genuine functional genetic variation from the broader genetic background, which represents a current challenge in human genetics.
Materials and methods
Sample description
A total of 249 unrelated individuals with paired data for WGS and EPIC DNA methylation profiling of peripheral blood collected at the time of enrollment were selected from the cohort collected via Pediatric Cardiac Genomic Consortium (PCGC) [41]. An extensive description of PCGC cases as well as further details about sample collection can be found in a summary publications released by the PCGC [41,42]. Briefly, our cohort comprises 249 individuals, aged from newborn to 47 years (mean 8.2 years) and diagnosed with a range of congenital heart defects, with conotruncal and left-sided obstructive lesions being the two most common diagnoses (S7 Table).
Ethics statement
This research was approved by the Institutional Review Board of the Icahn School of Medicine at Mount Sinai under IRB #17–01980. All samples were collected after obtaining written informed consent from each participant or their parents for broad genomic studies approved by the Institutional Review Boards of Boston’s Children’s Hospital, Brigham and Women’s Hospital, Children’s Hospital of Los Angeles, Children’s Hospital of Philadelphia, Columbia University Medical Center, Great Ormond Street Hospital, Icahn School of Medicine at Mount Sinai, Rochester School of Medicine and Dentistry, Steven and Alexandra Cohen Children’s Medical Center of New York and Yale School of Medicine. Further details of ethical issues regarding PCGC samples, including patient privacy, data storage, and return of results to participants, can be found in a summary publication released by the PCGC [41].
Illumina sequencing, quality control and variant filtering
After TruSeq DNA PCR-free (Illumina Inc., San Diego, CA, USA) library preparation, WGS was performed on genomic DNA isolated from peripheral blood at Baylor College of Medicine (Houston, TX, USA) using Illumina HiSeq instrument to obtain an average of 36-fold genome coverage (range from 25- to 39-fold). Paired-end reads were aligned to human reference genome (GRCh37/hg19) using BWA-mem [43]. SNVs were called from aligned reads using GATK (v2.7) as previously described [44–46].
After applying filters for variant quality (GATK VQSR, ≥10 total reads, genotype quality score ≥30 and alternate allele fraction between 0.2–0.8 and ≥0.9 for heterozygous and homozygous SNPs, respectively), we retained for further analysis those SNVs with MAF ≤1% in the 1000 Genomes Project (1000G) [6] and Genome Aggregation Database (GnomAD) [7] (http://gnomad.broadinstitute.org/about) databases, as well as in our cohort (n = 247).
Illumina Infinium MethylationEPIC BeadChip array profiling
Genome-wide DNA methylation profiling was performed on genomic DNA isolated from peripheral blood using the EPIC array (Illumina, Inc., San Diego, CA, USA) according to Infinium HD Methylation Assay Protocol (IIlumina). Methylation levels are given as β-values, representing the ratio of intensities between the methylated and unmethylated signals. Resulting β-values range from 0 for completely unmethylated to 1 for completely methylated.
For every sample, raw data files with β-values, intensity values per channel, and detection p-values per probe were obtained from the New York Genome Center. Before processing, QC of raw data was performed comprising a gender check comparing X and Y-chromosome data against reported sample gender, principal component analysis (PCA) plots and density plots of M values. Based on these steps, two samples were excluded from downstream analysis, as they were clear outliers on PCA plots. To prevent technical biases that could influence methylation measurements, we initially excluded EPIC probes with (i) internal common SNPs (MAF ≥5% in 1000G database) within the last 5 bases of the 3΄ end of the probe or (ii) non-unique mapping to the bisulfite-converted genome. Due to gender differences on the sex chromosomes, only β-values of autosomal probes were utilized in downstream processing. Raw β-values for 822,016 EPIC probes for the remaining samples (n = 247) were background and color corrected and quantile normalized using lumi and methylumi R packages [47]. Finally, Infinium probe design biases were corrected using BMIQ [48].
Since DNA methylation patterns are often tissue-specific and blood consists of multiple different cell-types [12], resulting β-values can simply reflect differences in cell-type composition among samples. To prevent this potential confounder, we first estimated the fraction of six blood cell types (CD4+, CD8+, B and natural killer cells, granulocytes and monocytes) in our samples applying Houseman algorithm [49] and correlated these estimates with β-values of each probe. We observed that the probes most significantly correlated with cell-type fraction (Bonferroni corrected P-value for number of probes <0.05) showed larger methylation differences, suggesting that heterogeneity in white blood cells across our samples had the potential to confound our DNA methylation measurements (S1 Fig, right panel). Subsequently, this bias was removed by excluding most highly correlated probes with every cell fraction up to 5% of the total number of processed probes on the EPIC array (S1 Fig, left panel). Following this approach, a total of 76,776 potentially cell-type confounded EPIC probes were excluded from our analysis. Furthermore, we also excluded EPIC probes with rare SNV-TFBSs present within the probe, at the interrogated CpG or its adjacent base.
Enrichment of extreme DNA methylation values due to the disruption of TFBS by rare SNVs
For this analysis, only autosomal SNVs derived from WGS data from 247 unrelated with minor allele frequency (MAF) ≤1% in GnomAD and 1000G databases as well as in our cohort (n = 247) (Number of rare SNVs per individual, mean = 84,245, max = 320,786, min = 47,391) were included. Furthermore, to capture rare events, SNVs that were present in >2 individuals were also excluded. Based on previous observations about the potential role of TFs in shaping DNA methylation profiles and how this process can be affected by genetic variation at their binding sites [15–17], we focused on rare SNVs overlapping TFBSs. We used a catalogue of 2,298,872 motifs for 133 human TFBSs identified from ChIP-Seq peaks generated by the ENCODE project [18,19,50]. Finally, to avoid redundancy in our analysis, we merged overlapping TFBSs for the same TF that share the same SNV, respectively. This resulted in a total of 131,357 SNVs that overlap 188,226 TFBSs.
After performing standard QC and processing (See above), we extracted β-values corresponding to EPIC probes located within the proximity of the SNV-disrupted TFBS, i.e., the SNV-disrupted TFBSs and their ±1-kb-extended flanks, and ranked these from lowest (1) to highest (247). As a result of intersecting our genotype and DNA methylation datasets, we identified a total of 91,356 rare SNV across 120,096 different binding motifs that were informative for DNA methylation, i.e., EPIC probes mapping within 1-kb-extended flanks of SNV-targeted TFBS (217,578 different EPIC probes). Bedtools (v2.27) was used to intersect rare SNVs derived from WGS with 2,298,872 motifs instances corresponding to 133 canonical TFBSs, obtained from the UCSC Genome Browser track “Transcription Factor ChIP-seq (161 factors) from ENCODE with Factorbook Motifs”. These motifs (size 8–24 bp) have been predicted using data obtained from ChIP-Seq experiments on 119 human TFs in multiple cell lines generated by the ENCODE Consortium [19]. (http://hgdownload.soe.ucsc.edu/goldenPath/hg19/database/factorbookMotifPos.txt.gz, March 16th 2014 release) [18,19].
To avoid redundancy in our analysis, overlapping motif instances of the same TFBS with a single SNV were merged into a single motif using bedtools groupby tool (v2.27). CpGs of interest, i.e. those within 1 kb flanking each TFBS, were identified using bedtools (v2.27). Next, the selected CpGs were ranked by increasing β-value and annotated for our extreme CpG methylation criteria, defined as within the 5% tails of rank distribution and a minimum absolute difference of 0.05 in β-value between SNV-TFBS carriers and non-carriers. Individuals were also annotated according to the presence (SNV carrier) or absence (SNV non-carrier) of a given SNV-TFBS. Enrichment for extreme methylation values was calculated as the ratio of the fraction of extreme CpG methylation values in the group of SNV-carriers over the fraction of extreme methylation values in the group of SNV-non-carriers. In order to assess statistically the ability of SNV-TFBSs to modulate local DNA methylation profiles, we performed 10,000 permutations where we randomized our data and re-calculated the ratio of extreme methylation values between the groups of SNV carriers and non-carriers. P-values were calculated as the fraction of permutations where this ratio was equal to or exceeded the observed data. We then applied a Bonferroni correction for the number of different TFBS tested (n = 133).
A similar approach was carried out to determine the effect of the disruption of TFBSs by rare SNVs on DNA methylation over distances >1 kb. However, here we first binned CpGs according to separation from the TFBS, and permutation testing was performed separately for each TFBS and bin, as described above. For each bin, resulting p-values from the permutation test were Bonferroni corrected for the number of different tested TFBSs (n = 133).
Identification of differentially methylated regions
We screened for DMRs in each sample using a sliding window algorithm, similar to that described in [16]. Briefly, for each individual sample, this algorithm searches for genomic intervals where, within a 1 Kb window:
at least 3 probes each have β values above the 95th percentile of the cohort distribution for that probe, and are ≥0.1 above the cohort mean (Hypermethylation).
at least 3 probes each have β values below the 5th percentile of the cohort distribution for that probe, and are ≤0.1 below the cohort mean (Hypomethylation).
Because our algorithm is only able to detect DMRs at genomic intervals where ≥3 probes within 1 kb are assayed by the EPIC array, we intersected our DMR list with the set of SNV-TFBS that had at least one sample with an overlapping mutation associated with outlier methylation, and which had 3 or more CpGs located within 1kb (n = 16,708).
Enrichment analysis for meQTL and HSMs
As the meQTL reported by Do et al. [21] are based on CpG sites present on the Infinium Human Methylation 450K BeadChip (450K), we first selected 145,041 CpG probes that are present within the proximity of mutated TFBS (±1Kb either side) and are present on both the 450K and EPIC array platforms. After categorizing these according to their presence in the list of meQTL provided by Do et al. [21] and their methylation profile (extreme and non-extreme), we estimated enrichment of these probes for meQTL according to their methylation profile.
Similarly, to assess the potential impact of HSMs on our results, we first categorized the set of 217,578 EPIC probes included in our study according to their presence within HSMs and their methylation profile (extreme and non-extreme). Subsequently, we estimated the enrichment of them for these genomic intervals according to their methylation profile. Statistical significance was tested by two-tailed Fisher’s Exact test.
Positional impact of genetic variants within TF motifs
The potential effect of each base substitution on predicted DNA-binding was estimated by computing the absolute difference in Position Weight Matrix (ΔPWM) scores of the two alleles at the variant site using available PWMs (http://hgdownload.soe.ucsc.edu/goldenPath/hg19/database/factorbookMotifPwm.txt.gz) (March 16th 2014 release). This file includes PWMs for 131 different canonical TFBSs generated from ChIP-Seq data generated by the ENCODE project [19].
We analyzed 45 TFBSs, at which SNVs were associated with extreme DNA methylation changes and for which PWMs were available. Prior to computing ΔPWM scores, for each tested TFBS we categorized rare SNV-TFBSs into two groups according to their association to extreme methylation values. One group included those SNV-TFBSs associated with local extreme methylation values and another consisting of the remaining SNV-TFBSs. In cases where multiple EPIC probes overlap the same 1-kb-extended TFBS interval, we considered rare SNV-TFBSs as associated with local extreme methylation when at least 1/3 of the CpGs fulfilled our criteria for extreme CpG methylation. Furthermore, to avoid ambiguity in our PWM-based approach, we only selected those variants within isolated TFBSs, i.e., that has not been previously merged with any overlapping TFBS. After this selection, a total of 62,400 (66.1%) rare SNV-TFBSs were included in this analysis. We next calculated the ΔPWM score between the alleles for the selected SNV and compared the median of the ΔPWM score between the two groups. For each TFBS, the median of the absolute ΔPWM scores was calculated for each SNV-TFBS group and statistically compared using a Wilcoxon matched-pairs signed rank test.
Functional constraint and effect size on DNA methylation
Each CpG was annotated with the absolute ΔPWM score based on the cis-linked SNV-TFBS. CpGs were binned according to their association with extreme CpG methylation values, and the magnitude of the methylation change between SNV-carriers and non-carriers. Pearson correlation between the median values of absolute ΔPWM scores and magnitude of methylation change was computed across each bin.
Rare TFBS vs non-TFBS rare SNVs effect on DNA methylation at nearby CpG sites
Rare SNVs that were present within 1-kb flanks of TFBS but did not overlap an annotated TFBS were identified by bedtools (v2.27). After ranking CpGs based on their β-values and annotating for extreme CpG methylation, we calculated enrichment for extreme methylation in individuals carrying SNVs outside TFBSs. A paired t-test was used for matched chromosomal enrichment for extreme CpG methylation values between SNV–TFBS and SNV outside.
Functional consequences of rare genetic variation at TFBS in heart tissues
Selection of rare SNVs
Similar to our previous analysis in blood, only rare SNVs derived from WGS with MAF ≤1% in the 1000G [6] and GnomAD [7] databases as well as in the PCGC cohort were retained for further analysis. These SNVs were intersected with motifs of 133 canonical TFBSs using bedtools (v2.27). Furthermore, overlapping motif instances for the same TFBS that are targeted for the same SNV were merged using bedtools groupby tool (v2.27).
RNA-Sequencing
Twenty four transcriptome profiles derived from RNA-Seq of ventricle tissue were collected from Richter F. et al.[51]. Briefly, RNA extracted from cardiac tissues that were either snap-frozen or conserved in RNA-later from 327 participants who underwent cardiac surgery were sequenced on the Illumina Hi-Seq X Ten or NextSeq instrument using 50-bp paired-end sequencing. After sequence reads were aligned to the hg19 reference genome using Subread [52], read counts were calculated per gene using featureCounts [53]. Genes showing a mean of ≥1 read per kb per million mapped reads (RPKM) across the samples were considered as expressed in cardiac tissue. Genes below this threshold were excluded from further analysis. For the remaining genes, expression levels were corrected for known covariates (tissue, library preparation, sequencing platform, tissue storage, age, and gender) and subsequently, converted to Z-scores. Expression levels >2 Z-scores from the mean were considered as outliers.
Although PCA plots for gene expression profiles did not show high variance across our samples, four samples were excluded from this analysis as they display a significantly higher number of genes showing outlier expression as compared to the remaining samples (S2 Fig).
Cardiac methylation profiling
Genome-wide DNA methylation profiling was performed on genomic DNA isolated from ventricular tissues obtained from 25 individuals who underwent cardiac surgery. One sample was excluded from this analysis as there was a mismatch between array-inferred gender utilizing data from the sex chromosomes and the reported gender. Raw data files with β-values, intensity values per channel, and detection p-values per probe were processed as previously described for blood DNA methylation profiling, without controlling for leukocyte distribution.
A total of 20 individuals with available WGS data, cardiac RNA sequencing (RNA-Seq), and EPIC DNA methylation profiles were selected. Processed β-values for CpGs present within the targeted TFBS by SNV or their 1-kb flanks were selected for downstream analysis. Unlike our previous analysis in blood samples, due to the smaller size of our cohort we considered extreme CpG methylation values as those at the top (1) or bottom (20) of the rank distribution and showing an absolute minimum β-value difference of 0.05 between SNV carriers and non-carriers. Subsequently, SNV-TFBSs were categorized according to the extreme CpG methylation content within its flank. In cases where there were multiple CpGs within 1 kb of a TFBS, we considered this SNV-TFBS associated with extreme methylation values when at least one third of the CpGs fulfilled our criteria for extreme methylation.
After DNA methylation and gene expression outlier annotation, SNV-TFBSs were annotated with the closest RefSeq gene using bedtools (v2.27). Only those SNV-TFBSs that fell within gene promoter (TSS±2kb) of the closest gene were selected. Expression and DNA methylation data were intersected using the dplyr R package and, subsequently, an enrichment analysis for expression outliers between SNV-TFBSs associated with extreme methylation and SNV-TFBSs that were not associated with extreme methylation was performed. Statistical significance was tested by two-tailed Fisher’s exact test.
Supporting information
S1 Fig [eps]
Blood cell-type specific composition influences DNA methylation.
S2 Fig [a]
Principal component analysis (PCA).
S1 Table [xlsx]
TFBS included in association study.
S2 Table [xlsx]
List of differentially methylated regions due to the disruption of TFBS by rare SNVs.
S3 Table [xlsx]
Permutation test for individual TFBSs to test local effect of the disruption of TFBSs by rare SNVs on DNA methylation.
S4 Table [xlsx]
List of EPIC probes associated with extreme DNA methylation upon disruption of TFBSs by rare SNVs.
S5 Table [xlsx]
Wider effect on DNA methylation profile of rare SNVs that disrupt TFBSs.
S6 Table [xlsx]
Genes showing outlier expression associated to local extreme DNA methylation profile due to the disruption of TFBSs by rare SNVs in ventricular tissues.
S7 Table [xlsx]
Phenotype of individuals included in the study.
Zdroje
1. McVicker G, van de Geijn B, Degner JF, Cain CE, Banovich NE, Raj A, et al. Identification of Genetic Variants That Affect Histone Modifications in Human Cells. Science (80-). 2013;342: 747–749. doi: 10.1126/science.1242429 24136359
2. Degner JF, Pai AA, Pique-Regi R, Veyrieras JB, Gaffney DJ, Pickrell JK, et al. DNase-I sensitivity QTLs are a major determinant of human expression variation. Nature. 2012;482: 390–394. doi: 10.1038/nature10808 22307276
3. Banovich NE, Lan X, McVicker G, van de Geijn B, Degner JF, Blischak JD, et al. Methylation QTLs Are Associated with Coordinated Changes in Transcription Factor Binding, Histone Modifications, and Gene Expression Levels. PLoS Genet. 2014;10: 1–12. doi: 10.1371/journal.pgen.1004663 25233095
4. Gibbs JR, van der Brug MP, Hernandez DG, Traynor BJ, Nalls MA, Lai S-L, et al. Abundant quantitative trait loci exist for DNA methylation and gene expression in human brain. PLoS Genet. 2010;6: e1000952. doi: 10.1371/journal.pgen.1000952 20485568
5. Drong AW, Nicholson G, Hedman AK, Meduri E, Grundberg E, Small KS, et al. The presence of methylation quantitative trait loci indicates a direct genetic influence on the level of DNA methylation in adipose tissue. PLoS One. 2013;8: e55923. doi: 10.1371/journal.pone.0055923 23431366
6. McVean GA, Altshuler (Co-Chair) DM, Durbin (Co-Chair) RM, Abecasis GR, Bentley DR, Chakravarti A, et al. An integrated map of genetic variation from 1,092 human genomes. Nature. 2012;491: 56–65. doi: 10.1038/nature11632 23128226
7. Lek M, Karczewski KJ, Minikel E V., Samocha KE, Banks E, Fennell T, et al. Analysis of protein-coding genetic variation in 60,706 humans. Nature. 2016;536: 285–291. doi: 10.1038/nature19057 27535533
8. Ardlie KG, Deluca DS, Segre A V., Sullivan TJ, Young TR, Gelfand ET, et al. The Genotype-Tissue Expression (GTEx) pilot analysis: Multitissue gene regulation in humans. Science (80-). 2015;348: 648–660. doi: 10.1126/science.1262110 25954001
9. Li X, Kim Y, Tsang EK, Davis JR, Damani FN, Chiang C, et al. The impact of rare variation on gene expression across tissues. Nature. 2017;550: 239–243. doi: 10.1038/nature24267 29022581
10. Montgomery SB, Lappalainen T, Gutierrez-Arcelus M, Dermitzakis ET. Rare and common regulatory variation in population-scale sequenced human genomes. PLoS Genet. 2011;7: e1002144. doi: 10.1371/journal.pgen.1002144 21811411
11. Zhao J, Akinsanmi I, Arafat D, Cradick TJ, Lee CM, Banskota S, et al. A Burden of Rare Variants Associated with Extremes of Gene Expression in Human Peripheral Blood. Am J Hum Genet. 2016;98: 299–309. doi: 10.1016/j.ajhg.2015.12.023 26849112
12. Robertson KD. DNA methylation and human disease. Nat Rev Genet. 2005;6: 597–610. doi: 10.1038/nrg1655 16136652
13. Beard C, Li E, Jaenisch R. Loss of methylation activates Xist in somatic but not in embryonic cells. Genes Dev. 1995;9: 2325–34. Available from: http://www.ncbi.nlm.nih.gov/pubmed/7557385. doi: 10.1101/gad.9.19.2325 7557385
14. Li E, Beard C, Jaenisch R. Role for DNA methylation in genomic imprinting. Nature. 1993;366: 362–365. doi: 10.1038/366362a0 8247133
15. Lienert F, Wirbelauer C, Som I, Dean A, Mohn F, Schübeler D. Identification of genetic elements that autonomously determine DNA methylation states. Nature Genetics. 2011. pp. 1091–1097. doi: 10.1038/ng.946 21964573
16. Barbosa M, Joshi RS, Garg P, Martin-Trujillo A, Patel N, Jadhav B, et al. Identification of rare de novo epigenetic variations in congenital disorders. Nat Commun. 2018;9: 2064. doi: 10.1038/s41467-018-04540-x 29802345
17. Stadler MB, Murr R, Burger L, Ivanek R, Lienert F, Schöler A, et al. DNA-binding factors shape the mouse methylome at distal regulatory regions. Nature. 2011;480: 490–495. doi: 10.1038/nature10716 22170606
18. Wang J, Zhuang J, Iyer S, Lin XY, Greven MC, Kim BH, et al. Factorbook.org: A Wiki-based database for transcription factor-binding data generated by the ENCODE consortium. Nucleic Acids Res. 2013;41: 171–176. doi: 10.1093/nar/gks1221 23203885
19. Dunham I, Kundaje A, Aldred SF, Collins PJ, Davis C a., Doyle F, et al. An integrated encyclopedia of DNA elements in the human genome. Nature. 2012;489: 57–74. doi: 10.1038/nature11247 22955616
20. Blattler A, Farnham PJ. Cross-talk between site-specific transcription factors and DNA methylation states. J Biol Chem. 2013;288: 34287–94. doi: 10.1074/jbc.R113.512517 24151070
21. Do C, Lang CF, Lin J, Darbary H, Krupska I, Gaba A, et al. Mechanisms and Disease Associations of Haplotype-Dependent Allele-Specific DNA Methylation. Am J Hum Genet. 2016;98: 934–955. doi: 10.1016/j.ajhg.2016.03.027 27153397
22. Bell CG, Gao F, Yuan W, Roos L, Acton RJ, Xia Y, et al. Obligatory and facilitative allelic variation in the DNA methylome within common disease-associated loci. Nat Commun. 2018;9: 8. doi: 10.1038/s41467-017-01586-1 29295990
23. Lambert SA, Jolma A, Campitelli LF, Das PK, Yin Y, Albu M, et al. The Human Transcription Factors. Cell. 2018;172: 650–665. doi: 10.1016/j.cell.2018.01.029 29425488
24. Onuchic V, Lurie E, Carrero I, Pawliczek P, Patel RY, Rozowsky J, et al. Allele-specific epigenome maps reveal sequence-dependent stochastic switching at regulatory loci. Science (80-). 2018;361. doi: 10.1126/science.aar3146 30139913
25. Krivega I, Dale RK, Dean A. Role of LDB1 in the transition from chromatin looping to transcription activation. Genes Dev. 2014;28: 1278–90. doi: 10.1101/gad.239749.114 24874989
26. Lee J, Krivega I, Dale RK, Dean A. The LDB1 Complex Co-opts CTCF for Erythroid Lineage-Specific Long-Range Enhancer Interactions. Cell Rep. 2017;19: 2490–2502. doi: 10.1016/j.celrep.2017.05.072 28636938
27. Dixon JR, Selvaraj S, Yue F, Kim A, Li Y, Shen Y, et al. Topological domains in mammalian genomes identified by analysis of chromatin interactions. Nature. 2012;485: 376–80. doi: 10.1038/nature11082 22495300
28. Phanstiel DH, Van Bortle K, Spacek D, Hess GT, Shamim MS, Machol I, et al. Static and Dynamic DNA Loops form AP-1-Bound Activation Hubs during Macrophage Development. Mol Cell. 2017;67: 1037–1048.e6. doi: 10.1016/j.molcel.2017.08.006 28890333
29. Jolma A, Yan J, Whitington T, Toivonen J, Nitta KR, Rastas P, et al. DNA-binding specificities of human transcription factors. Cell. 2013;152: 327–39. doi: 10.1016/j.cell.2012.12.009 23332764
30. Shi W, Fornes O, Mathelier A, Wasserman WW. Evaluating the impact of single nucleotide variants on transcription factor binding. Nucleic Acids Res. 2016;44: 10106–10116. doi: 10.1093/nar/gkw691 27492288
31. Behera V, Evans P, Face CJ, Hamagami N, Sankaranarayanan L, Keller CA, et al. Exploiting genetic variation to uncover rules of transcription factor binding and chromatin accessibility. Nat Commun. 2018;9: 782. doi: 10.1038/s41467-018-03082-6 29472540
32. Timothy E, Reddy TE, Gertz J, Pauli F, Kucera KS, Varley KE, et al. The effects of genome sequence on differential allelic transcription factor occupancy and gene expression. Genome Res. 2012;22: 860–869. doi: 10.1101/gr.131201.111 22300769
33. Deplancke B, Alpern D, Gardeux V. The Genetics of Transcription Factor DNA Binding Variation. Cell. 2016;166: 538–554. doi: 10.1016/j.cell.2016.07.012 27471964
34. Fujiki K, Shinoda A, Kano F, Sato R, Shirahige K, Murata M. PPARγ-induced PARylation promotes local DNA demethylation by production of 5-hydroxymethylcytosine. Nat Commun. 2013;4: 2262. doi: 10.1038/ncomms3262 23912449
35. Suzuki T, Shimizu Y, Furuhata E, Maeda S, Kishima M, Nishimura H, et al. RUNX1 regulates site specificity of DNA demethylation by recruitment of DNA demethylation machineries in hematopoietic cells. Blood Adv. 2017;1: 1699–1711. doi: 10.1182/bloodadvances.2017005710 29296817
36. Hervouet E, Peixoto P, Delage-Mourroux R, Boyer-Guittaut M, Cartron P-F. Specific or not specific recruitment of DNMTs for DNA methylation, an epigenetic dilemma. Clin Epigenetics. 2018;10: 17. doi: 10.1186/s13148-018-0450-y 29449903
37. Sato N, Kondo M, Arai K. The orphan nuclear receptor GCNF recruits DNA methyltransferase for Oct-3/4 silencing. Biochem Biophys Res Commun. 2006;344: 845–51. doi: 10.1016/j.bbrc.2006.04.007 16631596
38. Schoenherr CJ, Levorse JM, Tilghman SM. CTCF maintains differential methylation at the Igf2/H19 locus. Nat Genet. 2003;33: 66–9. doi: 10.1038/ng1057 12461525
39. Sugden K, Hannon EJ, Arseneault L, Belsky DW, Corcoran DL, Fisher HL, et al. Patterns of Reliability: Assessing the Reproducibility and Integrity of DNA Methylation Measurement. Patterns. 2020;1: 100014. doi: 10.1016/j.patter.2020.100014 32885222
40. Pidsley R, Zotenko E, Peters TJ, Lawrence MG, Risbridger GP, Molloy P, et al. Critical evaluation of the Illumina MethylationEPIC BeadChip microarray for whole-genome DNA methylation profiling. Genome Biol. 2016;17: 208. doi: 10.1186/s13059-016-1066-1 27717381
41. Gelb B, Brueckner M, Chung W, Goldmuntz E, Kaltman J, Pablo Kaski J, et al. The Congenital Heart Disease Genetic Network Study: Rationale, Design, and Early Results. Circ Res. 2013;112: 698–706. doi: 10.1161/CIRCRESAHA.111.300297 23410879
42. Hoang TT, Goldmuntz E, Roberts AE, Chung WK, Kline JK, Deanfield JE, et al. The Congenital Heart Disease Genetic Network Study: Cohort description. PLoS One. 2018;13: e0191319. doi: 10.1371/journal.pone.0191319 29351346
43. Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009;25: 1754–60. doi: 10.1093/bioinformatics/btp324 19451168
44. McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, et al. The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010;20: 1297–303. doi: 10.1101/gr.107524.110 20644199
45. DePristo MA, Banks E, Poplin R, Garimella K V, Maguire JR, Hartl C, et al. A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat Genet. 2011;43: 491–8. doi: 10.1038/ng.806 21478889
46. Auwera GA Van Der, Carneiro MO, Hartl C, Poplin R, Levy-moonshine A, Jordan T, et al. From FastQ data to high confidence varant calls: the Genonme Analysis Toolkit best practices pipeline. Curr Protoc Bioinformatics. 2014. doi: 10.1002/0471250953.bi1110s43.From
47. Du P, Kibbe W a., Lin SM. lumi: A pipeline for processing Illumina microarray. Bioinformatics. 2008;24: 1547–1548. doi: 10.1093/bioinformatics/btn224 18467348
48. Teschendorff AE, Marabita F, Lechner M, Bartlett T, Tegner J, Gomez-Cabrero D, et al. A beta-mixture quantile normalization method for correcting probe design bias in Illumina Infinium 450 k DNA methylation data. Bioinformatics. 2013;29: 189–96. doi: 10.1093/bioinformatics/bts680 23175756
49. Houseman EA, Accomando WP, Koestler DC, Christensen BC, Marsit CJ, Nelson HH, et al. DNA methylation arrays as surrogate measures of cell mixture distribution. BMC Bioinformatics. 2012;13: 86. doi: 10.1186/1471-2105-13-86 22568884
50. Wang J, Zhuang J, Iyer S, Jie Wang A, Lin X, Whitfield TW, et al. Sequence features and chromatin structure around the genomic regions bound by 119 human transcription factors Repository Citation Sequence features and chromatin structure around the genomic regions bound by 119 human transcription factors. Genome Res. 2012;9: 1798–1812. doi: 10.1101/gr.139105.112 22955990
51. Richter F, Hoffman GE, Manheimer KB, Patel N, Sharp AJ, McKean D, et al. ORE Identifies Extreme Expression Effects Enriched for Rare Variants. Bioinformatics. 2019. doi: 10.1093/bioinformatics/btz202 30903145
52. Liao Y, Smyth GK, Shi W. The Subread aligner: fast, accurate and scalable read mapping by seed-and-vote. Nucleic Acids Res. 2013;41: e108. doi: 10.1093/nar/gkt214 23558742
53. Liao Y, Smyth GK, Shi W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. 2014;30: 923–30. doi: 10.1093/bioinformatics/btt656 24227677
Článek vyšel v časopise
PLOS Genetics
2020 Číslo 11
- Prof. Libor Grubhoffer: Nová varianta mpox je příležitostí pokročit v poznání původce varioly
- Co obsahují nové pokyny WHO pro dozor nad antibiotickou rezistencí?
- AI místo lékaře? Výzkum zkoumal ochotu Čechů konzultovat své zdraví s umělou inteligencí
- Prof. Jan Škrha: Metformin je bezpečný, ale je třeba jej bezpečně užívat a léčbu kontrolovat
- Ibuprofen jako alternativa antibiotik při léčbě infekcí močových cest
Nejčtenější v tomto čísle
- Stability of SARS-CoV-2 phylogenies
- Formal commentary
- No association between SCN9A and monogenic human epilepsy disorders
- Oxidative stress antagonizes fluoroquinolone drug sensitivity via the SoxR-SUF Fe-S cluster homeostatic axis