scholarly journals Haplotype-based eQTL mapping finds evidence for complex gene regulatory regions poorly tagged by marginal SNPs

2018 ◽  
Author(s):  
Robert Brown ◽  
Sriram Sankararaman ◽  
Bogdan Pasaniuc

AbstractMotivationExpression quantitative trait loci (eQTLs), variations in the genome that impact gene expression, are identified through eQTL studies that test for a relationship between single nucleotide polymorphisms (SNPs) and gene expression levels. These studies typically assume an underlying additive model. Non-additive tests have been proposed, but are limited due to the increase in the multiple testing burden and are potentially biased by filtering criteria that relies on marginal association data. Here we propose using combinations of short haplotypes instead of SNPs as predictors for gene expression. Essentially, this method looks for genomic regions where haplotypes have different effect sizes. The differences in effect can be due to multiple genetic architectures such as a single SNP, a burden of rare SNPs, multiple SNPs with independent effect or multiple SNPs with an interaction effect occurring on the same haplotype.ResultsSimulations show that when haplotypes, rather than SNPs, are assigned non-zero effect sizes, our method has increased power compared to the marginal SNP method. In the GEUVADIS gene expression data, our method finds 101 more eGenes than the marginal method (5,202 vs. 5,101). The methods do not have full overlap in the eGenes that they find. Of the 5,202 eGenes found by our method, 707 are not found by the marginal method—even though it has a lower significance threshold. This indicates that many genes have regulatory architectures that are not well tagged by marginal SNPs and demonstrates the need to better model alternative archi-tectures.

2017 ◽  
Vol 2017 ◽  
pp. 1-5 ◽  
Author(s):  
Lijun Wu ◽  
Liwang Gao ◽  
Xiaoyuan Zhao ◽  
Meixian Zhang ◽  
Jianxin Wu ◽  
...  

Purpose. Genome-wide association studies have found two obesity-related single-nucleotide polymorphisms (SNPs), rs17782313 near the melanocortin-4 receptor (MC4R) gene and rs6265 near the brain-derived neurotrophic factor (BDNF) gene, but the associations of both SNPs with other obesity-related traits are not fully described, especially in children. The aim of the present study is to investigate the associations between the SNPs and adiponectin that has a regulatory role in glucose and lipid metabolism. Methods. We examined the associations of the SNPs with adiponectin in Beijing Child and Adolescent Metabolic Syndrome (BCAMS) study. A total of 3503 children participated in the study. Results. The SNP rs6265 was significantly associated with adiponectin under an additive model (P=0.02 and 0.024, resp.) after adjustment for age, gender, and BMI or obesity statuses. The SNP rs17782313 was significantly associated with low adiponectin under a recessive model. No statistical significance was found between the two SNPs and low adiponectin after correction for multiple testing. Conclusion. We demonstrate for the first time that the SNP rs17782313 near MC4R and the SNP rs6265 near BDNF are associated with adiponectin in Chinese children. These novel findings provide important evidence that adiponectin possibly mediates MC4R and BDNF involved in obesity.


2020 ◽  
Vol 117 (26) ◽  
pp. 15028-15035 ◽  
Author(s):  
Ronald Yurko ◽  
Max G’Sell ◽  
Kathryn Roeder ◽  
Bernie Devlin

To correct for a large number of hypothesis tests, most researchers rely on simple multiple testing corrections. Yet, new methodologies of selective inference could potentially improve power while retaining statistical guarantees, especially those that enable exploration of test statistics using auxiliary information (covariates) to weight hypothesis tests for association. We explore one such method, adaptiveP-value thresholding (AdaPT), in the framework of genome-wide association studies (GWAS) and gene expression/coexpression studies, with particular emphasis on schizophrenia (SCZ). Selected SCZ GWAS associationPvalues play the role of the primary data for AdaPT; single-nucleotide polymorphisms (SNPs) are selected because they are gene expression quantitative trait loci (eQTLs). This natural pairing of SNPs and genes allow us to map the following covariate values to these pairs: GWAS statistics from genetically correlated bipolar disorder, the effect size of SNP genotypes on gene expression, and gene–gene coexpression, captured by subnetwork (module) membership. In all, 24 covariates per SNP/gene pair were included in the AdaPT analysis using flexible gradient boosted trees. We demonstrate a substantial increase in power to detect SCZ associations using gene expression information from the developing human prefrontal cortex. We interpret these results in light of recent theories about the polygenic nature of SCZ. Importantly, our entire process for identifying enrichment and creating features with independent complementary data sources can be implemented in many different high-throughput settings to ultimately improve power.


