Evidence for Non-additive Influence of Single Nucleotide Polymorphisms within the Apolipoprotein E Gene

Quantitative trait loci, G×E and G×G for glycemic traits: response to metformin and placebo in the Diabetes Prevention Program (DPP)

The complex genetic architecture of type-2-diabetes (T2D) includes gene-by-environment (G×E) and gene-by-gene (G×G) interactions. To identify G×E and G×G, we screened markers for patterns indicative of interactions (relationship loci [rQTL] and variance heterogeneity loci [vQTL]). rQTL exist when the correlation between multiple traits varies by genotype and vQTL occur when the variance of a trait differs by genotype (potentially flagging G×G and G×E). In the metformin and placebo arms of the DPP (n = 1762) we screened 280,965 exomic and intergenic SNPs, for rQTL and vQTL patterns in association with year one changes from baseline in glycemia and related traits (insulinogenic index [IGI], insulin sensitivity index [ISI], fasting glucose and fasting insulin). Significant (p < 1.8 × 10^-7) rQTL and vQTL generated a priori hypotheses of individual G×E tests for a SNP × metformin treatment interaction and secondarily for G×G screens. Several rQTL and vQTL identified led to 6 nominally significant (p < 0.05) metformin treatment × SNP interactions (4 for IGI, one insulin, and one glucose) and 12G×G interactions (all IGI) that exceeded experiment-wide significance (p < 4.1 × 10^-9). Some loci are directly associated with incident diabetes, and others are rQTL and modify a trait's relationship with diabetes (2 diabetes/glucose, 2 diabetes/insulin, 1 diabetes/IGI). rs3197999, an ISI/insulin rQTL, is a possible gene damaging missense mutation in MST1, is associated with ulcerative colitis, sclerosing cholangitis, Crohn's disease, BMI and coronary artery disease. This study demonstrates evidence for context-dependent effects (G×G & G×E) and the complexity of these T2D-related traits.

Genome-wide association study for variants that modulate relationships between cerebrospinal fluid amyloid-beta 42, tau, and p-tau levels

Background

A relationship quantitative trait locus exists when the correlation between multiple traits varies by genotype for that locus. Relationship quantitative trait loci (rQTL) are often involved in gene-by-gene (G×G) interactions or gene-by-environmental interactions, making them a powerful tool for detecting G×G.

Methods

We performed genome-wide association studies to identify rQTL between tau and Aβ42 and ptau and Aβ42 with over 3000 individuals using age, gender, series, APOE ε2, APOE ε4, and two principal components for population structure as covariates. Each significant rQTL was separately screened for interactions with other loci for each trait in the rQTL model. Parametric bootstrapping was used to assess significance.

Results

We found four significant tau/Aβ42 rQTL from three unique locations and six ptau/Aβ42 rQTL from five unique locations. G×G screens with these rQTL produced four significant G×G interactions (one Aβ42, two ptau, and one tau) with four rQTL where each second locus was from a unique location. On follow-up, rs1036819 and rs74025622 were associated with Alzheimer’s disease (AD) case/control status; rs15205 and rs79099429 were associated with rate of decline.

Conclusions

The two most significant rQTL (rs8027714 and rs1036819) for ptau/Aβ42 are on different chromosomes and both are strong hits for pelvic organ prolapse. While diseases of the nervous system can cause pelvic organ prolapse, it is unlikely related to the ptau/Aβ42 relationship but may suggest that these two loci share a pathway. In addition to a ptau/Aβ42 rQTL and association with AD case/control status, rs1036819 is a strong rQTL for case/control status/Aβ42 and for tau/Aβ42. It resides in the ZFAT gene, which is related to autoimmune thyroid disease. For tau, rs9817620 interacts with the tau/Aβ42 rQTL rs74025622. It is in the CHL1 gene, which is a neural cell adhesion molecule and may be involved in signal transduction pathways. CHL1 is related to BACE1, which is a β-secretase enzyme that initiates production of the β-amyloid peptide involved in AD and is a primary drug target. Overall, there are numerous loci that affect the relationship between these important AD endophenotypes and some are due to interactions with other loci. Some affect the risk of AD and/or rate of progression.

Electronic supplementary material

The online version of this article (10.1186/s13195-018-0410-y) contains supplementary material, which is available to authorized users.

Genetic effects at pleiotropic loci are context-dependent with consequences for the maintenance of genetic variation in populations

