Using linkage analysis of large pedigrees to guide association analyses

To date, genome-wide association studies have yielded discoveries of common variants that partly explain familial aggregation of diseases and traits. Researchers are now turning their attention to less common variants because the price of sequencing has dropped drastically. However, because sequencing of the whole genome in large samples is costly, great care must be taken to prioritize which samples and which genomic regions are selected for sequencing. We are interested in identifying genomic regions for deep sequencing using large multiplex families collected as part of earlier linkage studies. We incorporate linkage analysis into our search for Q1-associated alleles. Overall, we found that power was low for both whole-exome and linkage-guided sequencing analysis. By restricting sequencing to regions with high LOD peaks, we found fewer associated single-nucleotide polymorphisms than by using whole-exome sequencing. However, incorporating linkage analysis enabled us to detect more than half of the associated susceptibility loci (52%) that would have been identified by whole-exome sequencing while examining only 2.5% of the exome. This result suggests that incorporating linkage results from large multiplex families might greatly increase the efficiency of sequencing to detect trait-associated alleles in complex disease

Associations of NINJ2 sequence variants with incident ischemic stroke in the Cohorts for Heart and Aging in Genomic Epidemiology (CHARGE) consortium

Background<p></p> Stroke, the leading neurologic cause of death and disability, has a substantial genetic component. We previously conducted a genome-wide association study (GWAS) in four prospective studies from the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) consortium and demonstrated that sequence variants near the NINJ2 gene are associated with incident ischemic stroke. Here, we sought to fine-map functional variants in the region and evaluate the contribution of rare variants to ischemic stroke risk.<p></p> Methods and Results<p></p> We sequenced 196 kb around NINJ2 on chromosome 12p13 among 3,986 European ancestry participants, including 475 ischemic stroke cases, from the Atherosclerosis Risk in Communities Study, Cardiovascular Health Study, and Framingham Heart Study. Meta-analyses of single-variant tests for 425 common variants (minor allele frequency [MAF] ≥ 1%) confirmed the original GWAS results and identified an independent intronic variant, rs34166160 (MAF = 0.012), most significantly associated with incident ischemic stroke (HR = 1.80, p = 0.0003). Aggregating 278 putatively-functional variants with MAF≤ 1% using count statistics, we observed a nominally statistically significant association, with the burden of rare NINJ2 variants contributing to decreased ischemic stroke incidence (HR = 0.81; p = 0.026).<p></p> Conclusion<p></p> Common and rare variants in the NINJ2 region were nominally associated with incident ischemic stroke among a subset of CHARGE participants. Allelic heterogeneity at this locus, caused by multiple rare, low frequency, and common variants with disparate effects on risk, may explain the difficulties in replicating the original GWAS results. Additional studies that take into account the complex allelic architecture at this locus are needed to confirm these findings

Atrial fibrillation genetic risk differentiates cardioembolic stroke from other stroke subtypes

Objective: We sought to assess whether genetic risk factors for atrial fibrillation (AF) can explain cardioembolic stroke risk. Methods: We evaluated genetic correlations between a previous genetic study of AF and AF in the presence of cardioembolic stroke using genome-wide genotypes from the Stroke Genetics Network (N = 3,190 AF cases, 3,000 cardioembolic stroke cases, and 28,026 referents). We tested whether a previously validated AF polygenic risk score (PRS) associated with cardioembolic and other stroke subtypes after accounting for AF clinical risk factors. Results: We observed a strong correlation between previously reported genetic risk for AF, AF in the presence of stroke, and cardioembolic stroke (Pearson r = 0.77 and 0.76, respectively, across SNPs with p < 4.4 × 10−4 in the previous AF meta-analysis). An AF PRS, adjusted for clinical AF risk factors, was associated with cardioembolic stroke (odds ratio [OR] per SD = 1.40, p = 1.45 × 10−48), explaining ∼20% of the heritable component of cardioembolic stroke risk. The AF PRS was also associated with stroke of undetermined cause (OR per SD = 1.07, p = 0.004), but no other primary stroke subtypes (all p > 0.1). Conclusions: Genetic risk of AF is associated with cardioembolic stroke, independent of clinical risk factors. Studies are warranted to determine whether AF genetic risk can serve as a biomarker for strokes caused by AF

Deep learning enables genetic analysis of the human thoracic aorta

