Family-based Association Analysis Research Articles

Single nucleotide polymorphisms in JAZF1 and BCL11A gene are nominally associated with type 2 diabetes in African-American families from the GENNID study.

Prior type 2 diabetes (T2D) genome-wide association studies (GWASs) have generated a list of well-replicated susceptibility loci in populations of European and Asian ancestry. In order to validate the trans-ethnic contribution of these GWAS SNPs, we performed a family-based association analysis of 32 selected GWAS SNPs in a cohort of 1496 African American (AA) subjects from the GENNID study. Functional roles of these SNPs were evaluated by screening cis-eQTLs in transformed lymphoblast cell lines available for a sub-group of GENNID families from Arkansas. Only three of the 32 GWAS-derived SNPs showed nominally significant association with T2D in our AA cohort. Among the replicated SNPs rs864745 in JAZF1 and rs10490072 in BCL11A gene (p= 0.006 and 0.03, respectively, after adjustment for BMI) were within the 1-lod drop support interval of T2D linkage peaks reported in these families. Genotyping of 19 Tag-SNPs in these two loci revealed no further common SNPs or haplotypes that may be a stronger predictor of T2D susceptibility than the index SNPs. Six T2D GWAS SNPs (rs6698181, rs9472138, rs730497, rs10811661, rs11037909, and rs1153188) were associated with nearby transcript expression in transformed lymphoblast cell lines of GENNID AA subjects. Thus, our study indicates a nominal role for JAZF1 and BCL11A variants in T2D susceptibility in AAs and suggested little overlap in known susceptibility to T2D between European and African derived populations when considering GWAS SNPs alone.

On the Validity of the Likelihood Ratio Test and Consistency of Resulting Parameter Estimates in Joint Linkage and Linkage Disequilibrium Analysis under Improperly Specified Parametric Models

It has been shown that parametric analysis of linkage disequilibrium conditional on linkage using an overly deterministic model can be optimal for family-based association analysis. However, if one applies this strategy carelessly, there is a risk of false inference. We analyse properties of such likelihood ratio tests when the assumed disease mode of inheritance is inaccurate. Under some conditions, problems result if one is not careful to consider what null hypothesis is being tested. We show that: (a) tests for which the null hypothesis assumes the absence of both linkage and association are independent of the true mode of inheritance; (b) likelihood ratio tests assuming either linkage or association under the null hypothesis may depend on the true mode of inheritance, leading to inconsistent parameter estimates, in particular under extremely deterministic models; (c) this problem cannot be eliminated by increasing sample size or adding population controls--as sample size increases, the chance of false positive inference goes to 100%; (d) this issue can lead to systematic false positive inference of association in regions of linkage. This is important because highly deterministic models are often used intentionally in model-based analyses because they can have more power than the true model, and are implicit in many model-free analysis methods.

Siblings With Ischemic Stroke Study

Ischemic stroke has a strong familial component to risk. The Siblings With Ischemic Stroke Study (SWISS) is a genome-wide, family-based analysis that included use of imputed genotypes. The Siblings With Ischemic Stroke Study was conducted to examine the associations between single-nucleotide polymorphisms (SNPs) and risk of stroke and stroke subtypes within pairs. The Siblings With Ischemic Stroke Study enrolled 312 probands with ischemic stroke from 70 US and Canadian centers. Affected siblings were ascertained by centers and confirmed by central record review; unaffected siblings were ascertained by telephone contact. Ischemic stroke was subtyped according to Trial of Org 10172 in Acute Stroke Treatment criteria. Genotyping was performed with an Illumina 610 quad array (probands) and an Illumina linkage V array (affected siblings). SNPs were imputed by using 1000 Genomes Project data and MACH software. Family-based association analyses were conducted by using the sibling transmission-disequilibrium test. For all pairs, the correlation of age at stroke within pairs of affected siblings was r=0.83 (95% CI, 0.78-0.86; P<2.2×10(-16)). The correlation did not differ substantially by subtype. The concordance of stroke subtypes among affected pairs was 33.8% (kappa=0.13; P=5.06×10(-4)) and did not differ by age at stroke in the proband. Although no SNP achieved genome-wide significance for risk of ischemic stroke, there was clustering of the most associated SNPs on chromosomes 3p (neuronal nitric oxide synthase) and 6p. Stroke subtype and age at stroke in affected sibling pairs exhibit significant clustering. No individual SNP reached genome-wide significance. However, 2 promising candidate loci were identified, including 1 that contains neuronal nitric oxide synthase, although these risk loci warrant further examination in larger sample collections.