Context-dependent genetic effects, including genotype-by-environment and genotype-by-sex interactions, are a potential mechanism by which genetic variation of complex traits is maintained in populations. Pleiotropic genetic effects are also thought to play an important role in evolution, reflecting functional and developmental relationships among traits. We examine context-dependent genetic effects at pleiotropic loci associated with normal variation in multiple metabolic syndrome (MetS) components (obesity, dyslipidemia, and diabetes-related traits). MetS prevalence is increasing in Western societies and, while environmental in origin, presents substantial variation in individual response. We identify 23 pleiotropic MetS quantitative trait loci (QTL) in an F₁₆ advanced intercross between the LG/J and SM/J inbred mouse strains (Wustl:LG,SM-G16; n = 1002). Half of each family was fed a high-fat diet and half fed a low-fat diet; and additive, dominance, and parent-of-origin imprinting genotypic effects were examined in animals partitioned into sex, diet, and sex-by-diet cohorts. We examine the context-dependency of the underlying additive, dominance, and imprinting genetic effects of the traits associated with these pleiotropic QTL. Further, we examine sequence polymorphisms (SNPs) between LG/J and SM/J as well as differential expression of positional candidate genes in these regions. We show that genetic associations are different in different sex, diet, and sex-by-diet settings. We also show that over- or underdominance and ecological cross-over interactions for single phenotypes may not be common, however multidimensional synthetic phenotypes at loci with pleiotropic effects can produce situations that favor the maintenance of genetic variation in populations. Our findings have important implications for evolution and the notion of personalized medicine.

We look at gene-by-diet and gene-by-sex interactions underlying natural variation in multiple metabolic traits mapping to the same regions of the genome in a mouse model. We find that the underlying genetic architecture of these traits is different in different sex and diet contexts. We further use expression data and whole-genome polymorphism data to identify compelling candidates for experimental follow-up. We use these results to examine theoretical evolutionary predictions about how variation in populations can be maintained. There has been much discussion of late on how to use evolutionary theory to inform medical genomics. Mouse models may be especially appropriate for bridging the divide between evolutionary and biomedical research, because they allow the study of the effects of natural alleles on normal variation and because human-mouse homology is well defined. Our study is unique in examining quantitative trait loci from both evolutionary and biomedical perspectives, and we highlight the complex connections of the traits comprising the metabolic syndrome and the evolutionary implications of their underlying genetic architecture. This is important for understanding disease etiology and is relevant to personalized medicine.

Analysis of gene-gene interactions

The goal of this unit is to introduce gene-gene interactions (epistasis) as a significant complicating factor in the search for disease susceptibility genes. This unit begins with an overview of gene-gene interactions and why they are likely to be common. Then, it reviews several statistical and computational methods for detecting and characterizing genes with effects that are dependent on other genes. The focus of this unit is genetic association studies of discrete and quantitative traits because most of the methods for detecting gene-gene interactions have been developed specifically for these study designs.

ATHENA: A knowledge-based hybrid backpropagation-grammatical evolution neural network algorithm for discovering epistasis among quantitative trait Loci

BACKGROUND

Growing interest and burgeoning technology for discovering genetic mechanisms that influence disease processes have ushered in a flood of genetic association studies over the last decade, yet little heritability in highly studied complex traits has been explained by genetic variation. Non-additive gene-gene interactions, which are not often explored, are thought to be one source of this "missing" heritability.

METHODS

Stochastic methods employing evolutionary algorithms have demonstrated promise in being able to detect and model gene-gene and gene-environment interactions that influence human traits. Here we demonstrate modifications to a neural network algorithm in ATHENA (the Analysis Tool for Heritable and Environmental Network Associations) resulting in clear performance improvements for discovering gene-gene interactions that influence human traits. We employed an alternative tree-based crossover, backpropagation for locally fitting neural network weights, and incorporation of domain knowledge obtainable from publicly accessible biological databases for initializing the search for gene-gene interactions. We tested these modifications in silico using simulated datasets.

RESULTS

We show that the alternative tree-based crossover modification resulted in a modest increase in the sensitivity of the ATHENA algorithm for discovering gene-gene interactions. The performance increase was highly statistically significant when backpropagation was used to locally fit NN weights. We also demonstrate that using domain knowledge to initialize the search for gene-gene interactions results in a large performance increase, especially when the search space is larger than the search coverage.

CONCLUSIONS

We show that a hybrid optimization procedure, alternative crossover strategies, and incorporation of domain knowledge from publicly available biological databases can result in marked increases in sensitivity and performance of the ATHENA algorithm for detecting and modelling gene-gene interactions that influence a complex human trait.

Analysis of gene-gene interactions

The goal of this unit is to introduce gene-gene interactions (epistasis) as a significant complicating factor in the search for disease susceptibility genes. This unit begins with an overview of gene-gene interactions and why they are likely to be common. Then, it reviews several statistical and computational methods for detecting and characterizing genes with effects that are dependent on other genes. The focus of this unit is genetic association studies of discrete and quantitative traits because most of the methods for detecting gene-gene interactions have been developed specifically for these study designs.

