Family-based Association Test Research Articles

Investigation of Novel Variants in Familial Multiple Sclerosis Through Family-based Association Tests and Segregation Analyses (S17.007)

Investigation of Novel Variants in Familial Multiple Sclerosis Through Family-based Association Tests and Segregation Analyses (S17.007)

Association between IRF6, TP63, GREM1 Gene Polymorphisms and Non-Syndromic Orofacial Cleft Phenotypes in Vietnamese Population: A Case-Control and Family-Based Study.

This study aims to identify potential variants in the TP63-IRF6 pathway and GREM1 for the etiology of non-syndromic orofacial cleft (NSOFC) among the Vietnamese population. By collecting 527 case-parent trios and 527 control samples, we conducted a stratified analysis based on different NSOFC phenotypes, using allelic, dominant, recessive and over-dominant models for case-control analyses, and family-based association tests for case-parent trios. Haplotype and linkage disequilibrium analyses were also conducted. IRF6 rs2235375 showed a significant association with an increased risk for non-syndromic cleft lip and palate (NSCLP) and cleft lip with or without cleft palate (NSCL/P) in the G allele, with pallele values of 0.0018 and 0.0003, respectively. Due to the recessive model (p = 0.0011) for the NSCL/P group, the reduced frequency of the GG genotype of rs2235375 was associated with a protective effect against NSCL/P. Additionally, offspring who inherited the G allele at rs2235375 had a 1.34-fold increased risk of NSCL/P compared to the C allele holders. IRF6 rs846810 and a G-G haplotype at rs2235375-rs846810 of IRF6 impacted NSCL/P, with p-values of 0.0015 and 0.0003, respectively. In conclusion, our study provided additional evidence for the association of IRF6 rs2235375 with NSCLP and NSCL/P. We also identified IRF6 rs846810 as a novel marker associated with NSCL/P, and haplotypes G-G and C-A at rs2235375-rs846810 of IRF6 associated with NSOFC.

A novel rare variants association test for binary traits in family-based designs via copulas.

With the cost-effectiveness technology in whole-genome sequencing, more sophisticated statistical methods for testing genetic association with both rare and common variants are being investigated to identify the genetic variation between individuals. Several methods which group variants, also called gene-based approaches, are developed. For instance, advanced extensions of the sequence kernel association test, which is a widely used variant-set test, have been proposed for unrelated samples and extended for family data. Family data have been shown to be powerful when analyzing rare variants. However, most of such methods capture familial relatedness using a random effect component within the generalized linear mixed model framework. Therefore, there is a need to develop unified and flexible methods to study the association between a set of genetic variants and a trait, especially for a binary outcome. Copulas are multivariate distribution functions with uniform margins on the interval and they provide suitable models to capture familial dependence structure. In this work, we propose a flexible family-based association test for both rare and common variants in the presence of binary traits. The method, termed novel rare variant association test (NRVAT), uses a marginal logistic model and a Gaussian Copula. The latter is employed to model the dependence between relatives. An analytic score-type test is derived. Through simulations, we show that our method can achieve greater power than existing approaches. The proposed model is applied to investigate the association between schizophrenia and bipolar disorder in a family-based cohort consisting of 17 extended families from Eastern Quebec.

A Family-based Association Test of the VDR Gene in Proximal Spinal Muscular Atrophy

Different factors may be involved in the clinical heterogeneity of spinal muscular atrophy disease. The vitamin D receptor (VDR) might be a candidate gene for this disease. Our study aimed to assess the preferential transmission of VDR polymorphisms from parents to SMA children. We genotyped 261 subjects (87 SMA nuclear families) for VDR FokI, BsmI, ApaI, and TaqI polymorphisms. The transmission of the genetic marker was estimated with Plink and FBAT software. It detected a preferential transmission of the rs731236 and rs7975232 variants to SMA1 patients and of rs1544410 variants to SMA2 patients. The variants of rs2228570 were preferentially transmitted to parents of all SMA patients. Haplotype analysis identified that haplotypes C-C-G-A and T-A-A-G seem to be involved in the booth type of SMA whereas the impact of T-A-A-A seems to be limited only to SMA2. Strong linkage disequilibrium (LD) between rs7975232 and rs1544410 was detected in samples from parents. Even though we investigated a small number of nuclear families the results suggest a potential link between VDR polymorphisms and SMA disease.