Genome-wide association analyses identify variants associated with thoracic aortic diameter. A polygenic score for ascending aortic diameter was associated with a diagnosis of thoracic aortic aneurysm in independent samples. Enlargement or aneurysm of the aorta predisposes to dissection, an important cause of sudden death. We trained a deep learning model to evaluate the dimensions of the ascending and descending thoracic aorta in 4.6 million cardiac magnetic resonance images from the UK Biobank. We then conducted genome-wide association studies in 39,688 individuals, identifying 82 loci associated with ascending and 47 with descending thoracic aortic diameter, of which 14 loci overlapped. Transcriptome-wide analyses, rare-variant burden tests and human aortic single nucleus RNA sequencing prioritized genes including SVIL, which was strongly associated with descending aortic diameter. A polygenic score for ascending aortic diameter was associated with thoracic aortic aneurysm in 385,621 UK Biobank participants (hazard ratio = 1.43 per s.d., confidence interval 1.32-1.54, P = 3.3 x 10(-20)). Our results illustrate the potential for rapidly defining quantitative traits with deep learning, an approach that can be broadly applied to biomedical images

Polygenic overlap between C-reactive protein, plasma lipids, and Alzheimer disease

Background—Epidemiological findings suggest a relationship between Alzheimer disease (AD), inflammation, and dyslipidemia, although the nature of this relationship is not well understood. We investigated whether this phenotypic association arises from a shared genetic basis. Methods and Results—Using summary statistics (P values and odds ratios) from genome-wide association studies of >200 000 individuals, we investigated overlap in single-nucleotide polymorphisms associated with clinically diagnosed AD and C-reactive protein (CRP), triglycerides, and high- and low-density lipoprotein levels. We found up to 50-fold enrichment of AD single-nucleotide polymorphisms for different levels of association with C-reactive protein, low-density lipoprotein, high-density lipoprotein, and triglyceride single-nucleotide polymorphisms using a false discovery rate threshold <0.05. By conditioning on polymorphisms associated with the 4 phenotypes, we identified 55 loci associated with increased AD risk. We then conducted a meta-analysis of these 55 variants across 4 independent AD cohorts (total: n=29 054 AD cases and 114 824 healthy controls) and discovered 2 genome-wide significant variants on chromosome 4 (rs13113697; closest gene, HS3ST1; odds ratio=1.07; 95% confidence interval=1.05–1.11; P=2.86×10−8) and chromosome 10 (rs7920721; closest gene, ECHDC3; odds ratio=1.07; 95% confidence interval=1.04–1.11; P=3.38×10−8). We also found that gene expression of HS3ST1 and ECHDC3 was altered in AD brains compared with control brains. Conclusions—We demonstrate genetic overlap between AD, C-reactive protein, and plasma lipids. By conditioning on the genetic association with the cardiovascular phenotypes, we identify novel AD susceptibility loci, including 2 genome-wide significant variants conferring increased risk for AD.acceptedVersio

Identification of additional risk loci for stroke and small vessel disease: a meta-analysis of genome-wide association studies

BACKGROUND: Genetic determinants of stroke, the leading neurological cause of death and disability, are poorly understood and have seldom been explored in the general population. Our aim was to identify additional loci for stroke by doing a meta-analysis of genome-wide association studies. METHODS: For the discovery sample, we did a genome-wide analysis of common genetic variants associated with incident stroke risk in 18 population-based cohorts comprising 84 961 participants, of whom 4348 had stroke. Stroke diagnosis was ascertained and validated by the study investigators. Mean age at stroke ranged from 45·8 years to 76·4 years, and data collection in the studies took place between 1948 and 2013. We did validation analyses for variants yielding a significant association (at p<5 × 10(-6)) with all-stroke, ischaemic stroke, cardioembolic ischaemic stroke, or non-cardioembolic ischaemic stroke in the largest available cross-sectional studies (70 804 participants, of whom 19 816 had stroke). Summary-level results of discovery and follow-up stages were combined using inverse-variance weighted fixed-effects meta-analysis, and in-silico lookups were done in stroke subtypes. For genome-wide significant findings (at p<5 × 10(-8)), we explored associations with additional cerebrovascular phenotypes and did functional experiments using conditional (inducible) deletion of the probable causal gene in mice. We also studied the expression of orthologs of this probable causal gene and its effects on cerebral vasculature in zebrafish mutants. FINDINGS: We replicated seven of eight known loci associated with risk for ischaemic stroke, and identified a novel locus at chromosome 6p25 (rs12204590, near FOXF2) associated with risk of all-stroke (odds ratio [OR] 1·08, 95% CI 1·05-1·12, p=1·48 × 10(-8); minor allele frequency 21%). The rs12204590 stroke risk allele was also associated with increased MRI-defined burden of white matter hyperintensity-a marker of cerebral small vessel disease-in stroke-free adults (n=21 079; p=0·0025). Consistently, young patients (aged 2-32 years) with segmental deletions of FOXF2 showed an extensive burden of white matter hyperintensity. Deletion of Foxf2 in adult mice resulted in cerebral infarction, reactive gliosis, and microhaemorrhage. The orthologs of FOXF2 in zebrafish (foxf2b and foxf2a) are expressed in brain pericytes and mutant foxf2b(-/-) cerebral vessels show decreased smooth muscle cell and pericyte coverage. INTERPRETATION: We identified common variants near FOXF2 that are associated with increased stroke susceptibility. Epidemiological and experimental data suggest that FOXF2 mediates this association, potentially via differentiation defects of cerebral vascular mural cells. Further expression studies in appropriate human tissues, and further functional experiments with long follow-up periods are needed to fully understand the underlying mechanisms