No association of polymorphisms in the CDK5, NDEL1, and LIS1 with autism in Chinese Han population

Autism is a pervasive neurodevelopmental disorder. CDK5 (cyclin-dependent kinase 5) and its interacting molecules are involved in neurodevelopment. We performed a family-based association analysis between CDK5, NDEL1, and LIS1 polymorphisms and autism in a Chinese Han population. Our study did not detect a significant association. It indicated that common genetic variations in these genes might not play a role in the genetic predisposition to autism.

Genetic and Functional Evidence Implicating DLL1 as the Gene That Influences Susceptibility to Visceral Leishmaniasis at Chromosome 6q27

Visceral leishmaniasis (VL) is caused by Leishmania donovani and Leishmania infantum chagasi. Genome-wide linkage studies from Sudan and Brazil identified a putative susceptibility locus on chromosome 6q27. Twenty-two single-nucleotide polymorphisms (SNPs) at genes PHF10, C6orf70, DLL1, FAM120B, PSMB1, and TBP were genotyped in 193 VL cases from 85 Sudanese families, and 8 SNPs at genes PHF10, C6orf70, DLL1, PSMB1, and TBP were genotyped in 194 VL cases from 80 Brazilian families. Family-based association, haplotype, and linkage disequilibrium analyses were performed. Multispecies comparative sequence analysis was used to identify conserved noncoding sequences carrying putative regulatory elements. Quantitative reverse-transcription polymerase chain reaction measured expression of candidate genes in splenic aspirates from Indian patients with VL compared with that in the control spleen sample. Positive associations were observed at PHF10, C6orf70, DLL1, PSMB1, and TBP in Sudan, but only at DLL1 in Brazil (combined P = 3 × 10(-4) at DLL1 across Sudan and Brazil). No functional coding region variants were observed in resequencing of 22 Sudanese VL cases. DLL1 expression was significantly (P = 2 × 10(-7)) reduced (mean fold change, 3.5 [SEM, 0.7]) in splenic aspirates from patients with VL, whereas other 6q27 genes showed higher levels (1.27 × 10(-6) < P < .01) than did the control spleen sample. A cluster of conserved noncoding sequences with putative regulatory variants was identified in the distal promoter of DLL1. DLL1, which encodes Delta-like 1, the ligand for Notch3, is strongly implicated as the chromosome 6q27 VL susceptibility gene.

Association BetweenMSX1Variants and Oral Clefts in Han Chinese in Western China

Previous animal and association studies have shown that the MSX1 gene is associated with oral clefts. Our aim was to investigate association between variants in the MSX1 gene and oral clefts in a Han Chinese population. Our study group consisted of 206 nonsyndromic oral cleft (NSOC) nuclear families (including the patients and their parents) and 224 controls. The three variants evaluated in this study were single-nucleotide polymorphisms rs3821949 and rs12532 and a missense mutation P147Q. Polymerase chain reaction-restriction fragment length polymorphism was used to genotype the three markers. Case-control and family-based association analyses were carried out. In the case-control analysis, no significant differences in genotypic or allelic frequencies were observed in any of the two single-nucleotide polymorphisms between patients and controls. Although the homozygous T allele for P147Q was not detected in any sample in this population, heterozygotes were more prevalent in NSOC (1.2%) when compared with the controls (0%). The analyses for family-based association did not suggest association between any of the three variants and NSOC. No significant association was found between NSOC and rs3821949 or rs12532 in MSX1 gene, whereas an association was observed between the P147Q variant and cleft lip with cleft palate in the case-control analysis.

Perspectives on genome‐wide multi‐stage family‐based association studies

