Family-based Genome-wide Association Research Articles

Causal interpretations of family GWAS in the presence of heterogeneous effects

Family-based genome-wide association studies (GWASs) are often claimed to provide an unbiased estimate of the average causal effects (or average treatment effects; ATEs) of alleles, on the basis of an analogy between the random transmission of alleles from parents to children and a randomized controlled trial. We show that this claim does not hold in general. Because Mendelian segregation only randomizes alleles among children of heterozygotes, the effects of alleles in the children of homozygotes are not observable. This feature will matter if an allele has different average effects in the children of homozygotes and heterozygotes, as can arise in the presence of gene-by-environment interactions, gene-by-gene interactions, or differences in linkage disequilibrium patterns. At a single locus, family-based GWAS can be thought of as providing an unbiased estimate of the average effect in the children of heterozygotes (i.e., a local average treatment effect; LATE). This interpretation does not extend to polygenic scores (PGSs), however, because different sets of SNPs are heterozygous in each family. Therefore, other than under specific conditions, the within-family regression slope of a PGS cannot be assumed to provide an unbiased estimate of the LATE for any subset or weighted average of families. In practice, the potential biases of a family-based GWAS are likely smaller than those that can arise from confounding in a standard, population-based GWAS, and so family studies remain important for the dissection of genetic contributions to phenotypic variation. Nonetheless, their causal interpretation is less straightforward than has been widely appreciated.

Family-based genome-wide association analysis of novelty seeking in a Korean schizophrenic population: A pilot study.

Schizophrenia (SPR) is the most devastating mental illness that causes severe deterioration in social and occupational functioning, but, the etiology remains unknown. The objective of this study is to explore the genetic underpinnings of novelty seeking behavior in schizophrenic family within the Korean population. By conducting a family-based genome-wide association study, we aim to identify potential genetic markers and variations associated with novelty seeking traits in the context of SPR. We have recruited 27 probands (with SPR) with their parents and siblings whenever possible. DNA was extracted from blood sampling of 58 individuals in 27 families and analyzed in an Illumina core exome single nucleotide polymorphism (SNP) array. A family-based association test (qFAM) was used to derive SNP association values across all chromosomes. Although none of the final 800,000 SNPs reached the genome-wide significant threshold of 8.45 × 10-7, the most significant 4 SNPs were within the 10-5 to 10-7. This study identifies genetic associations between novelty seeking behavior and SPR within families. RAPGEF5 emerges as a significant gene, along with other neuropsychiatric-related genes. Noteworthy genes like DRD4 and COMT did not show associations, possibly due to the focus on schizophrenic family. While shedding light on this complex relationship, larger studies are needed for robust conclusions and deeper mechanistic insights.

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.

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.

FGF20 and PGM2 variants are associated with childhood asthma in family-based whole-genome sequencing studies.

Asthma is a heterogeneous common respiratory disease that remains poorly understood. The established genetic associations fail to explain the high estimated heritability, and the prevalence of asthma differs between populations and geographic regions. Robust association analyses incorporating different genetic ancestries and whole-genome sequencing data may identify novel genetic associations. We performed family-based genome-wide association analyses of childhood-onset asthma based on whole-genome sequencing (WGS) data for the 'The Genetic Epidemiology of Asthma in Costa Rica' study (GACRS) and the Childhood Asthma Management Program (CAMP). Based on parent-child trios with children diagnosed with asthma, we performed a single variant analysis using an additive and a recessive genetic model and a region-based association analysis of low-frequency and rare variants. Based on 1180 asthmatic trios (894 GACRS trios and 286 CAMP trios, a total of 3540 samples with WGS data), we identified three novel genetic loci associated with childhood-onset asthma: rs4832738 on 4p14 ($P=1.72\ast{10}^{-9}$, recessive model), rs1581479 on 8p22 ($P=1.47\ast{10}^{-8}$, additive model) and rs73367537 on 10q26 ($P=1.21\ast{10}^{-8}$, additive model in GACRS only). Integrative analyses suggested potential novel candidate genes underlying these associations: PGM2 on 4p14 and FGF20 on 8p22. Our family-based whole-genome sequencing analysis identified three novel genetic loci for childhood-onset asthma. Gene expression data and integrative analyses point to PGM2 on 4p14 and FGF20 on 8p22 as linked genes. Furthermore, region-based analyses suggest independent potential low-frequency/rare variant associations on 8p22. Follow-up analyses are needed to understand the functional mechanisms and generalizability of these associations.

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.

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.

Assessing thermal adaptation using family-based association and FST outlier tests in a threatened trout species.

