The impact of phenotypic and genetic heterogeneity on results of genome wide association studies of complex diseases

Optimization of self- or parent-reported psychiatric phenotypes in longitudinal studies

BACKGROUND

The Adolescent Brain Cognitive Development (ABCD) study is a longitudinal study of US adolescents with a wide breadth of psychiatric, neuroimaging and genetic data that can be leveraged to better understand psychiatric diseases. The reliability and validity of the psychiatric data collected have not yet been examined. This study aims to explore and optimize the reliability/validity of psychiatric diagnostic constructs in the ABCD study.

METHODS

Parent-and-child-reported psychiatric data for 11,876 children (aged 9.5 ± 0.5 at first assessment) were examined over 4 years to derive specific constructs for psychiatric diagnoses using longitudinal information. Rates of psychiatric disorders were calculated and compared to those reported in the epidemiological literature.

RESULTS

The rates of self-reported psychiatric disorders at any single time point (broad diagnostic construct) were higher than indicated by epidemiological studies. Narrow diagnostic constructs, which required the endorsement of psychiatric disorders at a majority of longitudinal assessments, demonstrated a better rate approximation of literature-reported prevalences for most disorders (e.g. the prevalence of broad obsessive-compulsive disorder (OCD) was 13.3% compared to narrow OCD at 2.6% and a literature-reported prevalence of 2.3%). Analysis of comorbidity, using OCD as a representative example, also showed a better approximation of literature-reported comorbidity rates using the narrow construct, with some exceptions.

CONCLUSIONS

Self- or parent-report-based assessments tend to overestimate prevalences of psychiatric disorders in the ABCD Study, particularly when longitudinal data are summed to create lifetime prevalences. Such assessments should be accompanied by more in-depth assessments or clinician-administered structured interviews if using data where accurate disorder classifications are paramount.

Genome-wide discovery for biomarkers using quantile regression at biobank scale

Genome-wide association studies (GWAS) for biomarkers important for clinical phenotypes can lead to clinically relevant discoveries. Conventional GWAS for quantitative traits are based on simplified regression models modeling the conditional mean of a phenotype as a linear function of genotype. We draw attention here to an alternative, lesser known approach, namely quantile regression that naturally extends linear regression to the analysis of the entire conditional distribution of a phenotype of interest. Quantile regression can be applied efficiently at biobank scale, while having some unique advantages such as (1) identifying variants with heterogeneous effects across quantiles of the phenotype distribution; (2) accommodating a wide range of phenotype distributions including non-normal distributions, with invariance of results to trait transformations; and (3) providing more detailed information about genotype-phenotype associations even for those associations identified by conventional GWAS. We show in simulations that quantile regression is powerful across both homogeneous and various heterogeneous models. Applications to 39 quantitative traits in the UK Biobank demonstrate that quantile regression can be a helpful complement to linear regression in GWAS and can identify variants with larger effects on high-risk subgroups of individuals but with lower or no contribution overall.

Unsupervised clustering identified clinically relevant metabolic syndrome endotypes in UK and Taiwan Biobanks

Metabolic syndrome (MetS) is a collection of cardiovascular risk factors; however, the high prevalence and heterogeneity impede effective clinical management. We conducted unsupervised clustering on individuals from UK Biobank to reveal endotypes. Five MetS subgroups were identified: Cluster 1 (C1): non-descriptive, Cluster 2 (C2): hypertensive, Cluster 3 (C3): obese, Cluster 4 (C4): lipodystrophy-like, and Cluster 5 (C5): hyperglycemic. For all of the endotypes, we identified the corresponding cardiometabolic traits and their associations with clinical outcomes. Genome-wide association studies (GWASs) were conducted to identify associated genotypic traits. We then determined endotype-specific genotypic traits and constructed polygenic risk score (PRS) models specific to each endotype. GWAS of each MetS clusters revealed different genotypic traits. C1 GWAS revealed novel findings of TRIM63, MYBPC3, MYLPF, and RAPSN. Intriguingly, C1, C3, and C4 were associated with genes highly expressed in brain tissues. MetS clusters with comparable phenotypic and genotypic traits were identified in Taiwan Biobank.

Using Alternative Definitions of Controls to Increase Statistical Power in GWAS

Genome-wide association studies (GWAS) are often underpowered due to small effect sizes of common single nucleotide polymorphisms (SNPs) on phenotypes and extreme multiple testing thresholds. The most common approach for increasing statistical power is to increase sample size. We propose an alternative strategy of redefining case-control outcomes into ordinal case-subthreshold-asymptomatic variables. While maintaining the clinical case threshold, we subdivide controls into two groups: individuals who are symptomatic but do not meet the clinical criteria for diagnosis (subthreshold) and individuals who are effectively asymptomatic. We conducted a simulation study to examine the impact of effect size, minor allele frequency, population prevalence, and the prevalence of the subthreshold group on statistical power to detect genetic associations in three scenarios: a standard case-control, an ordinal, and a case-asymptomatic control analysis. Our results suggest the ordinal model consistently provides the greatest statistical power while the case-control model the least. Power in the case-asymptomatic control model reflects the case-control or ordinal model depending on the population prevalence and size of the subthreshold category. We then analyzed a major depression phenotype from the UK Biobank to corroborate our simulation results. Overall, the ordinal model improves statistical power in GWAS consistent with increasing the sample size by approximately 10%.

Genetics may affect the risk of undergoing surgery for rhizarthrosis

Osteoarthritis is a prevalent and severe disease. Involvement of the trapeziometacarpal joint is common and can lead to both pain and disability. Genetics are known to affect the risk of osteoarthritis, but it remains unclear how genetics affect disease trajectories. In this study, we investigated whether the genetic associations of trapeziometacarpal osteoarthritis (rhizarthrosis) vary with the need for surgical treatment. The study was conducted as a case-control genome-wide association study using individuals from the Copenhagen Hospital Biobank pain and degenerative musculoskeletal disease study and the Danish Blood Donor Study (N = 208,342). We identified patients diagnosed with rhizarthrosis and grouped them by treatment status, resulting in two case groups: surgical (N = 1083) and nonsurgical (N = 1888). The case groups were tested against osteoarthritis-free controls in two genome-wide association studies. We then compared variants suggestive of association (p < 10-6) in either of these analyses directly between the treatment groups (surgical vs. nonsurgical rhizarthrosis). We identified 10 variants suggestive of association with either surgical (seven variants) or nonsurgical (three variants) rhizarthrosis. None of the variants reached nominal significance in the opposite treatment group (p ≥ 0.14), and all 10 variants were significantly different between the treatment groups at a false discovery rate of 5%. These results suggest possible differences in the genetic associations of rhizarthrosis depending on surgical treatment. Clinical significance: Uncovering genetic differences between clinically distinct patient groups can reveal biological determinants of disease trajectories.

Longitudinal studies of bipolar patients and their families: translating findings to advance individualized risk prediction, treatment and research

BACKGROUND

Bipolar disorder is a broad diagnostic construct associated with significant phenotypic and genetic heterogeneity challenging progress in clinical practice and discovery research. Prospective studies of well-characterized patients and their family members have identified lithium responsive (LiR) and lithium non-responsive (LiNR) subtypes that hold promise for advancement.

METHOD

In this narrative review, relevant observations from published longitudinal studies of well-characterized bipolar patients and their families spanning six decades are highlighted. DSM diagnoses based on SADS-L interviews were decided in blind consensus reviews by expert clinicians. Genetic, neurobiological, and psychosocial factors were investigated in subsets of well-characterized probands and adult relatives. Systematic maintenance trials of lithium, antipsychotics, and lamotrigine were carried out. Clinical profiles that included detailed histories of the clinical course, symptom sets and disorders segregating in families were documented. Offspring of LiR and LiNR families were repeatedly assessed up to 20 years using KSADS-PL format interviews and DSM diagnoses and sub-threshold symptoms were decided by expert clinicians in blind consensus reviews using all available clinical and research data.

RESULTS

A characteristic clinical profile differentiated bipolar patients who responded to lithium stabilization from those who did not. The LiR subtype was characterized by a recurrent fully remitting course predominated by depressive episodes and a positive family history of episodic remitting mood disorders, and not schizophrenia. Response to lithium clustered in families and the characteristic clinical profile predicted lithium response, with the episodic remitting course being a strong correlate. There is accumulating evidence that genetic and neurobiological markers differ between LiR and LiNR subtypes. Further, offspring of bipolar parents subdivided by lithium response differed in developmental history, clinical antecedents and early course of mood disorders. Moreover, the nature of the emergent course bred true from parent to offspring, independent of the nature of emergent psychopathology.

CONCLUSIONS

Bipolar disorders are heterogeneous and response to long-term lithium is associated with a familial subtype with characteristic course, treatment response, family history and likely pathogenesis. Incorporating distinctive clinical profiles that index valid bipolar subtypes into routine practice and research will improve patient outcomes and advance the development and translation of novel treatment targets to improve prevention and early intervention.

Interpretation of 10 years of Alzheimer's disease genetic findings in the perspective of statistical heterogeneity

Common genetic variants and susceptibility loci associated with Alzheimer's disease (AD) have been discovered through large-scale genome-wide association studies (GWAS), GWAS by proxy (GWAX) and meta-analysis of GWAS and GWAX (GWAS+GWAX). However, due to the very low repeatability of AD susceptibility loci and the low heritability of AD, these AD genetic findings have been questioned. We summarize AD genetic findings from the past 10 years and provide a new interpretation of these findings in the context of statistical heterogeneity. We discovered that only 17% of AD risk loci demonstrated reproducibility with a genome-wide significance of P < 5.00E-08 across all AD GWAS and GWAS+GWAX datasets. We highlighted that the AD GWAS+GWAX with the largest sample size failed to identify the most significant signals, the maximum number of genome-wide significant genetic variants or maximum heritability. Additionally, we identified widespread statistical heterogeneity in AD GWAS+GWAX datasets, but not in AD GWAS datasets. We consider that statistical heterogeneity may have attenuated the statistical power in AD GWAS+GWAX and may contribute to explaining the low repeatability (17%) of genome-wide significant AD susceptibility loci and the decreased AD heritability (40-2%) as the sample size increased. Importantly, evidence supports the idea that a decrease in statistical heterogeneity facilitates the identification of genome-wide significant genetic loci and contributes to an increase in AD heritability. Collectively, current AD GWAX and GWAS+GWAX findings should be meticulously assessed and warrant additional investigation, and AD GWAS+GWAX should employ multiple meta-analysis methods, such as random-effects inverse variance-weighted meta-analysis, which is designed specifically for statistical heterogeneity.

Using Alternative Definitions of Controls to Increase Statistical Power in GWAS

Genome-wide association studies (GWAS) are underpowered due to small effect sizes of single nucleotide polymorphisms (SNPs) on phenotypes and extreme multiple testing thresholds. The most common approach for increasing statistical power is to increase sample size. We propose an alternative strategy of redefining case-control outcomes into ordinal case-subthreshold-asymptomatic variables. While maintaining the clinical case threshold, we subdivide controls into two groups: individuals who are symptomatic but do not meet the clinical criteria for diagnosis (subthreshold) and individuals who are effectively asymptomatic. We conducted a simulation study to examine the impact of effect size, minor allele frequency, population prevalence, and the prevalence of the subthreshold group on statistical power to detect genetic associations in three scenarios: a standard case-control, an ordinal, and a case-asymptomatic control analysis. Our results suggest the ordinal model consistently provides the most statistical power while the case-control model the least. Power in the case-asymptomatic control model reflects the case-control or ordinal model depending on the population prevalence and size of the subthreshold category. We then analyzed a major depression phenotype from the UK Biobank to corroborate our simulation results. Overall, the ordinal model improves statistical power in GWAS consistent with increasing the sample size by approximately 10%.

Omada: robust clustering of transcriptomes through multiple testing

BACKGROUND

Cohort studies increasingly collect biosamples for molecular profiling and are observing molecular heterogeneity. High-throughput RNA sequencing is providing large datasets capable of reflecting disease mechanisms. Clustering approaches have produced a number of tools to help dissect complex heterogeneous datasets, but selecting the appropriate method and parameters to perform exploratory clustering analysis of transcriptomic data requires deep understanding of machine learning and extensive computational experimentation. Tools that assist with such decisions without prior field knowledge are nonexistent. To address this, we have developed Omada, a suite of tools aiming to automate these processes and make robust unsupervised clustering of transcriptomic data more accessible through automated machine learning-based functions.

FINDINGS

The efficiency of each tool was tested with 7 datasets characterized by different expression signal strengths to capture a wide spectrum of RNA expression datasets. Our toolkit's decisions reflected the real number of stable partitions in datasets where the subgroups are discernible. Within datasets with less clear biological distinctions, our tools either formed stable subgroups with different expression profiles and robust clinical associations or revealed signs of problematic data such as biased measurements.

CONCLUSIONS

In conclusion, Omada successfully automates the robust unsupervised clustering of transcriptomic data, making advanced analysis accessible and reliable even for those without extensive machine learning expertise. Implementation of Omada is available at http://bioconductor.org/packages/omada/.

East Asian-specific and cross-ancestry genome-wide meta-analyses provide mechanistic insights into peptic ulcer disease

Peptic ulcer disease (PUD) refers to acid-induced injury of the digestive tract, occurring mainly in the stomach (gastric ulcer (GU)) or duodenum (duodenal ulcer (DU)). In the present study, we conducted a large-scale, cross-ancestry meta-analysis of PUD combining genome-wide association studies with Japanese and European studies (52,032 cases and 905,344 controls), and discovered 25 new loci highly concordant across ancestries. An examination of GU and DU genetic architecture demonstrated that GUs shared the same risk loci as DUs, although with smaller genetic effect sizes and higher polygenicity than DUs, indicating higher heterogeneity of GUs. Helicobacter pylori (HP)-stratified analysis found an HP-related host genetic locus. Integrative analyses using bulk and single-cell transcriptome profiles highlighted the genetic factors of PUD being enriched in the highly expressed genes in stomach tissues, especially in somatostatin-producing D cells. Our results provide genetic evidence that gastrointestinal cell differentiations and hormone regulations are critical in PUD etiology.

Microbiome-based risk prediction in incident heart failure: a community challenge

Heart failure (HF) is a major public health problem. Early identification of at-risk individuals could allow for interventions that reduce morbidity or mortality. The community-based FINRISK Microbiome DREAM challenge (synapse.org/finrisk) evaluated the use of machine learning approaches on shotgun metagenomics data obtained from fecal samples to predict incident HF risk over 15 years in a population cohort of 7231 Finnish adults (FINRISK 2002, n=559 incident HF cases). Challenge participants used synthetic data for model training and testing. Final models submitted by seven teams were evaluated in the real data. The two highest-scoring models were both based on Cox regression but used different feature selection approaches. We aggregated their predictions to create an ensemble model. Additionally, we refined the models after the DREAM challenge by eliminating phylum information. Models were also evaluated at intermediate timepoints and they predicted 10-year incident HF more accurately than models for 5- or 15-year incidence. We found that bacterial species, especially those linked to inflammation, are predictive of incident HF. This highlights the role of the gut microbiome as a potential driver of inflammation in HF pathophysiology. Our results provide insights into potential modeling strategies of microbiome data in prospective cohort studies. Overall, this study provides evidence that incorporating microbiome information into incident risk models can provide important biological insights into the pathogenesis of HF.

Bidirectional Mendelian Randomization Study of Personality Traits Reveals a Positive Feedback Loop Between Neuroticism and Back Pain

We conducted a bidirectional Mendelian randomization study to examine the causal effects of six personality traits (anxiety, neuroticism, extraversion, openness to experience, agreeableness, and conscientiousness) on back pain associated with health care use and the causal effect of back pain on the same risk factors. Genetic instruments for the personality traits and back pain were obtained from the largest published genome-wide association studies conducted in individuals of European ancestry. We used inverse weighted variance meta-analysis and Causal Analysis Using Summary Effect for primary analyses and sensitivity analyses to examine evidence for causal associations. We interpreted exposure-outcome associations as being consistent with a causal relationship if results of at least one primary analysis were statistically significant after accounting for multiple statistical testing (P-value < .0042), and the direction and magnitude of effect estimates were concordant between primary and sensitivity analyses. We found evidence for statistically significant bidirectional causal associations between neuroticism and back pain, with odds ratio 1.51 (95% confidence interval 1.37; 1.67) of back pain per neuroticism sum score standard deviation, P-value = 7.80e-16; and beta = .12, se = .04 of neuroticism sum score standard deviation per log odds of back pain, P-value = 2.48e-03. Other relationships did not meet our predefined criteria for causal association. PERSPECTIVE: The significant positive feedback loop between neuroticism and back pain highlights the importance of considering neuroticism in the management of patients with back pain.

SLCO5A1 and synaptic assembly genes contribute to impulsivity in juvenile myoclonic epilepsy

Elevated impulsivity is a key component of attention-deficit hyperactivity disorder (ADHD), bipolar disorder and juvenile myoclonic epilepsy (JME). We performed a genome-wide association, colocalization, polygenic risk score, and pathway analysis of impulsivity in JME (n = 381). Results were followed up with functional characterisation using a drosophila model. We identified genome-wide associated SNPs at 8q13.3 (P = 7.5 × 10-9) and 10p11.21 (P = 3.6 × 10-8). The 8q13.3 locus colocalizes with SLCO5A1 expression quantitative trait loci in cerebral cortex (P = 9.5 × 10-3). SLCO5A1 codes for an organic anion transporter and upregulates synapse assembly/organisation genes. Pathway analysis demonstrates 12.7-fold enrichment for presynaptic membrane assembly genes (P = 0.0005) and 14.3-fold enrichment for presynaptic organisation genes (P = 0.0005) including NLGN1 and PTPRD. RNAi knockdown of Oatp30B, the Drosophila polypeptide with the highest homology to SLCO5A1, causes over-reactive startling behaviour (P = 8.7 × 10-3) and increased seizure-like events (P = 6.8 × 10-7). Polygenic risk score for ADHD genetically correlates with impulsivity scores in JME (P = 1.60 × 10-3). SLCO5A1 loss-of-function represents an impulsivity and seizure mechanism. Synaptic assembly genes may inform the aetiology of impulsivity in health and disease.

The design of mapping populations: Impacts of geographic scale on genetic architecture and mapping efficacy for defense and immunity

Genome-wide association studies (GWAS) have yielded tremendous insight into the genetic architecture of trait variation. However, the collections of loci they uncover are far from exhaustive. As many of the complicating factors that confound or limit the efficacy of GWAS are exaggerated over broad geographic scales, a shift toward more analyses using mapping panels sampled from narrow geographic localities ("local" populations) could provide novel, complementary insights. Here, we present an overview of the major complicating factors, review mounting evidence from genomic analyses that these factors are pervasive, and synthesize theoretical and empirical evidence for the power of GWAS in local populations.

Prevalence and outcomes of rapid cycling bipolar disorder: Mixed method systematic meta-review

Rapid-cycling in bipolar disorder (RC-BD) is associated with greater illness morbidity and inferior treatment response but many aspects remain unclear, prompting this systematic review of its definitions, prevalence, and clinical characteristics. We searched multiple literature databases through April 2022 for systematic reviews or meta-analyses on RC-BD and extracted associated definitions, prevalence, risk-factors, and clinical outcomes. We assessed study quality (NIH Quality Assessment Tool) and levels of evidence (Oxford criteria). Of 146 identified reviews, 22 fulfilling selection criteria were included, yielding 30 studies involving 13,698 BD patients, of whom 3777 (27.6% [CI: 26.8-28.3]) were considered RC-BD, as defined in 14 reports by ≥4 recurrences/year within the past 12 months or in any year, without considering responsiveness to treatment. Random-effects meta-analytically pooled one-year prevalence was 22.3% [CI: 14.4-32.9] in 12 reports and lifetime prevalence was 35.5% [27.6-44.3] in 18 heterogenous reports. Meta-regression indicated greater lifetime prevalence of RC-BD among women than men (p=0.003). Association of RC-BD with suicide attempts, and unsatisfactory response to mood-stabilizers was supported by strong evidence (Level 1); associations with childhood maltreatment, mixed-features, female sex, and type-II BD had moderate evidence (Level 2). Other factors: genetic predisposition, metabolic disturbances or hypothyroidism, antidepressant exposure, predominant depressive polarity (Level 3), along with greater illness duration and immune-inflammatory dysfunction (Level 4) require further study. RC-BD was consistently recognized as having high prevalence (22.3%-35.5% of BD cases) and inferior treatment response. Identified associated factors can inform clinical practice. Long-term illness-course, metabolic factors, and optimal treatment require further investigation.

Quantile Regression for biomarkers in the UK Biobank

Genome-wide association studies (GWAS) for biomarkers important for clinical phenotypes can lead to clinically relevant discoveries. GWAS for quantitative traits are based on simplified regression models modeling the conditional mean of a phenotype as a linear function of genotype. An alternative and easy to apply approach is quantile regression that naturally extends linear regression to the analysis of the entire conditional distribution of a phenotype of interest by modeling conditional quantiles within a regression framework. Quantile regression can be applied efficiently at biobank scale using standard statistical packages in much the same way as linear regression, while having some unique advantages such as identifying variants with heterogeneous effects across different quantiles, including non-additive effects and variants involved in gene-environment interactions; accommodating a wide range of phenotype distributions with invariance to trait transformation; and overall providing more detailed information about the underlying genotype-phenotype associations. Here, we demonstrate the value of quantile regression in the context of GWAS by applying it to 39 quantitative traits in the UK Biobank (n > 300 , 000 individuals). Across these 39 traits we identify 7,297 significant loci, including 259 loci only detected by quantile regression. We show that quantile regression can help uncover replicable but unmodelled gene-environment interactions, and can provide additional key insights into poorly understood genotype-phenotype correlations for clinically relevant biomarkers at minimal additional cost.

Genome-wide association meta-analysis of knee and hip osteoarthritis uncovers genetic differences between patients treated with joint replacement and patients without joint replacement

OBJECTIVES

Osteoarthritis is a common and severe, multifactorial disease with a well-established genetic component. However, little is known about how genetics affect disease progression, and thereby the need for joint placement. Therefore, we aimed to investigate whether the genetic associations of knee and hip osteoarthritis differ between patients treated with joint replacement and patients without joint replacement.

METHODS

We included knee and hip osteoarthritis cases along with healthy controls, altogether counting >700 000 individuals. The cases were divided into two groups based on joint replacement status (surgical vs non-surgical) and included in four genome-wide association meta-analyses: surgical knee osteoarthritis (N = 22 525), non-surgical knee osteoarthritis (N = 38 626), surgical hip osteoarthritis (N = 20 221) and non-surgical hip osteoarthritis (N = 17 847). In addition, we tested for genetic correlation between the osteoarthritis groups and the pain phenotypes intervertebral disc disorder, dorsalgia, fibromyalgia, migraine and joint pain.

RESULTS

We identified 52 sequence variants associated with knee osteoarthritis (surgical: 17, non-surgical: 3) or hip osteoarthritis (surgical: 34, non-surgical: 1). For the surgical phenotypes, we identified 10 novel variants, including genes involved in autophagy (rs2447606 in ATG7) and mechanotransduction (rs202127176 in PIEZO1). One variant, rs13107325 in SLC39A8, associated more strongly with non-surgical knee osteoarthritis than surgical knee osteoarthritis. For all other variants, significance and effect sizes were higher for the surgical phenotypes. In contrast, genetic correlations with pain phenotypes tended to be stronger in the non-surgical groups.

CONCLUSIONS

Our results indicate differences in genetic associations between knee and hip osteoarthritis depending on joint replacement status.

An observational and Mendelian randomisation study on iron status and sepsis

Iron deficiency is associated with a substantial burden of morbidity. However, supplementation of iron has been linked to increased rates of serious infection in randomised trials of children in sub-Saharan Africa. Randomised trials in other settings have been inconclusive and it is unknown if changes in levels of iron biomarkers are linked to sepsis in these other settings. We used genetic variants associated with levels of iron biomarkers as instrumental variables in a Mendelian randomisation (MR) analysis to test the hypothesis that increasing levels of iron biomarkers increase the risk of sepsis. In observational and MR analyses we found that increases in iron biomarkers increase the odds of sepsis. In stratified analyses, we show that this risk may be larger in those with iron deficiency and/or anaemia. Taken together, results here suggest a required caution in supplementation of iron and underline the role of iron homeostasis in severe infection.

Familial risk of epithelial ovarian cancer after accounting for gynaecological surgery: a population-based study

BACKGROUND

Uptake of risk-reducing surgery has increased among women at high risk of epithelial ovarian cancer. We sought to characterise familial risk of epithelial ovarian cancer histotypes in a population-based study after accounting for gynaecological surgeries, including bilateral oophorectomy.

METHODS

We compared risk of epithelial ovarian cancer in relatives of 3536 epithelial ovarian cancer cases diagnosed in 1966-2016 and relatives of 35 326 matched controls. We used Cox competing risk models, incorporating bilateral oophorectomy as a competing risk, to estimate the relative risk of ovarian cancer in first-degree (FDR), second-degree (SDR) and third-degree (TDR) relatives from 1966 to 2016. We also estimated relative risks in time periods before (1966-1994, 1995-2004) and after (2005-2016) formal recommendations were made for prophylactic oophorectomy among women with pathogenic variants in BRCA1/2.

RESULTS

The relative risks of epithelial ovarian cancer in FDRs, SDRs and TDRs of cases versus controls were 1.68 (95% CI 1.39 to 2.04), 1.51 (95% CI 1.30 to 1.75) and 1.34 (95% CI 1.20 to 1.48), respectively. Relative risks were greatest for high-grade serous, mucinous and 'other epithelial' histotypes. Relative risks were attenuated for case FDRs, but not for SDRs or TDRs, from 2005 onwards, consistent with the timing of recommendations for prophylactic surgery.

CONCLUSION

Familial risk of epithelial ovarian cancer extends to TDRs, especially for high-grade serous and mucinous histotypes. Distant relatives share genes but minimal environment, highlighting the importance of germline inherited genetics in ovarian cancer aetiology. Increased ovarian cancer risk in distant relatives has implications for counselling and recommendations for prophylactic surgeries that, from our data, appear only to reach FDRs.

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.

Phenotypic subtyping via contrastive learning

Defining and accounting for subphenotypic structure has the potential to increase statistical power and provide a deeper understanding of the heterogeneity in the molecular basis of complex disease. Existing phenotype subtyping methods primarily rely on clinically observed heterogeneity or metadata clustering. However, they generally tend to capture the dominant sources of variation in the data, which often originate from variation that is not descriptive of the mechanistic heterogeneity of the phenotype of interest; in fact, such dominant sources of variation, such as population structure or technical variation, are, in general, expected to be independent of subphenotypic structure. We instead aim to find a subspace with signal that is unique to a group of samples for which we believe that subphenotypic variation exists (e.g., cases of a disease). To that end, we introduce Phenotype Aware Components Analysis (PACA), a contrastive learning approach leveraging canonical correlation analysis to robustly capture weak sources of subphenotypic variation. In the context of disease, PACA learns a gradient of variation unique to cases in a given dataset, while leveraging control samples for accounting for variation and imbalances of biological and technical confounders between cases and controls. We evaluated PACA using an extensive simulation study, as well as on various subtyping tasks using genotypes, transcriptomics, and DNA methylation data. Our results provide multiple strong evidence that PACA allows us to robustly capture weak unknown variation of interest while being calibrated and well-powered, far superseding the performance of alternative methods. This renders PACA as a state-of-the-art tool for defining de novo subtypes that are more likely to reflect molecular heterogeneity, especially in challenging cases where the phenotypic heterogeneity may be masked by a myriad of strong unrelated effects in the data.

E96V Mutation in the Kdelr3 Gene Is Associated with Type 2 Diabetes Susceptibility in Obese NZO Mice

Type 2 diabetes (T2D) represents a multifactorial metabolic disease with a strong genetic predisposition. Despite elaborate efforts in identifying the genetic variants determining individual susceptibility towards T2D, the majority of genetic factors driving disease development remain poorly understood. With the aim to identify novel T2D risk genes we previously generated an N2 outcross population using the two inbred mouse strains New Zealand obese (NZO) and C3HeB/FeJ (C3H). A linkage study performed in this population led to the identification of the novel T2D-associated quantitative trait locus (QTL) Nbg15 (NZO blood glucose on chromosome 15, Logarithm of odds (LOD) 6.6). In this study we used a combined approach of positional cloning, gene expression analyses and in silico predictions of DNA polymorphism on gene/protein function to dissect the genetic variants linking Nbg15 to the development of T2D. Moreover, we have generated congenic strains that associated the distal sublocus of Nbg15 to mechanisms altering pancreatic beta cell function. In this sublocus, Cbx6, Fam135b and Kdelr3 were nominated as potential causative genes associated with the Nbg15 driven effects. Moreover, a putative mutation in the Kdelr3 gene from NZO was identified, negatively influencing adaptive responses associated with pancreatic beta cell death and induction of endoplasmic reticulum stress. Importantly, knockdown of Kdelr3 in cultured Min6 beta cells altered insulin granules maturation and pro-insulin levels, pointing towards a crucial role of this gene in islets function and T2D susceptibility.

Melancholic features and typical neurovegetative symptoms of major depressive disorder show specific polygenic patterns

BACKGROUND

Major depressive disorder (MDD) is a highly prevalent psychiatric condition characterised by a heterogeneous clinical presentation and an estimated twin-based heritability of ~40-50 %. Different clinical MDD subtypes might partly reflect distinctive underlying genetics. This study aims to investigate if polygenic risk scores (PRSs) for different psychiatric disorders, personality traits, and substance use-related traits may be associated with different clinical subtypes of MDD (i.e., MDD with melancholic or psychotic features), higher symptom severity, or different clusters of depressive symptoms (i.e., sadness symptoms, typical neurovegetative symptoms, detachment symptoms, and negative thoughts).

METHODS

The target sample included 1149 patients with MDD, recruited by the European Group for the Study of Resistant Depression. PRSs for 25 psychiatric disorders and traits were computed based on the most recent publicly available summary statistics of the largest genome-wide association studies. PRSs were then used as predictors in regression models, adjusting for age, sex, population stratification, and recruitment sites.

RESULTS

Patients with MDD having higher PRS for MDD and loneliness were more likely to exhibit melancholic features of MDD (p = 0.0009 and p = 0.005, respectively). Moreover, patients with higher PRS for alcohol intake and post-traumatic stress disorder were more likely to experience greater typical neurovegetative symptoms (p = 0.0012 and p = 0.0045, respectively).

LIMITATIONS

The proportion of phenotypic variance explained by the PRSs was limited.

CONCLUSIONS

This study suggests that melancholic features and typical neurovegetative symptoms of MDD may show distinctive underlying genetics. Our findings provide a new contribution to the understanding of the genetic heterogeneity of MDD.

Genetic heterogeneity: Challenges, impacts, and methods through an associative lens

Genetic heterogeneity describes the occurrence of the same or similar phenotypes through different genetic mechanisms in different individuals. Robustly characterizing and accounting for genetic heterogeneity is crucial to pursuing the goals of precision medicine, for discovering novel disease biomarkers, and for identifying targets for treatments. Failure to account for genetic heterogeneity may lead to missed associations and incorrect inferences. Thus, it is critical to review the impact of genetic heterogeneity on the design and analysis of population level genetic studies, aspects that are often overlooked in the literature. In this review, we first contextualize our approach to genetic heterogeneity by proposing a high-level categorization of heterogeneity into "feature," "outcome," and "associative" heterogeneity, drawing on perspectives from epidemiology and machine learning to illustrate distinctions between them. We highlight the unique nature of genetic heterogeneity as a heterogeneous pattern of association that warrants specific methodological considerations. We then focus on the challenges that preclude effective detection and characterization of genetic heterogeneity across a variety of epidemiological contexts. Finally, we discuss systems heterogeneity as an integrated approach to using genetic and other high-dimensional multi-omic data in complex disease research.

Deconstructing a Syndrome: Genomic Insights Into PCOS Causal Mechanisms and Classification

Polycystic ovary syndrome (PCOS) is among the most common disorders in women of reproductive age, affecting up to 15% worldwide, depending on the diagnostic criteria. PCOS is characterized by a constellation of interrelated reproductive abnormalities, including disordered gonadotropin secretion, increased androgen production, chronic anovulation, and polycystic ovarian morphology. It is frequently associated with insulin resistance and obesity. These reproductive and metabolic derangements cause major morbidities across the lifespan, including anovulatory infertility and type 2 diabetes (T2D). Despite decades of investigative effort, the etiology of PCOS remains unknown. Familial clustering of PCOS cases has indicated a genetic contribution to PCOS. There are rare Mendelian forms of PCOS associated with extreme phenotypes, but PCOS typically follows a non-Mendelian pattern of inheritance consistent with a complex genetic architecture, analogous to T2D and obesity, that reflects the interaction of susceptibility genes and environmental factors. Genomic studies of PCOS have provided important insights into disease pathways and have indicated that current diagnostic criteria do not capture underlying differences in biology associated with different forms of PCOS. We provide a state-of-the-science review of genetic analyses of PCOS, including an overview of genomic methodologies aimed at a general audience of non-geneticists and clinicians. Applications in PCOS will be discussed, including strengths and limitations of each study. The contributions of environmental factors, including developmental origins, will be reviewed. Insights into the pathogenesis and genetic architecture of PCOS will be summarized. Future directions for PCOS genetic studies will be outlined.

Overlap between genetic variants associated with schizophrenia spectrum disorders and intelligence quotient: a systematic review

BACKGROUND

To study whether there is genetic overlap underlying the risk for schizophrenia spectrum disorders (SSDs) and low intelligence quotient (IQ), we reviewed and summarized the evidence on genetic variants associated with both traits.

METHODS

We performed this review in accordance with the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) and preregistered it in PROSPERO. We searched the Medline databases via PubMed, PsycInfo, Web of Science and Scopus. We included studies in adults with a diagnosis of SSD that explored genetic variants (single nucleotide polymorphisms [SNPs], copy number variants [CNVs], genomic insertions or genomic deletions), estimated IQ and studied the relationship between genetic variability and both traits (SSD and IQ). We synthesized the results and assessed risk of bias using the Quality of Genetic Association Studies (Q-Genie) tool.

RESULTS

Fifty-five studies met the inclusion criteria (45 case-control, 9 cross-sectional, 1 cohort), of which 55% reported significant associations for genetic variants involved in IQ and SSD. The SNPs more frequently explored through candidate gene studies were in COMT, DTNBP1, BDNF and TCF4. Through genome-wide association studies, 2 SNPs in CHD7 and GATAD2A were associated with IQ in patients with SSD. The studies on CNVs suggested significant associations between structural variants and low IQ in patients with SSD.

LIMITATIONS

Overall, primary studies used heterogeneous IQ measurement tools and had small samples. Grey literature was not screened.

CONCLUSION

Genetic overlap between SSD and IQ supports the neurodevelopmental hypothesis of schizophrenia. Most of the risk polymorphisms identified were in genes relevant to brain development, neural proliferation and differentiation, and synaptic plasticity.

Association Analysis of CMYA5 rs4704591 Polymorphism with Rheumatoid Arthritis in Caucasians

Einleitung Einzel nukleotid polymorphismen (SNPs) im Cardiomyopathy-Associated Protein 5 (CMYA5)-Gen wurden in genomweiten Assoziationsstudien mit rheumatoider Arthritis (RA) in Verbindung gebracht. In dieser Studie wollten wir die Assoziation zwischen CMYA5-Genpolymorphismen und RA in unabhängigen kaukasischen Fall-Kontroll-Kohorten replizieren und eine Metaanalyse durchführen, um die Assoziation von CMYA5-Genpolymorphismen mit RA in kaukasischen Populationen zu untersuchen.

Methoden Wir analysierten 2731 RA-Fälle und 1783 gesunde Kontrollen in vier unabhängigen kaukasischen Probensätzen. rs4704591 im CMYA5-Gen wurden unter Verwendung des TaqMan SNP-Genotypisierungsassays genotypisiert. Die Metaanalyse wurde über kaukasische Kohorten hinweg durchgeführt.

Ergebnisse Unsere Analyse ergab keine Hinweise auf eine Assoziation von rs4704591 mit RA in den Replikationsprobensätzen (P=0,941, OR=0,997). Die Metaanalyse zeigte eine schwache Assoziation zwischen dem kleinen Allel der CMYA5-Variante rs4704591 (C) und RA in den kombinierten RA-Kohorten (P=0,023, OR=0,938) unter Verwendung des logistischen Regressionsmodells in der Matched-Case-Control-Studie.

Schlussfolgerung Unsere Studie war nicht erfolgreich darin, die Assoziation der CMYA5-Variante rs4704591 mit RA zu replizieren. Daher können wir die Assoziation zwischen CMYA5-Genpolymorphismen und RA in der kaukasischen Bevölkerung nicht bestätigen.

Individual Genetic Heterogeneity

Genetic variation has been widely covered in literature, however, not from the perspective of an individual in any species. Here, a synthesis of genetic concepts and variations relevant for individual genetic constitution is provided. All the different levels of genetic information and variation are covered, ranging from whether an organism is unmixed or hybrid, has variations in genome, chromosomes, and more locally in DNA regions, to epigenetic variants or alterations in selfish genetic elements. Genetic constitution and heterogeneity of microbiota are highly relevant for health and wellbeing of an individual. Mutation rates vary widely for variation types, e.g., due to the sequence context. Genetic information guides numerous aspects in organisms. Types of inheritance, whether Mendelian or non-Mendelian, zygosity, sexual reproduction, and sex determination are covered. Functions of DNA and functional effects of variations are introduced, along with mechanism that reduce and modulate functional effects, including TARAR countermeasures and intraindividual genetic conflict. TARAR countermeasures for tolerance, avoidance, repair, attenuation, and resistance are essential for life, integrity of genetic information, and gene expression. The genetic composition, effects of variations, and their expression are considered also in diseases and personalized medicine. The text synthesizes knowledge and insight on individual genetic heterogeneity and organizes and systematizes the central concepts.

Genetic Influences on Cognitive Dysfunction in Schizophrenia

Schizophrenia is a severe and debilitating psychotic disorder that is highly heritable and relatively common in the population. The clinical heterogeneity associated with schizophrenia is substantial, with patients exhibiting a broad range of deficits and symptom severity. Large-scale genomic studies employing a case-control design have begun to provide some biological insight. However, this strategy combines individuals with clinically diverse symptoms and ignores the genetic risk that is carried by many clinically unaffected individuals. Consequently, the majority of the genetic architecture underlying schizophrenia remains unexplained, and the pathways by which the implicated variants contribute to the clinically observable signs and symptoms are still largely unknown. Parsing the complex, clinical phenotype of schizophrenia into biologically relevant components may have utility in research aimed at understanding the genetic basis of liability. Cognitive dysfunction is a hallmark symptom of schizophrenia that is associated with impaired quality of life and poor functional outcome. Here, we examine the value of quantitative measures of cognitive dysfunction to objectively target the underlying neurobiological pathways and identify genetic variants and gene networks contributing to schizophrenia risk. For a complex disorder, quantitative measures are also more efficient than diagnosis, allowing for the identification of associated genetic variants with fewer subjects. Such a strategy supplements traditional analyses of schizophrenia diagnosis, providing the necessary biological insight to help translate genetic findings into actionable treatment targets. Understanding the genetic basis of cognitive dysfunction in schizophrenia may thus facilitate the development of novel pharmacological and procognitive interventions to improve real-world functioning.

Stratification of rheumatoid arthritis cohort using Ayurveda based deep phenotyping approach identifies novel genes in a GWAS

Background and aim

Genome wide association studies have scaled up both in terms of sample size and range of complex disorders investigated, but these have explained relatively little phenotypic variance. Of the several reasons, phenotypic heterogeneity seems to be a likely contributor for missing out genetic associations of large effects. Ayurveda, the traditional Indian system of medicine is one such tool which adopts a holistic deep phenotyping approach and classifies individuals based on their body constitution/prakriti. We hypothesized that Ayurveda based phenotypic stratification of healthy and diseased individuals will allow us to achieve much desired homogeneous cohorts which would facilitate detection of genetic association of large effects. In this proof of concept study, we performed a genome wide association testing of clinically diagnosed rheumatoid arthritis patients and healthy controls, who were re-phenotyped into Vata, Pitta and Kapha predominant prakriti sub-groups.

Experimental procedure

Genotypes of rheumatoid arthritis cases (Vata = 49; Pitta = 117; Kapha = 78) and controls (Vata = 33; Pitta = 175; Kapha = 85) were retrieved from the total genotype data, used in a recent genome-wide association study performed in our laboratory. A total of 528461 SNPs were included after quality control. Prakriti-wise genome-wide association analysis was employed.

Results and conclusion

This study identified (i) prakriti-specific novel disease risk genes of high effect sizes; (ii) putative candidates of novel therapeutic potential; and (iii) a good correlation between genetic findings and clinical knowledge in Ayurveda. Adopting Ayurveda based deep phenotyping may facilitate explaining hitherto undiscovered heritability in complex traits and may propel much needed progress in personalized medicine.

Deep phenotyping for precision medicine in Parkinson's disease

A major challenge in medical genomics is to understand why individuals with the same disorder have different clinical symptoms and why those who carry the same mutation may be affected by different disorders. In every complex disorder, identifying the contribution of different genetic and non-genetic risk factors is a key obstacle to understanding disease mechanisms. Genetic studies rely on precise phenotypes and are unable to uncover the genetic contributions to a disorder when phenotypes are imprecise. To address this challenge, deeply phenotyped cohorts have been developed for which detailed, fine-grained data have been collected. These cohorts help us to investigate the underlying biological pathways and risk factors to identify treatment targets, and thus to advance precision medicine. The neurodegenerative disorder Parkinson's disease has a diverse phenotypical presentation and modest heritability, and its underlying disease mechanisms are still being debated. As such, considerable efforts have been made to develop deeply phenotyped cohorts for this disorder. Here, we focus on Parkinson's disease and explore how deep phenotyping can help address the challenges raised by genetic and phenotypic heterogeneity. We also discuss recent methods for data collection and computation, as well as methodological challenges that have to be overcome.

Melatonin and aggressive behavior: A systematic review of the literature on preclinical and clinical evidence

The melatonin system and circadian disruption have well-established links with aggressive behaviors; however, the biological underpinnings have not been thoroughly investigated. Here, we aimed at examining the current knowledge regarding the neurobiological and psychopharmacological involvement of the melatonin system in aggressive/violent behaviors. To this end, we performed a systematic review on Embase and Pubmed/MEDLINE of preclinical and clinical evidence linking the melatonin system, melatonin, and melatoninergic drugs with aggressive/violent behaviors. Two blinded raters performed an independent screening of the relevant literature. Overall, this review included 38 papers distributed between clinical and preclinical models. Eleven papers specifically addressed the existing evidence in rodent models, five in fish models, and 21 in humans. The data indicate that depending on the species, model, and timing of administration, melatonin may exert a complex influence on aggressive/violent behaviors. Particularly, the apparent contrasting findings on the link between the melatonin system and aggression/violence (with either increased, no, or decreased effect) shown in preclinical models underscore the need for further research to develop more accurate and fruitful translational models. Likewise, the significant heterogeneity found in the results of clinical studies does not allow yet to draw any firm conclusion on the efficacy of melatonin or melatonergic drugs on aggressive/violent behaviors. However, findings in children and in traits associated with aggressive/violent behavior, including irritability and anger, are emerging and deserve empirical attention given the low toxicity of melatonin and melatonergic drugs.

Digital Phenotyping in Clinical Neurology

Internet-connected devices, including personal computers, smartphones, smartwatches, and voice assistants, have evolved into powerful multisensor technologies that billions of people interact with daily to connect with friends and colleagues, access and share information, purchase goods, play games, and navigate their environment. Digital phenotyping taps into the data streams captured by these devices to characterize and understand health and disease. The purpose of this article is to summarize opportunities for digital phenotyping in neurology, review studies using everyday technologies to obtain motor and cognitive information, and provide a perspective on how neurologists can embrace and accelerate progress in this emerging field.

Genetic Regulatory Networks of Apolipoproteins and Associated Medical Risks

Apolipoproteins (APO proteins) are the lipoprotein family proteins that play key roles in transporting lipoproteins all over the body. There are nearly more than twenty members reported in the APO protein family, among which the A, B, C, E, and L play major roles in contributing genetic risks to several disorders. Among these genetic risks, the single nucleotide polymorphisms (SNPs), involving the variation of single nucleotide base pairs, and their contributing polymorphisms play crucial roles in the apolipoprotein family and its concordant disease heterogeneity that have predominantly recurred through the years. In this review, we have contributed a handful of information on such genetic polymorphisms that include APOE, ApoA1/B ratio, and A1/C3/A4/A5 gene cluster-based population genetic studies carried throughout the world, to elaborately discuss the effects of various genetic polymorphisms in imparting various medical conditions, such as obesity, cardiovascular, stroke, Alzheimer's disease, diabetes, vascular complications, and other associated risks.

Hyperfocus or flow? Attentional strengths in autism spectrum disorder

The comorbidity of autism spectrum disorder (ASD) and attention-deficit/hyperactivity disorder (ADHD) diagnoses is well established. An ASD diagnosis is associated with elevated ADHD traits and symptoms, as well as strengths in attention. In the ASD literature, attentional strengths have been described as maladaptive (e.g., hyperfocus), in contrast with positive portrayals in the typically developing population (e.g., flow). The objective of this study was to (1) compare profiles of attentional strengths and weaknesses in ASD and ADHD and (2) determine whether attentional strengths in ASD are associated with impairment, poorer cognitive flexibility, and perseveration/perfectionism. In a community sample of 5,744 children and youth, 131 children were reported as having a diagnosis of ASD (mean age 10.3 years) and 346 children were reported as having a diagnosis of ADHD (mean age 10.7 years). We used the Strengths and Weaknesses of Attention-Deficit/Hyperactivity-symptoms and Normal-behaviors (SWAN) rating scale to calculate attentional and hyperactive/impulse control strength and weakness counts and scores. The Autism-Spectrum Quotient Switching factor served as a measure of cognitive flexibility. Impairment was assessed with the Columbia Impairment Scale. We used the symmetry/ordering factor on the Toronto Obsessive-Compulsive Scale as a measure of perseveration/perfectionism. No differences were found between the ADHD and ASD groups in SWAN weakness scores, symptoms, or hyperactive/impulse control strengths; however, autistic children had higher rates of attentional strengths [odds ratio: 5.7, 95% CI (2.8, 11.6), p < 0.0001]. Post-hoc pairwise testing identified four attentional strengths with significantly higher rates in ASD than in ADHD. Attentional strength scores were not associated with impairment or poor cognitive flexibility, but predicted levels of perseveration/perfectionism. The effect of attentional strengths on impairment and cognitive flexibility did not differ between autistic and Control children, but the higher perseveration/perfectionism scores seen in ASD were not found in Control children. ASD is associated with a pattern of attentional strengths that is not found in ADHD Characterizing the full range of attentional abilities in autistic children may explain variability in outcomes such as quality-of-life indicators and identify protective factors, providing targets for strength-based behavioral interventions. The clinical and etiological implications of the subgroup of autistic children with attentional strengths require further investigation.

Comparison of symptom-based versus self-reported diagnostic measures of anxiety and depression disorders in the GLAD and COPING cohorts

BACKGROUND

Understanding and improving outcomes for people with anxiety or depression often requires large sample sizes. To increase participation and reduce costs, such research is typically unable to utilise "gold-standard" methods to ascertain diagnoses, instead relying on remote, self-report measures.

AIMS

Assess the comparability of remote diagnostic methods for anxiety and depression disorders commonly used in research.

METHOD

Participants from the UK-based GLAD and COPING NBR cohorts (N = 58,400) completed an online questionnaire between 2018 and 2020. Responses to detailed symptom reports were compared to DSM-5 criteria to generate symptom-based diagnoses of major depressive disorder (MDD), generalised anxiety disorder (GAD), specific phobia, social anxiety disorder, panic disorder, and agoraphobia. Participants also self-reported any prior diagnoses from health professionals, termed self-reported diagnoses. "Any anxiety" included participants with at least one anxiety disorder. Agreement was assessed by calculating accuracy, Cohen's kappa, McNemar's chi-squared, sensitivity, and specificity.

RESULTS

Agreement between diagnoses was moderate for MDD, any anxiety, and GAD, but varied by cohort. Agreement was slight to fair for the phobic disorders. Many participants with self-reported GAD did not receive a symptom-based diagnosis. In contrast, symptom-based diagnoses of the phobic disorders were more common than self-reported diagnoses.

CONCLUSIONS

Agreement for MDD, any anxiety, and GAD was higher for cases in the case-enriched GLAD cohort and for controls in the general population COPING NBR cohort. For anxiety disorders, self-reported diagnoses classified most participants as having GAD, whereas symptom-based diagnoses distributed participants more evenly across the anxiety disorders. Further validation against gold standard measures is required.

Pathfinder: a gamified measure to integrate general cognitive ability into the biological, medical, and behavioural sciences

Genome-wide association (GWA) studies have uncovered DNA variants associated with individual differences in general cognitive ability (g), but these are far from capturing heritability estimates obtained from twin studies. A major barrier to finding more of this 'missing heritability' is assessment--the use of diverse measures across GWA studies as well as time and the cost of assessment. In a series of four studies, we created a 15-min (40-item), online, gamified measure of g that is highly reliable (alpha = 0.78; two-week test-retest reliability = 0.88), psychometrically valid and scalable; we called this new measure Pathfinder. In a fifth study, we administered this measure to 4,751 young adults from the Twins Early Development Study. This novel g measure, which also yields reliable verbal and nonverbal scores, correlated substantially with standard measures of g collected at previous ages (r ranging from 0.42 at age 7 to 0.57 at age 16). Pathfinder showed substantial twin heritability (0.57, 95% CIs = 0.43, 0.68) and SNP heritability (0.37, 95% CIs = 0.04, 0.70). A polygenic score computed from GWA studies of five cognitive and educational traits accounted for 12% of the variation in g, the strongest DNA-based prediction of g to date. Widespread use of this engaging new measure will advance research not only in genomics but throughout the biological, medical, and behavioural sciences.

Beyond the average patient: how neuroimaging models can address heterogeneity in dementia

Dementia is a highly heterogeneous condition, with pronounced individual differences in age of onset, clinical presentation, progression rates and neuropathological hallmarks, even within a specific diagnostic group. However, the most common statistical designs used in dementia research studies and clinical trials overlook this heterogeneity, instead relying on comparisons of group average differences (e.g. patient versus control or treatment versus placebo), implicitly assuming within-group homogeneity. This one-size-fits-all approach potentially limits our understanding of dementia aetiology, hindering the identification of effective treatments. Neuroimaging has enabled the characterization of the average neuroanatomical substrates of dementias; however, the increasing availability of large open neuroimaging datasets provides the opportunity to examine patterns of neuroanatomical variability in individual patients. In this update, we outline the causes and consequences of heterogeneity in dementia and discuss recent research that aims to tackle heterogeneity directly, rather than assuming that dementia affects everyone in the same way. We introduce spatial normative modelling as an emerging data-driven technique, which can be applied to dementia data to model neuroanatomical variation, capturing individualized neurobiological 'fingerprints'. Such methods have the potential to detect clinically relevant subtypes, track an individual's disease progression or evaluate treatment responses, with the goal of moving towards precision medicine for dementia.

Polygenic scores differentially predict developmental trajectories of subtypes of social withdrawal in childhood

BACKGROUND

Children who consistently withdraw from social situations face increased risk for later socioemotional difficulties. Twin studies indicate that genetic factors substantially account for the persistence of social withdrawal over time. However, the molecular genetic etiology of chronic courses of social wariness and preference for solitude, two dimensions of social withdrawal, remains undocumented. The objectives of the present study were (a) to identify high-risk trajectories for social wariness and preference for solitude in childhood and (b) to examine whether falling on these high-risk trajectories can be predicted by specific polygenic scores for mental health traits and disorders and by a general polygenic predisposition to these traits.

METHODS

Teachers evaluated 971 genotyped children at five occasions (age 6 to 12 years) from two prospective longitudinal studies, the Quebec Newborn Twin Study and the Quebec Longitudinal Study of Child Development. Developmental trajectories for social wariness and preference for solitude were identified. We tested whether polygenic scores for attention deficit hyperactivity disorder, autism spectrum disorder, depression, loneliness, and subjective well-being, as well as a general mental health genetic risk score derived across these traits, were associated with the developmental trajectories.

RESULTS

Polygenic scores differentially predicted social wariness and preference for solitude. Only the loneliness polygenic score significantly predicted the high trajectory for social wariness. By contrast, the general mental health genetic risk score factor was associated with the trajectory depicting high-chronic preference for solitude.

CONCLUSIONS

Distinct associations were uncovered between the polygenic scores, social wariness, and preference for solitude.

Meta-Analysis of Joint Test of SNP and SNP-Environment Interaction with Heterogeneity

Many complex diseases are caused by single nucleotide polymorphisms (SNPs), environmental factors, and the interaction between SNPs and environment. Joint tests of the SNP and SNP-environment interaction effects (JMA) and meta-regression (MR) are commonly used to evaluate these SNP-environment interactions. However, these two methods do not consider genetic heterogeneity. We previously presented a random-effect MR, which provided higher power than the MR in datasets with high heterogeneity. However, this method requires group-level data, which sometimes are not available. Given this, we designed this study to evaluate the introduction of the random effects of SNP and SNP-environment interaction into the JMA, and then extended this to the random effect model. Likelihood ratio statistic is applied to test the JMA and the new method we proposed in this paper. We evaluated the null distributions of these tests, and the powers for this method. This method was verified by simulation and was shown to provide similar powers to the random effect meta-regression method (RMR). However, this method only requires study-level data which relaxed the condition of the RMR. Our study suggests that this method is more suitable for finding the association between SNP and diseases in the absence of group-level data.

Genetics of anorexia nervosa: an overview of genome-wide association studies and emerging biological links

Anorexia nervosa (AN) is a complex disorder with a strong genetic component. Comorbidities are frequent and there is substantial overlap with other disorders. The lack of understanding of the molecular and neuroanatomical causes has made it difficult to develop effective treatments and it is often difficult to treat in clinical practice. Recent advances in genetics have changed our understanding of polygenic diseases, increasing the possibility of understanding better how molecular pathways are intertwined. This review synthetizes the current state of genetic research providing an overview of genome-wide association studies (GWAS) findings in AN as well as overlap with other disorders, traits, pathways, and imaging results. This paper also discusses the different putative global pathways that are contributing to the disease including the evidence for metabolic and psychiatric origin of the disease.

Deciphering osteoarthritis genetics across 826,690 individuals from 9 populations

Osteoarthritis affects over 300 million people worldwide. Here, we conduct a genome-wide association study meta-analysis across 826,690 individuals (177,517 with osteoarthritis) and identify 100 independently associated risk variants across 11 osteoarthritis phenotypes, 52 of which have not been associated with the disease before. We report thumb and spine osteoarthritis risk variants and identify differences in genetic effects between weight-bearing and non-weight-bearing joints. We identify sex-specific and early age-at-onset osteoarthritis risk loci. We integrate functional genomics data from primary patient tissues (including articular cartilage, subchondral bone, and osteophytic cartilage) and identify high-confidence effector genes. We provide evidence for genetic correlation with phenotypes related to pain, the main disease symptom, and identify likely causal genes linked to neuronal processes. Our results provide insights into key molecular players in disease processes and highlight attractive drug targets to accelerate translation.

Infertility due to defective sperm flagella caused by an intronic deletion in DNAH17 that perturbs splicing

Artificial insemination in pig (Sus scrofa domesticus) breeding involves the evaluation of the semen quality of breeding boars. Ejaculates that fulfill predefined quality requirements are processed, diluted and used for inseminations. Within short time, eight Swiss Large White boars producing immotile sperm that had multiple morphological abnormalities of the sperm flagella were noticed at a semen collection center. The eight boars were inbred on a common ancestor suggesting that the novel sperm flagella defect is a recessive trait. Transmission electron microscopy cross-sections revealed that the immotile sperm had disorganized flagellar axonemes. Haplotype-based association testing involving microarray-derived genotypes at 41,094 SNPs of six affected and 100 fertile boars yielded strong association (P = 4.22 × 10-15) at chromosome 12. Autozygosity mapping enabled us to pinpoint the causal mutation on a 1.11 Mb haplotype located between 3,473,632 and 4,587,759 bp. The haplotype carries an intronic 13-bp deletion (Chr12:3,556,401-3,556,414 bp) that is compatible with recessive inheritance. The 13-bp deletion excises the polypyrimidine tract upstream exon 56 of DNAH17 (XM_021066525.1: c.8510-17_8510-5del) encoding dynein axonemal heavy chain 17. Transcriptome analysis of the testis of two affected boars revealed that the loss of the polypyrimidine tract causes exon skipping which results in the in-frame loss of 89 amino acids from DNAH17. Disruption of DNAH17 impairs the assembly of the flagellar axoneme and manifests in multiple morphological abnormalities of the sperm flagella. Direct gene testing may now be implemented to monitor the defective allele in the Swiss Large White population and prevent the frequent manifestation of a sterilizing sperm tail disorder in breeding boars.

Genome-wide association studies of low back pain and lumbar spinal disorders using electronic health record data identify a locus associated with lumbar spinal stenosis

ABSTRACT

Identifying genetic risk factors for lumbar spine disorders may lead to knowledge regarding underlying mechanisms and the development of new treatments. We conducted a genome-wide association study involving 100,811 participants with genotypes and longitudinal electronic health record data from the Electronic Medical Records and Genomics Network and Geisinger Health. Cases and controls were defined using validated algorithms and clinical diagnostic codes. Electronic health record-defined phenotypes included low back pain requiring healthcare utilization (LBP-HC), lumbosacral radicular syndrome (LSRS), and lumbar spinal stenosis (LSS). Genome-wide association study used logistic regression with additive genetic effects adjusting for age, sex, site-specific factors, and ancestry (principal components). A fixed-effect inverse-variance weighted meta-analysis was conducted. Genetic variants of genome-wide significance (P < 5 × 10-8) were carried forward for replication in an independent sample from UK Biobank. Phenotype prevalence was 48.8% for LBP-HC, 19.8% for LSRS, and 7.9% for LSS. No variants were significantly associated with LBP-HC. One locus was associated with LSRS (lead variant rs146153280:C>G, odds ratio [OR] = 1.17 for G, P = 2.1 × 10-9), but was not replicated. Another locus on chromosome 2 spanning GFPT1, NFU1, and AAK1 was associated with LSS (lead variant rs13427243:G>A, OR = 1.10 for A, P = 4.3 × 10-8) and replicated in UK Biobank (OR = 1.11, P = 5.4 × 10-5). This was the first genome-wide association study meta-analysis of lumbar spinal disorders using electronic health record data. We identified 2 novel associations with LSRS and LSS; the latter was replicated in an independent sample.

Depression and interleukin-6 signaling: A Mendelian Randomization study

BACKGROUND

A large body of research has reported associations between depression and elevated interleukin-6 (IL-6), a cytokine with several roles including pro-inflammatory signaling. The nature and directionality of this relationship are not yet clear. In this study we use Mendelian Randomization to examine the possibility of a causal relationship between IL-6 and depressive symptoms, and to explore multiple signaling pathways that could serve as mechanisms for this relationship.

METHODS

This study uses a two-sample Mendelian Randomization design. Data come from the UK Biobank (n = 89,119) and published summary statistics from six existing GWAS analyses. The primary analysis focuses on the soluble interleukin-6 receptor (sIL-6R), which is involved in multiple signaling pathways. Exploratory analyses use C-reactive protein (CRP) and soluble glycoprotein 130 (sgp130) to further examine potential underlying mechanisms.

RESULTS

Results are consistent with a causal effect of sIL-6R on depression (PCA-IVW Odds Ratio: 1.023 (95% Confidence Interval: 1.006-1.039), p = 0.006). Exploratory analyses demonstrate that the relationship could be consistent with either decreased classical signaling or increased trans signaling as the underlying mechanism.

DISCUSSION

These results strengthen the body evidence implicating IL-6 signaling in depression. When compared with existing observational and animal findings, the direction of these results suggests involvement of IL-6 trans signaling. Further study is needed to examine whether IL6R genetic variants might influence IL-6 trans signaling in the brain, as well as to explore other potential pathways linking depression and inflammation.

A Type 2 Diabetes Subtype Responsive to ACCORD Intensive Glycemia Treatment

OBJECTIVE

Current type 2 diabetes (T2D) management contraindicates intensive glycemia treatment in patients with high cardiovascular disease (CVD) risk and is partially motivated by evidence of harms in the Action to Control Cardiovascular Risk in Diabetes (ACCORD) trial. Heterogeneity in response to intensive glycemia treatment has been observed, suggesting potential benefit for some individuals.

RESEARCH DESIGN AND METHODS

ACCORD was a randomized controlled trial that investigated whether intensively treating glycemia in individuals with T2D would reduce CVD outcomes. Using a novel approach to cluster HbA1c trajectories, we identified groups in the intensive glycemia arm with modified CVD risk. Genome-wide analysis and polygenic score (PS) were developed to predict group membership. Mendelian randomization was performed to infer causality.

RESULTS

We identified four clinical groupings in the intensive glycemia arm, and clinical group 4 (C4) displayed fewer CVD (hazard ratio [HR] 0.34; P = 2.01 × 10-3) and microvascular outcomes (HR 0.86; P = 0.015) than those receiving standard treatment. A single-nucleotide polymorphism, rs220721, in MAS1 reached suggestive significance in C4 (P = 4.34 × 10-7). PS predicted C4 with high accuracy (area under the receiver operating characteristic curve 0.98), and this predicted C4 displayed reduced CVD risk with intensive versus standard glycemia treatment (HR 0.53; P = 4.02 × 10-6), but not reduced risk of microvascular outcomes (P < 0.05). Mendelian randomization indicated causality between PS, on-trial HbA1c, and reduction in CVD outcomes (P < 0.05).

CONCLUSIONS

We found evidence of a T2D clinical group in ACCORD that benefited from intensive glycemia treatment, and membership in this group could be predicted using genetic variants. This study generates new hypotheses with implications for precision medicine in T2D and represents an important development in this landmark clinical trial warranting further investigation.

Genome-wide association mapping in maize: status and prospects

Genome-wide association study (GWAS) provides a robust and potent tool to retrieve complex phenotypic traits back to their underlying genetics. Maize is an excellent crop for performing GWAS due to diverse genetic variability, rapid decay of linkage disequilibrium, availability of distinct sub-populations and abundant SNP information. The application of GWAS in maize has resulted in successful identification of thousands of genomic regions associated with many abiotic and biotic stresses. Many agronomic and quality traits of maize are severely affected by such stresses and, significantly affecting its growth and productivity. To improve productivity of maize crop in countries like India which contribute only 2% to the world's total production in 2019-2020, it is essential to understand genetic complexity of underlying traits. Various DNA markers and trait associations have been revealed using conventional linkage mapping methods. However, it has achieved limited success in improving polygenic complex traits due to lower resolution of trait mapping. The present review explores the prospects of GWAS in improving yield, quality and stress tolerance in maize besides, strengths and challenges of using GWAS for molecular breeding and genomic selection. The information gathered will facilitate elucidation of genetic mechanisms of complex traits and improve efficiency of marker-assisted selection in maize breeding.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-02799-4.