KnockoffTrio: A knockoff framework for the identification of putative causal variants in genome-wide association studies with trio design

Family-based designs can eliminate confounding due to population substructure and can distinguish direct from indirect genetic effects, but these designs are underpowered due to limited sample sizes. Here, we propose KnockoffTrio, a statistical method to identify putative causal genetic variants for father-mother-child trio design built upon a recently developed knockoff framework in statistics. KnockoffTrio controls the false discovery rate (FDR) in the presence of arbitrary correlations among tests and is less conservative and thus more powerful than the conventional methods that control the family-wise error rate via Bonferroni correction. Furthermore, KnockoffTrio is not restricted to family-based association tests and can be used in conjunction with more powerful, potentially nonlinear models to improve the power of standard family-based tests. We show, using empirical simulations, that KnockoffTrio can prioritize causal variants over associations due to linkage disequilibrium and can provide protection against confounding due to population stratification. In applications to 14,200 trios from three study cohorts for autism spectrum disorders (ASDs), including AGP, SPARK, and SSC, we show that KnockoffTrio can identify multiple significant associations that are missed by conventional tests applied to the same data. In particular, we replicate known ASD association signals with variants in several genes such as MACROD2, NRXN1, PRKAR1B, CADM2, PCDH9, and DOCK4 and identify additional associations with variants in other genes including ARHGEF10, SLC28A1, ZNF589, and HINT1 at FDR 10%.

Family-based association tests for rare variants

Next-generation sequencing has revolutionized family-based association tests for rare variants. As the lower power of genome wide association study for detecting casual rare variants, methods aggregating effects of multiple variants have been proposed, such as burden tests and variance component tests. This paper summarizes the methods of rare variants association test that can be applied for family data, introduces their principles, characteristics and applicable conditions and discusses the shortcomings and the improvement of the present methods.

Evidence of a casual relationship between vitamin D deficiency and hypertension: a family-based study.

An association between vitamin D deficiency and hypertension has been observed in numerous studies. However, blood pressure improvements resulting from supplementation with vitamin D have been inconsistent. The causal relationship between vitamin D deficiency and hypertension is still unclear and was investigated in this family-based study. A total of 1370 individuals from both vitamin D deficiency and hypertension families were included. First, the heritability of vitamin D deficiency was estimated by the Falconer method. Second, SNPs (single nucleotide polymorphisms) of vitamin D metabolic and functional pathway genes associated with vitamin D deficiency were screened by a family-based association test, and the findings were further verified in nuclear families with vitamin D deficiency. Finally, a family-based association test was applied to investigate the association between selected SNPs associated with vitamin D deficiency and hypertension. The heritability of vitamin D deficiency was 50.4% in this family-based study. Allele C of rs3847987 was a risk factor for vitamin D deficiency (OR: 1.639, 95% CI: 1.170-2.297, P = 0.004). Furthermore, a family-based association of rs3847987 with hypertension was found in both additive and recessive models (P < 0.05). In addition, vitamin D deficiency was associated with hypertension (OR: 1.317, 95% CI: 1.022-1.698, P = 0.033). In conclusion, rs3847987 in the VDR gene was associated with both vitamin D deficiency and hypertension. Therefore, vitamin D deficiency may be a causal factor for hypertension.

Cluster of Differentiation 24 Polymorphism Has No Significant Association with Chronic Hepatitis B Virus Infection in the Chinese Han Population: A Family-Based Association Study.

ObjectiveStudies have shown that cluster of differentiation (CD) 24 gene polymorphism is associated with several diseases. Among these, chronic hepatitis B (CHB) infection has not been investigated. This study aimed to assess the function of CD24 in CHB.MethodsThe study included 478 cases of CHB and 318 cases without CHB from 230 families that underwent genotyping. Polymerase chain reaction-restriction fragment length polymorphism was performed to assess the single nucleotide polymorphism (SNP) P170 of the CD24 gene. The detected genotypes were TT, CT, and CC. Then, family based-association analysis was carried out to investigate the association between CD24 gene polymorphism and susceptibility to CHB.ResultsIn the 478 patients with CHB, the frequencies of CD24 P170 T and C alleles were 35.5% and 64.5%, respectively, and the frequencies of CD24 P170 CC, CT, and TT genotypes were 39.3%, 50.4% and 10.3%, respectively. In a CD24 single-locus analysis by a family-based association test of P170 polymorphisms, T and C were not significantly associated with CHB in the additive (Z = 0.169, P = 0.866; Z = −0.169, P = 0.866, respectively), dominant (Z = 0.522, P = 0.602; Z = 0.428, P = 0.669, respectively), or recessive (Z = −0.428, P = 0.669; Z = −0.522, P = 0.602, respectively) models. Transmission-disequilibrium (TD) and sib-transmission disequilibrium (STD) tests revealed no excess of T or C alleles from heterozygous parents to their children with the disease or higher frequencies of these alleles in patients compared with their normal siblings (χ2 = 0.06, P = 0.897).ConclusionThe study findings suggest that the SNP P170 of CD24 has no significant association with susceptibility to the HB virus and related phenotypes in Chinese patients.

