Exploring causality via identification of SNPs or haplotypes responsible for a linkage signal

SNP prioritization in targeted sequencing data associated with humoral immune responses in chicken

Our study aimed to identify single nucleotide polymorphisms (SNPs) with a significant impact on the innate immunity represented by antibody response against lipopolysaccharide (LPS) and lipoteichoid acid (LTA) and the adaptive immune response represented toward keyhole limpet hemocyanin (KLH) using the SNP prioritization method. Data set consisted of 288 F2 experimental individuals, created by crossing Green-legged Partridgelike and White Leghorn. The analyzed SNPs were located within 24 short genomic regions of GGA1, GGA2, GGA3, GGA4, GGA9, GGA10, GGA14, GGA18, and GGZ, pre-targeted based on literature references and database information. For the specific antibody response toward KLH at d 0 the most highly prioritized SNP for additive and dominance effects were located on GGA2 in the 3’UTR of MYD88. For the response at d 7, the most highly prioritized SNP pointed at the 3’UTR of MYD88, but potential causal additive variants were located within ADIPOQ and one in PROCR. The highest priority for additive and dominance effects in the antibody response toward lipoteichoic acid at d 0 was attributed to the same SNP, located on GGA2 in the 3’UTR region of MYD88. Two SNPs among the top-10 for additive effect were located in the exon of NOCT. SNPs selected for their additive effect on antibody response toward lipopolysaccharide at d 0 marked 3 genes – NOCT, MYD88, and SNX8, while SNPs selected for their dominance effect marked – NOCT, ADIPOQ, and MYD88. The top-10 variants identified in our study were located in different functional parts of the genome. In the context of causality three groups can be distinguished: variants located in exons of protein coding genes (ADIPOQ, NOCT, PROCR, SNX8), variants within exons of non-coding transcripts, and variants located in genes’ UTR regions. Variants from the first group influence protein structure and variants from both latter groups’ exhibit regulatory roles on DNA (UTR) or RNA (lncRNA).

DNA methylation and breast cancer-associated variants

BACKGROUND

A breast cancer polygenic risk score (PRS) comprising 313 common variants reliably predicts disease risk. We examined possible relationships between genetic variation, regulation, and expression to clarify the molecular alterations associated with these variants.

METHODS

Genome-wide methylomic variation was quantified (MethylationEPIC) in Asian breast cancer patients (1152 buffy coats from peripheral whole blood). DNA methylation (DNAm) quantitative trait loci (mQTL) mapping was performed for 235 of the 313 variants with minor allele frequencies > 5%. Stability of identified mQTLs (p < 5e-8) across lifetime was examined using a public mQTL database. Identified mQTLs were also mapped to expression quantitative trait loci (eQTLs) in the Genotype-Tissue Expression Project and the eQTLGen Consortium.

RESULTS

Breast cancer PRS was not associated with DNAm. A higher proportion of significant cis-mQTLs were observed. Of 822 significant cis-mQTLs (179 unique variants) identified in our dataset, 141 (59 unique variants) were significant (p < 5e-8) in a public mQTL database. Eighty-six percent (121/141) of the matched mQTLs were consistent at multiple time points (birth, childhood, adolescence, pregnancy, middle age, post-diagnosis, or treatment). Ninety-three variants associated with DNAm were also cis-eQTLs (35 variants not genome-wide significant). Multiple loci in the breast cancer PRS are associated with DNAm, contributing to the polygenic nature of the disease. These mQTLs are mostly stable over time.

CONCLUSIONS

Consistent results from DNAm and expression data may reveal new candidate genes not previously associated with breast cancer.

Genome-wide linkage on chromosome 10q26 for a dimensional scale of major depression

Major depressive disorder (MDD) is a common and potentially life-threatening mood disorder. Identifying genetic markers for depression might provide reliable indicators of depression risk, which would, in turn, substantially improve detection, enabling earlier and more effective treatment. The aim of this study was to identify rare variants for depression, modeled as a continuous trait, using linkage and post-hoc association analysis. The sample comprised 1221 Mexican-American individuals from extended pedigrees. A single dimensional scale of MDD was derived using confirmatory factor analysis applied to all items from the Past Major Depressive Episode section of the Mini-International Neuropsychiatric Interview. Scores on this scale of depression were subjected to linkage analysis followed by QTL region-specific association analysis. Linkage analysis revealed a single genome-wide significant QTL (LOD=3.43) on 10q26.13, QTL-specific association analysis conducted in the entire sample revealed a suggestive variant within an intron of the gene LHPP (rs11245316, p=7.8×10(-04); LD-adjusted Bonferroni-corrected p=8.6×10(-05)). This region of the genome has previously been implicated in the etiology of MDD; the present study extends our understanding of the involvement of this region by highlighting a putative gene of interest (LHPP).

Combined linkage and family-based association analysis improves candidate gene detection in Genetic Analysis Workshop 18 simulation data

