Family-based Genome-wide Association Studies Research Articles

Genome-wide association studies have problems due to confounding: Are family-based designs the answer?

Genome-wide association studies (GWASs) can be affected by confounding. Family-based GWAS uses random, within-family genetic variation to avoid this. A study in PLOS Biology details how different sources of confounding affect GWAS and whether family-based designs offer a solution.

Interpreting population- and family-based genome-wide association studies in the presence of confounding.

A central aim of genome-wide association studies (GWASs) is to estimate direct genetic effects: the causal effects on an individual's phenotype of the alleles that they carry. However, estimates of direct effects can be subject to genetic and environmental confounding and can also absorb the "indirect" genetic effects of relatives' genotypes. Recently, an important development in controlling for these confounds has been the use of within-family GWASs, which, because of the randomness of mendelian segregation within pedigrees, are often interpreted as producing unbiased estimates of direct effects. Here, we present a general theoretical analysis of the influence of confounding in standard population-based and within-family GWASs. We show that, contrary to common interpretation, family-based estimates of direct effects can be biased by genetic confounding. In humans, such biases will often be small per-locus, but can be compounded when effect-size estimates are used in polygenic scores (PGSs). We illustrate the influence of genetic confounding on population- and family-based estimates of direct effects using models of assortative mating, population stratification, and stabilizing selection on GWAS traits. We further show how family-based estimates of indirect genetic effects, based on comparisons of parentally transmitted and untransmitted alleles, can suffer substantial genetic confounding. We conclude that, while family-based studies have placed GWAS estimation on a more rigorous footing, they carry subtle issues of interpretation that arise from confounding.

Polygenic score: An anchor holding the whole life course.

Polygenic score: An anchor holding the whole life course.

A disanalogy with RCTs and its implications for second-generation causal knowledge.

We are less optimistic than Madole & Harden that family-based genome-wide association studies (GWASs) will lead to significant second-generation causal knowledge. Despite bearing some similarities, family-based GWASs and randomised controlled trials (RCTs) are not identical. Most RCTs assess a relatively homogenous causal stimulus as a treatment, whereas GWASs assess highly heterogeneous causal stimuli. Thus, GWAS results will not translate so easily into second-generation causal knowledge.

Within-sibship genome-wide association analyses decrease bias in estimates of direct genetic effects

Estimates from genome-wide association studies (GWAS) of unrelated individuals capture effects of inherited variation (direct effects), demography (population stratification, assortative mating) and relatives (indirect genetic effects). Family-based GWAS designs can control for demographic and indirect genetic effects, but large-scale family datasets have been lacking. We combined data from 178,086 siblings from 19 cohorts to generate population (between-family) and within-sibship (within-family) GWAS estimates for 25 phenotypes. Within-sibship GWAS estimates were smaller than population estimates for height, educational attainment, age at first birth, number of children, cognitive ability, depressive symptoms and smoking. Some differences were observed in downstream SNP heritability, genetic correlations and Mendelian randomization analyses. For example, the within-sibship genetic correlation between educational attainment and body mass index attenuated towards zero. In contrast, analyses of most molecular phenotypes (for example, low-density lipoprotein-cholesterol) were generally consistent. We also found within-sibship evidence of polygenic adaptation on taller height. Here, we illustrate the importance of family-based GWAS data for phenotypes influenced by demographic and indirect genetic effects.

Genetic association studies of alterations in protein function expose recessive effects on cancer predisposition

The characterization of germline genetic variation affecting cancer risk, known as cancer predisposition, is fundamental to preventive and personalized medicine. Studies of genetic cancer predisposition typically identify significant genomic regions based on family-based cohorts or genome-wide association studies (GWAS). However, the results of such studies rarely provide biological insight or functional interpretation. In this study, we conducted a comprehensive analysis of cancer predisposition in the UK Biobank cohort using a new gene-based method for detecting protein-coding genes that are functionally interpretable. Specifically, we conducted proteome-wide association studies (PWAS) to identify genetic associations mediated by alterations to protein function. With PWAS, we identified 110 significant gene-cancer associations in 70 unique genomic regions across nine cancer types and pan-cancer. In 48 of the 110 PWAS associations (44%), estimated gene damage is associated with reduced rather than elevated cancer risk, suggesting a protective effect. Together with standard GWAS, we implicated 145 unique genomic loci with cancer risk. While most of these genomic regions are supported by external evidence, our results also highlight many novel loci. Based on the capacity of PWAS to detect non-additive genetic effects, we found that 46% of the PWAS-significant cancer regions exhibited exclusive recessive inheritance. These results highlight the importance of recessive genetic effects, without relying on familial studies. Finally, we show that many of the detected genes exert substantial cancer risk in the studied cohort determined by a quantitative functional description, suggesting their relevance for diagnosis and genetic consulting.

