Non-Traditional Approaches to Interrogating Genome-Wide SNP Data

The field of human molecular genetics has undergone a substantial technological transformation in the past decade, allowing researchers to identify and analyze genetic variation across the human genome with unprecedented depth and precision. A central goal in utilizing this technological advancement is to discover the genetic variation that underlies complex heritable traits and disorders. In recent years, much of the focus has been on large-scale genome- wide association studies (GWAS) in an attempt to identify effects of common single nucleotide polymorphisms (SNPs) on a phenotype. Progress in this arena, however, has been limited, as validated findings for most phenotypes represent only a small fraction of the variance attributed to genetics known from family and twin studies, leaving a large proportion of heritability to be explained. My research is in large part motivated by the issues surrounding GWAS, as additional methodological techniques and population genetic theory can help explain phenotypic variance unaccounted for by the traditional genetic association design. In my first study, I look at a case/control sample of bipolar disorder, examining how prior information from linkage studies can inform GWAS signals. The primary aim is to ease the burden of multiple-testing correction applied in GWAS using empirically informed weighting to tease out true signals supported by prior genetic evidence. In my second study, I determine the best practices to detect signatures of distant inbreeding, via runs of homozygosity, in genome-wide SNP data. The motivation for studying this phenomenon is due to the extensive evidence of inbreeding depression on fitness related traits, where the effects of recessive or partially recessive alleles are expressed. Using an extensive simulation design, I test multiple programs to determine the optimal method to identify runs of homozygosity caused by distant inbreeding. In my third and final study, I apply my work from the second study to a comprehensive dataset of cognitive measures to understand the extent to which distant inbreeding affects variation in general cognitive ability

Joint contributions of rare copy number variants and common SNPs to risk for schizophrenia

Objective: Both rare copy number variants (CNVs) and common single-nucleotide polymorphisms (SNPs) contribute to liability to schizophrenia, but their etiological relationship has not been fully elucidated. The authors evaluated an additive model whereby risk of schizophrenia requires less contribution from common SNPs in the presence of a rare CNV, and tested for interactions. Method: Genetic data from 21,094 case subjects with schizophrenia and 20,227 control subjects from the Psychiatric Genomics Consortium were examined. Three classes of rare CNVs were assessed: CNVs previously associated with schizophrenia, CNVs with large deletions ≥500 kb, and total CNV burden. The mean polygenic risk scores (PRSs) between study subjects with and without rare CNVs were compared, and joint effects of PRS and CNVs on schizophrenia liability were modeled by using logistic regression. Results: Schizophrenia case subjects carrying risk CNVs had a lower polygenic risk than case subjects without risk CNVs but a higher risk than control subjects. For case subjects carrying known risk CNVs, the PRS was diminished in proportion to the effect size of the CNV. The strongly associated 22q11.2 deletion required little added PRS to produce schizophrenia. Large deletions and increased CNV burden were also associated with lower polygenic risk in schizophrenia case subjects but not in control subjects or after removal of known risk CNV carriers. Conclusions: The authors found evidence for interactive effects of PRS and previously associated CNVs for risk for schizophrenia, and the results for large deletions and total CNV burden support an additive model. These findings offer insights into the genetic architecture of schizophrenia by illuminating how different established genetic risk factors act and interact to influence liability to schizophrenia

Runs of Homozygosity Implicate Autozygosity as a Schizophrenia Risk Factor

Autozygosity occurs when two chromosomal segments that are identical from a common ancestor are inherited from each parent. This occurs at high rates in the offspring of mates who are closely related (inbreeding), but also occurs at lower levels among the offspring of distantly related mates. Here, we use runs of homozygosity in genome-wide SNP data to estimate the proportion of the autosome that exists in autozygous tracts in 9,388 cases with schizophrenia and 12,456 controls. We estimate that the odds of schizophrenia increase by ∼17% for every 1% increase in genome-wide autozygosity. This association is not due to one or a few regions, but results from many autozygous segments spread throughout the genome, and is consistent with a role for multiple recessive or partially recessive alleles in the etiology of schizophrenia. Such a bias towards recessivity suggests that alleles that increase the risk of schizophrenia have been selected against over evolutionary time

Detecting autozygosity through runs of homozygosity: A comparison of three autozygosity detection algorithms

<p>Abstract</p> <p>Background</p> <p>A central aim for studying runs of homozygosity (ROHs) in genome-wide SNP data is to detect the effects of autozygosity (stretches of the two homologous chromosomes within the same individual that are identical by descent) on phenotypes. However, it is unknown which current ROH detection program, and which set of parameters within a given program, is optimal for differentiating ROHs that are truly autozygous from ROHs that are homozygous at the marker level but vary at unmeasured variants between the markers.</p> <p>Method</p> <p>We simulated 120 Mb of sequence data in order to know the true state of autozygosity. We then extracted common variants from this sequence to mimic the properties of SNP platforms and performed ROH analyses using three popular ROH detection programs, PLINK, GERMLINE, and BEAGLE. We varied detection thresholds for each program (e.g., prior probabilities, lengths of ROHs) to understand their effects on detecting known autozygosity.</p> <p>Results</p> <p>Within the optimal thresholds for each program, PLINK outperformed GERMLINE and BEAGLE in detecting autozygosity from distant common ancestors. PLINK's sliding window algorithm worked best when using SNP data pruned for linkage disequilibrium (LD).</p> <p>Conclusion</p> <p>Our results provide both general and specific recommendations for maximizing autozygosity detection in genome-wide SNP data, and should apply equally well to research on whole-genome autozygosity burden or to research on whether specific autozygous regions are predictive using association mapping methods.</p