Multi-ethnic genome-wide association study for atrial fibrillation

Atrial fibrillation (AF) affects more than 33 million individuals worldwide and has a complex heritability. We conducted the largest meta-analysis of genome-wide association studies (GWAS) for AF to date, consisting of more than half a million individuals, including 65,446 with AF. In total, we identified 97 loci significantly associated with AF, including 67 that were novel in a combined-ancestry analysis, and 3 that were novel in a European-specific analysis. We sought to identify AF-associated genes at the GWAS loci by performing RNA-sequencing and expression quantitative trait locus analyses in 101 left atrial samples, the most relevant tissue for AF. We also performed transcriptome-wide analyses that identified 57 AF-associated genes, 42 of which overlap with GWAS loci. The identified loci implicate genes enriched within cardiac developmental, electrophysiological, contractile and structural pathways. These results extend our understanding of the biological pathways underlying AF and may facilitate the development of therapeutics for AF

Multi-ancestry GWAS of the electrocardiographic PR interval identifies 202 loci underlying cardiac conduction

The electrocardiographic PR interval reflects atrioventricular conduction, and is associated with conduction abnormalities, pacemaker implantation, atrial fibrillation (AF), and cardiovascular mortality. Here we report a multi-ancestry (N=293,051) genome-wide association meta-analysis for the PR interval, discovering 202 loci of which 141 have not previously been reported. Variants at identified loci increase the percentage of heritability explained, from 33.5% to 62.6%. We observe enrichment for cardiac muscle developmental/contractile and cytoskeletal genes, highlighting key regulation processes for atrioventricular conduction. Additionally, 8 loci not previously reported harbor genes underlying inherited arrhythmic syndromes and/or cardiomyopathies suggesting a role for these genes in cardiovascular pathology in the general population. We show that polygenic predisposition to PR interval duration is an endophenotype for cardiovascular disease, including distal conduction disease, AF, and atrioventricular pre-excitation. These findings advance our understanding of the polygenic basis of cardiac conduction, and the genetic relationship between PR interval duration and cardiovascular disease. On the electrocardiogram, the PR interval reflects conduction from the atria to ventricles and also serves as risk indicator of cardiovascular morbidity and mortality. Here, the authors perform genome-wide meta-analyses for PR interval in multiple ancestries and identify 141 previously unreported genetic loci.Peer reviewe

Rare coding variants in PLCG2, ABI3, and TREM2 implicate microglial-mediated innate immunity in Alzheimer's disease

We identified rare coding variants associated with Alzheimer’s disease (AD) in a 3-stage case-control study of 85,133 subjects. In stage 1, 34,174 samples were genotyped using a whole-exome microarray. In stage 2, we tested associated variants (P<1×10-4) in 35,962 independent samples using de novo genotyping and imputed genotypes. In stage 3, an additional 14,997 samples were used to test the most significant stage 2 associations (P<5×10-8) using imputed genotypes. We observed 3 novel genome-wide significant (GWS) AD associated non-synonymous variants; a protective variant in PLCG2 (rs72824905/p.P522R, P=5.38×10-10, OR=0.68, MAFcases=0.0059, MAFcontrols=0.0093), a risk variant in ABI3 (rs616338/p.S209F, P=4.56×10-10, OR=1.43, MAFcases=0.011, MAFcontrols=0.008), and a novel GWS variant in TREM2 (rs143332484/p.R62H, P=1.55×10-14, OR=1.67, MAFcases=0.0143, MAFcontrols=0.0089), a known AD susceptibility gene. These protein-coding changes are in genes highly expressed in microglia and highlight an immune-related protein-protein interaction network enriched for previously identified AD risk genes. These genetic findings provide additional evidence that the microglia-mediated innate immune response contributes directly to AD development