With the establishment of large consortiums of researchers, genome-wide association (GWA) studies have become increasingly popular and feasible. Although most of these association studies focus on unrelated individuals, a lot of advantages can be exploited by including families in the analysis as well. To overcome the additional genotyping cost, multi-stage designs are particularly useful. In this article, I offer a perspective view on genome-wide family-based association analyses, both within a model-based and model-free paradigm. I highlight how multi-stage designs and analysis techniques, which are quite popular in clinical epidemiology, can enter GWA settings. I furthermore discuss how they have proven successful in reducing analysis complexity, and in overcoming one of the most cumbersome statistical hurdles in the genome-wide context, namely controlling increased false positives due to multiple testing.

Polymorphisms in ABLIM1 are Associated with Personality Traits and Alcohol Dependence

Personality traits like novelty seeking (NS), harm avoidance (HA), and reward dependence (RD) are known to be moderately heritable (30-60%). These personality traits and their comorbidities, such as alcohol dependence (AD), may share genetic components. We examined 11,120 single nucleotide polymorphisms (SNPs) genotyped in 292 nuclear families from the Genetic Analysis Workshop 14, a subset from the Collaborative Study on the Genetics of Alcoholism (COGA). A family-based association analysis was performed using the FBAT program. NS, HA, and RD were treated as quantitative traits and AD as a binary trait. Based on a multivariate association test of three quantitative traits in FBAT, we observed 20 SNPs with p < 10(-3). Interestingly, several genes (TESK2, TIPARP, THEMIS, ABLIM1, RFX4, STON2 and LILRA1) are associated with three personality traits with p < 10(-3) using single trait analysis and AD. Especially, SNP rs727532 within ABLIM1 gene at 10q25 showed the most significant association (p = 6.4 × 10(-5)) in the multivariate test and strong associations with NS, HA, RD, and AD (p = 4.48 × 10(-4), 1.2 × 10(-5), 5.6 × 10(-5), 3.12 × 10(-4), respectively) in the COGA sample. In addition, the association of rs727532 with AD was confirmed in a replication study. This study reports some newly recognized associations between several genetic loci and both AD and three personality traits.

Evidence of an Involvement of TFAP2A Gene in Nonsyndromic Cleft LIP with or without Cleft Palate: An Italian Study

Unraveling of factors involved in multifactorial diseases is a great challenge. Different approaches can be contemplate and applied to a variety of congenital malformations. In the present investigation TFAP2A has been considered a good candidate gene for nonsyndromic cleft lip with or without cleft palate (NSCLP) aetiology, basing on a sum of considerations. TFAP2A has been seen involved in orofacial development in mice; it is located in the NSCLP candidate region 6p24; it codes for a transcription factor which regulates expression of IRF6, a gene implied in NSCLP; finally, it is embroiled in the branchiooculofacial syndrome, that includes clefting as feature. A family based association analysis was performed with a sample study of 405 NSCLP triads. Evidence of association was obtained with both single marker and haplotype analyses, thus providing a support for TFAP2A in NSCLP aetiology.

Study of ABCB1 Multidrug Resistance Protein in a Common Orofacial Malformation

The onset of embryonic malformations is greatly determined by the intrauterine environment, conditioned by maternal lifestyle, diet, drugs and medication intake, in addition to both foetal and maternal genotypes. Maternal C677T MTHFR genotype has been identified as important factor in cleft lip with or without cleft palate (CL/P) etiology. In the present study we evaluated the possible interaction between maternal methylenetetrahydrofolate reductase (MTHFR) and foetal ABCB1 genotypes. ABCB1 gene codes for a drug-transport pump in charge to protect the cell by extruding a variety of harmful exogens, but with a reduced activity in a folate-restricted condition. Maternal 677T genotype is translated in a reduced folate availability for the developing embryo who consequently may becomes more exposed to external insults. A family based association analysis was performed to test the effect of ABCB1 polymorphisms in clefting, in the whole sample and in the stratified sample accordingly to maternal MTHFR genotype. No evidence of association between ABCB1 polymorphisms and CL/P was detected. This suggests that ABCB1 or ABCB1-MTHFR feto-maternal interaction could have no effect in orofacial clefting or could play a role in a limited number of cases.

Inducible t cell costimulator ligand (ICOSLG) influences crohn's disease susceptibility in the scottish paediatric ibd population

IntroductionInducible T cell costimulator ligand (ICOSLG) is intimately involved in the proliferation and differentiation of T lymphocytes. A locus on Chr. 21q22 harbouring the ICOSLG gene has been shown to...