Family-Based Genome-Wide Association Study of Autism Spectrum Disorder in Middle Eastern Families.

Autism spectrum disorder (ASD) is a neurodevelopmental disease characterized by abnormalities in language and social communication with substantial clinical heterogeneity. Genetic factors play an important role in ASD with heritability estimated between 70% to 80%. Genome-wide association studies (GWAS) have identified multiple loci associated with ASD. However, most studies were performed on European populations and little is known about the genetic architecture of ASD in Middle Eastern populations. Here, we report the first GWAS of ASD in the Middle eastern population of Qatar. We analyzed 171 families with ASD, using linear mixed models adjusting for relatedness and other confounders. Results showed that common single nucleotide polymorphisms (SNP) in seven loci are associated with ASD (p < 1 × 10−5). Although the identified loci did not reach genome-wide significance, many of the top associated SNPs are located within or near genes that have been implicated in ASD or related neurodevelopmental disorders. These include GORASP2, GABBR2, ANKS6, THSD4, ERCC6L, ARHGEF6, and HDAC8. Additionally, three of the top associated SNPs were significantly associated with gene expression. We also found evidence of association signals in two previously reported ASD-susceptibility loci (rs10099100 and rs4299400). Our results warrant further functional studies and replication to provide further insights into the genetic architecture of ASD.

Heritability and family-based GWAS analyses of the N-acyl ethanolamine and ceramide plasma lipidome.

Signalling lipids of the N-acyl ethanolamine (NAE) and ceramide (CER) classes have emerged as potential biomarkers of cardiovascular disease (CVD). We sought to establish the heritability of plasma NAEs (including the endocannabinoid anandamide) and CERs, to identify common DNA variants influencing the circulating concentrations of the heritable lipids, and assess causality of these lipids in CVD using 2-sample Mendelian randomization (2SMR). Nine NAEs and 16 CERs were analyzed in plasma samples from 999 members of 196 British Caucasian families, using targeted ultra-performance liquid chromatography with tandem mass spectrometry. All lipids were significantly heritable (h2 = 36–62%). A missense variant (rs324420) in the gene encoding the enzyme fatty acid amide hydrolase (FAAH), which degrades NAEs, associated at genome-wide association study (GWAS) significance (P < 5 × 10−8) with four NAEs (DHEA, PEA, LEA and VEA). For CERs, rs680379 in the SPTLC3 gene, which encodes a subunit of the rate-limiting enzyme in CER biosynthesis, associated with a range of species (e.g. CER[N(24)S(19)]; P = 4.82 × 10−27). We observed three novel associations between SNPs at the CD83, SGPP1 and DEGS1 loci, and plasma CER traits (P < 5 × 10−8). 2SMR in the CARDIoGRAMplusC4D cohorts (60 801 cases; 123 504 controls) and in the DIAGRAM cohort (26 488 cases; 83 964 controls), using the genetic instruments from our family-based GWAS, did not reveal association between genetically determined differences in CER levels and CVD or diabetes. Two of the novel GWAS loci, SGPP1 and DEGS1, suggested a casual association between CERs and a range of haematological phenotypes, through 2SMR in the UK Biobank, INTERVAL and UKBiLEVE cohorts (n = 110 000–350 000).

ADAMTS12, a new candidate gene for pediatric stroke.

