Pleiotropy and genetically inferred causality linking multisite chronic pain to substance use disorders.

Individuals suffering from chronic pain develop substance use disorders (SUDs) more often than others. Understanding the shared genetic influences underlying the comorbidity between chronic pain and SUDs will lead to a greater understanding of their biology. Genome-wide association statistics were obtained from the UK Biobank for multisite chronic pain (MCP, Neffective = 387,649) and from the Million Veteran Program and the Psychiatric Genomics Consortium meta-analyses for alcohol use disorder (AUD, Neffective = 296,974), cannabis use disorder (CanUD, Neffective = 161,053), opioid use disorder (OUD, Neffective = 57,120), and problematic tobacco use (PTU, Neffective = 270,120). SNP-based heritability was estimated for each of the traits and genetic correlation (rg) analyses were performed to assess MCP-SUD pleiotropy. Bidirectional Mendelian Randomization analyses evaluated possible causal relationships. Finally, to identify and characterize individual loci, we performed a genome-wide pleiotropy analysis and a brain-wide analysis using imaging phenotypes available from the UK Biobank. MCP was positively genetically correlated with AUD (rg = 0.26, p = 7.55 × 10-18), CanUD (rg = 0.37, p = 8.21 × 10-37), OUD (rg = 0.20, p = 1.50 × 10-3), and PTU (rg = 0.29, p = 8.53 × 10-12). Although the MR analyses supported bi-directional relationships, MCP had larger effects on AUD (pain-exposure: beta = 0.18, p = 8.21 × 10-4; pain-outcome: beta = 0.07, p = 0.018), CanUD (pain-exposure: beta = 0.58, p = 2.70 × 10-6; pain-outcome: beta = 0.05, p = 0.014) and PTU (pain-exposure: beta = 0.43, p = 4.16 × 10-8; pain-outcome: beta = 0.09, p = 3.05 × 10-6) than the reverse. The genome-wide analysis identified two SNPs pleiotropic between MCP and all SUD investigated: IHO1 rs7652746 (ppleiotropy = 2.69 × 10-8), and CADM2 rs1248857 (ppleiotropy = 1.98 × 10-5). In the brain-wide analysis, rs7652746 was associated with multiple cerebellum and amygdala imaging phenotypes. When analyzing MCP pleiotropy with each SUD separately, we found 25, 22, and 4 pleiotropic variants for AUD, CanUD, and OUD, respectively. To our knowledge, this is the first large-scale study to provide evidence of potential causal relationships and shared genetic mechanisms underlying MCP-SUD comorbidity.

Genome-wide significant risk loci for mood disorders in the Old Order Amish founder population.

Genome-wide association studies (GWAS) of mood disorders in large case-control cohorts have identified numerous risk loci, yet pathophysiological mechanisms remain elusive, primarily due to the very small effects of common variants. We sought to discover risk variants with larger effects by conducting a genome-wide association study of mood disorders in a founder population, the Old Order Amish (OOA, n = 1,672). Our analysis revealed four genome-wide significant risk loci, all of which were associated with >2-fold relative risk. Quantitative behavioral and neurocognitive assessments (n = 314) revealed effects of risk variants on sub-clinical depressive symptoms and information processing speed. Network analysis suggested that OOA-specific risk loci harbor novel risk-associated genes that interact with known neuropsychiatry-associated genes via gene interaction networks. Annotation of the variants at these risk loci revealed population-enriched, non-synonymous variants in two genes encoding neurodevelopmental transcription factors, CUX1 and CNOT1. Our findings provide insight into the genetic architecture of mood disorders and a substrate for mechanistic and clinical studies.

Identifying neurofunctional domains across substance use disorders.