Autoantibodies to Zinc Transporter 8 and SLC30A8 Genotype in Type 1 Diabetes Childhood: A Pioneering Study in North Africa.

Background Type 1 diabetes (T1D) occurs as a result of insulin deficiency due to destructive lesions of pancreatic β cells. In addition to classical autoantibodies (Abs) to islet cell antigens, antizinc transporter 8 Abs (ZnT8-Ab) have been recently described in T1D. Objective As no data on ZnT8-Ab in Tunisian patients has been reported, we aim to evaluate the relationships between ZnT8-Ab, ZnT8 coding gene (SLC30A8) promoter polymorphism, and T1D risk in newly diagnosed children. Methods ZnT8-Ab were measured in the serum of T1D newly affected children (n = 156) who were admitted to the pediatric department of the Hedi Chaker University Hospital of Sfax. Rs13266634 was genotyped in T1D children and 79 of their first-degree parents. The SPSS software was used to analyze the serological data. Allelic association analysis was conducted with family-based association tests implemented in the FBAT program v1.5.1. Results ZnT8-Ab was detected in 66/156 (42.3%) of T1D newly diagnosed children. Among them, 6 (9%) presented ZnT8-Ab as the only humoral marker. The inclusion of ZnT8-Ab increased the number of Ab-positive patients to 90% and reduced the negative ones by 27%. There was no evidence of any overtransmission of any allele of the rs13266634 C/T polymorphism from parents to affected T1D children, nor of any correlation with any clinical or serological parameter. After the T1D disease onset age adjustment, a significant association was observed with the C allele suggesting that it could have a susceptibility role. Conclusion ZnT8-Ab appears as a relevant diagnostic marker for T1D in Tunisian children, especially at the onset of the disease as teenagers.

Machine Learning Prediction of ADHD Severity: Association and Linkage to ADGRL3, DRD4, and SNAP25

Objective: To investigate whether single nucleotide polymorphisms (SNPs) in the ADGRL3, DRD4, and SNAP25 genes are associated with and predict ADHD severity in families from a Caribbean community. Method: ADHD severity was derived using latent class cluster analysis of DSM-IV symptomatology. Family-based association tests were conducted to detect associations between SNPs and ADHD severity latent phenotypes. Machine learning algorithms were used to build predictive models of ADHD severity based on demographic and genetic data. Results: Individuals with ADHD exhibited two seemingly independent latent class severity configurations. SNPs harbored in DRD4, SNAP25, and ADGRL3 showed evidence of linkage and association to symptoms severity and a potential pleiotropic effect on distinct domains of ADHD severity. Predictive models discriminate severe from non-severe ADHD in specific symptom domains. Conclusion: This study supports the role of DRD4, SNAP25, and ADGRL3 genes in outlining ADHD severity, and a new prediction framework with potential clinical use.

Adjusting for Treatment Effect when Estimating or Testing Genetic Effect is of Main Interest

It is known that “standard methods for estimating the causal effect of a time-varying treatment on the mean of a repeated measures outcome (for example, GEE regression) may be biased when there are time-dependent variables that are simultaneously confounders of the effect of interest and are predicted by previous treatment” (Hern´an et al. 2002). Inverse-probability of treatment weighted (IPTW) methods are developed in the literature of causal inference. In genetic studies, however, the main interest is to estimate or test the genetic effect rather than the treatment effect. In this work, we describe an IPTW method that provides unbiased estimate for the genetic effect, and discuss how to develop a family-based association test using IPTW for family-based studies. We apply the developed methods to systolic blood pressure data in Framingham Heart Study, where some subjects took antihypertensive treatment during the course of study.

