Genome-wide Association Studies in Ancestrally Diverse Populations: Opportunities, Methods, Pitfalls, and Recommendations.

Genome-wide association studies (GWASs) have focused primarily on populations of European descent, but it is essential that diverse populations become better represented. Increasing diversity among study participants will advance our understanding of genetic architecture in all populations and ensure that genetic research is broadly applicable. To facilitate and promote research in multi-ancestry and admixed cohorts, we outline key methodological considerations and highlight opportunities, challenges, solutions, and areas in need of development. Despite the perception that analyzing genetic data from diverse populations is difficult, it is scientifically and ethically imperative, and there is an expanding analytical toolbox to do it well.

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, P < 2.14 × 10-6) and 32 genes at a false discovery rate of <5%. These genes have the greatest expression in central nervous system neurons and have diverse molecular functions that include the formation, structure and function of the synapse. The associations of the NMDA (N-methyl-D-aspartate) receptor subunit GRIN2A and AMPA (α-amino-3-hydroxy-5-methyl-4-isoxazole propionic acid) receptor subunit GRIA3 provide support for dysfunction of the glutamatergic system as a mechanistic hypothesis in the pathogenesis of schizophrenia. We observe an overlap of rare variant risk among schizophrenia, autism spectrum disorders1, epilepsy and severe neurodevelopmental disorders2, although different mutation types are implicated in some shared genes. Most genes described here, however, are not implicated in neurodevelopment. We demonstrate that genes prioritized from common variant analyses of schizophrenia are enriched in rare variant risk3, suggesting that common and rare genetic risk factors converge at least partially on the same underlying pathogenic biological processes. Even after excluding significantly associated genes, schizophrenia cases still carry a substantial excess of URVs, which indicates that more risk genes await discovery using this approach.

Mapping genomic loci implicates genes and synaptic biology in schizophrenia.

Schizophrenia has a heritability of 60-80%1, much of which is attributable to common risk alleles. Here, in a two-stage genome-wide association study of up to 76,755 individuals with schizophrenia and 243,649 control individuals, we report common variant associations at 287 distinct genomic loci. Associations were concentrated in genes that are expressed in excitatory and inhibitory neurons of the central nervous system, but not in other tissues or cell types. Using fine-mapping and functional genomic data, we identify 120 genes (106 protein-coding) that are likely to underpin associations at some of these loci, including 16 genes with credible causal non-synonymous or untranslated region variation. We also implicate fundamental processes related to neuronal function, including synaptic organization, differentiation and transmission. Fine-mapped candidates were enriched for genes associated with rare disruptive coding variants in people with schizophrenia, including the glutamate receptor subunit GRIN2A and transcription factor SP4, and were also enriched for genes implicated by such variants in neurodevelopmental disorders. We identify biological processes relevant to schizophrenia pathophysiology; show convergence of common and rare variant associations in schizophrenia and neurodevelopmental disorders; and provide a resource of prioritized genes and variants to advance mechanistic studies.

BCFtools/liftover: an accurate and comprehensive tool to convert genetic variants across genome assemblies.

MOTIVATION: Many genetics studies report results tied to genomic coordinates of a legacy genome assembly. However, as assemblies are updated and improved, researchers are faced with either realigning raw sequence data using the updated coordinate system or converting legacy datasets to the updated coordinate system to be able to combine results with newer datasets. Currently available tools to perform the conversion of genetic variants have numerous shortcomings, including poor support for indels and multi-allelic variants, that lead to a higher rate of variants being dropped or incorrectly converted. As a result, many researchers continue to work with and publish using legacy genomic coordinates. RESULTS: Here we present BCFtools/liftover, a tool to convert genomic coordinates across genome assemblies for variants encoded in the variant call format with improved support for indels represented by different reference alleles across genome assemblies and full support for multi-allelic variants. It further supports variant annotation fields updates whenever the reference allele changes across genome assemblies. The tool has the lowest rate of variants being dropped with an order of magnitude less indels dropped or incorrectly converted and is an order of magnitude faster than other tools typically used for the same task. It is particularly suited for converting variant callsets from large cohorts to novel telomere-to-telomere assemblies as well as summary statistics from genome-wide association studies tied to legacy genome assemblies. AVAILABILITY AND IMPLEMENTATION: The tool is written in C and freely available under the MIT open source license as a BCFtools plugin available at http://github.com/freeseek/score.