Background: Substance use disorders (SUDs) are heterogeneous across multiple functional domains. Various frameworks posit that domains (e.g., executive function) contribute to the persistence of SUDs; however, the domains identified in different studies vary.Objectives: We used factor analysis to identify the underlying latent domains present in a large sample (N = 5,244, 55.8% male) with a variety of SUDs to yield findings more generalizable than studies with a narrower focus.Method: Participants (1,384 controls and 3,860 participants with one or more SUDs including alcohol, cocaine, cannabis, and/or opioid use disorders) completed the Semi-Structured Assessment for Drug Dependence and Alcoholism, the NEO Personality Inventory, and the Wisconsin Card Sorting Test. Exploratory factor analysis (EFA) and fit indices (root mean-squared error of approximation (RMSEA), Comparative Fit Index (CFI), and Tucker-Lewis Index (TLI)) were used to examine different latent variable models. A multiple indicators, multiple causes (MIMIC) approach-tested associations of the latent variables with sociodemographics, substance use, and a history of abuse/neglect.Results: A six-factor model (predominant alcohol, predominant cocaine, predominant opioid, externalizing, personality, and executive function) provided the best fit [RMSEA = 0.063 (90% CI 0.060, 0.066), CFI = 0.98, TLI = 0.96]. All factors were moderately correlated (coefficient = 0.25-0.55, p < .05) with the exception of executive function. MIMIC analysis revealed different patterns of associations (all p < .0001) with sociodemographics, substance use, and a history of abuse/neglect among the factors.Conclusions: The domains identified, particularly executive function, were parallel to those observed previously. These factors underscore the heterogeneous nature of SUDs and may be useful in developing more targeted clinical interventions.

Investigation of convergent and divergent genetic influences underlying schizophrenia and alcohol use disorder.

BACKGROUND: Alcohol use disorder (AUD) and schizophrenia (SCZ) frequently co-occur, and large-scale genome-wide association studies (GWAS) have identified significant genetic correlations between these disorders. METHODS: We used the largest published GWAS for AUD (total cases = 77 822) and SCZ (total cases = 46 827) to identify genetic variants that influence both disorders (with either the same or opposite direction of effect) and those that are disorder specific. RESULTS: We identified 55 independent genome-wide significant single nucleotide polymorphisms with the same direction of effect on AUD and SCZ, 8 with robust effects in opposite directions, and 98 with disorder-specific effects. We also found evidence for 12 genes whose pleiotropic associations with AUD and SCZ are consistent with mediation via gene expression in the prefrontal cortex. The genetic covariance between AUD and SCZ was concentrated in genomic regions functional in brain tissues (p = 0.001). CONCLUSIONS: Our findings provide further evidence that SCZ shares meaningful genetic overlap with AUD.

Genome-wide association study in individuals of European and African ancestry and multi-trait analysis of opioid use disorder identifies 19 independent genome-wide significant risk loci.

Despite the large toll of opioid use disorder (OUD), genome-wide association studies (GWAS) of OUD to date have yielded few susceptibility loci. We performed a large-scale GWAS of OUD in individuals of European (EUR) and African (AFR) ancestry, optimizing genetic informativeness by performing MTAG (Multi-trait analysis of GWAS) with genetically correlated substance use disorders (SUDs). Meta-analysis included seven cohorts: the Million Veteran Program, Psychiatric Genomics Consortium, iPSYCH, FinnGen, Partners Biobank, BioVU, and Yale-Penn 3, resulting in a total N = 639,063 (Ncases = 20,686;Neffective = 77,026) across ancestries. OUD cases were defined as having a lifetime OUD diagnosis, and controls as anyone not known to meet OUD criteria. We estimated SNP-heritability (h2SNP) and genetic correlations (rg). Based on genetic correlation, we performed MTAG on OUD, alcohol use disorder (AUD), and cannabis use disorder (CanUD). A leave-one-out polygenic risk score (PRS) analysis was performed to compare OUD and OUD-MTAG PRS as predictors of OUD case status in Yale-Penn 3. The EUR meta-analysis identified three genome-wide significant (GWS; p ≤ 5 × 10-8) lead SNPs-one at FURIN (rs11372849; p = 9.54 × 10-10) and two OPRM1 variants (rs1799971, p = 4.92 × 10-09; rs79704991, p = 1.11 × 10-08; r2 = 0.02). Rs1799971 (p = 4.91 × 10-08) and another OPRM1 variant (rs9478500; p = 1.95 × 10-08; r2 = 0.03) were identified in the cross-ancestry meta-analysis. Estimated h2SNP was 12.75%, with strong rg with CanUD (rg = 0.82; p = 1.14 × 10-47) and AUD (rg = 0.77; p = 6.36 × 10-78). The OUD-MTAG resulted in a GWAS Nequivalent = 128,748 and 18 independent GWS loci, some mapping to genes or gene regions that have previously been associated with psychiatric or addiction phenotypes. The OUD-MTAG PRS accounted for 3.81% of OUD variance (beta = 0.61;s.e. = 0.066; p = 2.00 × 10-16) compared to 2.41% (beta = 0.45; s.e. = 0.058; p = 2.90 × 10-13) explained by the OUD PRS. The current study identified OUD variant associations at OPRM1, single variant associations with FURIN, and 18 GWS associations in the OUD-MTAG. The genetic architecture of OUD is likely influenced by both OUD-specific loci and loci shared across SUDs.

