Adaptive and Innate Immunity Are Key Drivers of Age at Onset of Multiple Sclerosis.

Background and Objectives: Multiple sclerosis (MS) age at onset (AAO) is a clinical predictor of long-term disease outcomes, independent of disease duration. Little is known about the genetic and biological mechanisms underlying age of first symptoms. We conducted a genome-wide association study (GWAS) to investigate associations between individual genetic variation and the MS AAO phenotype. Methods: The study population was comprised participants with MS in 6 clinical trials: ADVANCE (N = 655; relapsing-remitting [RR] MS), ASCEND (N = 555; secondary-progressive [SP] MS), DECIDE (N = 1,017; RRMS), OPERA1 (N = 581; RRMS), OPERA2 (N = 577; RRMS), and ORATORIO (N = 529; primary-progressive [PP] MS). Altogether, 3,905 persons with MS of European ancestry were analyzed. GWAS were conducted for MS AAO in each trial using linear additive models controlling for sex and 10 principal components. Resultant summary statistics across the 6 trials were then meta-analyzed, for a total of 8.3 × 10-6 single nucleotide polymorphisms (SNPs) across all trials after quality control and filtering for heterogeneity. Gene-based tests of associations, pathway enrichment analyses, and Mendelian randomization analyses for select exposures were also performed. Results: Four lead SNPs within 2 loci were identified (p < 5 × 10-8), including a) 3 SNPs in the major histocompatibility complex and their effects were independent of HLA-DRB1*15:01 and b) a LOC105375167 variant on chromosome 7. At the gene level, the top association was HLA-C (p = 1.2 × 10-7), which plays an important role in antiviral immunity. Functional annotation revealed the enrichment of pathways related to T-cell receptor signaling, autoimmunity, and the complement cascade. Mendelian randomization analyses suggested a link between both earlier age at puberty and shorter telomere length and earlier AAO, while there was no evidence for a role for either body mass index or vitamin D levels. Discussion: Two genetic loci associated with MS AAO were identified, and functional annotation demonstrated an enrichment of genes involved in adaptive and complement immunity. There was also evidence supporting a link with age at puberty and telomere length. The findings suggest that AAO in MS is multifactorial, and the factors driving onset of symptoms overlap with those influencing MS risk.

Cross-ancestry atlas of gene, isoform, and splicing regulation in the developing human brain.

Neuropsychiatric genome-wide association studies (GWASs), including those for autism spectrum disorder and schizophrenia, show strong enrichment for regulatory elements in the developing brain. However, prioritizing risk genes and mechanisms is challenging without a unified regulatory atlas. Across 672 diverse developing human brains, we identified 15,752 genes harboring gene, isoform, and/or splicing quantitative trait loci, mapping 3739 to cellular contexts. Gene expression heritability drops during development, likely reflecting both increasing cellular heterogeneity and the intrinsic properties of neuronal maturation. Isoform-level regulation, particularly in the second trimester, mediated the largest proportion of GWAS heritability. Through colocalization, we prioritized mechanisms for about 60% of GWAS loci across five disorders, exceeding adult brain findings. Finally, we contextualized results within gene and isoform coexpression networks, revealing the comprehensive landscape of transcriptome regulation in development and disease.

Whole-exome sequencing in UK Biobank reveals rare genetic architecture for depression.

Nearly two hundred common-variant depression risk loci have been identified by genome-wide association studies (GWAS). However, the impact of rare coding variants on depression remains poorly understood. Here, we present whole-exome sequencing analyses of depression with seven different definitions based on survey, questionnaire, and electronic health records in 320,356 UK Biobank participants. We showed that the burden of rare damaging coding variants in loss-of-function intolerant genes is significantly associated with risk of depression with various definitions. We compared the rare and common genetic architecture across depression definitions by genetic correlation and showed different genetic relationships between definitions across common and rare variants. In addition, we demonstrated that the effects of rare damaging coding variant burden and polygenic risk score on depression risk are additive. The gene set burden analyses revealed overlapping rare genetic variant components with developmental disorder, autism, and schizophrenia. Our study provides insights into the contribution of rare coding variants, separately and in conjunction with common variants, on depression with various definitions and their genetic relationships with neurodevelopmental disorders.