Because the genotype-phenotype correlation information is investigated differently by linkage and association analyses, various efforts have been made to model linkage and association jointly. However, joint modeling methods are usually computationally intensive; hence they cannot currently accommodate large pedigrees with dense markers. This article proposes a simple method to combine the linkage and association evidence with the aim of improving the detection power of disease susceptibility genes. Our detection power comparisons show that the combined linkage-association p values can improve remarkably the causal gene detection power in Genetic Analysis Workshop 18 simulation data.

Complex pedigrees in the sequencing era: to track transmissions or decorrelate?

Next-generation sequencing (NGS) studies are becoming commonplace, and the NGS field is continuing to develop rapidly. Analytic methods aimed at testing for the various roles that genetic susceptibility plays in disease are also rapidly being developed and optimized. Studies that incorporate large, complex pedigrees are of particular importance because they provide detailed information about inheritance patterns and can be analyzed in a variety of complementary ways. The nine contributions from our Genetic Analysis Workshop 18 working group on family-based tests of association for rare variants using simulated data examined analytic methods for testing genetic association using whole-genome sequencing data from 20 large pedigrees with 200 phenotype simulation replicates. What distinguishes the approaches explored is how the complexities of analyzing familial genetic data were handled. Here, we explore the methods that either harness inheritance patterns and transmission information or attempt to adjust for the correlation between family members in order to utilize computationally and conceptually simpler statistical testing procedures. Although directly comparing these two classes of approaches across contributions is difficult, we note that the two classes balance robustness to population stratification and computational complexity (the transmission-based approaches) with simplicity and increased power, assuming no population stratification or proper adjustment for it (decorrelation approaches).

A novel method, the Variant Impact On Linkage Effect Test (VIOLET), leads to improved identification of causal variants in linkage regions

The Human Genome Project was expected to individualize medicine by rapidly advancing knowledge of common complex disease through discovery of disease-causing genetic variants. However, this has proved challenging. Although linkage analysis has identified replicated chromosomal regions, subsequent detection of causal variants for complex traits has been limited. One explanation for this difficulty is that utilization of association to follow up linkage is problematic given that linkage and association are not required to co-occur. Indeed, co-occurrence is likely to occur only in special circumstances, such as Mendelian inheritance, but cannot be universally expected. To overcome this problem, we propose a novel method, the Variant Impact On Linkage Effect Test (VIOLET), which differs from other quantitative methods in that it is designed to follow up linkage by identifying variants that influence the variance explained by a quantitative trait locus. VIOLET's performance was compared with measured genotype and combined linkage association in two data sets with quantitative traits. Using simulated data, VIOLET had high power to detect the causal variant and reduced false positives compared with standard methods. Using real data, VIOLET identified a single variant, which explained 24% of linkage; this variant exhibited only nominal association (P=0.04) using measured genotype and was not identified by combined linkage association. These results demonstrate that VIOLET is highly specific while retaining low false-negative results. In summary, VIOLET overcomes a barrier to gene discovery and thus may be broadly applicable to identify underlying genetic etiology for traits exhibiting linkage.

A method to detect single-nucleotide polymorphisms accounting for a linkage signal using covariate-based affected relative pair linkage analysis

We evaluate an approach to detect single-nucleotide polymorphisms (SNPs) that account for a linkage signal with covariate-based affected relative pair linkage analysis in a conditional-logistic model framework using all 200 replicates of the Genetic Analysis Workshop 17 family data set. We begin by combining the multiple known covariate values into a single variable, a propensity score. We also use each SNP as a covariate, using an additive coding based on the number of minor alleles. We evaluate the distribution of the difference between LOD scores with the propensity score covariate only and LOD scores with the propensity score covariate and a SNP covariate. The inclusion of causal SNPs in causal genes increases LOD scores more than the inclusion of noncausal SNPs either within causal genes or outside causal genes. We compare the results from this method to results from a family-based association analysis and conclude that it is possible to identify SNPs that account for the linkage signals from genes using a SNP-covariate-based affected relative pair linkage approach.

Incorporation of covariates in simultaneous localization of two linked loci using affected relative pairs

Background

Many dichotomous traits for complex diseases are often involved more than one locus and/or associated with quantitative biomarkers or environmental factors. Incorporating these quantitative variables into linkage analysis as well as localizing two linked disease loci simultaneously could therefore improve the efficiency in mapping genes. We extended the robust multipoint Identity-by-Descent (IBD) approach with incorporation of covariates developed previously to simultaneously estimate two linked loci using different types of affected relative pairs (ARPs).

Results

We showed that the efficiency was enhanced by incorporating a quantitative covariate parametrically or non-parametrically while localizing two disease loci using ARPs. In addition to its help in identifying factors associated with the disease and in improving the efficiency in estimating disease loci, this extension also allows investigators to account for heterogeneity in risk-ratios for different ARPs. Data released from the collaborative study on the genetics of alcoholism (COGA) for Genetic Analysis Workshop 14 (GAW 14) were used to illustrate the application of this extended method.