Performance of polygenic risk scores for cancer prediction in a racially diverse academic biobank.

PURPOSE: Genome-wide association studies have identified hundreds of single nucleotide variations (formerly single nucleotide polymorphisms) associated with several cancers, but the predictive ability of polygenic risk scores (PRSs) is unclear, especially among non-Whites. METHODS: PRSs were derived from genome-wide significant single-nucleotide variations for 15 cancers in 20,079 individuals in an academic biobank. We evaluated the improvement in discriminatory accuracy by including cancer-specific PRS in patients of genetically-determined African and European ancestry. RESULTS: Among the individuals of European genetic ancestry, PRSs for breast, colon, melanoma, and prostate were significantly associated with their respective cancers. Among the individuals of African genetic ancestry, PRSs for breast, colon, prostate, and thyroid were significantly associated with their respective cancers. The area under the curve of the model consisting of age, sex, and principal components was 0.621 to 0.710, and it increased by 1% to 4% with the inclusion of PRS in individuals of European genetic ancestry. In individuals of African genetic ancestry, area under the curve was overall higher in the model without the PRS (0.723-0.810) but increased by <1% with the inclusion of PRS for most cancers. CONCLUSION: PRS moderately increased the ability to discriminate the cancer status in individuals of European but not African ancestry. Further large-scale studies are needed to identify ancestry-specific genetic factors in non-White populations to incorporate PRS into cancer risk assessment.

Ancestral haplotype reconstruction in endogamous populations using identity-by-descent.

In this work we develop a novel algorithm for reconstructing the genomes of ancestral individuals, given genotype or sequence data from contemporary individuals and an extended pedigree of family relationships. A pedigree with complete genomes for every individual enables the study of allele frequency dynamics and haplotype diversity across generations, including deviations from neutrality such as transmission distortion. When studying heritable diseases, ancestral haplotypes can be used to augment genome-wide association studies and track disease inheritance patterns. The building blocks of our reconstruction algorithm are segments of Identity-By-Descent (IBD) shared between two or more genotyped individuals. The method alternates between identifying a source for each IBD segment and assembling IBD segments placed within each ancestral individual. Unlike previous approaches, our method is able to accommodate complex pedigree structures with hundreds of individuals genotyped at millions of SNPs. We apply our method to an Old Order Amish pedigree from Lancaster, Pennsylvania, whose founders came to North America from Europe during the early 18th century. The pedigree includes 1338 individuals from the past 12 generations, 394 with genotype data. The motivation for reconstruction is to understand the genetic basis of diseases segregating in the family through tracking haplotype transmission over time. Using our algorithm thread, we are able to reconstruct an average of 224 ancestral individuals per chromosome. For these ancestral individuals, on average we reconstruct 79% of their haplotypes. We also identify a region on chromosome 16 that is difficult to reconstruct-we find that this region harbors a short Amish-specific copy number variation and the gene HYDIN. thread was developed for endogamous populations, but can be applied to any extensive pedigree with the recent generations genotyped. We anticipate that this type of practical ancestral reconstruction will become more common and necessary to understand rare and complex heritable diseases in extended families.

Genetic liability for substance use associated with medical comorbidities in electronic health records of African- and European-ancestry individuals.