We recently reported a family-based genome wide association study (GWAS) for pediatric stroke pointing our attention to two significantly associated genes of the ADAMTS (a disintegrin and metalloproteinase with thrombospondin motifs) gene family ADAMTS2 (rs469568, p = 8x10-6) and ADAMTS12 (rs1364044, p = 2.9x10-6). To further investigate these candidate genes, we applied a targeted resequencing approach on 48 discordant sib-pairs for pediatric stroke followed by genotyping of the detected non-synonymous variants in the full cohort of 270 offspring trios and subsequent fine mapping analysis. We identified eight non-synonymous SNPs in ADAMTS2 and six in ADAMTS12 potentially influencing the respective protein function. These variants were genotyped within a cohort of 270 affected offspring trios, association analysis revealed the ADAMTS12 variant rs77581578 to be significantly under-transmitted (p = 6.26x10-3) to pediatric stroke patients. The finding was validated in a pediatric venous thromboembolism (VTE) cohort of 189 affected trios. Subsequent haplotype analysis of ADAMTS12 detected a significantly associated haplotype comprising the originally identified GWAS variant. Several ADAMTS genes such as ADAMTS13 are involved in thromboembolic disease process. Here, we provide further evidence for ADAMTS12 to likely play a role in pediatric stroke. Further functional studies are warranted to assess the functional role of ADAMTS12 in the pathogenesis of stroke.

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.

A Family-Based Genome Wide Association Study of Externalizing Behaviors.

Shared genetic factors contribute to the high degree of comorbidity among externalizing problems (e.g. substance use and antisocial behavior). We leverage this common genetic etiology to identify genetic influences externalizing problems in participants from the Collaborative Study on the Genetics of Alcoholism (European ancestry = 7568; African ancestry = 3274). We performed a family-based genome-wide association study (GWAS) on externalizing scores derived from criterion counts of five DSM disorders (alcohol dependence, alcohol abuse, illicit drug dependence, illicit drug abuse, and either antisocial personality disorder or conduct disorder). We meta analyzed these results with a similar measure of externalizing in an independent sample, Spit for Science (combined sample N = 15,112). We did not discover any robust genome-wide significant signals. Polygenic scores derived from the ancestry-specific GWAS summary statistics predicted externalizing problems in an independent European ancestry sample, but not in those of African ancestry. However, these PRS were no longer significant after adjusting for multiple testing. Larger samples with deep phenotyping are necessary for the discovery of SNPs related to externalizing problems.

Family-based genome wide association analysis for salmon lice (Lepeophtheirus salmonis) resistance in North American Atlantic salmon using a 50 K SNP array

The salmon louse, Lepeophtheirus salmonis, is an economically important parasite on aquaculture stocks of Atlantic salmon grown in sea cages. The aim of this study was to perform a family-based genome wide association study to identify the genomic architecture of salmon louse resistance in a commercial population of Atlantic salmon using five year classes of data. A total of 1756 fish were genotyped from 248 families, which were experimentally infected with L. salmonis copepodids, using a 50 K SNP array designed specifically for North American Atlantic salmon. A within family-based association test for quantitative traits was run using QFAM in PLINK with salmon lice counts pre-adjusted for contemporary group as the phenotype. Salmon lice counts varied significantly among families and salmon lice resistance had an estimated heritability of 0.21 ± 0.03. After correcting for multiple testing, seven SNPs from three chromosomes surpassed a suggestive chromosome wide level. Ssa04 had the most suggestive SNP, and four other suggestive SNPs within a close range on the chromosome. A single SNP was suggestive on each of Ssa14 and Ssa20. The most suggestive SNP, and two suggestive tag SNPs, were located near a predicted gene on Ssa04 that codes for ceramide-1-phophate transfer protein which is involved with pathophysiological functions such as cell survival and inflammation. No loci with major effects were detected suggesting a polygenic trait architecture. However, several loci were detected which could be included directly in models used to estimate genomic breeding values. Therefore, genomic selection should continue to be used to breed Atlantic salmon naturally more resistant to L. salmonis infection.

A rare missense variant of CASP7 is associated with familial late-onset Alzheimer's disease

IntroductionThe genetic architecture of Alzheimer's disease (AD) is only partially understood. MethodsWe conducted an association study for AD using whole sequence data from 507 genetically enriched AD cases (i.e., cases having close relatives affected by AD) and 4917 cognitively healthy controls of European ancestry (EA) and 172 enriched cases and 179 controls of Caribbean Hispanic ancestry. Confirmation of top findings from stage 1 was sought in two family-based genome-wide association study data sets and in a whole genome–sequencing data set comprising members from 42 EA and 115 Caribbean Hispanic families. ResultsWe identified associations in EAs with variants in 12 novel loci. The most robust finding is a rare CASP7 missense variant (rs116437863; P = 2.44 × 10−10) which improved when combined with results from stage 2 data sets (P = 1.92 × 10−10). DiscussionOur study demonstrated that an enriched case design can strengthen genetic signals, thus allowing detection of associations that would otherwise be missed in a traditional case-control study.