Exome sequencing in bipolar disorder identifies AKAP11 as a risk gene shared with schizophrenia

We report results from the Bipolar Exome (BipEx) collaboration analysis of whole-exome sequencing of 13,933 patients with bipolar disorder (BD) matched with 14,422 controls. We find an excess of ultra-rare protein-truncating variants (PTVs) in patients with BD among genes under strong evolutionary constraint in both major BD subtypes. We find enrichment of ultra-rare PTVs within genes implicated from a recent schizophrenia exome meta-analysis (SCHEMA; 24,248 cases and 97,322 controls) and among binding targets of CHD8. Genes implicated from genome-wide association studies (GWASs) of BD, however, are not significantly enriched for ultra-rare PTVs. Combining gene-level results with SCHEMA, AKAP11 emerges as a definitive risk gene (odds ratio (OR) = 7.06, P = 2.83 × 10-9). At the protein level, AKAP-11 interacts with GSK3B, the hypothesized target of lithium, a primary treatment for BD. Our results lend support to BD's polygenicity, demonstrating a role for rare coding variation as a significant risk factor in BD etiology

Analysis of protein-coding genetic variation in 60,706 humans

Large-scale reference data sets of human genetic variation are critical for the medical and functional interpretation of DNA sequence changes. We describe the aggregation and analysis of high-quality exome (protein-coding region) sequence data for 60,706 individuals of diverse ethnicities generated as part of the Exome Aggregation Consortium (ExAC). This catalogue of human genetic diversity contains an average of one variant every eight bases of the exome, and provides direct evidence for the presence of widespread mutational recurrence. We have used this catalogue to calculate objective metrics of pathogenicity for sequence variants, and to identify genes subject to strong selection against various classes of mutation; identifying 3,230 genes with near-complete depletion of truncating variants with 72% having no currently established human disease phenotype. Finally, we demonstrate that these data can be used for the efficient filtering of candidate disease-causing variants, and for the discovery of human “knockout” variants in protein-coding genes

Rare coding variants in ten genes confer substantial risk for schizophrenia

Rare coding variation has historically provided the most direct connections between gene function and disease pathogenesis. By meta-analysing the whole exomes of 24,248 schizophrenia cases and 97,322 controls, we implicate ultra-rare coding variants (URVs) in 10 genes as conferring substantial risk for schizophrenia (odds ratios of 3-50, PPeer reviewe

No Reliable Association between Runs of Homozygosity and Schizophrenia in a Well-Powered Replication Study

It is well known that inbreeding increases the risk of recessive monogenic diseases, but it is less certain whether it contributes to the etiology of complex diseases such as schizophrenia. One way to estimate the effects of inbreeding is to examine the association between disease diagnosis and genome-wide autozygosity estimated using runs of homozygosity (ROH) in genome-wide single nucleotide polymorphism arrays. Using data for schizophrenia from the Psychiatric Genomics Consortium (n = 21,868), Keller et al. (2012) estimated that the odds of developing schizophrenia increased by approximately 17% for every additional percent of the genome that is autozygous (β = 16.1, CI(β) = [6.93, 25.7], Z = 3.44, p = 0.0006). Here we describe replication results from 22 independent schizophrenia case-control datasets from the Psychiatric Genomics Consortium (n = 39,830). Using the same ROH calling thresholds and procedures as Keller et al. (2012), we were unable to replicate the significant association between ROH burden and schizophrenia in the independent PGC phase II data, although the effect was in the predicted direction, and the combined (original + replication) dataset yielded an attenuated but significant relationship between Froh and schizophrenia (β = 4.86,CI(β) = [0.90,8.83],Z = 2.40,p = 0.02). Since Keller et al. (2012), several studies reported inconsistent association of ROH burden with complex traits, particularly in case-control data. These conflicting results might suggest that the effects of autozygosity are confounded by various factors, such as socioeconomic status, education, urbanicity, and religiosity, which may be associated with both real inbreeding and the outcome measures of interest

Contribution of copy number variants to schizophrenia from a genome-wide study of 41,321 subjects

Copy number variants (CNVs) have been strongly implicated in the genetic etiology of schizophrenia (SCZ). However, genome-wide investigation of the contribution of CNV to risk has been hampered by limited sample sizes. We sought to address this obstacle by applying a centralized analysis pipeline to a SCZ cohort of 21,094 cases and 20,227 controls. A global enrichment of CNV burden was observed in cases (OR=1.11, P=5.7×10−15), which persisted after excluding loci implicated in previous studies (OR=1.07, P=1.7 ×10−6). CNV burden was enriched for genes associated with synaptic function (OR = 1.68, P = 2.8 ×10−11) and neurobehavioral phenotypes in mouse (OR = 1.18, P= 7.3 ×10−5). Genome-wide significant evidence was obtained for eight loci, including 1q21.1, 2p16.3 (NRXN1), 3q29, 7q11.2, 15q13.3, distal 16p11.2, proximal 16p11.2 and 22q11.2. Suggestive support was found for eight additional candidate susceptibility and protective loci, which consisted predominantly of CNVs mediated by non-allelic homologous recombination