Polygenic risk scores (PRS) represent an individual's summed genetic risk for a trait and can serve as biomarkers for disease. Less is known about the utility of PRS as a means to quantify genetic risk for substance use disorders (SUDs) than for many other traits. Nonetheless, the growth of large, electronic health record-based biobanks makes it possible to evaluate the association of SUD PRS with other traits. We calculated PRS for smoking initiation, alcohol use disorder (AUD), and opioid use disorder (OUD) using summary statistics from the Million Veteran Program sample. We then tested the association of each PRS with its primary phenotype in the Penn Medicine BioBank (PMBB) using all available genotyped participants of African or European ancestry (AFR and EUR, respectively) (N = 18,612). Finally, we conducted phenome-wide association analyses (PheWAS) separately by ancestry and sex to test for associations across disease categories. Tobacco use disorder was the most common SUD in the PMBB, followed by AUD and OUD, consistent with the population prevalence of these disorders. All PRS were associated with their primary phenotype in both ancestry groups. PheWAS results yielded cross-trait associations across multiple domains, including psychiatric disorders and medical conditions. SUD PRS were associated with their primary phenotypes; however, they are not yet predictive enough to be useful diagnostically. The cross-trait associations of the SUD PRS are indicative of a broader genetic liability. Future work should extend findings to additional population groups and for other substances of abuse.

Association of APOL1 Risk Alleles With Cardiovascular Disease in Blacks in the Million Veteran Program.