A coalescent simulation of marker selection strategy for candidate gene association studies

Recent efforts have focused on the challenges of finding alleles that contribute to health-related phenotypes in genome-wide association studies. However, in candidate gene studies, where the genomic region of interest is small and recombination is limited, factors that affect the ability to detect disease-susceptibility alleles remain poorly understood. In particular, it is unclear how varying the number of markers on a haplotype, the type of marker (e.g., single nucleotide polymorphism (SNP), short tandem repeat (STR)), including the causative site (cs) as a genetic marker, or population demographics influences the power to detect a candidate gene. We evaluated the power of association tests using coalescent-modeled computer simulations. Results show that an effective number of markers on a haplotype is dependent on whether the cs is included as a marker. When the analyses include the cs, highest power is achieved with a single-marker association test. However, when the cs is excluded from analyses, the addition of more nonfunctional SNPs on the haplotype increases power to a certain point under most scenarios. We find a rapidly expanding population always has lower power compared to a population of constant size; although utilizing markers with a frequency of at least 5% improves the chance of detecting an association. Comparing the mutational properties of a nonfunctional SNP versus an STR, multi-allelic STRs provide more or comparable power than a bi-allelic SNP unless SNP frequencies are constrained to 10% or more. Similarly, including an STR with SNPs on a haplotype improves power unless SNP frequencies are 5% or more.

Generalized genomic distance-based regression methodology for multilocus association analysis

Large-scale, multilocus genetic association studies require powerful and appropriate statistical-analysis tools that are designed to relate genotype and haplotype information to phenotypes of interest. Many analysis approaches consider relating allelic, haplotypic, or genotypic information to a trait through use of extensions of traditional analysis techniques, such as contingency-table analysis, regression methods, and analysis-of-variance techniques. In this work, we consider a complementary approach that involves the characterization and measurement of the similarity and dissimilarity of the allelic composition of a set of individuals' diploid genomes at multiple loci in the regions of interest. We describe a regression method that can be used to relate variation in the measure of genomic dissimilarity (or "distance") among a set of individuals to variation in their trait values. Weighting factors associated with functional or evolutionary conservation information of the loci can be used in the assessment of similarity. The proposed method is very flexible and is easily extended to complex multilocus-analysis settings involving covariates. In addition, the proposed method actually encompasses both single-locus and haplotype-phylogeny analysis methods, which are two of the most widely used approaches in genetic association analysis. We showcase the method with data described in the literature. Ultimately, our method is appropriate for high-dimensional genomic data and anticipates an era when cost-effective exhaustive DNA sequence data can be obtained for a large number of individuals, over and above genotype information focused on a few well-chosen loci.

The contribution of individual and pairwise combinations of SNPs in the APOA1 and APOC3 genes to interindividual HDL-C variability

Apolipoproteins (apo) A-I and C-III are components of high-density lipoprotein-cholesterol (HDL-C), a quantitative trait negatively correlated with risk of cardiovascular disease (CVD). We analyzed the contribution of individual and pairwise combinations of single nucleotide polymorphisms (SNPs) in the APOA1/APOC3 genes to HDL-C variability to evaluate (1) consistency of published single-SNP studies with our single-SNP analyses; (2) consistency of single-SNP and two-SNP phenotype-genotype relationships across race-, gender-, and geographical location-dependent contexts; and (3) the contribution of single SNPs and pairs of SNPs to variability beyond that explained by plasma apo A-I concentration. We analyzed 45 SNPs in 3,831 young African-American (N=1,858) and European-American (N=1,973) females and males ascertained by the Coronary Artery Risk Development in Young Adults (CARDIA) study. We found three SNPs that significantly impact HDL-C variability in both the literature and the CARDIA sample. Single-SNP analyses identified only one of five significant HDL-C SNP genotype relationships in the CARDIA study that was consistent across all race-, gender-, and geographical location-dependent contexts. The other four were consistent across geographical locations for a particular race-gender context. The portion of total phenotypic variance explained by single-SNP genotypes and genotypes defined by pairs of SNPs was less than 3%, an amount that is miniscule compared to the contribution explained by variability in plasma apo A-I concentration. Our findings illustrate the impact of context-dependence on SNP selection for prediction of CVD risk factor variability.

Evidence for consistent intragenic and intergenic interactions between SNP effects in the APOA1/C3/A4/A5 gene cluster

OBJECTIVE

Evaluate the consistency of the contribution of interactions between single nucleotide polymorphism (SNP) genotype effects to variation in measures of lipid metabolism across ethnic strata within gender.

METHODS AND RESULTS