PICALO: principal interaction component analysis for the identification of discrete technical, cell-type, and environmental factors that mediate eQTLs.

Expression quantitative trait loci (eQTL) offer insights into the regulatory mechanisms of trait-associated variants, but their effects often rely on contexts that are unknown or unmeasured. We introduce PICALO, a method for hidden variable inference of eQTL contexts. PICALO identifies and disentangles technical from biological context in heterogeneous blood and brain bulk eQTL datasets. These contexts are biologically informative and reproducible, outperforming cell counts or expression-based principal components. Furthermore, we show that RNA quality and cell type proportions interact with thousands of eQTLs. Knowledge of hidden eQTL contexts may aid in the inference of functional mechanisms underlying disease variants.

Rare genetic variants impact muscle strength.

Muscle strength is highly heritable and predictive for multiple adverse health outcomes including mortality. Here, we present a rare protein-coding variant association study in 340,319 individuals for hand grip strength, a proxy measure of muscle strength. We show that the exome-wide burden of rare protein-truncating and damaging missense variants is associated with a reduction in hand grip strength. We identify six significant hand grip strength genes, KDM5B, OBSCN, GIGYF1, TTN, RB1CC1, and EIF3J. In the example of the titin (TTN) locus we demonstrate a convergence of rare with common variant association signals and uncover genetic relationships between reduced hand grip strength and disease. Finally, we identify shared mechanisms between brain and muscle function and uncover additive effects between rare and common genetic variation on muscle strength.

The impact of rare protein coding genetic variation on adult cognitive function.

Compelling evidence suggests that human cognitive function is strongly influenced by genetics. Here, we conduct a large-scale exome study to examine whether rare protein-coding variants impact cognitive function in the adult population (n = 485,930). We identify eight genes (ADGRB2, KDM5B, GIGYF1, ANKRD12, SLC8A1, RC3H2, CACNA1A and BCAS3) that are associated with adult cognitive function through rare coding variants with large effects. Rare genetic architecture for cognitive function partially overlaps with that of neurodevelopmental disorders. In the case of KDM5B we show how the genetic dosage of one of these genes may determine the variability of cognitive, behavioral and molecular traits in mice and humans. We further provide evidence that rare and common variants overlap in association signals and contribute additively to cognitive function. Our study introduces the relevance of rare coding variants for cognitive function and unveils high-impact monogenic contributions to how cognitive function is distributed in the normal adult population.

Cross-ancestry, cell-type-informed atlas of gene, isoform, and splicing regulation in the developing human brain.

Genomic regulatory elements active in the developing human brain are notably enriched in genetic risk for neuropsychiatric disorders, including autism spectrum disorder (ASD), schizophrenia, and bipolar disorder. However, prioritizing the specific risk genes and candidate molecular mechanisms underlying these genetic enrichments has been hindered by the lack of a single unified large-scale gene regulatory atlas of human brain development. Here, we uniformly process and systematically characterize gene, isoform, and splicing quantitative trait loci (xQTLs) in 672 fetal brain samples from unique subjects across multiple ancestral populations. We identify 15,752 genes harboring a significant xQTL and map 3,739 eQTLs to a specific cellular context. We observe a striking drop in gene expression and splicing heritability as the human brain develops. Isoform-level regulation, particularly in the second trimester, mediates the greatest proportion of heritability across multiple psychiatric GWAS, compared with eQTLs. Via colocalization and TWAS, we prioritize biological mechanisms for ~60% of GWAS loci across five neuropsychiatric disorders, nearly two-fold that observed in the adult brain. Finally, we build a comprehensive set of developmentally regulated gene and isoform co-expression networks capturing unique genetic enrichments across disorders. Together, this work provides a comprehensive view of genetic regulation across human brain development as well as the stage-and cell type-informed mechanistic underpinnings of neuropsychiatric disorders.