Family-based association analysis of alcohol dependence in the COGA sample and replication in the Australian twin-family study

Family, twin, and adoption studies have indicated that genetic and environmental factors contribute to the development of alcohol dependence (AD). We conducted a low-density genome-wide association analysis to identify genetic variants influencing AD. We used 11,120 SNPs from the Affymetrix 10K Genechips genotyped in 116 Caucasian pedigrees (272 nuclear families) from Genetic Analysis Workshop 14, a subset from the Collaborative Study on the Genetics of Alcoholism (COGA). Family-based association analyses for AD were performed by the PBAT program for autosomal SNPs and by the FBAT program for X-chromosome SNPs. We identified 37 SNPs associated with AD (P<10(-3)), thirteen of which were located in known genes. The most significant association with AD was observed with SNP rs1986644 (P=8.51×10(-6)) at 13q22 near EDNRB gene. The next best signal was at 1q41 in USH2A (rs532342, P=1.07×10(-5)) and the third region was at 3q25.31 in TIPARP (rs1367311, P=2.31×10(-5)). Furthermore, we found support for association of MAOA gene (P=4.14×10(-4) for rs979606). Six of the 37 AD associated SNPs were confirmed to be associated with AD in Australian twin-family study sample (P<0.05). Interestingly, four SNPs in DSCAML1 at 11q23 reached the genome-wide significance (the top SNP is rs10892169 with P=5.31×10(-9)), while rs637547 in NKAIN2 at 6q21 showed strong association with AD (P=5.11×10(-7)) in the replication sample. These findings offer the potential for new insights into the pathogenesis of AD and will serve as a resource for replication in other populations to elucidate the potential role of these genetic variants in AD.

Comprehensive family‐based association study of the glutamate transporter gene SLC1A1 in obsessive‐compulsive disorder

SLC1A1 encodes a neuronal glutamate transporter and is a promising candidate gene for obsessive-compulsive disorder (OCD). Several independent research groups have reported significant associations between OCD and single nucleotide polymorphisms (SNPs) in this gene. Previously, we evaluated 13 SNPs in, or near, SLC1A1 and reported a strong association signal with rs301443, a SNP 7.5 kb downstream of the gene [Shugart et al. (2009); Am J Med Genet Part B 150B:886–892]. The aims of the current study were first, to further investigate this finding by saturating the region around rs301443; and second, to explore the entire gene more thoroughly with a dense panel of SNP markers. We genotyped an additional 111 SNPs in or near SLC1A1, covering from 9 kb upstream to 84 kb downstream of the gene at average spacing of 1.7 kb per SNP, and conducted family-based association analyses in 1,576 participants in 377 families.We found that none of the surrounding markers were in linkage disequilibrium with rs301443, nor were any associated with OCD. We also found that SNP rs4740788, located about 8.8 kb upstream of the gene, was associated with OCD in all families (P = 0.003) and in families with male affecteds (P = 0.002). A three-SNP haplotype (rs4740788–rs10491734–rs10491733) was associated with OCD in the total sample (P = 0.00015) and in families with male affecteds (P = 0.0007). Although of nominal statistical significance considering the number of comparisons, these findings provide further support for the involvement of SLC1A1 in the pathogenesis of OCD.

HLA-DRB1*04 is a novel fetal susceptibility allele in congenital heart block

ObjectiveCongenital heart block may develop in the fetus of Ro52 autoantibody positive women. The reported recurrence rate for autoantibody associated congenital heart block is however only 10–25%, indicating that other...

Association of the arginine vasopressin receptor 1A (AVPR1A) haplotypes with listening to music