We considered 80 SNPs within the apolipoprotein (APO) A1/C3/A4/A5 gene cluster using an over-parameterized general linear model to identify SNPs whose genotype effects combine non-additively to influence plasma levels of high density lipoprotein cholesterol (HDL-C), total cholesterol (TC) and triglycerides (TG) in a consistent manner across ethnic strata. We analyzed population-based samples of unrelated 18 to 30 year old African-Americans (n = 1,858) and European-Americans (n = 1,973) ascertained without regard to health at four field centers (Birmingham, Ala.; Chicago, Ill.; Minneapolis, Minn. and Oakland, Calif., USA) by the Coronary Artery Risk Development in Young Adults (CARDIA) study. To identify which SNP genotype effects combine non-additively we used a two-tier analysis strategy. We first required that pairs of SNPs show statistically significant non-additivity in both ethnic strata within a gender, where experiment-wise significance was evaluated using a permutation test to determine the probability of observing the number of tests significant in both ethnic strata by chance alone. Second, we required no significant evidence of heterogeneity of the relationship between the phenotype and the two SNP genotypes across ethnic strata and across field centers within each ethnic group. From this strategy we identified ten pairs of SNPs, involving thirteen SNPs, that displayed statistically significant non-additivity of SNP genotype effects on TC. Only one of these thirteen SNPs had statistically significant genotype effects that were consistent across samples.

CONCLUSION

Our analyses suggest that ignoring the contribution of interactions between SNP genotype effects when modeling multi-SNP genotype-phenotype relationships may result in an underestimate of the contribution of genetic variation to variation in quantitative cardiovascular disease (CVD) risk factor traits.

TagSNP analyses of the PON gene cluster: effects on PON1 activity, LDL oxidative susceptibility, and vascular disease

Paraoxonase 1 (PON1) activity is consistently predictive of vascular disease, although the genotype at four functional PON1 polymorphisms is not. To address this inconsistency, we investigated the role of all common PON1 genetic variability, as measured by tagging single-nucleotide polymorphisms (tagSNPs), in predicting PON1 activity for phenylacetate hydrolysis, LDL susceptibility to oxidation ex vivo, plasma homocysteine (Hcy) levels, and carotid artery disease (CAAD) status. The biological goal was to establish whether additional common genetic variation beyond consideration of the four known functional SNPs improves prediction of these phenotypes. PON2 and PON3 tagSNPs were secondarily evaluated. Expanded analysis of an additional 26 tagSNPs found evidence of previously undescribed common PON1 polymorphisms that affect PON1 activity independently of the four known functional SNPs. PON1 activity was not significantly correlated with LDL oxidative susceptibility, but genotypes at the PON1(-108) promoter polymorphism and several other PON1 SNPs were. Neither PON1 activity nor PON1 genotype was significantly correlated with plasma Hcy levels. This study revealed previously undetected common functional PON1 polymorphisms that explain 4% of PON1 activity and a high rate of recombination in PON1, but the sum of the common PON1 locus variation does not explain the relationship between PON1 activity and CAAD.

Haplotype Effect of the Matrix Metalloproteinase-1 Gene on Risk of Myocardial Infarction

Myocardial infarction (MI) is commonly caused by atherosclerotic plaque rupture following excessive degradation of collagen fibers in the atherosclerotic lesion. We investigated whether interindividual variability in risk of MI was related to polymorphisms in the gene encoding matrix metalloproteinase (MMP)-1, a key fibrillar collagen–degrading enzyme. Several single nucleotide polymorphisms in the MMP1 gene promoter were identified following sequencing DNA samples from 30 individuals. An analysis of the polymorphisms in a cohort of British whites with coronary atherosclerosis, including 639 patients with MI and 538 non-MI subjects, revealed a haplotype effect of the −519A>G and −340T>C polymorphisms on risk of MI, with the A −519 -C −340 and G −519 -T −340 haplotypes being protective (odds ratio=0.70 [0.57 to 0.86]; P =0.0007), whereas the G −519 -C −340 haplotype increased MI risk (odds ratio=1.94 [1.15 to 3.28]; P =0.013). This finding was replicated in a subsequent analysis of 387 Swedish MI patients and 387 healthy controls (odds ratio=0.70 [0.55 to 0.89], P =0.003, for A −519 -C −340 and G −519 -T −340 ; odds ratio=1.54 [0.97 to 2.46], P =0.07, for G −519 -C −340 ). In vitro assays showed that compared with the A −519 -T −340 haplotype, the A −519 -C −340 and G −519 -T −340 haplotypes had lower promoter activity, whereas the G −519 -C −340 haplotype had greater promoter strength, in driving gene expression in human macrophages. Haplotype-specific differences in MMP1 mRNA level in atherosclerotic tissues were also detected. The data indicate that MMP1 gene variation is a genetic factor contributing to interindividual differences in MI risk.