An exonic splicing enhancer mutation in DUOX2 causes aberrant alternative splicing and severe congenital hypothyroidism in Bama pigs.

ABSTRACTPigs share many similarities with humans in terms of anatomy, physiology and genetics, and have long been recognized as important experimental animals in biomedical research. Using an N-ethyl-N-nitrosourea (ENU) mutagenesis screen, we previously identified a large number of pig mutants, which could be further established as human disease models. However, the identification of causative mutations in large animals with great heterogeneity remains a challenging endeavor. Here, we select one pig mutant, showing congenital nude skin and thyroid deficiency in a recessive inheritance pattern. We were able to efficiently map the causative mutation using family-based genome-wide association studies combined with whole-exome sequencing and a small sample size. A loss-of-function variant (c.1226 A>G) that resulted in a highly conserved amino acid substitution (D409G) was identified in the DUOX2 gene. This mutation, located within an exonic splicing enhancer motif, caused aberrant splicing of DUOX2 transcripts and resulted in lower H2O2 production, which might cause a severe defect in thyroid hormone production. Our findings suggest that exome sequencing is an efficient way to map causative mutations and that DUOX2D409G/D409G mutant pigs could be a potential large animal model for human congenital hypothyroidism.

Family-Based Genome-Wide Association Study of South Indian Pedigrees Supports WNT7B as a Central Corneal Thickness Locus.

PurposeTo identify genetic risk factors contributing to central corneal thickness (CCT) in individuals from South India, a population with a high prevalence of ocular disorders.MethodsOne hundred ninety-five individuals from 15 large South Indian pedigrees were genotyped using the Omni2.5 bead array. Family-based association for CCT was conducted using the score test in MERLIN.ResultsGenome-wide association study (GWAS) identified strongest association for single nucleotide polymorphisms (SNPs) in the first intron of WNT7B and CCT (top SNP rs9330813; β = −0.57, 95% confidence interval [CI]: −0.78 to −0.36; P = 1.7 × 10−7). We further investigated rs9330813 in a Latino cohort and four independent European cohorts. A meta-analysis of these data sets demonstrated statistically significant association between rs9330813 and CCT (β = −3.94, 95% CI: −5.23 to −2.66; P = 1.7 × 10−9). WNT7B SNPs located in the same genomic region that includes rs9330813 have previously been associated with CCT in Latinos but with other ocular quantitative traits related to myopia (corneal curvature and axial length) in a Japanese population (rs10453441 and rs200329677). To evaluate the specificity of the observed WNT7B association with CCT in the South Indian families, we completed an ocular phenome-wide association study (PheWAS) for the top WNT7B SNPs using 45 ocular traits measured in these same families including corneal curvature and axial length. The ocular PheWAS results indicate that in the South Indian families WNT7B SNPs are primarily associated with CCT.ConclusionsThe results indicate robust evidence for association between WNT7B SNPs and CCT in South Indian pedigrees, and suggest that WNT7B SNPs can have population-specific effects on ocular quantitative traits.

Genome-Wide Analysis Reveals Four Novel Loci for Attention-Deficit Hyperactivity Disorder in Korean Youths.

Objectives:The molecular mechanisms underlying attention-deficit hyperactivity disorder (ADHD) remain unclear. Therefore, this study aimed to identify the genetic susceptibility loci for ADHD in Korean children with ADHD. We performed a case-control and a family-based genome-wide association study (GWAS), as well as genome-wide quantitative trait locus (QTL) analyses, for two symptom traits.Methods:A total of 135 subjects (71 cases and 64 controls), for the case-control analysis, and 54 subjects (27 probands and 27 unaffected siblings), for the family-based analysis, were included.Results:The genome-wide QTL analysis identified four single nucleotide polymorphisms (SNPs) (rs7684645 near APELA, rs12538843 near YAE1D1 and POU6F2, rs11074258 near MCTP2, and rs34396552 near CIDEA) that were significantly associated with the number of inattention symptoms in ADHD. These SNPs showed possible association with ADHD in the family-based GWAS, and with hyperactivity-impulsivity in genome-wide QTL analyses. Moreover, association signals in the family-based QTL analysis for the number of inattention symptoms were clustered near genes IL10, IL19, SCL5A9, and SKINTL.Conclusion:We have identified four QTLs with genome-wide significance and several promising candidates that could potentially be associated with ADHD (CXCR4, UPF1, SETD5, NALCN-AS1, ERC1, SOX2-OT, FGFR2, ANO4, and TBL1XR1). Further replication studies with larger sample sizes are needed.