Discovering genetic markers associated with phenotypic or ecological characteristics can improve our understanding of adaptation and guide conservation of key evolutionary traits. The Lahontan cutthroat trout (Oncorhynchus clarkii henshawi) of the northern Great Basin Desert, USA, demonstrated exceptional tolerance to high temperatures in the desert lakes where it resided historically. This trait is central to a conservation hatchery effort to protect the genetic legacy of the nearly extinct lake ecotype. We genotyped full-sibling families from this conservation broodstock and samples from the only two remaining, thermally distinct, native lake populations at 4,644 new single nucleotide polymorphisms (SNPs). Family-based genome-wide association testing of the broodstock identified nine and 26 SNPs associated with thermal tolerance (p<0.05 and p<0.1), measured in a previous thermal challenge experiment. Genes near the associated SNPs had complex functions related to immunity, growth, metabolism and ion homeostasis. Principal component analysis using the thermotolerance-related SNPs showed unexpected divergence between the conservation broodstock and the native lake populations at these loci. FST outlier tests on the native lake populations identified 18 loci shared between two or more of the tests, with two SNPs identified by all three tests (p<0.01); none overlapped with loci identified by association testing in the broodstock. A recent history of isolation and the complex genetic and demographic backgrounds of Lahontan cutthroat trout probably limited our ability to find shared thermal tolerance loci. Our study extends the still relatively rare application of genomic tools testing for markers associated with important phenotypic or environmental characteristics in species of conservation concern.

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 of inflammation biomarkers and fenofibrate treatment response in the GOLDN study

In this paper we analyzed whole-genome genetic information provided by GAW20 from the Genetics of Lipid Lowering Drugs and Diet Network (GOLDN) study for family data. Lipid levels such as triglycerides (TGs) and high-density lipoprotein (HDL) are measured at different time points before and after administration of an anti-inflammatory drug fenofibrate. Apart from that, the data contain some covariates and whole-genome genotype information. We propose 2 novel approaches based on Henderson’s iterative mixed model to identify associated loci corresponding to (a) inflammatory biomarkers like TGs and HDLs together over time, and (b) the response to fenofibrate treatment. We developed a mixed-model approach using both TG and HDL phenotypes at all 4 time points for a genetic association study whereas we used TGs only to study genetic association with response to the drug. We expect that use of complete family data in a longitudinal framework under a single model involving the appropriate correlation structures would be able to exploit the maximum possible information contained in the sample. Our analysis of whole-genome single nucleotide polymorphisms (SNPs) and genomic regions corresponding to drug treatment finds no significant locus after multiple correction. Arguably, the moderately small sample size of the data set, as compared to the sample size usually used in genome-wide association studies (GWAS), could be a reason for such a result. Nevertheless, we report the top 20 SNPs associated with the phenotypes, and the top 20 SNPs and genomic regions associated with a response to fenofibrate treatment. Application of our methods to larger GWAS and further functional validation of the reported top SNPs and genomic regions might provide important biological insight into the genetic constitution of the trait.

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.

The Central Role of KNG1 Gene as a Genetic Determinant of Coagulation Pathway-Related Traits: Exploring Metaphenotypes.

Traditional genetic studies of single traits may be unable to detect the pleiotropic effects involved in complex diseases. To detect the correlation that exists between several phenotypes involved in the same biological process, we introduce an original methodology to analyze sets of correlated phenotypes involved in the coagulation cascade in genome-wide association studies. The methodology consists of a two-stage process. First, we define new phenotypic meta-variables (linear combinations of the original phenotypes), named metaphenotypes, by applying Independent Component Analysis for the multivariate analysis of correlated phenotypes (i.e. the levels of coagulation pathway–related proteins). The resulting metaphenotypes integrate the information regarding the underlying biological process (i.e. thrombus/clot formation). Secondly, we take advantage of a family based Genome Wide Association Study to identify genetic elements influencing these metaphenotypes and consequently thrombosis risk. Our study utilized data from the GAIT Project (Genetic Analysis of Idiopathic Thrombophilia). We obtained 15 metaphenotypes, which showed significant heritabilities, ranging from 0.2 to 0.7. These results indicate the importance of genetic factors in the variability of these traits. We found 4 metaphenotypes that showed significant associations with SNPs. The most relevant were those mapped in a region near the HRG, FETUB and KNG1 genes. Our results are provocative since they show that the KNG1 locus plays a central role as a genetic determinant of the entire coagulation pathway and thrombus/clot formation. Integrating data from multiple correlated measurements through metaphenotypes is a promising approach to elucidate the hidden genetic mechanisms underlying complex diseases.

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