Dataset Information

A linkage disequilibrium-based approach to selecting disease-associated rare variants.

ABSTRACT: Rare variants have increasingly been cited as major contributors in the disease etiology of several complex disorders. Recently, several approaches have been proposed for analyzing the association of rare variants with disease. These approaches include collapsing rare variants, summing rare variant test statistics within a particular locus to improve power, and selecting a subset of rare variants for association testing, e.g., the step-up approach. We found that (a) if the variants being pooled are in linkage disequilibrium, the standard step-up method of selecting the best subset of variants results in loss of power compared to a model that pools all rare variants and (b) if the variants are in linkage equilibrium, performing a subset selection using step-based selection methods results in a gain of power of association compared to a model that pools all rare variants. Therefore, we propose an approach to selecting the best subset of variants to include in the model that is based on the linkage disequilibrium pattern among the rare variants. The proposed linkage disequilibrium-based variant selection model is flexible and borrows strength from the model that pools all rare variants when the rare variants are in linkage disequilibrium and from step-based selection methods when the variants are in linkage equilibrium. We performed simulations under three different realistic scenarios based on: (1) the HapMap3 dataset of the DRD2 gene, and CHRNA3/A5/B4 gene cluster (2) the block structure of linkage disequilibrium, and (3) linkage equilibrium. We proposed a permutation-based approach to control the type 1 error rate. The power comparisons after controlling the type 1 error show that the proposed linkage disequilibrium-based subset selection approach is an attractive alternative method for subset selection of rare variants.

SUBMITTER: Talluri R

PROVIDER: S-EPMC3708889 | biostudies-literature | 2013

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

A linkage disequilibrium-based approach to selecting disease-associated rare variants.

Talluri Rajesh R Shete Sanjay S

PloS one 20130711 7

Rare variants have increasingly been cited as major contributors in the disease etiology of several complex disorders. Recently, several approaches have been proposed for analyzing the association of rare variants with disease. These approaches include collapsing rare variants, summing rare variant test statistics within a particular locus to improve power, and selecting a subset of rare variants for association testing, e.g., the step-up approach. We found that (a) if the variants being pooled ...[more]

PMID: 23874919

Similar Datasets

Project description:The Lundehund is an old dog breed with remarkable anatomical features including polydactyly in all four limbs and extraordinary flexibility of the spine. We genotyped 28 Lundehund using the canine Illumina high density beadchip to estimate the effective population size (Ne) and inbreeding coefficients as well as to identify potential regions of positive selection. The decay of linkage disequilibrium was slow with r2=0.95 in 50 kb distance. The last 7-200 generations ago, Ne was at 10-13. An increase of Ne was noted in the very recent generations with a peak value of 19 for Ne at generation 4. The FROH estimated for 50-, 65- and 358-SNP windows were 0.87, 087 and 0.81, respectively. The most likely estimates for FROH after removing identical-by-state segments due to linkage disequilibria were at 0.80-0.81. The extreme loss of heterozygosity has been accumulated through continued inbreeding over 200 generations within a probably closed population with a small effective population size. The mean inbreeding coefficient based on pedigree data for the last 11 generations (FPed=0.10) was strongly biased downwards due to the unknown coancestry of the founders in this pedigree data. The long-range haplotype test identified regions with genes involved in processes of immunity, olfaction, woundhealing and neuronal development as potential targets of selection. The genes QSOX2, BMPR1B and PRRX2 as well as MYOM1 are candidates for selection on the Lundehund characteristics small body size, increased number of digits per paw and extraordinary mobility, respectively. The objectives of the present study were to genotype 28 Lundehund using the canine Illumina high density beadchip (Illumina, San Diego, CA, USA) to estimate the effective population size (Ne) from data on linkage disequilibria (LD). We identified runs of homozygosity (ROH) as regions with a local loss of genetic variation. Inbreeding coefficients were calculated on the basis of pedigree data (FPed) and genotype information (FROH, FIS) to compare the results of these different methods for this population. We applied the long-range haplotype test to search for potential selective sweeps.

Project description:BACKGROUND: Hirschsprung disease (HSCR) is a complex disorder with traditional germline mutations in RET in up to 30% of familial cases and in 3% of sporadic cases in a population-based series. We have previously demonstrated that an ancestral haplotype at the 5' end of RET (haplotype 0) was strongly associated with a large subset of isolated HSCR cases and that a putative low penetrance susceptibility locus was encompassed within this ancestral haplotype, anchored by exon 2 SNP A45A. OBJECTIVE: To determine the 5' extent of the HSCR-associated ancestral haplotype by defining the linkage disequilibrium breakpoint in search for the low penetrance susceptibility locus. METHODS: Systematic screening of the region upstream of the anchoring A45A SNP, comprising RET intron 1, exon 1, and promoter in 117 population-based HSCR cases and 100 controls. Dual luciferase assay to determine differential activities between SNP combinations near a transcription start site. RESULTS: New SNP's were found which formed upstream haplotypes, anchored by A45A, in linkage disequilibrium with HSCR (2 = 76.96, p<0.00000001). Linkage disequilibrium appeared to break at the -1249C/T SNP. Further, the HSCR-associated genotype (00) was found in >60% of HSCR but only 2% of controls. Because only 2 variants, -200A>G and -196C>A, lie within the promoter region and are in proximity to the transcriptional start site (at -195), we modelled these combinations into constructs for luciferase reporter assay. The HSCR-associated SNP combination showed the lowest activity and the control-associated combination, the highest. CONCLUSIONS: Our observations seem to discard the existence of a HSCR-causing mutation as it is conceived in the traditional sense, but strengthen the idea of a specific combination of variants conferring susceptibility to the disease in a low penetrance fashion. The data derived from our functional "in vitro" studies would suggest that the HSCR-associated haplotype 0 may result in a lower level of expression of the RET gene [corrected]

Dataset Information

A linkage disequilibrium-based approach to selecting disease-associated rare variants.

Publications

A linkage disequilibrium-based approach to selecting disease-associated rare variants.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets