Genome-wide association study of esophageal squamous cell cancer identifies shared and distinct risk variants in African and Chinese populations.

Esophageal squamous cell carcinoma (ESCC) has a high disease burden in sub-Saharan Africa and has a very poor prognosis. Genome-wide association studies (GWASs) of ESCC in predominantly East Asian populations indicate a substantial genetic contribution to its etiology, but no genome-wide studies have been done in populations of African ancestry. Here, we report a GWAS in 1,686 African individuals with ESCC and 3,217 population-matched control individuals to investigate its genetic etiology. We identified a genome-wide-significant risk locus on chromosome 9 upstream of FAM120A (rs12379660, p = 4.58 × 10-8, odds ratio = 1.28, 95% confidence interval = 1.22-1.34), as well as a potential African-specific risk locus on chromosome 2 (rs142741123, p = 5.49 × 10-8) within MYO1B. FAM120A is a component of oxidative stress-induced survival signals, and the associated variants at the FAM120A locus co-localized with highly significant cis-eQTLs in FAM120AOS in both esophageal mucosa and esophageal muscularis tissue. A trans-ethnic meta-analysis was then performed with the African ESCC study and a Chinese ESCC study in a combined total of 3,699 ESCC-affected individuals and 5,918 control individuals, which identified three genome-wide-significant loci on chromosome 9 at FAM120A (rs12379660, pmeta = 9.36 × 10-10), chromosome 10 at PLCE1 (rs7099485, pmeta = 1.48 × 10-8), and chromosome 22 at CHEK2 (rs1033667, pmeta = 1.47 × 10-9). This indicates the existence of both shared and distinct genetic risk loci for ESCC in African and Asian populations. Our GWAS of ESCC conducted in a population of African ancestry indicates a substantial genetic contribution to ESCC risk in Africa.

Stratified genome-wide association analysis of type 2 diabetes reveals subgroups with genetic and environmental heterogeneity.

Type 2 diabetes (T2D) is a heterogeneous illness caused by genetic and environmental factors. Previous genome-wide association studies (GWAS) have identified many genetic variants associated with T2D and found evidence of differing genetic profiles by age-at-onset. This study seeks to explore further the genetic and environmental drivers of T2D by analyzing subgroups on the basis of age-at-onset of diabetes and body mass index (BMI). In the UK Biobank, 36 494 T2D cases were stratified into three subgroups, and GWAS was performed for all T2D cases and for each subgroup relative to 421 021 controls. Altogether, 18 single nucleotide polymorphisms were significantly associated with T2D genome-wide in one or more subgroups and also showed evidence of heterogeneity between the subgroups (Cochrane's Q P < 0.01), with two SNPs remaining significant after multiple testing (in CDKN2B and CYTIP). Combined risk scores, on the basis of genetic profile, BMI and age, resulted in excellent diabetes prediction [area under the ROC curve (AUC) = 0.92]. A modest improvement in prediction (AUC = 0.93) was seen when the contribution of genetic and environmental factors was evaluated separately for each subgroup. Increasing sample sizes of genetic studies enables us to stratify disease cases into subgroups, which have sufficient power to highlight areas of genetic heterogeneity. Despite some evidence that optimizing combined risk scores by subgroup improves prediction, larger sample sizes are likely needed for prediction when using a stratification approach.

Methylome-wide association study of early life stressors and adult mental health.

The environment and events that we are exposed to in utero, during birth and in early childhood influence our future physical and mental health. The underlying mechanisms that lead to these outcomes are unclear, but long-term changes in epigenetic marks, such as DNA methylation, could act as a mediating factor or biomarker. DNA methylation data were assayed at 713 522 CpG sites from 9537 participants of the Generation Scotland: Scottish Family Health Study, a family-based cohort with extensive genetic, medical, family history and lifestyle information. Methylome-wide association studies of eight early life environment phenotypes and two adult mental health phenotypes (major depressive disorder and brief resilience scale) were conducted using DNA methylation data collected from adult whole blood samples. Two genes involved with different developmental pathways (PRICKLE2, Prickle Planar Cell Polarity Protein 2 and ABI1, Abl-Interactor-1) were annotated to CpG sites associated with preterm birth (P < 1.27 × 10-9). A further two genes important to the development of sensory pathways (SOBP, Sine Oculis Binding Protein Homolog and RPGRIP1, Retinitis Pigmentosa GTPase Regulator Interacting Protein) were annotated to sites associated with low birth weight (P < 4.35 × 10-8). The examination of methylation profile scores and genes and gene-sets annotated from associated CpGs sites found no evidence of overlap between the early life environment and mental health conditions. Birth date was associated with a significant difference in estimated lymphocyte and neutrophil counts. Previous studies have shown that early life environments influence the risk of developing mental health disorders later in life; however, this study found no evidence that this is mediated by stable changes to the methylome detectable in peripheral blood.

The impact of major depressive disorder on glycaemic control in type 2 diabetes: a longitudinal cohort study using UK Biobank primary care records.

BACKGROUND: This study evaluates longitudinal associations between glycaemic control, measured by mean and within-patient variability of glycated haemaglobin (HbA1c) levels, and major depressive disorder (MDD) in individuals with type 2 diabetes (T2D), focusing on the timings of these diagnoses. METHODS: In UK Biobank, T2D was defined using self-report and linked health outcome data, then validated using polygenic scores. Repeated HbA1c measurements (mmol/mol) over the 10 years following T2D diagnosis were outcomes in mixed effects models, with disease duration included using restricted cubic splines. Four MDD exposures were considered: MDD diagnosis prior to T2D diagnosis (pre-T2D MDD), time between pre-T2D MDD diagnosis and T2D, new MDD diagnosis during follow-up (post-T2D MDD) and time since post-T2D MDD diagnosis. Models with and without covariate adjustment were considered. RESULTS: T2D diagnostic criteria were robustly associated with T2D polygenic scores. In 11,837 T2D cases (6.9 years median follow-up), pre-T2D MDD was associated with a 0.92 increase in HbA1c (95% CI: [0.00, 1.84]), but earlier pre-T2D MDD diagnosis correlated with lower HbA1c. These pre-T2D MDD effects became non-significant after covariate adjustment. Post-T2D MDD individuals demonstrated increasing HbA1c with years since MDD diagnosis ( ß = 0.51 , 95% CI: [0.17, 0.86]). Retrospectively, across study follow-up, within-patient variability in HbA1c was 1.16 (95% CI: 1.13-1.19) times higher in post-T2D MDD individuals. CONCLUSIONS: The timing of MDD diagnosis is important for understanding glycaemic control in T2D. Poorer control was observed in MDD diagnosed post-T2D, highlighting the importance of depression screening in T2D, and closer monitoring for individuals who develop MDD after T2D.

Estimating the impact of transmitted and non-transmitted psychiatric and neurodevelopmental polygenic scores on youth emotional problems.

Anxiety and depression (emotional disorders) are familial and heritable, especially when onset is early. However, other cross-generational studies suggest transmission of youth emotional problems is explained by mainly environmental risks. We set out to test the contribution of parental non-transmitted genetic liability, as indexed by psychiatric/neurodevelopmental common polygenic liability, to youth emotional problems using a UK population-based cohort: the Millennium Cohort Study. European (N = 6328) and South Asian (N = 814) ancestries were included, as well as a subset with genomic data from both parents (European: N = 2809; South Asian: N = 254). We examined the association of transmitted (PGST) and non-transmitted polygenic scores (PGSNT) for anxiety, depression, bipolar disorder and neurodevelopmental disorders (attention-deficit/hyperactivity disorder [ADHD], autism spectrum disorder [ASD], schizophrenia) with youth emotional disorder and symptom scores, measured using the parent- and self-reported Strengths and Difficulties Questionnaire emotional subscale at 6 timepoints between ages 3-17 years. In the European sample, PGST for anxiety and depression, but not bipolar disorder, were associated with emotional disorder and symptom scores across all ages, except age 3, with strongest association in adolescence. ADHD and ASD PGST also showed association across ages 11-17 years. In the South Asian sample, evidence for associations between all PGST and outcome measures were weaker. There was weak evidence of association between PGSNT for anxiety and depression and age 17 symptom scores in the South Asian sample, but not in the European sample for any outcome. Overall, PGST for depression, anxiety, ADHD and ASD contributed to youth emotional problems, with stronger associations in adolescence. There was limited support for non-transmitted genetic effects: these findings do not support the hypothesis that parental polygenic psychiatric/neurodevelopmental liability confer risk to offspring emotional problems through non-transmitted rearing/nurture effects.

Stratifying early-onset emotional disorders: using genetics to assess persistence in young people of European and South Asian ancestry.

BACKGROUND: Depression and anxiety are the most common mental health problems in young people. Currently, clinicians are advised to wait before initiating treatment for young people with these disorders as many spontaneously remit. However, others develop recurrent disorder but this subgroup cannot be identified at the outset. We examined whether psychiatric polygenic scores (PGS) could help inform stratification efforts to predict those at higher risk of recurrence. METHODS: Probable emotional disorder was examined in two UK population cohorts using the emotional symptoms subscale of the Strengths and Difficulties Questionnaire (SDQ). Those with emotional disorder at two or more time points between ages 5 and 25 years were classed as 'recurrent emotional disorder' (n = 1,643) and those with emotional disorder at one time point as having 'single episode emotional disorder' (n = 1,435, controls n = 8,715). We first examined the relationship between psychiatric PGS and emotional disorders in childhood and adolescence. Second, we tested whether psychiatric PGS added to predictor variables of known association with emotional disorder (neurodevelopmental comorbidity, special educational needs, family history of depression and socioeconomic status) when discriminating between single-episode and recurrent emotional disorder. Analyses were conducted separately in individuals of European and South Asian ancestry. RESULTS: Probable emotional disorder was associated with higher PGS for major depressive disorder (MDD), anxiety, broad depression, ADHD and autism spectrum disorder (ASD) in those of European ancestry. Higher MDD and broad depression PGS were associated with emotional disorder in people of South Asian ancestry. Recurrent, compared to single-episode, emotional disorder was associated with ASD and parental psychiatric history. PGS were not associated with episode recurrence, and PGS did not improve discrimination of recurrence when combined with clinical predictors. CONCLUSIONS: Our findings do not support the use of PGS as a tool to assess the likelihood of recurrence in young people experiencing their first episode of emotional disorder.

Evaluation of polygenic prediction methodology within a reference-standardized framework.

The predictive utility of polygenic scores is increasing, and many polygenic scoring methods are available, but it is unclear which method performs best. This study evaluates the predictive utility of polygenic scoring methods within a reference-standardized framework, which uses a common set of variants and reference-based estimates of linkage disequilibrium and allele frequencies to construct scores. Eight polygenic score methods were tested: p-value thresholding and clumping (pT+clump), SBLUP, lassosum, LDpred1, LDpred2, PRScs, DBSLMM and SBayesR, evaluating their performance to predict outcomes in UK Biobank and the Twins Early Development Study (TEDS). Strategies to identify optimal p-value thresholds and shrinkage parameters were compared, including 10-fold cross validation, pseudovalidation and infinitesimal models (with no validation sample), and multi-polygenic score elastic net models. LDpred2, lassosum and PRScs performed strongly using 10-fold cross-validation to identify the most predictive p-value threshold or shrinkage parameter, giving a relative improvement of 16-18% over pT+clump in the correlation between observed and predicted outcome values. Using pseudovalidation, the best methods were PRScs, DBSLMM and SBayesR. PRScs pseudovalidation was only 3% worse than the best polygenic score identified by 10-fold cross validation. Elastic net models containing polygenic scores based on a range of parameters consistently improved prediction over any single polygenic score. Within a reference-standardized framework, the best polygenic prediction was achieved using LDpred2, lassosum and PRScs, modeling multiple polygenic scores derived using multiple parameters. This study will help researchers performing polygenic score studies to select the most powerful and predictive analysis methods.

Clarifying the causes of consistent and inconsistent findings in genetics.

As research in genetics has advanced, some findings have been unexpected or shown to be inconsistent between studies or datasets. The reasons these inconsistencies arise are complex. Results from genetic studies can be affected by various factors including statistical power, linkage disequilibrium, quality control, confounding and selection bias, as well as real differences from interactions and effect modifiers, which may be informative about the mechanisms of traits and disease. Statistical artefacts can manifest as differences between results but they can also conceal underlying differences, which implies that their critical examination is important for understanding the underpinnings of traits. In this review, we examine these factors and outline how they can be identified and conceptualised with structural causal models. We explain the consequences they have on genetic estimates, such as genetic associations, polygenic scores, family- and genome-wide heritability, and describe methods to address them to aid in the estimation of true effects of genetic variation. Clarifying these factors can help researchers anticipate when results are likely to diverge and aid researchers' understanding of causal relationships between genes and complex traits.

Exploring polygenic-environment and residual-environment interactions for depressive symptoms within the UK Biobank.

Substantial advances have been made in identifying genetic contributions to depression, but little is known about how the effect of genes can be modulated by the environment, creating a gene-environment interaction. Using multivariate reaction norm models (MRNMs) within the UK Biobank (N = 61294-91644), we investigate whether the polygenic and residual variance components of depressive symptoms are modulated by 17 a priori selected covariate traits-12 environmental variables and 5 biomarkers. MRNMs, a mixed-effects modelling approach, provide unbiased polygenic-covariate interaction estimates for a quantitative trait by controlling for outcome-covariate correlations and residual-covariate interactions. A continuous depressive symptom variable was the outcome in 17 MRNMs-one for each covariate trait. Each MRNM had a fixed-effects model (fixed effects included the covariate trait, demographic variables, and principal components) and a random effects model (where polygenic-covariate and residual-covariate interactions are modelled). Of the 17 selected covariates, 11 significantly modulate deviations in depressive symptoms through the modelled interactions, but no single interaction explains a large proportion of phenotypic variation. Results are dominated by residual-covariate interactions, suggesting that covariate traits (including neuroticism, childhood trauma, and BMI) typically interact with unmodelled variables, rather than a genome-wide polygenic component, to influence depressive symptoms. Only average sleep duration has a polygenic-covariate interaction explaining a demonstrably nonzero proportion of the variability in depressive symptoms. This effect is small, accounting for only 1.22% (95% confidence interval: [0.54, 1.89]) of variation. The presence of an interaction highlights a specific focus for intervention, but the negative results here indicate a limited contribution from polygenic-environment interactions.

Kaposi sarcoma-associated herpesvirus, HIV-1 and Kaposi sarcoma risk in black South Africans diagnosed with cancer during antiretroviral treatment rollout.

Kaposi sarcoma-associated herpesvirus (KSHV) causes Kaposi sarcoma (KS). The risk of KS is amplified in HIV-immunosuppressed individuals and antiretroviral therapy (ART) reduces KS incidence. Reliable data on the relationship between these factors are lacking in Africa. We used questionnaires and serum from 7886 black South Africans (18-74 years) with incident cancer, recruited between 1995 and 2016. ART rollout started in 2004. We measured associations between KS, HIV-1 and KSHV before and after ART rollout. We measured seropositivity to HIV-1, KSHV latency-associated nuclear antigen (LANA) and glycoprotein (K8.1) and calculated case-control-adjusted odds ratios (ORadj ) and 95% confidence intervals (CI) in relation to KS and KSHV infection, before (1995-2004), early (2005-2009) and late (2010-2016) ART rollout periods. KSHV seropositivity among 1237 KS cases was 98%. Among 6649 controls, KSHV seropositivity was higher in males (ORadj  = 1.4 [95%CI 1.23-1.52]), in persons with HIV, (ORadj  = 4.2 [95%CI 3.74-4.73]) and lower in high school leavers (ORadj  = 0.7 [95%CI 0.59-0.83]). KSHV seropositivity declined over the three ART rollout periods (37%, 28% and 28%, Ptrend < .001) coinciding with increases in high school leavers over the same periods (46%, 58% and 67%, Ptrend < .001). HIV-1 seroprevalence increased from 10% in the pre-ART period to 22% in the late ART period (Ptrend < .001). Compared to HIV-1 and KSHV seronegatives, KSHV seropositives yielded an OR for KS of 26 (95%CI 11-62) in HIV-1 seronegative participants and an OR of 2501 (95%CI 1083-5776) in HIV-1 seropositive participants. HIV-1 increases the risk of KS in those infected with KSHV by 100-fold. Declines in KSHV seroprevalence coincide with ART rollout and with improvements in educational standards and general hygiene.

Imputed gene expression risk scores: a functionally informed component of polygenic risk.

Integration of functional genomic annotations when estimating polygenic risk scores (PRS) can provide insight into aetiology and improve risk prediction. This study explores the predictive utility of gene expression risk scores (GeRS), calculated using imputed gene expression and transcriptome-wide association study (TWAS) results. The predictive utility of GeRS was evaluated using 12 neuropsychiatric and anthropometric outcomes measured in two target samples: UK Biobank and the Twins Early Development Study. GeRS were calculated based on imputed gene expression levels and TWAS results, using 53 gene expression-genotype panels, termed single nucleotide polymorphism (SNP)-weight sets, capturing expression across a range of tissues. We compare the predictive utility of elastic net models containing GeRS within and across SNP-weight sets, and models containing both GeRS and PRS. We estimate the proportion of SNP-based heritability attributable to cis-regulated gene expression. GeRS significantly predicted a range of outcomes, with elastic net models combining GeRS across SNP-weight sets improving prediction. GeRS were less predictive than PRS, but models combining GeRS and PRS improved prediction for several outcomes, with relative improvements ranging from 0.3% for height (P = 0.023) to 4% for rheumatoid arthritis (P = 5.9 × 10-8). The proportion of SNP-based heritability attributable to cis-regulated expression was modest for most outcomes, even when restricting GeRS to colocalized genes. GeRS represent a component of PRS and could be useful for functional stratification of genetic risk. Only in specific circumstances can GeRS substantially improve prediction over PRS alone. Future research considering functional genomic annotations when estimating genetic risk is warranted.

Depressive symptoms and neuroticism-related traits are the main factors associated with wellbeing independent of the history of lifetime depression in the UK Biobank.

BACKGROUND: Wellbeing has a fundamental role in determining life expectancy and major depressive disorder (MDD) is one of the main modulating factors of wellbeing. This study evaluated the modulators of wellbeing in individuals with lifetime recurrent MDD (RMDD), single-episode MDD (SMDD) and no MDD in the UK Biobank. METHODS: Scores of happiness, meaningful life and satisfaction about functioning were condensed in a functioning-wellbeing score (FWS). We evaluated depression and anxiety characteristics, neuroticism-related traits, physical diseases, lifestyle and polygenic risk scores (PRSs) of psychiatric disorders. Other than individual predictors, we estimated the cumulative contribution to FWS of each group of predictors. We tested the indirect role of neuroticism on FWS through the modulation of depression manifestations using a mediation analysis. RESULTS: We identified 47 966, 21 117 and 207 423 individuals with lifetime RMDD, SMDD and no MDD, respectively. Depression symptoms and personality showed the largest impact on FWS (variance explained ~20%), particularly self-harm, worthlessness feelings during the worst depression, chronic depression, loneliness and neuroticism. Personality played a stronger role in SMDD. Anxiety characteristics showed a higher effect in SMDD and no MDD groups. Neuroticism played indirect effects through specific depressive symptoms that modulated FWS. Physical diseases and lifestyle explained only 4-5% of FWS variance. The PRS of MDD showed the largest effect on FWS compared to other PRSs. CONCLUSIONS: This was the first study to comprehensively evaluate the predictors of wellbeing in relation to the history of MDD. The identified variables are important to identify individuals at risk and promote wellbeing.

Anxiety disorders and age-related changes in physiology.

BACKGROUND: Anxiety disorders are leading contributors to the global disease burden, highly prevalent across the lifespan and associated with substantially increased morbidity and early mortality. AIMS: The aim of this study was to examine age-related changes across a wide range of physiological measures in middle-aged and older adults with a lifetime history of anxiety disorders compared with healthy controls. METHOD: The UK Biobank study recruited >500 000 adults, aged 37-73, between 2006 and 2010. We used generalised additive models to estimate non-linear associations between age and hand-grip strength, cardiovascular function, body composition, lung function and heel bone mineral density in a case group and in a control group. RESULTS: The main data-set included 332 078 adults (mean age 56.37 years; 52.65% females). In both sexes, individuals with anxiety disorders had a lower hand-grip strength and lower blood pressure, whereas their pulse rate and body composition measures were higher than in the healthy control group. Case-control group differences were larger when considering individuals with chronic and/or severe anxiety disorders, and differences in body composition were modulated by depression comorbidity status. Differences in age-related physiological changes between females in the anxiety disorder case group and healthy controls were most evident for blood pressure, pulse rate and body composition, whereas this was the case in males for hand-grip strength, blood pressure and body composition. Most differences in physiological measures between the case and control groups decreased with increasing age. CONCLUSIONS: Findings in individuals with a lifetime history of anxiety disorders differed from a healthy control group across multiple physiological measures, with some evidence of case-control group differences by age. The differences observed varied by chronicity/severity and depression comorbidity.

Latent subtypes of manic and/or irritable episode symptoms in two population-based cohorts.

BACKGROUND: Mood disorders are characterised by pronounced symptom heterogeneity, which presents a substantial challenge both to clinical practice and research. Identification of subgroups of individuals with homogeneous symptom profiles that cut across current diagnostic categories could provide insights in to the transdiagnostic relevance of individual symptoms, which current categorical diagnostic systems cannot impart. AIMS: To identify groups of people with homogeneous clinical characteristics, using symptoms of manic and/or irritable mood, and explore differences between groups in diagnoses, functional outcomes and genetic liability. METHOD: We used latent class analysis on eight binary self-reported symptoms of manic and irritable mood in the UK Biobank and PROTECT studies, to investigate how individuals formed latent subgroups. We tested associations between the latent classes and diagnoses of psychiatric disorders, sociodemographic characteristics and polygenic risk scores. RESULTS: Five latent classes were derived in UK Biobank (N = 42 183) and were replicated in the independent PROTECT cohort (N = 4445), including 'minimally affected', 'inactive restless', active restless', 'focused creative' and 'extensively affected' individuals. These classes differed in disorder risk, polygenic risk score and functional outcomes. One class that experienced disruptive episodes of mostly irritable mood largely comprised cases of depression/anxiety, and a class of individuals with increased confidence/creativity reported comparatively lower disruptiveness and functional impairment. CONCLUSIONS: Findings suggest that data-driven investigations of psychopathological symptoms that include sub-diagnostic threshold conditions can complement research of clinical diagnoses. Improved classification systems of psychopathology could investigate a weighted approach to symptoms, toward a more dimensional classification of mood disorders.

Using major depression polygenic risk scores to explore the depressive symptom continuum.

BACKGROUND: Major depression (MD) is often characterised as a categorical disorder; however, observational studies comparing sub-threshold and clinical depression suggest MD is continuous. Many of these studies do not explore the full continuum and are yet to consider genetics as a risk factor. This study sought to understand if polygenic risk for MD could provide insight into the continuous nature of depression. METHODS: Factor analysis on symptom-level data from the UK Biobank (N = 148 957) was used to derive continuous depression phenotypes which were tested for association with polygenic risk scores (PRS) for a categorical definition of MD (N = 119 692). RESULTS: Confirmatory factor analysis showed a five-factor hierarchical model, incorporating 15 of the original 18 items taken from the PHQ-9, GAD-7 and subjective well-being questionnaires, produced good fit to the observed covariance matrix (CFI = 0.992, TLI = 0.99, RMSEA = 0.038, SRMR = 0.031). MD PRS associated with each factor score (standardised ß range: 0.057-0.064) and the association remained when the sample was stratified into case- and control-only subsets. The case-only subset had an increased association compared to controls for all factors, shown via a significant interaction between lifetime MD diagnosis and MD PRS (p value range: 2.23 × 10-3-3.94 × 10-7). CONCLUSIONS: An association between MD PRS and a continuous phenotype of depressive symptoms in case- and control-only subsets provides support against a purely categorical phenotype; indicating further insights into MD can be obtained when this within-group variation is considered. The stronger association within cases suggests this variation may be of particular importance.

Depression with atypical neurovegetative symptoms shares genetic predisposition with immuno-metabolic traits and alcohol consumption.

BACKGROUND: Depression is a highly prevalent and heterogeneous disorder. This study aims to determine whether depression with atypical features shows different heritability and different degree of overlap with polygenic risk for psychiatric and immuno-metabolic traits than other depression subgroups. METHODS: Data included 30 069 European ancestry individuals from the UK Biobank who met criteria for lifetime major depression. Participants reporting both weight gain and hypersomnia were classified as ↑WS depression (N = 1854) and the others as non-↑WS depression (N = 28 215). Cases with non-↑WS depression were further classified as ↓WS depression (i.e. weight loss and insomnia; N = 10 142). Polygenic risk scores (PRS) for 22 traits were generated using genome-wide summary statistics (Bonferroni corrected p = 2.1 × 10-4). Single-nucleotide polymorphism (SNP)-based heritability of depression subgroups was estimated. RESULTS: ↑WS depression had a higher polygenic risk for BMI [OR = 1.20 (1.15-1.26), p = 2.37 × 10-14] and C-reactive protein [OR = 1.11 (1.06-1.17), p = 8.86 × 10-06] v. non-↑WS depression and ↓WS depression. Leptin PRS was close to the significance threshold (p = 2.99 × 10-04), but the effect disappeared when considering GWAS summary statistics of leptin adjusted for BMI. PRS for daily alcohol use was inversely associated with ↑WS depression [OR = 0.88 (0.83-0.93), p = 1.04 × 10-05] v. non-↑WS depression. SNP-based heritability was not significantly different between ↑WS depression and ↓WS depression (14.3% and 12.2%, respectively). CONCLUSIONS: ↑WS depression shows evidence of distinct genetic predisposition to immune-metabolic traits and alcohol consumption. These genetic signals suggest that biological targets including immune-cardio-metabolic pathways may be relevant to therapies in individuals with ↑WS depression.

Exploring the genetic heterogeneity in major depression across diagnostic criteria.

Major depressive disorder (MDD) is defined differently across genetic research studies and this may be a key source of heterogeneity. While previous literature highlights differences between minimal and strict phenotypes, the components contributing to this heterogeneity have not been identified. Using the cardinal symptoms (depressed mood/anhedonia) as a baseline, we build MDD phenotypes using five components-(1) five or more symptoms, (2) episode duration, (3) functional impairment, (4) episode persistence, and (5) episode recurrence-to determine the contributors to such heterogeneity. Thirty-two depression phenotypes which systematically incorporate different combinations of MDD components were created using the mental health questionnaire data within the UK Biobank. SNP-based heritabilities and genetic correlations with three previously defined major depression phenotypes were calculated (Psychiatric Genomics Consortium (PGC) defined depression, 23andMe self-reported depression and broad depression) and differences between estimates analysed. All phenotypes were heritable (h2SNP range: 0.102-0.162) and showed substantial genetic correlations with other major depression phenotypes (Rg range: 0.651-0.895 (PGC); 0.652-0.837 (23andMe); 0.699-0.900 (broad depression)). The strongest effect on SNP-based heritability was from the requirement for five or more symptoms (1.4% average increase) and for a long episode duration (2.7% average decrease). No significant differences were noted between genetic correlations. While there is some variation, the two cardinal symptoms largely reflect the genetic aetiology of phenotypes incorporating more MDD components. These components may index severity, however, their impact on heterogeneity in genetic results is likely to be limited.

Genetic and clinical characteristics of treatment-resistant depression using primary care records in two UK cohorts.

Treatment-resistant depression (TRD) is a major contributor to the disability caused by major depressive disorder (MDD). Primary care electronic health records provide an easily accessible approach to investigate TRD clinical and genetic characteristics. MDD defined from primary care records in UK Biobank (UKB) and EXCEED studies was compared with other measures of depression and tested for association with MDD polygenic risk score (PRS). Using prescribing records, TRD was defined from at least two switches between antidepressant drugs, each prescribed for at least 6 weeks. Clinical-demographic characteristics, SNP-based heritability (h2SNP) and genetic overlap with psychiatric and non-psychiatric traits were compared in TRD and non-TRD MDD cases. In 230,096 and 8926 UKB and EXCEED participants with primary care data, respectively, the prevalence of MDD was 8.7% and 14.2%, of which 13.2% and 13.5% was TRD, respectively. In both cohorts, MDD defined from primary care records was strongly associated with MDD PRS, and in UKB it showed overlap of 71-88% with other MDD definitions. In UKB, TRD vs healthy controls and non-TRD vs healthy controls h2SNP was comparable (0.25 [SE = 0.04] and 0.19 [SE = 0.02], respectively). TRD vs non-TRD was positively associated with the PRS of attention deficit hyperactivity disorder, with lower socio-economic status, obesity, higher neuroticism and other unfavourable clinical characteristics. This study demonstrated that MDD and TRD can be reliably defined using primary care records and provides the first large scale population assessment of the genetic, clinical and demographic characteristics of TRD.

Clinical, demographic, and genetic risk factors of treatment-attributed suicidality in >10,000 Australian adults taking antidepressants.

Emergence of suicidal symptoms has been reported as a potential antidepressant adverse drug reaction. Identifying risk factors associated could increase our understanding of this phenomenon and stratify individuals at higher risk. Logistic regressions were used to identify risk factors of self-reported treatment-attributed suicidal ideation (TASI). We then employed classifiers to test the predictive ability of the variables identified. A TASI GWAS, as well as SNP-based heritability estimation, were performed. GWAS replication was sought from an independent study. Significant associations were found for age and comorbid conditions, including bipolar and personality disorders. Participants reporting TASI from one antidepressant were more likely to report TASI from other antidepressants. No genetic loci associated with TAS I (p < 5e-8) were identified. Of 32 independent variants with suggestive association (p < 1e-5), 27 lead SNPs were available in a replication dataset from the GENDEP study. Only one variant showed a consistent effect and nominal association in the independent replication sample. Classifiers were able to stratify non-TASI from TASI participants (AUC = 0.77) and those reporting treatment-attributed suicide attempts (AUC = 0.85). The pattern of TASI co-occurrence across participants suggest nonspecific factors underlying its etiology. These findings provide insights into the underpinnings of TASI and serve as a proof-of-concept of the use of classifiers for risk stratification.