Music is listened in all cultures. We hypothesize that willingness to produce and perceive sound and music is social communication that needs musical aptitude. Here, listening to music was surveyed using a web-based questionnaire and musical aptitude using the auditory structuring ability test (Karma Music test) and Carl Seashores tests for pitch and for time. Three highly polymorphic microsatellite markers (RS3, RS1 and AVR) of the arginine vasopressin receptor 1A (AVPR1A) gene, previously associated with social communication and attachment, were genotyped and analyzed in 31 Finnish families (n=437 members) using family-based association analysis. A positive association between the AVPR1A haplotype (RS1 and AVR) and active current listening to music (permuted P=0.0019) was observed. Other AVPR1A haplotype (RS3 and AVR) showed association with lifelong active listening to music (permuted P=0.0022). In addition to AVPR1A, two polymorphisms (5-HTTLPR and variable number of tandem repeat) of human serotonin transporter gene (SLC6A4), a candidate gene for many neuropsychiatric disorders and previously associated with emotional processing, were analyzed. No association between listening to music and the polymorphisms of SLC6A4 were detected. The results suggest that willingness to listen to music is related to neurobiological pathways affecting social affiliation and communication.

Linkage Analysis of Autopsy-Confirmed Familial Alzheimer Disease Supports an Alzheimer Disease Locus in 8q24

Background/Aims: We have previously reported the results of an extended genome-wide scan of Swedish Alzheimer disease (AD)-affected families; in this paper, we analyzed a subset of these families with autopsy-confirmed AD. Methods: We report the fine-mapping, using both microsatellite markers and single-nucleotide polymorphisms (SNPs), in the observed maximum logarithm of the odds (LOD)-2 unit (LOD_max-2) region under the identified linkage peak, linkage analysis of the fine-mapping data with additionally analyzed pedigrees, and association analysis of SNPs selected from candidate genes in the linked interval. The subset was made on the criterion of at least one autopsy-confirmed AD case per family, resulting in 24 families. Results: Linkage analysis of a family subset having at least one autopsy-confirmed AD case showed a significant nonparametric single-point LOD score of 4.4 in 8q24. Fine-mapping under the linkage peak with 10 microsatellite markers yielded an increase in the multipoint (mpt) LOD score from 2.1 to 3.0. SNP genotyping was performed on 21 selected candidate transcripts of the LOD_max-2 region. Both family-based association and linkage analysis were performed on extended material from 30 families, resulting in a suggestive linkage at peak marker rs6577853 (mpt LOD score = 2.4). Conclusion: The 8q24 region has been implicated to be involved in AD etiology.

Role of gene–gene/gene–environment interaction in the etiology of eastern Indian ADHD probands

Associations between attention deficit hyperactivity disorder (ADHD) and genetic polymorphisms in the dopamine receptors, transporter and metabolizing enzymes have been reported in different ethnic groups. Gene variants may affect disease outcome by acting synergistically or antagonistically and thus their combined effect becomes an important aspect to study in the disease etiology. In the present investigation, interaction between ten functional polymorphisms in DRD4, DAT1, MAOA, COMT, and DBH genes were explored in the Indo-Caucasoid population. ADHD cases were recruited based on DSM-IV criteria. Peripheral blood samples were collected from ADHD probands (N=126), their parents (N=233) and controls (N=96) after obtaining informed written consent for participation. Genomic DNA was subjected to PCR based analysis of single nucleotide polymorphisms and variable number of tandem repeats (VNTRs). Data obtained was examined for population as well as family-based association analyses. While case-control analysis revealed higher occurrence of DAT1 intron 8 VNTR 5R allele (P=0.02) in cases, significant preferential transmission of the 7R-T (DRD4 exon3 VNTR-rs1800955) and 3R-T (MAOA-u VNTR-rs6323) haplotypes were noticed from parents to probands (P=0.02 and 0.002 respectively). Gene-gene interaction analysis revealed significant additive effect of DBH rs1108580 and DRD4 rs1800955 with significant main effects of DRD4 exon3 VNTR, DAT1 3'UTR and intron 8 VNTR, MAOA u-VNTR, rs6323, COMT rs4680, rs362204, DBH rs1611115 and rs1108580 thereby pointing towards a strong association of these markers with ADHD. Correlation between gene variants, high ADHD score and low DBH enzymatic activity was also noticed, especially in male probands. From these observations, an impact of the studied sites on the disease etiology could be speculated in this ethnic group.

Polymorphisms in programmed death-1 gene are not associated with chronic HBV infection in Chinese patients

