Identification of tissue-specific and common methylation quantitative trait loci in healthy individuals using MAGAR.
Ontology highlight
ABSTRACT: Background
Understanding the influence of genetic variants on DNA methylation is fundamental for the interpretation of epigenomic data in the context of disease. There is a need for systematic approaches not only for determining methylation quantitative trait loci (methQTL), but also for discriminating general from cell type-specific effects.Results
Here, we present a two-step computational framework MAGAR ( https://bioconductor.org/packages/MAGAR ), which fully supports the identification of methQTLs from matched genotyping and DNA methylation data, and additionally allows for illuminating cell type-specific methQTL effects. In a pilot analysis, we apply MAGAR on data in four tissues (ileum, rectum, T cells, B cells) from healthy individuals and demonstrate the discrimination of common from cell type-specific methQTLs. We experimentally validate both types of methQTLs in an independent data set comprising additional cell types and tissues. Finally, we validate selected methQTLs located in the PON1, ZNF155, and NRG2 genes by ultra-deep local sequencing. In line with previous reports, we find cell type-specific methQTLs to be preferentially located in enhancer elements.Conclusions
Our analysis demonstrates that a systematic analysis of methQTLs provides important new insights on the influences of genetic variants to cell type-specific epigenomic variation.
SUBMITTER: Scherer M
PROVIDER: S-EPMC8444396 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Interaction of DNA methylation and sequence variants that are methylation quantitative trait loci (mQTLs) may influence susceptibility to diseases such as alcohol dependence (AD). We used genome-wide genotype data from 268 African Americans (AAs: 129 AD cases and 139 controls) and 143 European Americans (EAs: 129 AD cases and 14 controls) to identify mQTLs that were associated with promoter CpGs in 82 AD risk genes. 282 significant mQTL-CpG pairs (9.9 × 10(-100) ? P(nominal) ? 7.7 × 10(-8)) in AAs and 313 significant mQTL-CpG pairs (2.7 × 10(-53) ? P(nominal) ? 9.9 × 10(-8)) in EAs were identified [i.e., mQTL-CpG associations survived multiple-testing correction, q values (false discovery rate) ? 0.05]. The most significant mQTL was rs1800759, which was strongly associated with CpG cg12011299 in both AAs (P(nominal) = 9.9 × 10(-100); q = 6.7 × 10(-91)) and EAs (P(nominal) = 2.7 × 10(-53); q = 1.4 × 10(-44)). Rs1800759 (previously known to be associated to AD) and CpG cg12011299 (distance: 37 bp) are both located in alcohol dehydrogenase (ADH) 4 gene (ADH4) promoter region. In general, the strength of association between mQTLs and CpGs was inversely correlated with the distance between them. Association was also influenced by race and AD. Additionally, 48.3 % of the mQTLs identified in AAs and 65.6 % of the mQTLs identified in EAs were predicted to be expression QTLs. Three mQTLs (rs2173201, rs4147542, and rs4147541 in ADH1B-AHD1C gene cluster region) found in AAs were previously identified by our genome-wide association studies as being significantly associated with AD in AAs. Thus, DNA methylation, which can be influenced by sequence variants and is implicated in gene expression regulation, appears to at least partially underlie the association of genetic variation with AD.
| S-EPMC4127343 | biostudies-literature
Project description:Individual genotypes at specific loci can result in different patterns of DNA methylation. These methylation quantitative trait loci (meQTLs) influence methylation across extended genomic regions and may underlie direct SNP associations or gene-environment interactions. We hypothesized that the detection of meQTLs varies with ancestral population, developmental stage, and tissue type. We explored this by analyzing seven datasets that varied by ancestry (African American vs. Caucasian), developmental stage (neonate vs. adult), and tissue type (blood vs. four regions of postmortem brain) with genome-wide DNA methylation and SNP data. We tested for meQTLs by constructing linear regression models of methylation levels at each CpG site on SNP genotypes within 50 kb under an additive model controlling for multiple tests.Most meQTLs mapped to intronic regions, although a limited number appeared to occur in synonymous or nonsynonymous coding SNPs. We saw significant overlap of meQTLs between ancestral groups, developmental stages, and tissue types, with the highest rates of overlap within the four brain regions. Compared with a random group of SNPs with comparable frequencies, meQTLs were more likely to be 1) represented among the most associated SNPs in the WTCCC bipolar disorder results and 2) located in microRNA binding sites.These data give us insight into how SNPs impact gene regulation and support the notion that peripheral blood may be a reliable correlate of physiological processes in other tissues.
| S-EPMC4028873 | biostudies-literature
Project description:We present an approach for quantitative trait locus (QTL) mapping, termed as 'lineage-specific QTL mapping', for inferring allelic changes of QTL evolution along with branches in a phylogeny. We describe and analyze the simplest case: by adding a third taxon into the normal procedure of QTL mapping between pairs of taxa, such inferences can be made along lineages to a presumed common ancestor. Although comparisons of QTL maps among species can identify homology of QTLs by apparent co-location, lineage-specific mapping of QTL can classify homology into (1) orthology (shared origin of QTL) versus (2) paralogy (independent origin of QTL within resolution of map distance). In this light, we present a graphical method that identifies six modes of QTL evolution in a three taxon comparison. We then apply our model to map lineage-specific QTLs for inbreeding among three taxa of yellow monkey-flower: Mimulus guttatus and two inbreeders M. platycalyx and M. micranthus, but critically assuming outcrossing was the ancestral state. The two most common modes of homology across traits were orthologous (shared ancestry of mutation for QTL alleles). The outbreeder M. guttatus had the fewest lineage-specific QTL, in accordance with the presumed ancestry of outbreeding. Extensions of lineage-specific QTL mapping to other types of data and crosses, and to inference of ancestral QTL state, are discussed.
| S-EPMC3716265 | biostudies-literature
Project description:DNA methylation may be regulated by genetic variants within a genomic region, referred to as methylation quantitative trait loci (mQTLs). The changes of methylation levels can further lead to alterations of gene expression, and influence the risk of various complex human diseases. Detecting mQTLs may provide insights into the underlying mechanism of how genotypic variations may influence the disease risk. In this article, we propose a methylation random field (MRF) method to detect mQTLs by testing the association between the methylation level of a CpG site and a set of genetic variants within a genomic region. The proposed MRF has two major advantages over existing approaches. First, it uses a beta distribution to characterize the bimodal and interval properties of the methylation trait at a CpG site. Second, it considers multiple common and rare genetic variants within a genomic region to identify mQTLs. Through simulations, we demonstrated that the MRF had improved power over other existing methods in detecting rare variants of relatively large effect, especially when the sample size is small. We further applied our method to a study of congenital heart defects with 83 cardiac tissue samples and identified two mQTL regions, MRPS10 and PSORS1C1, which were colocalized with expression QTL in cardiac tissue. In conclusion, the proposed MRF is a useful tool to identify novel mQTLs, especially for studies with limited sample sizes.
| S-EPMC8575051 | biostudies-literature
Project description:A genome-wide association study (GWAS) was conducted to examine expression quantitative trait loci (eQTLs) for histone genes. We examined common eQTLs for multiple histone genes in 373 European lymphoblastoid cell lines (LCLs). A linear regression model was employed to identify single-nucleotide polymorphisms (SNPs) associated with expression of the histone genes, and the number of eQTLs was determined by linkage disequilibrium analysis. Additional associations of the identified eQTLs with other genes were also examined. We identified 31 eQTLs for 29 histone genes through genome-wide analysis using 29 histone genes (P < 2.97 × 10-10). Among them, 12 eQTLs were associated with the expression of multiple histone genes. Transcriptome-wide association analysis using the identified eQTLs showed their associations with additional 80 genes (P < 4.75 × 10-6). In particular, expression of RPPH1, SCARNA2, and SCARNA7 genes was associated with 26, 25, and 23 eQTLs, respectively. This study suggests that histone genes shared 12 common eQTLs that might regulate cell cycle-dependent transcription of histone and other genes. Further investigations are needed to elucidate the transcriptional mechanisms of these genes.
| S-EPMC5591967 | biostudies-literature
Project description:Interpreting the susceptible loci documented by genome-wide association studies (GWASs) is of utmost importance in the post-GWAS era. Since most complex traits are contributed by multiple tissues, analyzing tissue-specific effects of expression quantitative trait loci (eQTLs) is a promising approach. Here we describe "opposite eQTL effects", i.e., gene expression effects of eQTLs that are in the opposite direction between different tissues, as the biologically meaningful annotations of genes and genetic variants for understanding the GWAS loci. The genes and single-nucleotide polymorphisms (SNPs) associated with the opposite eQTL effects (opp-multi-eQTL-Genes and opp-multi-eQTL-SNPs) were extracted from the largest eQTL database provided by the Genotype-Tissue Expression (GTEx) project (release version 7). The opposite eQTL effects were detected even between closely related tissues such as cerebellum and brain cortex, and a significant proportion of the genes having eQTLs were annotated as the opp-multi-eQTL-Genes (2,323 out of 31,212; 7.4%). The opp-multi-eQTL-SNPs showed locational enrichment at the transcription start site and also possible involvement of epigenetic regulation. The biological importance of the opposite eQTL effects was also assessed using the SNPs reported in GWASs (GWAS-SNPs), which demonstrated that a high proportion of the opp-multi-eQTL-SNPs are in linkage disequilibrium with the GWAS-SNPs (2,498 out of 9,290; 26.9%). Based on the results, the opposite eQTL effects can be a common phenomenon in the tissue-specific gene regulation with a possible contribution to the development of complex traits.
| S-EPMC6871526 | biostudies-literature
Project description:Variation in gene expression is heritable and has been mapped to the genome in humans and model organisms as expression quantitative trait loci (eQTLs). We applied integrated genome-wide expression profiling and linkage analysis to the regulation of gene expression in fat, kidney, adrenal, and heart tissues using the BXH/HXB panel of rat recombinant inbred strains. Here, we report the influence of heritability and allelic effect of the quantitative trait locus on detection of cis- and trans-acting eQTLs and discuss how these factors operate in a tissue-specific context. We identified several hundred major eQTLs in each tissue and found that cis-acting eQTLs are highly heritable and easier to detect than trans-eQTLs. The proportion of heritable expression traits was similar in all tissues; however, heritability alone was not a reliable predictor of whether an eQTL will be detected. We empirically show how the use of heritability as a filter reduces the ability to discover trans-eQTLs, particularly for eQTLs with small effects. Only 3% of cis- and trans-eQTLs exhibited large allelic effects, explaining more than 40% of the phenotypic variance, suggestive of a highly polygenic control of gene expression. Power calculations indicated that, across tissues, minor differences in genetic effects are expected to have a significant impact on detection of trans-eQTLs. Trans-eQTLs generally show smaller effects than cis-eQTLs and have a higher false discovery rate, particularly in more heterogeneous tissues, suggesting that small biological variability, likely relating to tissue composition, may influence detection of trans-eQTLs in this system. We delineate the effects of genetic architecture on variation in gene expression and show the sensitivity of this experimental design to tissue sampling variability in large-scale eQTL studies.
| S-EPMC1617131 | biostudies-literature
Project description:BackgroundKernel weight and morphology are important traits affecting cereal yields and quality. Dissecting the genetic basis of thousand kernel weight (TKW) and its related traits is an effective method to improve wheat yield.ResultsIn this study, we performed quantitative trait loci (QTL) analysis using recombinant inbred lines derived from the cross 'PuBing3228 × Gao8901' (PG-RIL) to dissect the genetic basis of kernel traits. A total of 17 stable QTLs related to kernel traits were identified, notably, two stable QTLs QTkw.cas-1A.2 and QTkw.cas-4A explained the largest portion of the phenotypic variance for TKW and kernel length (KL), and the other two stable QTLs QTkw.cas-6A.1 and QTkw.cas-7D.2 contributed more effects on kernel width (KW). Conditional QTL analysis revealed that the stable QTLs for TKW were mainly affected by KW. The QTLs QTkw.cas-7D.2 and QKw.cas-7D.1 associated with TKW and KW were delimited to the physical interval of approximately 3.82 Mb harboring 47 candidate genes. Among them, the candidate gene TaFT-D1 had a 1 bp insertions/deletion (InDel) within the third exon, which might be the reason for diversity in TKW and KW between the two parents. A Kompetitive Allele-Specific PCR (KASP) marker of TaFT-D1 allele was developed and verified by PG-RIL and a natural population consisted of 141 cultivar/lines. It was found that the favorable TaFT-D1 (G)-allele has been positively selected during Chinese wheat breeding. Thus, these results can be used for further positional cloning and marker-assisted selection in wheat breeding programs.ConclusionsSeventeen stable QTLs related to kernel traits were identified. The stable QTLs for thousand kernel weight were mainly affected by kernel width. TaFT-D1 could be the candidate gene for QTLs QTkw.cas-7D.2 and QKw.cas-7D.1.
| S-EPMC7682089 | biostudies-literature
Project description:DNA methylation (DNAm) is an epigenetic regulator of gene expression and a hallmark of gene-environment interaction. Using whole-genome bisulfite sequencing, we have surveyed DNAm in 344 samples of human postmortem brain tissue from neurotypical subjects and individuals with schizophrenia. We identify genetic influence on local methylation levels throughout the genome, both at CpG sites and CpH sites, with 86% of SNPs and 55% of CpGs being part of methylation quantitative trait loci (meQTLs). These associations can further be clustered into regions that are differentially methylated by a given SNP, highlighting the genes and regions with which these loci are epigenetically associated. These findings can be used to better characterize schizophrenia GWAS-identified variants as epigenetic risk variants. Regions differentially methylated by schizophrenia risk-SNPs explain much of the heritability associated with risk loci, despite covering only a fraction of the genomic space. We provide a comprehensive, single base resolution view of association between genetic variation and genomic methylation, and implicate schizophrenia GWAS-associated variants as influencing the epigenetic plasticity of the brain.
| S-EPMC8413445 | biostudies-literature
Project description:The analysis of gene sets is usually carried out based on gene ontology terms and known biological pathways. These approaches may not establish any formal relation between genotype and trait specific phenotype. In plant biology and breeding, analysis of gene sets with trait specific Quantitative Trait Loci (QTL) data are considered as great source for biological knowledge discovery. Therefore, we proposed an innovative statistical approach called Gene Set Analysis with QTLs (GSAQ) for interpreting gene expression data in context of gene sets with traits. The utility of GSAQ was studied on five different complex abiotic and biotic stress scenarios in rice, which yields specific trait/stress enriched gene sets. Further, the GSAQ approach was more innovative and effective in performing gene set analysis with underlying QTLs and identifying QTL candidate genes than the existing approach. The GSAQ approach also provided two potential biological relevant criteria for performance analysis of gene selection methods. Based on this proposed approach, an R package, i.e., GSAQ ( https://cran.r-project.org/web/packages/GSAQ ) has been developed. The GSAQ approach provides a valuable platform for integrating the gene expression data with genetically rich QTL data.
| S-EPMC5799309 | biostudies-literature