2014 ◽  
Vol 2014 ◽  
pp. 1-8 ◽  
Author(s):  
Patricia Vázquez ◽  
Otsanda Ruiz-Larrañaga ◽  
Joseba M. Garrido ◽  
Mikel Iriondo ◽  
Carmen Manzano ◽  
...  

A genetic susceptibility toMycobacterium aviumsubsp.paratuberculosis(MAP) infections in ruminants has been longtime suspected to exist. Recently, natural infections in cattle have been reclassified intolatentandpatentforms based on histopathological findings and their associations with immunological and microbiological variables. This study aims to explore whether these newly defined phenotypes are associated with twenty-four single-nucleotide polymorphisms (SNPs) in six bovine candidate genes:nucleotide-binding oligomerization domain 2 (NOD2), solute carrier family 11 member A1 (SLC11A1), nuclear body protein SP110 (SP110), toll-like receptors (TLRs) 2and4, andCD209(also known asDC-SIGN, dendritic cell-specific ICAM3-grabbing nonintegrin). SNPs were genotyped for 772 Holstein-Friesian animals (52.6%apparently free; 38.1%latent; 9.3%patent) by TaqMan OpenArray technology. Genotypic-phenotypic associations were assessed by logistic regression analysis adjusted for age at slaughter, under five models (codominant, dominant, recessive, overdominant, and log-additive), and corrected for multiple testing. The rs208222804 C allele (CD209gene) was found to be associated withlatentparatuberculosis (log-additive model:P<0.0034after permutation procedure; OR = 0.64, 95% CI = 0.48–0.86). No significant association was detected between any SNP and thepatentphenotype. Consequently,CD209gene may play a key role in the pathogenesis of bovine paratuberculosis.


Author(s):  
Daniel Berner ◽  
Valentin Amrhein

A paradigm shift away from null hypothesis significance testing seems in progress. Based on simulations, we illustrate some of the underlying motivations. First, P-values vary strongly from study to study, hence dichotomous inference using significance thresholds is usually unjustified. Second, statistically significant results have overestimated effect sizes, a bias declining with increasing statistical power. Third, statistically non-significant results have underestimated effect sizes, and this bias gets stronger with higher statistical power. Fourth, the tested statistical hypotheses generally lack biological justification and are often uninformative. Despite these problems, a screen of 48 papers from the 2020 volume of the Journal of Evolutionary Biology exemplifies that significance testing is still used almost universally in evolutionary biology. All screened studies tested the default null hypothesis of zero effect with the default significance threshold of p = 0.05, none presented a pre-planned alternative hypothesis, and none calculated statistical power and the probability of &lsquo;false negatives&rsquo; (beta error). The papers reported 49 significance tests on average. Of 41 papers that contained verbal descriptions of a &lsquo;statistically non-significant&rsquo; result, 26 (63%) falsely claimed the absence of an effect. We conclude that our studies in ecology and evolutionary biology are mostly exploratory and descriptive. We should thus shift from claiming to &ldquo;test&rdquo; specific hypotheses statistically to describing and discussing many hypotheses (effect sizes) that are most compatible with our data, given our statistical model. We already have the means for doing so, because we routinely present compatibility (&ldquo;confidence&rdquo;) intervals covering these hypotheses.


Nutrients ◽  
2019 ◽  
Vol 11 (4) ◽  
pp. 935 ◽  
Author(s):  
Veronika Fedirko ◽  
Mazda Jenab ◽  
Catherine Méplan ◽  
Jeb S. Jones ◽  
Wanzhe Zhu ◽  
...  

Selenoprotein genetic variations and suboptimal selenium (Se) levels may contribute to the risk of colorectal cancer (CRC) development. We examined the association between CRC risk and genotype for single nucleotide polymorphisms (SNPs) in selenoprotein and Se metabolic pathway genes. Illumina Goldengate assays were designed and resulted in the genotyping of 1040 variants in 154 genes from 1420 cases and 1421 controls within the European Prospective Investigation into Cancer and Nutrition (EPIC) study. Multivariable logistic regression revealed an association of 144 individual SNPs from 63 Se pathway genes with CRC risk. However, regarding the selenoprotein genes, only TXNRD1 rs11111979 retained borderline statistical significance after adjustment for correlated tests (PACT = 0.10; PACT significance threshold was P < 0.1). SNPs in Wingless/Integrated (Wnt) and Transforming growth factor (TGF) beta-signaling genes (FRZB, SMAD3, SMAD7) from pathways affected by Se intake were also associated with CRC risk after multiple testing adjustments. Interactions with Se status (using existing serum Se and Selenoprotein P data) were tested at the SNP, gene, and pathway levels. Pathway analyses using the modified Adaptive Rank Truncated Product method suggested that genes and gene x Se status interactions in antioxidant, apoptosis, and TGF-beta signaling pathways may be associated with CRC risk. This study suggests that SNPs in the Se pathway alone or in combination with suboptimal Se status may contribute to CRC development.