Brain expression quantitative trait locus and network analyses reveal downstream effects and putative drivers for brain-related diseases.

Identification of therapeutic targets from genome-wide association studies (GWAS) requires insights into downstream functional consequences. We harmonized 8,613 RNA-sequencing samples from 14 brain datasets to create the MetaBrain resource and performed cis- and trans-expression quantitative trait locus (eQTL) meta-analyses in multiple brain region- and ancestry-specific datasets (n ≤ 2,759). Many of the 16,169 cortex cis-eQTLs were tissue-dependent when compared with blood cis-eQTLs. We inferred brain cell types for 3,549 cis-eQTLs by interaction analysis. We prioritized 186 cis-eQTLs for 31 brain-related traits using Mendelian randomization and co-localization including 40 cis-eQTLs with an inferred cell type, such as a neuron-specific cis-eQTL (CYP24A1) for multiple sclerosis. We further describe 737 trans-eQTLs for 526 unique variants and 108 unique genes. We used brain-specific gene-co-regulation networks to link GWAS loci and prioritize additional genes for five central nervous system diseases. This study represents a valuable resource for post-GWAS research on central nervous system diseases.

Defining healthcare never events to effect system change: A protocol for systematic review.

INTRODUCTION: A never event is the most egregious of patient safety incidents. It refers to events that should theoretically never happen, such amputating the wrong limb. The term "never event" is used around the world by a variety of medical and patient safety organizations and is synonymous with sentinel events and serious reportable events. Unfortunately, there is little consensus about which events, in particular, are never events. These differing lists hinder potential collaboration or large-scale analyses. A recent systematic review by Hegarty et al. (2020) identified the need for a standardized definition for serious reportable events. The objective of our systematic review is to build on this by identifying which events are consistently or frequently identified as never events in order to isolate those which are core never events. MATERIALS AND METHODS: A systematic review will be conducted using Medline, Medline in Process, Scopus, PsychINFO, Embase via OVID, and CINAHL via EBSCO databases, as well as grey literature. We will include articles of any study design that discuss never events or one of its synonymous terms in the context of medical care. Four independent reviewers will conduct the title and abstract as well as the full-text screening, and 2 reviewers will abstract data. Data will be analyzed using narrative synthesis. Results will be categorized by year and geographic location, and by other factors determined during full-text screening. DISCUSSION AND CONCLUSION: The lack of consensus regarding never events hinders progress in reducing their occurrence. Differing data sources makes comparison challenging, and limits the ability for patient safety groups to work collaboratively and share learnings with others. Identifying a core set of never events will serve as a first step to focus our efforts to reduce these harmful incidents.

Patterns and trends among physicians-in-training named in civil legal cases: a retrospective analysis of Canadian Medical Protective Association data from 1993 to 2017.

BACKGROUND: Medico-legal data show opportunities to improve safe medical care; little is published on the experience of physicians-in-training with medical malpractice. The purpose of this study was to examine closed civil legal cases involving physicians-in-training over time and provide novel insights on case and physicians characteristics. METHODS: We conducted a retrospective descriptive study of closed civil legal cases at the Canadian Medical Protective Association (CMPA), a mutual medico-legal defence organization for more than 105 000 physicians, representing an estimated 95% of physicians in Canada. Eligible cases involved at least 1 physician-in-training and were closed between 1993 and 2017 (for time trends) or 2008 and 2017 (for descriptive analyses). We analyzed case rates over time using Poisson regression and the annualized change rate. Descriptive analyses addressed case duration, medico-legal outcome and patient harm. We explored physician specialties and practice characteristics in a subset of cases. RESULTS: Over a 25-year period (1993-2017), 4921 physicians-in-training were named in 2951 closed civil legal cases, and case rates decreased significantly (ß = -0.04, 95% confidence interval -0.05 to -0.03, where ß was the 1-year difference in log case rates). The annualized change rate was -1.1% per year. Between 2008 and 2017, 1901 (4.1%) of 45 967 physicians-in-training were named in 1107 civil legal cases. Cases with physicians-in-training generally involved more severe patient harm than cases without physicians-in-training. In a subgroup with available information (n = 951), surgical specialties were named most often (n = 531, 55.8%). INTERPRETATION: The rate of civil legal cases involving physicians-in-training has diminished over time, but more recent cases featured severe patient harm and death. Efforts to promote patient safety may enhance medical care and reduce the frequency and severity of malpractice issues for physicians-in-training.

Systematic single-variant and gene-based association testing of thousands of phenotypes in 394,841 UK Biobank exomes.

Genome-wide association studies have successfully discovered thousands of common variants associated with human diseases and traits, but the landscape of rare variations in human disease has not been explored at scale. Exome-sequencing studies of population biobanks provide an opportunity to systematically evaluate the impact of rare coding variations across a wide range of phenotypes to discover genes and allelic series relevant to human health and disease. Here, we present results from systematic association analyses of 4,529 phenotypes using single-variant and gene tests of 394,841 individuals in the UK Biobank with exome-sequence data. We find that the discovery of genetic associations is tightly linked to frequency and is correlated with metrics of deleteriousness and natural selection. We highlight biological findings elucidated by these data and release the dataset as a public resource alongside the Genebass browser for rapidly exploring rare-variant association results.

The burden of rare protein-truncating genetic variants on human lifespan.

Genetic predisposition has been shown to contribute substantially to the age at which we die. Genome-wide association studies (GWASs) have linked more than 20 loci to phenotypes related to human lifespan1. However, little is known about how lifespan is impacted by gene loss of function. Through whole-exome sequencing of 352,338 UK Biobank participants of European ancestry, we assessed the relevance of protein-truncating variant (PTV) gene burden on individual and parental survival. We identified four exome-wide significant (P < 4.2 × 10-7) human lifespan genes, BRCA1, BRCA2, ATM and TET2. Gene and gene-set, PTV-burden, phenome-wide association studies support known roles of these genes in cancer to impact lifespan at the population level. The TET2 PTV burden was associated with a lifespan through somatic mutation events presumably due to clonal hematopoiesis. The overlap between PTV burden and common variant-based lifespan GWASs was modest, underscoring the value of exome sequencing in well-powered biobank cohorts to complement GWASs for identifying genes underlying complex traits.

Pairwise effects between lipid GWAS genes modulate lipid plasma levels and cellular uptake.

Complex traits are characterized by multiple genes and variants acting simultaneously on a phenotype. However, studying the contribution of individual pairs of genes to complex traits has been challenging since human genetics necessitates very large population sizes, while findings from model systems do not always translate to humans. Here, we combine genetics with combinatorial RNAi (coRNAi) to systematically test for pairwise additive effects (AEs) and genetic interactions (GIs) between 30 lipid genome-wide association studies (GWAS) genes. Gene-based burden tests from 240,970 exomes show that in carriers with truncating mutations in both, APOB and either PCSK9 or LPL ("human double knock-outs") plasma lipid levels change additively. Genetics and coRNAi identify overlapping AEs for 12 additional gene pairs. Overlapping GIs are observed for TOMM40/APOE with SORT1 and NCAN. Our study identifies distinct gene pairs that modulate plasma and cellular lipid levels primarily via AEs and nominates putative drug target pairs for improved lipid-lowering combination therapies.

Advancing human genetics research and drug discovery through exome sequencing of the UK Biobank.