Chromosomal regions strongly associated with waist circumference and body mass index in metabolic syndrome in a family-based study

Obesity is the most crucial phenotype in metabolic syndrome (MetS), and waist circumference (WC) and body mass index (BMI) are two common indexes to define obesity. It is an accepted fact that genetic and environmental interaction influence obesity and MetS. Microsatellites are a subcategory of tandem repeats with a length of 1 to 10 nucleotides. Tandem repeats make up repetitive genomic regions. Differences in the number of tandem repeats or their variation (alleles) result in microsatellite polymorphisms. Thus, we attempted to find microsatellite variation associated with WC and BMI in a family-based study. Twelve microsatellite markers were selected to investigate possible genes or chromosomal regions in 91 families with at least one affected MetS. The cut-off values for BMI and WC were considered 25 kg/m2 and 90 cm, respectively. In all members of the families, the strongest association was observed between the marker D11S1304 (allele 1) with both WC and BMI, independently, by the biallelic model in the family-based association test analysis (P < 0.05). Besides, when we compared high- and low-level groups in members with MetS, the markers D8S1743 and D11S1304 (allele 1) showed a strong association with WC (P = 0.0080) and BMI (P = 0.0074), respectively. When the simultaneous detection of the high WC and MetS status was used as a trait, the strongest association was observed with the marker D8S1743 (P = 0.0034). Moreover, when BMI with the high MetS status was used as a trait, the strongest association was observed with the marker D8S1743 (allele 4) (P = 0.0034). The obtained results showed a relationship between obesity and MetS with markers on the selected regions on chromosomes 8 and 11, and to a lesser degree, on chromosome 12.

Sex-dependent complex association of TPH2 with multiple dimensions of ADHD

BackgroundTryptophan hydroxylase 2 (TPH2) is a key enzyme in the biosynthesis of serotonin in the brain. This study aims to investigate the role of a functional variant in TPH2 (rs17110747) in the pathophysiology of ADHD. This variant has been implicated in mood disorders in recent meta-analysis. This study uses a comprehensive approach that combines association testing and pharmaco-dynamic evaluation of behaviour, in a large sample of children with ADHD (n = 570). MethodsThe association between various ADHD relevant traits and rs17110747 was analyzed using family-based association tests (FBAT). Children were assessed by parents, teachers and research staff under three experimental conditions (EC): baseline, placebo, and methylphenidate using a double-blind placebo-controlled crossover trial. OutcomesFBAT analysis conducted in a sample stratified based on sex of the proband, showed that there was a highly significant overtransmission of the G allele from parents to affected girls. In addition, significant association with several behavioral and cognitive dimensions of ADHD was observed only when the proband was female. Further, girls with the G/G genotype (rs17110747) had greater response to placebo when evaluated by parents. InterpretationThese results suggest that there may be a complex association of TPH2 in the etiology of ADHD, with a sex-specific effect.

Whole-exome sequencing identifies multiple pathogenic variants in a large South Indian family with primary open-angle glaucoma.

Purpose:To identify the pathogenic variants associated with primary open-angle glaucoma (POAG) using whole-exome sequencing (WES) data of a large South Indian family.Methods:We recruited a large five-generation South Indian family (n = 84) with a positive family history of POAG (n = 19). All study participants had a comprehensive ocular evaluation. We performed WES for 16 samples (nine POAG and seven unaffected controls) since Sanger sequencing of the POAG candidate genes (MYOC, OPTN, and TBK1) showed no genetic variation. We used an in-house pipeline for prioritizing the pathogenic variants based on their segregation among the POAG individual.Results:We identified one novel and five low-frequency pathogenic variants with consistent co-segregation in all affected individuals. The variant c.G3719A in RPGR-interacting domain of RPGRIP1 that segregated heterozygously with the six POAG cases is distinct from variants causing photoreceptor dystrophies, reported affecting the RPGR protein complex signaling in primary cilia. The cilia in trabecular meshwork (TM) cells has been reported to mediate the intraocular pressure (IOP) sensation. Furthermore, we identified a novel c.A1295G variant in Rho guanine nucleotide exchange factors Gene 40 (ARHGEF40) and a likely pathogenic variant in the RPGR gene, suggesting that they may alter the RhoA activity essential for IOP regulation.Conclusion:Our study supports that low-frequency pathogenic variants in multiple genes and pathways probably affect Primary Open Angle Glaucoma’s pathogenesis in the large South Indian family. Furthermore, it requires larger case-controls to perform family-based association tests and to strengthen our analysis.