Neurology ◽  
2017 ◽  
Vol 89 (17) ◽  
pp. 1829-1839 ◽  
Author(s):  
Kristiina Rannikmäe ◽  
Vhinoth Sivakumaran ◽  
Henry Millar ◽  
Rainer Malik ◽  
Christopher D. Anderson ◽  
...  

Objective:To determine whether common variants in familial cerebral small vessel disease (SVD) genes confer risk of sporadic cerebral SVD.Methods:We meta-analyzed genotype data from individuals of European ancestry to determine associations of common single nucleotide polymorphisms (SNPs) in 6 familial cerebral SVD genes (COL4A1, COL4A2, NOTCH3, HTRA1, TREX1, and CECR1) with intracerebral hemorrhage (ICH) (deep, lobar, all; 1,878 cases, 2,830 controls) and ischemic stroke (IS) (lacunar, cardioembolic, large vessel disease, all; 19,569 cases, 37,853 controls). We applied data quality filters and set statistical significance thresholds accounting for linkage disequilibrium and multiple testing.Results:A locus in COL4A2 was associated (significance threshold p < 3.5 × 10−4) with both lacunar IS (lead SNP rs9515201: odds ratio [OR] 1.17, 95% confidence interval [CI] 1.11–1.24, p = 6.62 × 10−8) and deep ICH (lead SNP rs4771674: OR 1.28, 95% CI 1.13–1.44, p = 5.76 × 10−5). A SNP in HTRA1 was associated (significance threshold p < 5.5 × 10−4) with lacunar IS (rs79043147: OR 1.23, 95% CI 1.10–1.37, p = 1.90 × 10−4) and less robustly with deep ICH. There was no clear evidence for association of common variants in either COL4A2 or HTRA1 with non-SVD strokes or in any of the other genes with any stroke phenotype.Conclusions:These results provide evidence of shared genetic determinants and suggest common pathophysiologic mechanisms of distinct ischemic and hemorrhagic cerebral SVD stroke phenotypes, offering new insights into the causal mechanisms of cerebral SVD.


2021 ◽  
Vol 21 (1) ◽  
Author(s):  
Yu Zhang ◽  
Li Hua ◽  
Quan-Hua Liu ◽  
Shu-Yuan Chu ◽  
Yue-Xin Gan ◽  
...  

Abstract Background A number of studies have examined the association between mold exposure and childhood asthma. However, the conclusions were inconsistent, which might be partly attributable to the lack of consideration of gene function, especially the key genes affecting the pathogenesis of childhood asthma. Research on the interactions between genes and mold exposure on childhood asthma is still very limited. We therefore examined whether there is an interaction between inflammation-related genes and mold exposure on childhood asthma. Methods A case–control study with 645 asthmatic children and 910 non-asthmatic children aged 3–12 years old was conducted. Eight single nucleotide polymorphisms (SNPs) in inflammation-related genes were genotyped using MassARRAY assay. Mold exposure was defined as self-reported visible mold on the walls. Associations between visible mold exposure, SNPs and childhood asthma were evaluated using logistic regression models. In addition, crossover analyses were used to estimate the gene-environment interactions on childhood asthma on an additive scale. Results After excluding children without information on visible mold exposure or SNPs, 608 asthmatic and 839 non-asthmatic children were included in the analyses. Visible mold exposure was reported in 151 asthmatic (24.8%) and 119 non-asthmatic children (14.2%) (aOR 2.19, 95% CI 1.62–2.97). The rs7216389 SNP in gasdermin B gene (GSDMB) increased the risk of childhood asthma with each C to T substitution in a dose-dependent pattern (additive model, aOR 1.32, 95% CI 1.11–1.57). Children carrying the rs7216389 T allele and exposed to visible mold dramatically increased the risk of childhood asthma (aOR 3.21; 95% CI 1.77–5.99). The attributable proportion due to the interaction (AP: 0.47, 95% CI 0.03–0.90) and the relative excess risk due to the interaction (RERI: 1.49, 95% CI 0–2.99) were statistically significant. Conclusions In the present study, there was a significant additive interaction between visible mold exposure and rs7216389 SNP on childhood asthma. Future studies need to consider the gene-environment interactions when exploring the risk factors of childhood asthma.


2021 ◽  
Vol 21 (1) ◽  
Author(s):  
Erika Calvano Küchler ◽  
Agnes Schröder ◽  
Vinicius Broska Teodoro ◽  
Ute Nazet ◽  
Rafaela Scariot ◽  
...  