BACKGROUND: Approximately 13% of black individuals carry 2 copies of the apolipoprotein L1 (APOL1) risk alleles G1 or G2, which are associated with 1.5- to 2.5-fold increased risk of chronic kidney disease. There have been conflicting reports as to whether an association exists between APOL1 risk alleles and cardiovascular disease (CVD) that is independent of the effects of APOL1 on kidney disease. We sought to test the association of APOL1 G1/G2 alleles with coronary artery disease, peripheral artery disease, and stroke among black individuals in the Million Veteran Program. METHODS: We performed a time-to-event analysis of retrospective electronic health record data using Cox proportional hazard and competing-risks Fine and Gray subdistribution hazard models. The primary exposure was APOL1 risk allele status. The primary outcome was incident coronary artery disease among individuals without chronic kidney disease during the 12.5-year follow-up period. We separately analyzed the cross-sectional association of APOL1 risk allele status with lipid traits and 115 cardiovascular diseases using phenome-wide association. RESULTS: Among 30 903 black Million Veteran Program participants, 3941 (13%) carried the 2 APOL1 risk allele high-risk genotype. Individuals with normal kidney function at baseline with 2 risk alleles had slightly higher risk of developing coronary artery disease compared with those with no risk alleles (hazard ratio, 1.11 [95% CI, 1.01-1.21]; P=0.039). Similarly, modest associations were identified with incident stroke (hazard ratio, 1.20 [95% CI, 1.05-1.36; P=0.007) and peripheral artery disease (hazard ratio, 1.15 [95% CI, 1.01-1.29l; P=0.031). When both cardiovascular and renal outcomes were modeled, APOL1 was strongly associated with incident renal disease, whereas no significant association with the CVD end points could be detected. Cardiovascular phenome-wide association analyses did not identify additional significant associations with CVD subsets. CONCLUSIONS: APOL1 risk variants display a modest association with CVD, and this association is likely mediated by the known APOL1 association with chronic kidney disease.

Genomics-First Evaluation of Heart Disease Associated With Titin-Truncating Variants.

BACKGROUND: Truncating variants in the Titin gene (TTNtvs) are common in individuals with idiopathic dilated cardiomyopathy (DCM). However, a comprehensive genomics-first evaluation of the impact of TTNtvs in different clinical contexts, and the evaluation of modifiers such as genetic ancestry, has not been performed. METHODS: We reviewed whole exome sequence data for >71 000 individuals (61 040 from the Geisinger MyCode Community Health Initiative (2007 to present) and 10 273 from the PennMedicine BioBank (2013 to present) to identify anyone with TTNtvs. We further selected individuals with TTNtvs in exons highly expressed in the heart (proportion spliced in [PSI] >0.9). Using linked electronic health records, we evaluated associations of TTNtvs with diagnoses and quantitative echocardiographic measures, including subanalyses for individuals with and without DCM diagnoses. We also reviewed data from the Jackson Heart Study to validate specific analyses for individuals of African ancestry. RESULTS: Identified with a TTNtv in a highly expressed exon (hiPSI) were 1.2% individuals in PennMedicine BioBank and 0.6% at Geisinger. The presence of a hiPSI TTNtv was associated with increased odds of DCM in individuals of European ancestry (odds ratio [95% CI]: 18.7 [9.1-39.4] {PennMedicine BioBank} and 10.8 [7.0-16.0] {Geisinger}). hiPSI TTNtvs were not associated with DCM in individuals of African ancestry, despite a high DCM prevalence (odds ratio, 1.8 [0.2-13.7]; P=0.57). Among 244 individuals of European ancestry with DCM in PennMedicine BioBank, hiPSI TTNtv carriers had lower left ventricular ejection fraction (ß=-12%, P=3×10-7), and increased left ventricular diameter (ß=0.65 cm, P=9×10-3). In the Geisinger cohort, hiPSI TTNtv carriers without a cardiomyopathy diagnosis had more atrial fibrillation (odds ratio, 2.4 [1.6-3.6]) and heart failure (odds ratio, 3.8 [2.4-6.0]), and lower left ventricular ejection fraction (ß=-3.4%, P=1×10-7). CONCLUSIONS: Individuals of European ancestry with hiPSI TTNtv have an abnormal cardiac phenotype characterized by lower left ventricular ejection fraction, irrespective of the clinical manifestation of cardiomyopathy. Associations with arrhythmias, including atrial fibrillation, were observed even when controlling for cardiomyopathy diagnosis. In contrast, no association between hiPSI TTNtvs and DCM was discerned among individuals of African ancestry. Given these findings, clinical identification of hiPSI TTNtv carriers may alter clinical management strategies.

Phenome-wide association analysis suggests the APOL1 linked disease spectrum primarily drives kidney-specific pathways.

The relationship between commonly occurring genetic variants (G1 and G2) in the APOL1 gene in African Americans and different disease traits, such as kidney disease, cardiovascular disease, and pre-eclampsia, remains the subject of controversy. Here we took a genotype-first approach, a phenome-wide association study, to define the spectrum of phenotypes associated with APOL1 high-risk variants in 1,837 African American participants of Penn Medicine Biobank and 4,742 African American participants of Vanderbilt BioVU. In the Penn Medicine Biobank, outpatient creatinine measurement-based estimated glomerular filtration rate and multivariable regression models were used to evaluate the association between high-risk APOL1 status and renal outcomes. In meta-analysis of both cohorts, the strongest phenome-wide association study associations were for the high-risk APOL1 variants and diagnoses codes were highly significant for "kidney dialysis" (odds ratio 3.75) and "end stage kidney disease" (odds ratio 3.42). A number of phenotypes were associated with APOL1 high-risk genotypes in an analysis adjusted only for demographic variables. However, no associations were detected with non-renal phenotypes after controlling for chronic/end stage kidney disease status. Using calculated estimated glomerular filtration rate -based phenotype analysis in the Penn Medicine Biobank, APOL1 high-risk status was associated with prevalent chronic/end stage kidney disease /kidney transplant (odds ratio 2.27, 95% confidence interval 1.67-3.08). In high-risk participants, the estimated glomerular filtration rate was 15.4 mL/min/1.73m2; significantly lower than in low-risk participants. Thus, although APOL1 high-risk variants are associated with a range of phenotypes, the risks for other associated phenotypes appear much lower and in our dataset are driven by a primary effect on renal disease.

A genome-first approach to aggregating rare genetic variants in LMNA for association with electronic health record phenotypes.

PURPOSE: "Genome-first" approaches, in which genetic sequencing is agnostically linked to associated phenotypes, can enhance our understanding of rare variants' contributions to disease. Loss-of-function variants in LMNA cause a range of rare diseases, including cardiomyopathy. METHODS: We leveraged exome sequencing from 11,451 unselected individuals in the Penn Medicine Biobank to associate rare variants in LMNA with diverse electronic health record (EHR)-derived phenotypes. We used Rare Exome Variant Ensemble Learner (REVEL) to annotate rare missense variants, clustered predicted deleterious and loss-of-function variants into a "gene burden" (N = 72 individuals), and performed a phenome-wide association study (PheWAS). Major findings were replicated in DiscovEHR. RESULTS: The LMNA gene burden was significantly associated with primary cardiomyopathy (p = 1.78E-11) and cardiac conduction disorders (p = 5.27E-07). Most patients had not been clinically diagnosed with LMNA cardiomyopathy. We also noted an association with chronic kidney disease (p = 1.13E-06). Regression analyses on echocardiography and serum labs revealed that LMNA variant carriers had dilated cardiomyopathy and primary renal disease. CONCLUSION: Pathogenic LMNA variants are an underdiagnosed cause of cardiomyopathy. We also find that LMNA loss of function may be a primary cause of renal disease. Finally, we show the value of aggregating rare, annotated variants into a gene burden and using PheWAS to identify novel ontologies for pleiotropic human genes.

Differentiating Types of Self-Reported Alcohol Abstinence.

We contrast three types of abstinence: quit after alcohol associated problems (Q-AP), quit for other reasons (Q-OR), and lifetime abstainer (LTA). We summarized the characteristics of people living with HIV (PLWH), and matched uninfected individuals, by levels of alcohol use and types of abstinence. We then identified factors that differentiate abstinence and determined whether the association with an alcohol biomarker or a genetic polymorphism is improved by differentiating abstinence. Among abstainers, 34% of PLWH and 38% of uninfected were Q-AP; 53% and 53% were Q-OR; and 12% and 10% were LTA. Logistic regression models found smoking, alcohol, cocaine, and hepatitis C increased odds of Q-AP, whereas smoking and marijuana decreased odds of LTA. Differentiating types of abstinence improved association. Q-APs and LTAs can be readily differentiated by an alcohol biomarker and genetic polymorphism. Differentiating type of abstinence may enhance understanding of alcohol health effects.

Rare variants in the genetic background modulate cognitive and developmental phenotypes in individuals carrying disease-associated variants.

PURPOSE: To assess the contribution of rare variants in the genetic background toward variability of neurodevelopmental phenotypes in individuals with rare copy-number variants (CNVs) and gene-disruptive variants. METHODS: We analyzed quantitative clinical information, exome sequencing, and microarray data from 757 probands and 233 parents and siblings who carry disease-associated variants. RESULTS: The number of rare likely deleterious variants in functionally intolerant genes ("other hits") correlated with expression of neurodevelopmental phenotypes in probands with 16p12.1 deletion (n=23, p=0.004) and in autism probands carrying gene-disruptive variants (n=184, p=0.03) compared with their carrier family members. Probands with 16p12.1 deletion and a strong family history presented more severe clinical features (p=0.04) and higher burden of other hits compared with those with mild/no family history (p=0.001). The number of other hits also correlated with severity of cognitive impairment in probands carrying pathogenic CNVs (n=53) or de novo pathogenic variants in disease genes (n=290), and negatively correlated with head size among 80 probands with 16p11.2 deletion. These co-occurring hits involved known disease-associated genes such as SETD5, AUTS2, and NRXN1, and were enriched for cellular and developmental processes. CONCLUSION: Accurate genetic diagnosis of complex disorders will require complete evaluation of the genetic background even after a candidate disease-associated variant is identified.

Patterns and Correlates of Prescription Opioid Receipt Among US Veterans: A National, 18-Year Observational Cohort Study.

A better understanding of predisposition to transition to high-dose, long-term opioid therapy after initial opioid receipt could facilitate efforts to prevent opioid use disorder (OUD). We extracted data on 69,268 patients in the Veterans Aging Cohort Study who received any opioid prescription between 1998 and 2015. Using latent growth mixture modelling, we identified four distinguishable dose trajectories: low (53%), moderate (29%), escalating (13%), and rapidly escalating (5%). Compared to low dose trajectory, those in the rapidly escalating dose trajectory were proportionately more European-American (59% rapidly escalating vs. 38% low); had a higher prevalence of HIV (31% vs. 29%) and hepatitis C (18% vs. 12%); and during follow-up, had a higher incidence of OUD diagnoses (13% vs. 3%); were hospitalised more often [18.1/100 person-years (PYs) vs. 12.5/100 PY]; and had higher all-cause mortality (4.7/100 PY vs. 1.8/100 PY, all p < 0.0001). These measures can potentially be used in future prevention research, including genetic discovery.

Increased burden of deleterious variants in essential genes in autism spectrum disorder.

Autism spectrum disorder (ASD) is a heterogeneous, highly heritable neurodevelopmental syndrome characterized by impaired social interaction, communication, and repetitive behavior. It is estimated that hundreds of genes contribute to ASD. We asked if genes with a strong effect on survival and fitness contribute to ASD risk. Human orthologs of genes with an essential role in pre- and postnatal development in the mouse [essential genes (EGs)] are enriched for disease genes and under strong purifying selection relative to human orthologs of mouse genes with a known nonlethal phenotype [nonessential genes (NEGs)]. This intolerance to deleterious mutations, commonly observed haploinsufficiency, and the importance of EGs in development suggest a possible cumulative effect of deleterious variants in EGs on complex neurodevelopmental disorders. With a comprehensive catalog of 3,915 mammalian EGs, we provide compelling evidence for a stronger contribution of EGs to ASD risk compared with NEGs. By examining the exonic de novo and inherited variants from 1,781 ASD quartet families, we show a significantly higher burden of damaging mutations in EGs in ASD probands compared with their non-ASD siblings. The analysis of EGs in the developing brain identified clusters of coexpressed EGs implicated in ASD. Finally, we suggest a high-priority list of 29 EGs with potential ASD risk as targets for future functional and behavioral studies. Overall, we show that large-scale studies of gene function in model organisms provide a powerful approach for prioritization of genes and pathogenic variants identified by sequencing studies of human disease.

Association of the V122I Hereditary Transthyretin Amyloidosis Genetic Variant With Heart Failure Among Individuals of African or Hispanic/Latino Ancestry.

Importance: Hereditary transthyretin (TTR) amyloid cardiomyopathy (hATTR-CM) due to the TTR V122I variant is an autosomal-dominant disorder that causes heart failure in elderly individuals of African ancestry. The clinical associations of carrying the variant, its effect in other African ancestry populations including Hispanic/Latino individuals, and the rates of achieving a clinical diagnosis in carriers are unknown. Objective: To assess the association between the TTR V122I variant and heart failure and identify rates of hATTR-CM diagnosis among carriers with heart failure. Design, Setting, and Participants: Cross-sectional analysis of carriers and noncarriers of TTR V122I of African ancestry aged 50 years or older enrolled in the Penn Medicine Biobank between 2008 and 2017 using electronic health record data from 1996 to 2017. Case-control study in participants of African and Hispanic/Latino ancestry with and without heart failure in the Mount Sinai BioMe Biobank enrolled between 2007 and 2015 using electronic health record data from 2007 to 2018. Exposures: TTR V122I carrier status. Main Outcomes and Measures: The primary outcome was prevalent heart failure. The rate of diagnosis with hATTR-CM among TTR V122I carriers with heart failure was measured. Results: The cross-sectional cohort included 3724 individuals of African ancestry with a median age of 64 years (interquartile range, 57-71); 1755 (47%) were male, 2896 (78%) had a diagnosis of hypertension, and 753 (20%) had a history of myocardial infarction or coronary revascularization. There were 116 TTR V122I carriers (3.1%); 1121 participants (30%) had heart failure. The case-control study consisted of 2307 individuals of African ancestry and 3663 Hispanic/Latino individuals; the median age was 73 years (interquartile range, 68-80), 2271 (38%) were male, 4709 (79%) had a diagnosis of hypertension, and 1008 (17%) had a history of myocardial infarction or coronary revascularization. There were 1376 cases of heart failure. TTR V122I was associated with higher rates of heart failure (cross-sectional cohort: n = 51/116 TTR V122I carriers [44%], n = 1070/3608 noncarriers [30%], adjusted odds ratio, 1.7 [95% CI, 1.2-2.4], P = .006; case-control study: n = 36/1376 heart failure cases [2.6%], n = 82/4594 controls [1.8%], adjusted odds ratio, 1.8 [95% CI, 1.2-2.7], P = .008). Ten of 92 TTR V122I carriers with heart failure (11%) were diagnosed as having hATTR-CM; the median time from onset of symptoms to clinical diagnosis was 3 years. Conclusions and Relevance: Among individuals of African or Hispanic/Latino ancestry enrolled in 2 academic medical center-based biobanks, the TTR V122I genetic variant was significantly associated with heart failure.

Genomic view of bipolar disorder revealed by whole genome sequencing in a genetic isolate.

Bipolar disorder is a common, heritable mental illness characterized by recurrent episodes of mania and depression. Despite considerable effort to elucidate the genetic underpinnings of bipolar disorder, causative genetic risk factors remain elusive. We conducted a comprehensive genomic analysis of bipolar disorder in a large Old Order Amish pedigree. Microsatellite genotypes and high-density SNP-array genotypes of 388 family members were combined with whole genome sequence data for 50 of these subjects, comprising 18 parent-child trios. This study design permitted evaluation of candidate variants within the context of haplotype structure by resolving the phase in sequenced parent-child trios and by imputation of variants into multiple unsequenced siblings. Non-parametric and parametric linkage analysis of the entire pedigree as well as on smaller clusters of families identified several nominally significant linkage peaks, each of which included dozens of predicted deleterious variants. Close inspection of exonic and regulatory variants in genes under the linkage peaks using family-based association tests revealed additional credible candidate genes for functional studies and further replication in population-based cohorts. However, despite the in-depth genomic characterization of this unique, large and multigenerational pedigree from a genetic isolate, there was no convergence of evidence implicating a particular set of risk loci or common pathways. The striking haplotype and locus heterogeneity we observed has profound implications for the design of studies of bipolar and other related disorders.

Copy number variants encompassing Mendelian disease genes in a large multigenerational family segregating bipolar disorder.

BACKGROUND: Bipolar affective disorder (BP) is a common, highly heritable psychiatric disorder characterized by periods of depression and mania. Using dense SNP genotype data, we characterized CNVs in 388 members of an Old Order Amish Pedigree with bipolar disorder. We identified CNV regions arising from common ancestral mutations by utilizing the pedigree information. By combining this analysis with whole genome sequence data in the same individuals, we also explored the role of compound heterozygosity. RESULTS: Here we describe 541 inherited CNV regions, of which 268 are rare in a control population of European origin but present in a large number of Amish individuals. In addition, we highlight a set of CNVs found at higher frequencies in BP individuals, and within genes known to play a role in human development and disease. As in prior reports, we find no evidence for an increased burden of CNVs in BP individuals, but we report a trend towards a higher burden of CNVs in known Mendelian disease loci in bipolar individuals (BPI and BPII, p = 0.06). CONCLUSIONS: We conclude that CNVs may be contributing factors in the phenotypic presentation of mood disorders and co-morbid medical conditions in this family. These results reinforce the hypothesis of a complex genetic architecture underlying BP disorder, and suggest that the role of CNVs should continue to be investigated in BP data sets.

Application of polygenic scores to a deeply phenotyped sample enriched for substance use disorders reveals extensive pleiotropy with psychiatric and medical traits.

Co-occurring psychiatric, medical, and substance use disorders (SUDs) are common, but the complex pathways leading to such comorbidities are poorly understood. A greater understanding of genetic influences on this phenomenon could inform precision medicine efforts. We used the Yale-Penn dataset, a cross-sectional sample enriched for individuals with SUDs, to examine pleiotropic effects of genetic liability for psychiatric and medical traits. Participants completed an in-depth interview that provides information on demographics, environment, medical illnesses, and psychiatric and SUDs. Polygenic scores (PGS) for psychiatric disorders and medical traits were calculated in European-ancestry (EUR; n=5,691) participants and, when discovery datasets were available, for African-ancestry (AFR; n=4,918) participants. Phenome-wide association studies (PheWAS) were then conducted. In AFR participants, the only PGS with significant associations was bipolar disorder (BD), all of which were with substance use phenotypes. In EUR participants, PGS for major depressive disorder (MDD), generalized anxiety disorder (GAD), post-traumatic stress disorder (PTSD), schizophrenia (SCZ), body mass index (BMI), coronary artery disease (CAD), and type 2 diabetes (T2D) all showed significant associations, the majority of which were with phenotypes in the substance use categories. For instance, PGS MDD was associated with over 200 phenotypes, 15 of which were depression-related (e.g., depression criterion count), 55 of which were other psychiatric phenotypes, and 126 of which were substance use phenotypes; and PGS BMI was associated with 138 phenotypes, 105 of which were substance related. Genetic liability for psychiatric and medical traits is associated with numerous phenotypes across multiple categories, indicative of the broad genetic liability of these traits.