Correlates of suicidal behaviors and genetic risk among United States veterans with schizophrenia or bipolar I disorder.

Persons diagnosed with schizophrenia (SCZ) or bipolar I disorder (BPI) are at high risk for self-injurious behavior, suicidal ideation, and suicidal behaviors (SB). Characterizing associations between diagnosed health problems, prior pharmacological treatments, and polygenic scores (PGS) has potential to inform risk stratification. We examined self-reported SB and ideation using the Columbia Suicide Severity Rating Scale (C-SSRS) among 3,942 SCZ and 5,414 BPI patients receiving care within the Veterans Health Administration (VHA). These cross-sectional data were integrated with electronic health records (EHRs), and compared across lifetime diagnoses, treatment histories, follow-up screenings, and mortality data. PGS were constructed using available genomic data for related traits. Genome-wide association studies were performed to identify and prioritize specific loci. Only 20% of the veterans who reported SB had a corroborating ICD-9/10 EHR code. Among those without prior SB, more than 20% reported new-onset SB at follow-up. SB were associated with a range of additional clinical diagnoses, and with treatment with specific classes of psychotropic medications (e.g., antidepressants, antipsychotics, etc.). PGS for externalizing behaviors, smoking initiation, suicide attempt, and major depressive disorder were associated with SB. The GWAS for SB yielded no significant loci. Among individuals with a diagnosed mental illness, self-reported SB were strongly associated with clinical variables across several EHR domains. Analyses point to sequelae of substance-related and psychiatric comorbidities as strong correlates of prior and subsequent SB. Nonetheless, past SB was frequently not documented in health records, underscoring the value of regular screening with direct, in-person assessments, especially among high-risk individuals.

Modeling the longitudinal changes of ancestry diversity in the Million Veteran Program.