To investigate the association between the programmed death-1(PD-1) polymorphisms and genetic susceptibility of chronic hepatitis B virus (HBV) infection in Chinese patients. Two single nucleotide polymorphisms (SNPs), PD-1.1 G > A and PD-1.2 G > A, were genotyped in 539 patients with chronic HBV infection and 353 other family members (HbsAg-) from 256 nuclear families using polymerase chain reactiorestriction fragment length polymorphisms assay. The associations between PD-1 polymorphisms and genetic susceptibility of chronic HBV infection were analyzed usng the family-based association analysis method. No association or linkage was detected among 539 patients. Univariate (single-marker) family-based association tests demonstrated that PD-1 genotypes, alleles and transmitted haplotypes are not associated with chronic HBV infection (all with P value more than 0.05). Transmission/disequilibrium test and sibship disequilibrium test analysis showed no excess of the alleles from heterozygous parents to affected offspring (P = 0.688880, P = 1.000000 respectively). The data demonstrated that PD-1.1 and PD-1.2 polymorphisms are not associated with chronic HBV infection in Chinese patients.

Linkage and association on 8p21.2-p21.1 in schizophrenia

In the past decade, we and others have consistently reported linkage to a schizophrenia (SZ) susceptibility region on chromosome 8p21. Most recently, in the largest SZ linkage sample to date, a multi-site international collaboration performed a SNP-based linkage scan (~6,000 SNPs; 831 pedigrees; 121 from Johns Hopkins (JHU)), that showed the strongest evidence for linkage in a 1 Mb region of chr 8p21 from rs1561817 to rs9797 (Z(max) = 3.22, P = 0.0004) [Holmans et al. 2009. Mol Psychiatry]. We have investigated this 8p21 peak region further in two ways: first by linkage and family-based association in 106 8p-linked European-Caucasian (EUC) JHU pedigrees using 1,402 SNPs across a 4.4 Mb region surrounding the peak; second, by an independent case-control association study in the genetically more homogeneous Ashkenazim (AJ) (709 cases, 1,547 controls) using 970 SNPs in a further narrowed 2.8 Mb region. Family-based association analyses in EUC pedigrees and case-control analyses in AJ samples reveal significant associations for SNPs in and around DPYSL2 and ADRA1A, candidate genes previously associated with SZ in our work and others. Further, several independent gene expression studies have shown that DPYSL2 is differentially expressed in SZ brains [Beasley et al. 2006. Proteomics 6(11):3414–3425; Edgar et al. 2000. Mol Psychiatry 5(1):85–90; Johnston-Wilson et al. 2000. Mol Psychiatry 5(2):142–149] or in response to psychosis-inducing pharmaceuticals [Iwazaki et al. 2007. Proteomics 7(7):1131–1139; Paulson et al. 2004. Proteomics 4(3):819–825]. Taken together, this work further supports DPYSL2 and the surrounding genomic region as a susceptibility locus for SZ.

Fine mapping of candidate regions for bipolar disorder provides strong evidence for susceptibility loci on chromosomes 7q

Genomewide scans of bipolar disorder (BP) have not produced consistent linkage findings. Follow-up studies using enlarged samples and enhanced marker density can bolster or refute claims of linkage and pave the way to gene discovery. We conducted linkage and association analyses, using a ~3-cM density map of 10 candidate regions, in a large BP pedigree sample (865 individuals from 56 pedigrees). The candidate regions were identified in a previous 10-cM genome-wide scan using a subset of this sample (373 individuals from 40 pedigrees). The present sample consists of the expanded original pedigrees ("core" pedigrees) and 16 additional pedigrees. We obtained experiment-wide significant linkage on chromosome 7q34 (LOD score 3.53, P < 0.001), substantially stronger than that observed in the genome-wide scan. Support for linkage was sustained on chromosomes 2p13, 4q31, 8q13, 13q32, 14q21, and 17q11, though at a more modest level. Family-based association analysis was consistent with the linkage results at all regions with linkage evidence, except 4q an 8q, but the results fell short of statistical significance. Three of the previously implicated regions-9q31, 10q21 and 10q24-showed substantial reduction in evidence of linkage. Our results strongly support 7q34 as a region harboring susceptibility locus for BP. Somewhat lesser, yet notable support was obtained for 2p13, 4q31, 8q13, 13q32, 14q21, and 17q11. These regions could be considered prime candidates for future gene finding efforts.