A family-based genome-wide association study of chronic rhinosinusitis with nasal polyps implicates several genes in the disease pathogenesis.

BackgroundThe pathogenesis of chronic rhinosinusitis with nasal polyps is largely unknown. Previous studies have given valuable information about genetic variants associated with this disease but much is still unexplained. Our goal was to identify genetic markers and genes associated with susceptibility to chronic rhinosinusitis with nasal polyps using a family-based genome-wide association study.Methods427 patients (293 males and 134 females) with CRSwNP and 393 controls (175 males and 218 females) were recruited from several Swedish hospitals. SNP association values were generated using DFAM (implemented in PLINK) and Efficient Mixed Model Association eXpedited (EMMAX). Analyses of pathway enrichment, gene expression levels and expression quantitative trait loci were then performed in turn.ResultsNone of the analysed SNPs reached genome wide significant association of 5.0 x 10−8. Pathway analyses using our top 1000 markers with the most significant association p-values resulted in 138 target genes. A comparison between our target genes and gene expression data from the NCBI Gene Expression Omnibus database showed significant overlap for 36 of these genes. Comparisons with data from expression quantitative trait loci showed the most skewed allelic distributions in cases with chronic rhinosinusitis with nasal polyps compared with controls for the genes HLCS, HLA-DRA, BICD2, VSIR and SLC5A1.ConclusionOur study indicates that HLCS, HLA-DRA, BICD2, VSIR and SLC5A1 could be involved in the pathogenesis of chronic rhinosinusitis with nasal polyps. HLA-DRA has been associated with chronic rhinosinusitis with nasal polyps in previous studies and HLCS, BICD2, VSIR and SLC5A1 may be new targets for future research.

Genome-wide copy number variation analysis identified deletions in SFMBT1 associated with fasting plasma glucose in a Han Chinese population

BackgroundFasting glucose and fasting insulin are glycemic traits closely related to diabetes, and understanding the role of genetic factors in these traits can help reveal the etiology of type 2 diabetes. Although single nucleotide polymorphisms (SNPs) in several candidate genes have been found to be associated with fasting glucose and fasting insulin, copy number variations (CNVs), which have been reported to be associated with several complex traits, have not been reported for association with these two traits. We aimed to identify CNVs associated with fasting glucose and fasting insulin.ResultsWe conducted a genome-wide CNV association analysis for fasting plasma glucose (FPG) and fasting plasma insulin (FPI) using a family-based genome-wide association study sample from a Han Chinese population in Taiwan. A family-based CNV association test was developed in this study to identify common CNVs (i.e., CNVs with frequencies ≥ 5%), and a generalized estimating equation approach was used to test the associations between the traits and counts of global rare CNVs (i.e., CNVs with frequencies <5%). We found a significant genome-wide association for common deletions with a frequency of 5.2% in the Scm-like with four mbt domains 1 (SFMBT1) gene with FPG (association p-value = 2×10−4 and an adjusted p-value = 0.0478 for multiple testing). No significant association was observed between global rare CNVs and FPG or FPI. The deletions in 20 individuals with DNA samples available were successfully validated using PCR-based amplification. The association of the deletions in SFMBT1 with FPG was further evaluated using an independent population-based replication sample obtained from the Taiwan Biobank. An association p-value of 0.065, which was close to the significance level of 0.05, for FPG was obtained by testing 9 individuals with CNVs in the SFMBT1 gene region and 11,692 individuals with normal copies in the replication cohort.ConclusionsPrevious studies have found that SNPs in SFMBT1 are associated with blood pressure and serum urate concentration, suggesting that SFMBT1 may have functional implications in some metabolic-related traits.

A family-based genome-wide association study reveals an association of spondyloarthritis with MAPK14

ObjectiveMore than 40 loci have been associated with ankylosing spondylitis (AS), but less is known about genetic associations in spondyloarthritis (SpA) as a whole. We conducted a family-based genome-wide association...

Abstract IA03: Lung cancer risk in African Americans: Familial aggregation and genetic susceptibility

Abstract IA03: Lung cancer risk in African Americans: Familial aggregation and genetic susceptibility