Abstract Background This study aimed to investigate, if different physiological concentrations of vitamin D (25(OH)D3) and single nucleotide polymorphisms in vitamin D receptor (VDR) gene have an impact on gene expression in human periodontal ligament (hPDL) fibroblasts induced by simulated orthodontic compressive strain. Methods A pool of hPDL fibroblasts was treated in absence or presence of 25(OH)D3 in 3 different concentrations (10, 40 and 60 ng/ml). In order to evaluate the role of single nucleotide polymorphisms in the VDR gene, hPDL fibroblasts from 9 patients were used and treated in absence or presence of 40 ng/ml 25(OH)D3. Each experiment was performed with and without simulated orthodontic compressive strain. Real-time PCR was used for gene expression and allelic discrimination analysis. Relative expression of dehydrocholesterol reductase (DHCR7), Sec23 homolog A, amidohydrolase domain containing 1 (AMDHD1), vitamin D 25-hydroxylase (CYP2R1), Hydroxyvitamin D-1-α hydroxylase, receptor activator of nuclear factor-κB ligand (RANKL), osteoprotegerin (OPG), cyclooxygenase-2 (COX-2) and interleukin-6 (IL6) was assessed. Three single nucleotide polymorphisms in VDR were genotyped. Parametric or non-parametric tests were used with an alpha of 5%. Results RANKL, RANKL:OPG ratio, COX-2, IL-6, DHCR7, CYP2R1 and AMDHD1 were differentially expressed during simulated orthodontic compressive strain (p < 0.05). The RANKL:OPG ratio was downregulated by all concentrations (10 ng/ml, 40 ng/ml and 60 ng/ml) of 25(OH)D3 (mean = 0.96 ± 0.68, mean = 1.61 ± 0.66 and mean = 1.86 ± 0.78, respectively) in comparison to the control (mean 2.58 ± 1.16) (p < 0.05). CYP2R1 gene expression was statistically modulated by the different 25(OH)D3 concentrations applied (p = 0.008). Samples from individuals carrying the GG genotype in rs739837 presented lower VDR mRNA expression and samples from individuals carrying the CC genotype in rs7975232 presented higher VDR mRNA expression (p < 0.05). Conclusions Simulated orthodontic compressive strain and physiological concentrations of 25(OH)D3 seem to regulate the expression of orthodontic tooth movement and vitamin-D-related genes in periodontal ligament fibroblasts in the context of orthodontic compressive strain. Our study also suggests that single nucleotide polymorphisms in the VDR gene regulate VDR expression in periodontal ligament fibroblasts in the context of orthodontic compressive strain.


2010 ◽  
Vol 70 (4) ◽  
pp. 668-674 ◽  
Author(s):  
P Dieudé ◽  
M Guedj ◽  
J Wipff ◽  
B Ruiz ◽  
G Riemekasten ◽  
...  

BackgroundRecent evidence has highlighted a potential role of interleukin 1β (IL-1β) in systemic sclerosis (SSc). NLRP1 provides a scaffold for the assembly of the inflammasome that promotes the processing and maturation of pro-IL-1β. In addition, NLRP1 variants were found to confer susceptibility to autoimmune disorders.ObjectiveTo study a possible association of the NLRP1 rs6502867, rs2670660 and rs8182352, rs12150220 and rs4790797 with SSc in the European Caucasian population.MethodsNLRP1 single nucleotide polymorphisms were genotyped in 3227 individuals comprising a discovery set (870 SSc patients and 962 controls) and a replication set including individuals from Germany (532 SSc patients and 324 controls) and Italy (527 SSc patients and 301 controls), all individuals being of European Caucasian origin.ResultsConditional analyses revealed a significant association for the NLRP1 rs8182352 variant with both anti-topoisomerase-positive and SSc-related fibrosing alveolitis (FA) subsets under an additive model: p=0.0042, OR 1.23 (95% CI 1.07 to 1.41) and p=0.0065 OR 1.19 (95% CI 1.05 to 1.36), respectively. Logistic regression analysis showed an additive effect of IRF5 rs2004640, STAT4 rs7574865 and NLRP1 rs8182352 risk alleles on SSc-related FA.ConclusionsOur results establish NLRP1 as a new genetic susceptibility factor for SSc-related pulmonary fibrosis and anti-topoisomerase-positive SSc phenotypes. This provides new insights into the pathogenesis of SSc, underlining the potential role of innate immunity in particular in the FA-positive SSc subphenotype, which represents a severe subset of the disease.


Sign in / Sign up

Export Citation Format

Share Document