MF-TOWmuT: Testing an optimally weighted combination of common and rare variants with multiple traits using family data.

With rapid advancements of sequencing technologies and accumulations of electronic health records, a large number of genetic variants and multiple correlated human complex traits have become available in many genetic association studies. Thus, it becomes necessary and important to develop new methods that can jointly analyze the association between multiple genetic variants and multiple traits. Compared with methods that only use a single marker or trait, the joint analysis of multiple genetic variants and multiple traits is more powerful since such an analysis can fully incorporate the correlation structure of genetic variants and/or traits and their mutual dependence patterns. However, most of existing methods that simultaneously analyze multiple genetic variants and multiple traits are only applicable to unrelated samples. We develop a new method called MF-TOWmuT to detect association of multiple phenotypes and multiple genetic variants in a genomic region with family samples. MF-TOWmuT is based on an optimally weighted combination of variants. Our method can be applied to both rare and common variants and both qualitative and quantitative traits. Our simulation results show that (1) the type I error of MF-TOWmuT is preserved; (2) MF-TOWmuT outperforms two existing methods such as Multiple Family-based Quasi-Likelihood Score Test and Multivariate Family-based Rare Variant Association Testin terms of power. We also illustrate the usefulness of MF-TOWmuT by analyzing genotypic and phenotipic data from the Genetics of Kidneys in Diabetes study. R program is available at https://github.com/gaochengPRC/MF-TOWmuT.

Genetic Variants of the MTMR9 Gene Are Associated with Nonspecific Intellectual Disability: A Family-Based Association Study.

Background: Mutations within the myotubularin-related protein 9 gene (MTMR9) have been identified in several families with nonsyndromic intellectual disability (NSID), a generalized neurodevelopmental disorder; however, the relationship between MTMR9 and NSID needs to be verified using a larger sample size. Aim: To explore whether genetic variants in the MTMR9 gene are linked to susceptibility of NSID among the Chinese population. Materials and Methods: Seven single nucleotide polymorphisms (SNPs) of the MTMR9 gene (rs4559208, rs3824211, rs2164272, rs2164273, rs1897951, rs6991606, and rs7815802) were analyzed using family-based association testing among 258 Han Chinese NSID families. Results: Three SNPs of MTMR9 were significantly associated with NSID (z = 2.152, p = 0.031 for rs4559208; z = 2.403, p = 0.016 for rs2164273; and z = 2.758, p = 0.006 for rs7815802). Three alleles of these SNPs were more likely to be transferred from the carrier parents to the affected offspring. Haplotypes constructed using these SNPs also showed a similar transmitting trend (z = 2.505, p = 0.012, χ2(3) = 8.835, and global p = 0.032). Carriers with the G-G-C haplotype showed a higher risk of NSID (odds ratio = 1.46, 95% confidence interval [1.01-2.09], p = 0.04) than others. In silico functional predictions supported an etiological role for these three SNPs in NSID biology. Conclusions: This study provides additional insights into the association of NSID with specific alleles, and haplotypes within the MTMR9 gene. Genotypic analyses of the MTMR9 gene should be considered for patients presenting with NSID of unknown etiology.

A novel locus for exertional dyspnoea in childhood asthma.

Most children diagnosed with asthma have respiratory symptoms such as cough, dyspnoea and wheezing, which are also important markers of overall respiratory function. A decade of genome-wide association studies (GWAS) have investigated genetic susceptibility to asthma itself, but few have focused on important respiratory symptoms that characterise childhood asthma.Using whole-genome sequencing (WGS) data for 894 asthmatic trios from a Costa Rican cohort, we performed family-based association tests (FBATs) to assess the association between genetic variants and multiple asthma-relevant respiratory phenotypes: cough, phlegm, wheezing, exertional dyspnoea and exertional chest tightness. We tested whether genome-wide significant associations were replicated in two additional studies: 1) 286 asthmatic trios from the Childhood Asthma Management Program (CAMP), and 2) 2691 African American current or former smokers from the COPDGene study.In the 894 Costa Rican trios, we identified a genome-wide significant association (p=2.16×10-9) between exertional dyspnoea and the single nucleotide polymorphism (SNP) rs10165869, located on chromosome 2q37.3, that was replicated in the CAMP cohort (p=0.023) with the same direction of association (combined p=3.28×10-10). This association was not found in the African American participants from COPDGene. We also found suggestive evidence for an association between SNP rs10165869 and the atypical chemokine receptor 3 (ACKR3).Our finding encourages the secondary association analysis of a wider range of phenotypes that characterise respiratory symptoms in other airway diseases/studies.