The UK Biobank Exome Sequencing Consortium (UKB-ESC) is a private-public partnership between the UK Biobank (UKB) and eight biopharmaceutical companies that will complete the sequencing of exomes for all ~500,000 UKB participants. Here, we describe the early results from ~200,000 UKB participants and the features of this project that enabled its success. The biopharmaceutical industry has increasingly used human genetics to improve success in drug discovery. Recognizing the need for large-scale human genetics data, as well as the unique value of the data access and contribution terms of the UKB, the UKB-ESC was formed. As a result, exome data from 200,643 UKB enrollees are now available. These data include ~10 million exonic variants-a rich resource of rare coding variation that is particularly valuable for drug discovery. The UKB-ESC precompetitive collaboration has further strengthened academic and industry ties and has provided teams with an opportunity to interact with and learn from the wider research community.

The Parkinson's disease-associated gene ITPKB protects against α-synuclein aggregation by regulating ER-to-mitochondria calcium release.

Inositol-1,4,5-triphosphate (IP3) kinase B (ITPKB) is a ubiquitously expressed lipid kinase that inactivates IP3, a secondary messenger that stimulates calcium release from the endoplasmic reticulum (ER). Genome-wide association studies have identified common variants in the ITPKB gene locus associated with reduced risk of sporadic Parkinson's disease (PD). Here, we investigate whether ITPKB activity or expression level impacts PD phenotypes in cellular and animal models. In primary neurons, knockdown or pharmacological inhibition of ITPKB increased levels of phosphorylated, insoluble α-synuclein pathology following treatment with α-synuclein preformed fibrils (PFFs). Conversely, ITPKB overexpression reduced PFF-induced α-synuclein aggregation. We also demonstrate that ITPKB inhibition or knockdown increases intracellular calcium levels in neurons, leading to an accumulation of calcium in mitochondria that increases respiration and inhibits the initiation of autophagy, suggesting that ITPKB regulates α-synuclein pathology by inhibiting ER-to-mitochondria calcium transport. Furthermore, the effects of ITPKB on mitochondrial calcium and respiration were prevented by pretreatment with pharmacological inhibitors of the mitochondrial calcium uniporter complex, which was also sufficient to reduce α-synuclein pathology in PFF-treated neurons. Taken together, these results identify ITPKB as a negative regulator of α-synuclein aggregation and highlight modulation of ER-to-mitochondria calcium flux as a therapeutic strategy for the treatment of sporadic PD.

The medico-legal helpline: A content analysis of postgraduate medical trainee advice calls.

CONTEXT: Available literature exploring medical liability and postgraduate medical education consistently posits that postgraduate trainees worry about their exposure to medico-legal liability. This assumption has formed the basis for research and curriculum development. OBJECTIVES: The aim of this study was to describe the encounters that lead physicians-in-training to seek external medico-legal guidance. We sought to provide empirical evidence on trends and themes related to medico-legal advice requests from physicians-in-training. METHODS: Our primary dataset consisted of records of calls from physicians-in-training to the medico-legal helpline of the Canadian Medical Protective Association (CMPA), a national mutual defence organisation providing medico-legal advice and liability protection for over 95% of Canada's physicians. We conducted a trend analysis of the frequency of calls for advice over 10 years from physician-in-training compared with non-trainee physicians. Furthermore, we performed a content analysis of calls made over the most recent 2 years (2016-2017) to elucidate the concerns that led to trainees seeking medico-legal advice. RESULTS: The 10-year trend analysis revealed that the annual growth in the number of physician-in-training advice calls (8.8%) exceeded other CMPA physician groups and was in excess of trainee population growth over the same period. The content analysis identified four core themes: managing confidential information, complex care situations, academic matters and patient safety incidents. CONCLUSIONS: Our findings indicate that trainees are asking questions about their medico-legal liability with increasing frequency. This study contributes new evidence on the issues that lead to trainees seeking help. We believe that understanding trainees' medico-legal advice requests will support medical educators to tailor quality improvement education to learners' needs.