BACKGROUND: The Million Veteran Program (MVP) participants represent 100 years of US history, including significant social and demographic changes over time. Our study assessed two aspects of the MVP: (i) longitudinal changes in population diversity and (ii) how these changes can be accounted for in genome-wide association studies (GWAS). To investigate these aspects, we divided MVP participants into five birth cohorts (N-range = 123,888 [born from 1943 to 1947] to 136,699 [born from 1948 to 1953]). RESULTS: Ancestry groups were defined by (i) HARE (harmonized ancestry and race/ethnicity) and (ii) a random-forest clustering approach using the 1000 Genomes Project and the Human Genome Diversity Project (1kGP + HGDP) reference panels (77 world populations representing six continental groups). In these groups, we performed GWASs of height, a trait potentially affected by population stratification. Birth cohorts demonstrate important trends in ancestry diversity over time. More recent HARE-assigned Europeans, Africans, and Hispanics had lower European ancestry proportions than older birth cohorts (0.010 < Cohen's d < 0.259, p < 7.80 × 10-4). Conversely, HARE-assigned East Asians showed an increase in European ancestry proportion over time. In GWAS of height using HARE assignments, genomic inflation due to population stratification was prevalent across all birth cohorts (linkage disequilibrium score regression intercept = 1.08 ± 0.042). The 1kGP + HGDP-based ancestry assignment significantly reduced the population stratification (mean intercept reduction = 0.045 ± 0.007, p < 0.05) confounding in the GWAS statistics. CONCLUSIONS: This study provides a characterization of ancestry diversity of the MVP cohort over time and compares two strategies to infer genetically defined ancestry groups by assessing differences in controlling population stratification in genome-wide association studies.

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.

Latin American Trans-ancestry INitiative for OCD genomics (LATINO): Study protocol.

Obsessive-compulsive disorder (OCD) is a debilitating psychiatric disorder. Worldwide, its prevalence is ~2% and its etiology is mostly unknown. Identifying biological factors contributing to OCD will elucidate underlying mechanisms and might contribute to improved treatment outcomes. Genomic studies of OCD are beginning to reveal long-sought risk loci, but >95% of the cases currently in analysis are of homogenous European ancestry. If not addressed, this Eurocentric bias will result in OCD genomic findings being more accurate for individuals of European ancestry than other ancestries, thereby contributing to health disparities in potential future applications of genomics. In this study protocol paper, we describe the Latin American Trans-ancestry INitiative for OCD genomics (LATINO, https://www.latinostudy.org). LATINO is a new network of investigators from across Latin America, the United States, and Canada who have begun to collect DNA and clinical data from 5000 richly phenotyped OCD cases of Latin American ancestry in a culturally sensitive and ethical manner. In this project, we will utilize trans-ancestry genomic analyses to accelerate the identification of OCD risk loci, fine-map putative causal variants, and improve the performance of polygenic risk scores in diverse populations. We will also capitalize on rich clinical data to examine the genetics of treatment response, biologically plausible OCD subtypes, and symptom dimensions. Additionally, LATINO will help elucidate the diversity of the clinical presentations of OCD across cultures through various trainings developed and offered in collaboration with Latin American investigators. We believe this study will advance the important goal of global mental health discovery and equity.

Characterisation of age and polarity at onset in bipolar disorder.

BACKGROUND: Studying phenotypic and genetic characteristics of age at onset (AAO) and polarity at onset (PAO) in bipolar disorder can provide new insights into disease pathology and facilitate the development of screening tools. AIMS: To examine the genetic architecture of AAO and PAO and their association with bipolar disorder disease characteristics. METHOD: Genome-wide association studies (GWASs) and polygenic score (PGS) analyses of AAO (n = 12 977) and PAO (n = 6773) were conducted in patients with bipolar disorder from 34 cohorts and a replication sample (n = 2237). The association of onset with disease characteristics was investigated in two of these cohorts. RESULTS: Earlier AAO was associated with a higher probability of psychotic symptoms, suicidality, lower educational attainment, not living together and fewer episodes. Depressive onset correlated with suicidality and manic onset correlated with delusions and manic episodes. Systematic differences in AAO between cohorts and continents of origin were observed. This was also reflected in single-nucleotide variant-based heritability estimates, with higher heritabilities for stricter onset definitions. Increased PGS for autism spectrum disorder (ß = -0.34 years, s.e. = 0.08), major depression (ß = -0.34 years, s.e. = 0.08), schizophrenia (ß = -0.39 years, s.e. = 0.08), and educational attainment (ß = -0.31 years, s.e. = 0.08) were associated with an earlier AAO. The AAO GWAS identified one significant locus, but this finding did not replicate. Neither GWAS nor PGS analyses yielded significant associations with PAO. CONCLUSIONS: AAO and PAO are associated with indicators of bipolar disorder severity. Individuals with an earlier onset show an increased polygenic liability for a broad spectrum of psychiatric traits. Systematic differences in AAO across cohorts, continents and phenotype definitions introduce significant heterogeneity, affecting analyses.

Leveraging genome-wide data to investigate differences between opioid use vs. opioid dependence in 41,176 individuals from the Psychiatric Genomics Consortium.

To provide insights into the biology of opioid dependence (OD) and opioid use (i.e., exposure, OE), we completed a genome-wide analysis comparing 4503 OD cases, 4173 opioid-exposed controls, and 32,500 opioid-unexposed controls, including participants of European and African descent (EUR and AFR, respectively). Among the variants identified, rs9291211 was associated with OE (exposed vs. unexposed controls; EUR z = -5.39, p = 7.2 × 10-8). This variant regulates the transcriptomic profiles of SLC30A9 and BEND4 in multiple brain tissues and was previously associated with depression, alcohol consumption, and neuroticism. A phenome-wide scan of rs9291211 in the UK Biobank (N > 360,000) found association of this variant with propensity to use dietary supplements (p = 1.68 × 10-8). With respect to the same OE phenotype in the gene-based analysis, we identified SDCCAG8 (EUR + AFR z = 4.69, p = 10-6), which was previously associated with educational attainment, risk-taking behaviors, and schizophrenia. In addition, rs201123820 showed a genome-wide significant difference between OD cases and unexposed controls (AFR z = 5.55, p = 2.9 × 10-8) and a significant association with musculoskeletal disorders in the UK Biobank (p = 4.88 × 10-7). A polygenic risk score (PRS) based on a GWAS of risk-tolerance (n = 466,571) was positively associated with OD (OD vs. unexposed controls, p = 8.1 × 10-5; OD cases vs. exposed controls, p = 0.054) and OE (exposed vs. unexposed controls, p = 3.6 × 10-5). A PRS based on a GWAS of neuroticism (n = 390,278) was positively associated with OD (OD vs. unexposed controls, p = 3.2 × 10-5; OD vs. exposed controls, p = 0.002) but not with OE (p = 0.67). Our analyses highlight the difference between dependence and exposure and the importance of considering the definition of controls in studies of addiction.

Contributions of common genetic variants to risk of schizophrenia among individuals of African and Latino ancestry.

Schizophrenia is a common, chronic and debilitating neuropsychiatric syndrome affecting tens of millions of individuals worldwide. While rare genetic variants play a role in the etiology of schizophrenia, most of the currently explained liability is within common variation, suggesting that variation predating the human diaspora out of Africa harbors a large fraction of the common variant attributable heritability. However, common variant association studies in schizophrenia have concentrated mainly on cohorts of European descent. We describe genome-wide association studies of 6152 cases and 3918 controls of admixed African ancestry, and of 1234 cases and 3090 controls of Latino ancestry, representing the largest such study in these populations to date. Combining results from the samples with African ancestry with summary statistics from the Psychiatric Genomics Consortium (PGC) study of schizophrenia yielded seven newly genome-wide significant loci, and we identified an additional eight loci by incorporating the results from samples with Latino ancestry. Leveraging population differences in patterns of linkage disequilibrium, we achieve improved fine-mapping resolution at 22 previously reported and 4 newly significant loci. Polygenic risk score profiling revealed improved prediction based on trans-ancestry meta-analysis results for admixed African (Nagelkerke's R2 = 0.032; liability R2 = 0.017; P < 10-52), Latino (Nagelkerke's R2 = 0.089; liability R2 = 0.021; P < 10-58), and European individuals (Nagelkerke's R2 = 0.089; liability R2 = 0.037; P < 10-113), further highlighting the advantages of incorporating data from diverse human populations.

Evidence of shared familial factors influencing neurocognitive endophenotypes in adult- and childhood-onset schizophrenia.

BACKGROUND: The aggregation of neurocognitive deficits among the non-psychotic first-degree relatives of adult- and childhood-onset schizophrenia patients suggests that there may be a common etiology for these deficits in childhood- and adult-onset illness. However, there is considerable heterogeneity in the presentation of neurobiological abnormalities, and whether there are differences in the extent of familial transmission for specific domains of cognitive function has not been systematically addressed. METHODS: We employed variance components analysis, as implemented in SOLAR-Eclipse, to evaluate the evidence of familial transmission for empirically derived composite scores representing attention, working memory, verbal learning, verbal retention, and memory for faces. We contrast estimates for adult- and childhood-onset schizophrenia families and matched community control pedigrees, and compare our findings to previous reports based on analogous neurocognitive assessments. RESULTS: We observed varying degrees of familial transmission; attention and working memory yielded comparable, significant estimates for adult-onset and community control pedigrees; verbal learning was significant for childhood-onset and community control pedigrees; and facial memory demonstrated significant familial transmission only for childhood-onset schizophrenia. Model-fitting analyses indicated significant differences in familiality between adult- and childhood-onset schizophrenia for attention, working memory, and verbal learning. CONCLUSIONS: By comprehensively assessing a wide range of neurocognitive domains in adult- and childhood-onset schizophrenia families, we provide additional support for specific neurocognitive domains as schizophrenia endophenotypes. Whereas comparable estimates of familial transmission for certain dimensions of cognitive functioning support a shared etiology of adult- and childhood-onset neurocognitive function, observed differences may be taken as preliminary evidence of partially divergent multifactorial architectures.

Pathway-based polygene risk for severe depression implicates drug metabolism in CONVERGE.

BACKGROUND: The Psychiatric Genomics Consortium (PGC) has made major advances in the molecular etiology of MDD, confirming that MDD is highly polygenic. Pathway enrichment results from PGC meta-analyses can also be used to help inform molecular drug targets. Prior to any knowledge of molecular biomarkers for MDD, drugs targeting molecular pathways (MPs) proved successful in treating MDD. It is possible that examining polygenicity within specific MPs implicated in MDD can further refine molecular drug targets. METHODS: Using a large case-control GWAS based on low-coverage whole genome sequencing (N = 10 640) in Han Chinese women, we derived polygenic risk scores (PRS) for MDD and for MDD specific to each of over 300 MPs previously shown to be relevant to psychiatric diagnoses. We then identified sets of PRSs, accounting for critical covariates, significantly predictive of case status. RESULTS: Over and above global MDD polygenic risk, polygenic risk within the GO: 0017144 drug metabolism pathway significantly predicted recurrent depression after multiple testing correction. Secondary transcriptomic analysis suggests that among genes in this pathway, CYP2C19 (family of Cytochrome P450) and CBR1 (Carbonyl Reductase 1) might be most relevant to MDD. Within the cases, pathway-based risk was additionally associated with age at onset of MDD. CONCLUSIONS: Results indicate that pathway-based risk might inform etiology of recurrent major depression. Future research should examine whether polygenicity of the drug metabolism gene pathway has any association with clinical presentation or treatment response. We discuss limitations to the generalizability of these preliminary findings, and urge replication in future research.

Genome-wide association study of cognitive performance in U.S. veterans with schizophrenia or bipolar disorder.

Cognitive impairment is a frequent and serious problem in patients with various forms of severe mental illnesses (SMI), including schizophrenia (SZ) and bipolar disorder (BP). Recent research suggests genetic links to several cognitive phenotypes in both SMI and in the general population. Our goal in this study was to identify potential genomic signatures of cognitive functioning in veterans with severe mental illness and compare them to previous findings for cognition across different populations. Veterans Affairs (VA) Cooperative Studies Program (CSP) Study #572 evaluated cognitive and functional capacity measures among SZ and BP patients. In conjunction with the VA Million Veteran Program, 3,959 European American (1,095 SZ, 2,864 BP) and 2,601 African American (1,095 SZ, 2,864 BP) patients were genotyped using a custom Affymetrix Axiom Biobank array. We performed a genome-wide association study of global cognitive functioning, constructed polygenic scores for SZ and cognition in the general population, and examined genetic correlations with 2,626 UK Biobank traits. Although no single locus attained genome-wide significance, observed allelic effects were strongly consistent with previous studies. We observed robust associations between global cognitive functioning and polygenic scores for cognitive performance, intelligence, and SZ risk. We also identified significant genetic correlations with several cognition-related traits in UK Biobank. In a diverse cohort of U.S. veterans with SZ or BP, we demonstrate broad overlap of common genetic effects on cognition in the general population, and find that greater polygenic loading for SZ risk is associated with poorer cognitive performance.

Polygenic Risk Score Prediction of Alcohol Dependence Symptoms Across Population-Based and Clinically Ascertained Samples.

BACKGROUND: Despite consistent evidence of the heritability of alcohol use disorders (AUDs), few specific genes with an etiological role have been identified. It is likely that AUDs are highly polygenic; however, the etiological pathways and genetic variants involved may differ between populations. The aim of this study was thus to evaluate whether aggregate genetic risk for AUDs differed between clinically ascertained and population-based epidemiological samples. METHODS: Four independent samples were obtained: 2 from unselected birth cohorts (Avon Longitudinal Study of Parents and Children [ALSPAC], N = 4,304; FinnTwin12 [FT12], N = 1,135) and 2 from families densely affected with AUDs, identified from treatment-seeking patients (Collaborative Study on the Genetics of Alcoholism, N = 2,097; Irish Affected Sib Pair Study of Alcohol Dependence, N = 706). AUD symptoms were assessed with clinical interviews, and participants of European ancestry were genotyped. Genomewide association was conducted separately in each sample, and the resulting association weights were used to create polygenic risk scores in each of the other samples (12 total discovery-validation pairs), and from meta-analyses within sample type. We then tested how well these aggregate genetic scores predicted AUD outcomes within and across sample types. RESULTS: Polygenic scores derived from 1 population-based sample (ALSPAC) significantly predicted AUD symptoms in another population-based sample (FT12), but not in either clinically ascertained sample. Trend-level associations (uncorrected p < 0.05) were found for polygenic score predictions within sample types but no or negative predictions across sample types. Polygenic scores accounted for 0 to 1% of the variance in AUD symptoms. CONCLUSIONS: Though preliminary, these results provide suggestive evidence of differences in the genetic etiology of AUDs based on sample characteristics such as treatment-seeking status, which may index other important clinical or demographic factors that moderate genetic influences. Although the variance accounted for by genomewide polygenic scores remains low, future studies could improve gene identification efforts by amassing very large samples, or reducing genetic heterogeneity by informing analyses with other phenotypic information such as sample characteristics. Multiple complementary approaches may be needed to make progress in gene identification for this complex disorder.

Genomewide Association Study of Alcohol Dependence Identifies Risk Loci Altering Ethanol-Response Behaviors in Model Organisms.

BACKGROUND: Alcohol dependence (AD) shows evidence for genetic liability, but genes influencing risk remain largely unidentified. METHODS: We conducted a genomewide association study in 706 related AD cases and 1,748 unscreened population controls from Ireland. We sought replication in 15,496 samples of European descent. We used model organisms (MOs) to assess the role of orthologous genes in ethanol (EtOH)-response behaviors. We tested 1 primate-specific gene for expression differences in case/control postmortem brain tissue. RESULTS: We detected significant association in COL6A3 and suggestive association in 2 previously implicated loci, KLF12 and RYR3. None of these signals are significant in replication. A suggestive signal in the long noncoding RNA LOC339975 is significant in case:control meta-analysis, but not in a population sample. Knockdown of a COL6A3 ortholog in Caenorhabditis elegans reduced EtOH sensitivity. Col6a3 expression correlated with handling-induced convulsions in mice. Loss of function of the KLF12 ortholog in C. elegans impaired development of acute functional tolerance (AFT). Klf12 expression correlated with locomotor activation following EtOH injection in mice. Loss of function of the RYR3 ortholog reduced EtOH sensitivity in C. elegans and rapid tolerance in Drosophila. The ryanodine receptor antagonist dantrolene reduced motivation to self-administer EtOH in rats. Expression of LOC339975 does not differ between cases and controls but is reduced in carriers of the associated rs11726136 allele in nucleus accumbens (NAc). CONCLUSIONS: We detect association between AD and COL6A3, KLF12, RYR3, and LOC339975. Despite nonreplication of COL6A3, KLF12, and RYR3 signals, orthologs of these genes influence behavioral response to EtOH in MOs, suggesting potential involvement in human EtOH response and AD liability. The associated LOC339975 allele may influence gene expression in human NAc. Although the functions of long noncoding RNAs are poorly understood, there is mounting evidence implicating these genes in multiple brain functions and disorders.

CHRONICITY OF DEPRESSION AND MOLECULAR MARKERS IN A LARGE SAMPLE OF HAN CHINESE WOMEN.

BACKGROUND: Major depressive disorder (MDD) has been associated with changes in mean telomere length and mitochondrial DNA (mtDNA) copy number. This study investigates if clinical features of MDD differentially impact these molecular markers. METHODS: Data from a large, clinically ascertained sample of Han Chinese women with recurrent MDD were used to examine whether symptom presentation, severity, and comorbidity were related to salivary telomere length and/or mtDNA copy number (maximum N = 5,284 for both molecular and phenotypic data). RESULTS: Structural equation modeling revealed that duration of longest episode was positively associated with mtDNA copy number, while earlier age of onset of most severe episode and a history of dysthymia were associated with shorter telomeres. Other factors, such as symptom presentation, family history of depression, and other comorbid internalizing disorders, were not associated with these molecular markers. CONCLUSIONS: Chronicity of depressive symptoms is related to more pronounced telomere shortening and increased mtDNA copy number among individuals with a history of recurrent MDD. As these molecular markers have previously been implicated in physiological aging and morbidity, individuals who experience prolonged depressive symptoms are potentially at greater risk of adverse medical outcomes.

Genome-wide association study reveals greater polygenic loading for schizophrenia in cases with a family history of illness.

Genome-wide association studies (GWAS) of schizophrenia have yielded more than 100 common susceptibility variants, and strongly support a substantial polygenic contribution of a large number of small allelic effects. It has been hypothesized that familial schizophrenia is largely a consequence of inherited rather than environmental factors. We investigated the extent to which familiality of schizophrenia is associated with enrichment for common risk variants detectable in a large GWAS. We analyzed single nucleotide polymorphism (SNP) data for cases reporting a family history of psychotic illness (N = 978), cases reporting no such family history (N = 4,503), and unscreened controls (N = 8,285) from the Psychiatric Genomics Consortium (PGC1) study of schizophrenia. We used a multinomial logistic regression approach with model-fitting to detect allelic effects specific to either family history subgroup. We also considered a polygenic model, in which we tested whether family history positive subjects carried more schizophrenia risk alleles than family history negative subjects, on average. Several individual SNPs attained suggestive but not genome-wide significant association with either family history subgroup. Comparison of genome-wide polygenic risk scores based on GWAS summary statistics indicated a significant enrichment for SNP effects among family history positive compared to family history negative cases (Nagelkerke's R(2 ) = 0.0021; P = 0.00331; P-value threshold <0.4). Estimates of variability in disease liability attributable to the aggregate effect of genome-wide SNPs were significantly greater for family history positive compared to family history negative cases (0.32 and 0.22, respectively; P = 0.031). We found suggestive evidence of allelic effects detectable in large GWAS of schizophrenia that might be specific to particular family history subgroups. However, consideration of a polygenic risk score indicated a significant enrichment among family history positive cases for common allelic effects. Familial illness might, therefore, represent a more heritable form of schizophrenia, as suggested by previous epidemiological studies.

Correlates of Risk for Disinhibited Behaviors in the Million Veteran Program Cohort.

Importance: Many psychiatric outcomes share a common etiologic pathway reflecting behavioral disinhibition, generally referred to as externalizing (EXT) disorders. Recent genome-wide association studies (GWASs) have demonstrated the overlap between EXT disorders and important aspects of veterans' health, such as suicide-related behaviors and substance use disorders (SUDs). Objective: To explore correlates of risk for EXT disorders within the Veterans Health Administration (VA) Million Veteran Program (MVP). Design, Setting, and Participants: A series of phenome-wide association studies (PheWASs) of polygenic risk scores (PGSs) for EXT disorders was conducted using electronic health records. First, ancestry-specific PheWASs of EXT PGSs were conducted in the African, European, and Hispanic or Latin American ancestries. Next, a conditional PheWAS, covarying for PGSs of comorbid psychiatric problems (depression, schizophrenia, and suicide attempt; European ancestries only), was performed. Lastly, to adjust for unmeasured confounders, a within-family analysis of significant associations from the main PheWAS was performed in full siblings (European ancestries only). This study included the electronic health record data from US veterans from VA health care centers enrolled in MVP. Analyses took place from February 2022 to August 2023 covering a period from October 1999 to January 2020. Exposures: PGSs for EXT, depression, schizophrenia, and suicide attempt. Main Outcomes and Measures: Phecodes for diagnoses derived from the International Statistical Classification of Diseases, Ninth and Tenth Revisions, Clinical Modification, codes from electronic health records. Results: Within the MVP (560 824 patients; mean [SD] age, 67.9 [14.3] years; 512 593 male [91.4%]), the EXT PGS was associated with 619 outcomes, of which 188 were independent of risk for comorbid problems or PGSs (from odds ratio [OR], 1.02; 95% CI, 1.01-1.03 for overweight/obesity to OR, 1.44; 95% CI, 1.42-1.47 for viral hepatitis C). Of the significant outcomes, 73 (11.9%) were significant in the African results and 26 (4.5%) were significant in the Hispanic or Latin American results. Within-family analyses uncovered robust associations between EXT PGS and consequences of SUDs, including liver disease, chronic airway obstruction, and viral hepatitis C. Conclusions and Relevance: Results of this cohort study suggest a shared polygenic basis of EXT disorders, independent of risk for other psychiatric problems. In addition, this study found associations between EXT PGS and diagnoses related to SUDs and their sequelae. Overall, this study highlighted the potential negative consequences of EXT disorders for health and functioning in the US veteran population.

Detecting and Adjusting for Hidden Biases due to Phenotype Misclassification in Genome-Wide Association Studies.

With the advent of healthcare-based genotyped biobanks, genome-wide association studies (GWAS) leverage larger sample sizes, incorporate patients with diverse ancestries and introduce noisier phenotypic definitions. Yet the extent and impact of phenotypic misclassification on large-scale datasets is not currently well understood due to a lack of statistical methods to estimate relevant parameters from empirical data. Here, we develop a statistical method and scalable software, PheMED, Phenotypic Measurement of Effective Dilution, to quantify phenotypic misclassification across GWAS using only summary statistics. We illustrate how the parameters estimated by PheMED relate to the negative and positive predictive value of the labeled phenotype, compared to ground truth, and how misclassification of the phenotype yields diluted effect-sizes of variant-phenotype associations. Furthermore, we apply our methodology to detect multiple instances of statistically significant dilution in real-world data. We demonstrate how effective dilution biases downstream GWAS replication and heritability analyses despite utilizing current best practices, and provide a dilution-aware meta-analysis approach that outperforms existing methods. Consequently, we anticipate that PheMED will be a valuable tool for researchers to address phenotypic data quality issues both within and across cohorts.