Conclusions

The simulation studies and example illustrated that the efficiency in estimating disease loci was demonstratively enhanced by incorporating a quantitative covariate and by using all relative pairs while mapping two linked loci simultaneously.

Evaluation of approaches to identify associated SNPs that explain the linkage evidence in nuclear families with affected siblings

Linkage analysis is often followed by association mapping to localize disease variants. In this paper, we evaluate approaches to determine how much of the observed linkage evidence, namely the identity-by-descent (IBD) sharing at the linkage peak, is explained by associated SNPs. We study several methods: Homozygote Sharing Tests (HST), Genotype Identity-by-Descent Sharing Test (GIST), and a permutation approach. We also propose a new approach, HSTMLB, combining HST and the Maximum Likelihood Binomial (MLB) linkage statistic. These methods can identify SNPs partially explaining the linkage peak, but only HST and HSTMLB can identify SNPs that do not fully explain the linkage evidence and be applied to multiple-SNPs. We contrast these methods with the association tests implemented in the software LAMP. In our simulations, GIST is more powerful at finding SNPs that partially explain the linkage peak, while HST and HSTMLB are equally powerful at identifying SNPs that do not fully explain the linkage peak. When applied to the North American Rheumatoid Arthritis Consortium data, HST and HSTMLB identify marker pairs that may fully explain the linkage peak on chromosome 6. In conclusion, HST and HSTMLB provide simple and flexible tools to identify SNPs that explain the IBD sharing at the linkage peak.

A polygenic model for integration of linkage and pathway information

We introduce an approximate model for linkage curves which accommodates the polygenic structure of complex diseases and accounts for the simultaneous action of closely located genes. The model is extended so that information on biological pathways can be integrated. Using data on rheumatoid arthritis, we describe some of the many applications which the model allows: it can be used to test for residual linkage in the presence of already established loci, to derive a global test for linkage, to test for the relevance of a gene list in terms of linkage and to help in candidate gene prioritization by integration of gene-pathway annotation data.

A composite-likelihood approach for identifying polymorphisms that are potentially directly associated with disease

If a linkage signal can be fully accounted for by the association of a particular polymorphism with the disease, this polymorphism may be the sole causal variant in the region. On the other hand, if the linkage signal exceeds that explained by the association, different or additional directly associated loci must exist in the region. Several methods have been proposed for testing the hypothesis that association with a particular candidate single-nucleotide polymorphism (SNP) can explain an observed linkage signal. When several candidate SNPs exist, all of the existing methods test the hypothesis for each candidate SNP separately, by fitting the appropriate model for each individual candidate SNP. Here we propose a method that combines analyses of two or more candidate SNPs using a composite-likelihood approach. We use simulations to demonstrate that the proposed method can lead to substantial power increases over the earlier single SNP analyses.

Using linkage and association to identify and model genetic effects: summary of GAW15 Group 4

Group 4 at Genetic Analysis Workshop 15 focused on methods that exploited both linkage and association information to map disease loci. All contributions considered the dichotomous trait of rheumatoid arthritis, using either affected sibpairs and/or unrelated controls. While one contribution investigated linkage and association approaches separately in genome-wide analyses, the remaining others focused on joint linkage and association methods in specific genomic regions. The latter contributions proposed new methods and/or examined existing methods that addressed whether one or more polymorphisms partially or fully explained a linkage signal, particularly the methods proposed by Li et al. that are implemented in the computer program Linkage and Association Modeling in Pedigrees (LAMP). Using simulated SNP data under linkage peaks, several contributions found that existing family-based association approaches such as those of Martin et al. and Lake et al. had power similar to LAMP and to several methods proposed by the contributors for testing that a single nucleotide polymorphism partially explains a linkage peak. In evaluating methods for identifying if a polymorphism or a set of polymorphisms fully accounted for a linkage signal, several contributions found that it was important to understand that these methods may be subject to low power in some situations and thus, a non-significant result was not necessarily indicative of the polymorphism(s) being fully responsible for the linkage signal. Finally, modeling the disease using association evidence conditional on linkage may improve understanding of the etiology of disease.

Joint linkage and association analysis for identification of potentially causal polymorphisms in GAW15 data

In a small chromosomal region, a number of polymorphisms may be both linked to and associated with a disease. Potentially directly associated causal loci may be distinguished from indirectly associated loci by determining which associations can explain the observed linkage signal. We apply methods for testing whether association with a particular polymorphism or haplotype can explain an observed linkage signal to the Genetic Analysis Workshop 15 simulated (Problem 3) data, to try to identify potentially causal polymorphisms. We compare the power of several methods for testing the null hypothesis that association with a particular variant can explain the observed linkage signal, and discuss scenarios under which the various methods may be advantageous.