A family-based genome-wide association study of recurrent aphthous stomatitis.

The aetiology of recurrent aphthous stomatitis (RAS) remains unknown. Individuals may share features of genetic susceptibility, and there may also be a hereditary component. The aim was to identify patterns of association and segregation for genetic variants and to identify the genes and signalling pathways that determine the risk of developing RAS, through a family-based genome-wide association study (GWAS). DNA was extracted from buccal swabs of 91 individuals in 16 families and analysed in an Illumina core exome single nucleotide polymorphism (SNP) array. A family-based association test (dFAM) was used to derive SNP association values across all chromosomes. None of the final 288,452 SNPs reached the genome-wide significant threshold of 5×10-8 . The most significant pathways were the Ras and PI3K-Akt signalling pathways, pathways in cancer, circadian entrainment and the Rap 1 signalling pathway. This confirms that RAS is not monogenic but results as a consequence of interactions between multiple host genes and possibly also environmental factors. The present approach provides novel insights into the mechanisms underlying RAS and raises the possibility of identifying individuals at risk of acquiring this condition.

Identification of Novel Alzheimer\u2019s Disease Loci Using Sex-Specific Family-Based Association Analysis of Whole-Genome Sequence Data

With the advent of whole genome-sequencing (WGS) studies, family-based designs enable sex-specific analysis approaches that can be applied to only affected individuals; tests using family-based designs are attractive because they are completely robust against the effects of population substructure. These advantages make family-based association tests (FBATs) that use siblings as well as parents especially suited for the analysis of late-onset diseases such as Alzheimer’s Disease (AD). However, the application of FBATs to assess sex-specific effects can require additional filtering steps, as sensitivity to sequencing errors is amplified in this type of analysis. Here, we illustrate the implementation of robust analysis approaches and additional filtering steps that can minimize the chances of false positive-findings due to sex-specific sequencing errors. We apply this approach to two family-based AD datasets and identify four novel loci (GRID1, RIOK3, MCPH1, ZBTB7C) showing sex-specific association with AD risk. Following stringent quality control filtering, the strongest candidate is ZBTB7C (Pinter = 1.83 × 10−7), in which the minor allele of rs1944572 confers increased risk for AD in females and protection in males. ZBTB7C encodes the Zinc Finger and BTB Domain Containing 7C, a transcriptional repressor of membrane metalloproteases (MMP). Members of this MMP family were implicated in AD neuropathology.

Common genetic variation in the Angelman syndrome imprinting centre affects the imprinting of chromosome 15

Angelman syndrome (AS) is a rare neurogenetic imprinting disorder caused by the loss of function of UBE3A. In ~3–5% of AS patients, the disease is due to an imprinting defect (ID). These patients lack DNA methylation of the maternal SNRPN promotor so that a large SNRPN sense/UBE3A antisense transcript (SNHG14) is expressed, which silences UBE3A. In very rare cases, the ID is caused by a deletion of the AS imprinting centre (AS-IC). To search for sequence alterations, we sequenced this region in 168 patients without an AS-IC deletion, but did not detect any sequence alteration. However, the AS-IC harbours six common variants (five single nucleotide variants and one TATG insertion/deletion variant), which constitute five common haplotypes. To determine if any of these haplotypes is associated with an increased risk for an ID, we investigated 119 informative AS-ID trios with the transmission disequilibrium test, which is a family-based association test that measures the over-transmission of an allele or haplotype from heterozygous parents to affected offspring. By this we observed maternal over-transmission of haplotype H-AS3 (p = 0.0073). Interestingly, H-AS3 is the only haplotype that includes the TATG deletion allele. We conclude that this haplotype and possibly the TATG deletion, which removes a SOX2 binding site, increases the risk for a maternal ID and AS. Our data strengthen the notion that the AS-IC is important for establishing and/or maintaining DNA methylation at the SNRPN promotor and show that common genetic variation can affect genomic imprinting.