Estimating the relative frequency of leukodystrophies and recommendations for carrier screening in the era of next-generation sequencing.

Leukodystrophies are a heterogeneous group of heritable disorders characterized by abnormal brain white matter signal on magnetic resonance imaging (MRI) and primary involvement of the cellular components of myelin. Previous estimates suggest the incidence of leukodystrophies as a whole to be 1 in 7,000 individuals, however the frequency of specific diagnoses relative to others has not been described. Next generation sequencing approaches offer the opportunity to redefine our understanding of the relative frequency of different leukodystrophies. We assessed the relative frequency of all 30 leukodystrophies (associated with 55 genes) in more than 49,000 exomes. We identified a relatively high frequency of disorders previously thought of as very rare, including Aicardi Goutières Syndrome, TUBB4A-related leukodystrophy, Peroxisomal biogenesis disorders, POLR3-related Leukodystrophy, Vanishing White Matter, and Pelizaeus-Merzbacher Disease. Despite the relative frequency of these conditions, carrier-screening laboratories regularly test only 20 of the 55 leukodystrophy-related genes, and do not test at all, or test only one or a few, genes for some of the higher frequency disorders. Relative frequency of leukodystrophies previously considered very rare suggests these disorders may benefit from expanded carrier screening.

Exome Sequencing in Individuals with Isolated Biliary Atresia.

Biliary atresia (BA) is a severe pediatric liver disease resulting in necroinflammatory obliteration of the extrahepatic biliary tree. BA presents within the first few months of life as either an isolated finding or with additional syndromic features. The etiology of isolated BA is unknown, with evidence for infectious, environmental, and genetic risk factors described. However, to date, there are no definitive causal genes identified for isolated BA in humans, and the question of whether single gene defects play a major role remains open. We performed exome-sequencing in 101 North American patients of European descent with isolated BA (including 30 parent-child trios) and considered several experimental designs to identify potentially deleterious protein-altering variants that may be involved in the disease. In a case-only analysis, we did not identify genes with variants shared among more than two probands, and burden tests of rare variants using a case-case control design did not yield significant results. In the trio analysis of 30 simplex families (patient and parent trios), we identified 66 de novo variants in 66 genes including potentially deleterious variants in STIP1 and REV1. STIP1 is a co-chaperone for the heat-shock protein, HSP90, and has been shown to have diverse functions in yeast, flies and mammals, including stress-responses. REV1 is known to be a key player in DNA repair pathway and to interact with HSP90. In conclusion, our results do not support the hypothesis that a simple genetic model is responsible for the majority of cases of isolated BA. Our finding of de novo variants in genes linked to evolutionarily conserved stress responses (STIP1 and REV1) suggests that exploration of how genetic susceptibility and environmental exposure may interact to cause BA is warranted.

Characterizing reduced coverage regions through comparison of exome and genome sequencing data across 10 centers.

PURPOSE: As massively parallel sequencing is increasingly being used for clinical decision making, it has become critical to understand parameters that affect sequencing quality and to establish methods for measuring and reporting clinical sequencing standards. In this report, we propose a definition for reduced coverage regions and describe a set of standards for variant calling in clinical sequencing applications. METHODS: To enable sequencing centers to assess the regions of poor sequencing quality in their own data, we optimized and used a tool (ExCID) to identify reduced coverage loci within genes or regions of particular interest. We used this framework to examine sequencing data from 500 patients generated in 10 projects at sequencing centers in the National Human Genome Research Institute/National Cancer Institute Clinical Sequencing Exploratory Research Consortium. RESULTS: This approach identified reduced coverage regions in clinically relevant genes, including known clinically relevant loci that were uniquely missed at individual centers, in multiple centers, and in all centers. CONCLUSION: This report provides a process road map for clinical sequencing centers looking to perform similar analyses on their data.