New Enhancing MRI Lesions Associate with IL-17, Neutrophil Degranulation and Integrin Microparticles: Multi-Omics Combined with Frequent MRI in Multiple Sclerosis

Background: Blood–barrier (BBB) breakdown and active inflammation are hallmarks of relapsing multiple sclerosis (RMS), but the molecular events contributing to the development of new lesions are not well explored. Leaky endothelial junctions are associated with increased production of endothelial-derived extracellular microvesicles (EVs) and result in the entry of circulating immune cells into the brain. MRI with intravenous gadolinium (Gd) can visualize acute blood–barrier disruption as the initial event of the evolution of new lesions. Methods: Here, weekly MRI with Gd was combined with proteomics, multiplex immunoassay, and endothelial stress-optimized EV array to identify early markers related to BBB disruption. Five patients with RMS with no disease-modifying treatment were monitored weekly using high-resolution 3T MRI scanning with intravenous gadolinium (Gd) for 8 weeks. Patients were then divided into three groups (low, medium, or high MRI activity) defined by the number of new, total, and maximally enhancing Gd-enhancing lesions and the number of new FLAIR lesions. Plasma samples taken at each MRI were analyzed for protein biomarkers of inflammation by quantitative proteomics, and cytokines using multiplex immunoassays. EVs were characterized with an optimized endothelial stress EV array based on exosome surface protein markers for the detection of soluble secreted EVs. Results: Proteomics analysis of plasma yielded quantitative information on 208 proteins at each patient time point (n = 40). We observed the highest number of unique dysregulated proteins (DEPs) and the highest functional enrichment in the low vs. high MRI activity comparison. Complement activation and complement/coagulation cascade were also strongly overrepresented in the low vs. high MRI activity comparison. Activation of the alternative complement pathway, pathways of blood coagulation, extracellular matrix organization, and the regulation of TLR and IGF transport were unique for the low vs. high MRI activity comparison as well, with these pathways being overrepresented in the patient with high MRI activity. Principal component analysis indicated the individuality of plasma profiles in patients. IL-17 was upregulated at all time points during 8 weeks in patients with high vs. low MRI activity. Hierarchical clustering of soluble markers in the plasma indicated that all four MRI outcomes clustered together with IL-17, IL-12p70, and IL-1β. MRI outcomes also showed clustering with EV markers CD62E/P, MIC A/B, ICAM-1, and CD42A. The combined cluster of these cytokines, EV markers, and MRI outcomes clustered also with IL-12p40 and IL-7. All four MRI outcomes correlated positively with levels of IL-17 (p < 0.001, respectively), and EV-ICAM-1 (p < 0.0003, respectively). IL-1β levels positively correlated with the number of new Gd-enhancing lesions (p < 0.01), new FLAIR lesions (p < 0.001), and total number of Gd-enhancing lesions (p < 0.05). IL-6 levels positively correlated with the number of new FLAIR lesions (p < 0.05). Random Forests and linear mixed models identified IL-17, CCL17/TARC, CCL3/MIP-1α, and TNF-α as composite biomarkers predicting new lesion evolution. Conclusions: Combination of serial frequent MRI with proteome, neuroinflammation markers, and protein array data of EVs enabled assessment of temporal changes in inflammation and endothelial dysfunction in RMS related to the evolution of new and enhancing lesions. Particularly, the Th17 pathway and IL-1β clustered and correlated with new lesions and Gd enhancement, indicating their importance in BBB disruption and initiating acute brain inflammation in MS. In addition to the Th17 pathway, abundant protein changes between MRI activity groups suggested the role of EVs and the coagulation system along with innate immune responses including acute phase proteins, complement components, and neutrophil degranulation.Background: Blood-barrier (BBB) breakdown and active inflammation are hallmarks of relapsing multiple sclerosis (RMS), but the molecular events contributing to the development of new lesions are not well explored. Leaky endothelial junctions are associated with increased production of endothelial-derived extracellular microvesicles (EVs) and result in the entry of circulating immune cells into the brain. MRI with intravenous gadolinium (Gd) can visualize acute blood-barrier disruption as the initial event of the evolution of new lesions.Methods: Here, weekly MRI with Gd was combined with proteomics, multiplex immunoassay, and endothelial stress-optimized EV array to identify early markers related to BBB disruption. Five patients with RMS with no disease-modifying treatment were monitored weekly using high-resolution 3T MRI scanning with intravenous gadolinium (Gd) for 8 weeks. Patients were then divided into three groups (low, medium, or high MRI activity) defined by the number of new, total, and maximally enhancing Gd-enhancing lesions and the number of new FLAIR lesions. Plasma samples taken at each MRI were analyzed for protein biomarkers of inflammation by quantitative proteomics, and cytokines using multiplex immunoassays. EVs were characterized with an optimized endothelial stress EV array based on exosome surface protein markers for the detection of soluble secreted EVs.Results: Proteomics analysis of plasma yielded quantitative information on 208 proteins at each patient time point (n = 40). We observed the highest number of unique dysregulated proteins (DEPs) and the highest functional enrichment in the low vs. high MRI activity comparison. Complement activation and complement/coagulation cascade were also strongly overrepresented in the low vs. high MRI activity comparison. Activation of the alternative complement pathway, pathways of blood coagulation, extracellular matrix organization, and the regulation of TLR and IGF transport were unique for the low vs. high MRI activity comparison as well, with these pathways being overrepresented in the patient with high MRI activity. Principal component analysis indicated the individuality of plasma profiles in patients. IL-17 was upregulated at all time points during 8 weeks in patients with high vs. low MRI activity. Hierarchical clustering of soluble markers in the plasma indicated that all four MRI outcomes clustered together with IL-17, IL-12p70, and IL-1β. MRI outcomes also showed clustering with EV markers CD62E/P, MIC A/B, ICAM-1, and CD42A. The combined cluster of these cytokines, EV markers, and MRI outcomes clustered also with IL-12p40 and IL-7. All four MRI outcomes correlated positively with levels of IL-17 (p < 0.001, respectively), and EV-ICAM-1 (p < 0.0003, respectively). IL-1β levels positively correlated with the number of new Gd-enhancing lesions (p < 0.01), new FLAIR lesions (p < 0.001), and total number of Gd-enhancing lesions (p < 0.05). IL-6 levels positively correlated with the number of new FLAIR lesions (p < 0.05). Random Forests and linear mixed models identified IL-17, CCL17/TARC, CCL3/MIP-1α, and TNF-α as composite biomarkers predicting new lesion evolution.Conclusions: Combination of serial frequent MRI with proteome, neuroinflammation markers, and protein array data of EVs enabled assessment of temporal changes in inflammation and endothelial dysfunction in RMS related to the evolution of new and enhancing lesions. Particularly, the Th17 pathway and IL-1β clustered and correlated with new lesions and Gd enhancement, indicating their importance in BBB disruption and initiating acute brain inflammation in MS. In addition to the Th17 pathway, abundant protein changes between MRI activity groups suggested the role of EVs and the coagulation system along with innate immune responses including acute phase proteins, complement components, and neutrophil degranulation

Novel multiple sclerosis susceptibility loci implicated in epigenetic regulation

We conducted a genome-wide association study (GWAS) on multiple sclerosis (MS) susceptibility in German cohorts with 4888 cases and 10,395 controls. In addition to associations within the major histocompatibility complex (MHC) region, 15 non-MHC loci reached genome-wide significance. Four of these loci are novel MS susceptibility loci. They map to the genes L3MBTL3, MAZ, ERG, and SHMT1. The lead variant at SHMT1 was replicated in an independent Sardinian cohort. Products of the genes L3MBTL3, MAZ, and ERG play important roles in immune cell regulation. SHMT1 encodes a serine hydroxymethyltransferase catalyzing the transfer of a carbon unit to the folate cycle. This reaction is required for regulation of methylation homeostasis, which is important for establishment and maintenance of epigenetic signatures. Our GWAS approach in a defined population with limited genetic substructure detected associations not found in larger, more heterogeneous cohorts, thus providing new clues regarding MS pathogenesis

Genome-Wide Association Study with Targeted and Non-targeted NMR Metabolomics Identifies 15 Novel Loci of Urinary Human Metabolic Individuality

Genome-wide association studies with metabolic traits (mGWAS) uncovered many genetic variants that influence human metabolism. These genetically influenced metabotypes (GIMs) contribute to our metabolic individuality, our capacity to respond to environmental challenges, and our susceptibility to specific diseases. While metabolic homeostasis in blood is a well investigated topic in large mGWAS with over 150 known loci, metabolic detoxification through urinary excretion has only been addressed by few small mGWAS with only 11 associated loci so far. Here we report the largest mGWAS to date, combining targeted and non-targeted 1H NMR analysis of urine samples from 3,861 participants of the SHIP-0 cohort and 1,691 subjects of the KORA F4 cohort. We identified and replicated 22 loci with significant associations with urinary traits, 15 of which are new (HIBCH, CPS1, AGXT, XYLB, TKT, ETNPPL, SLC6A19, DMGDH, SLC36A2, GLDC, SLC6A13, ACSM3, SLC5A11, PNMT, SLC13A3). Two-thirds of the urinary loci also have a metabolite association in blood. For all but one of the 6 loci where significant associations target the same metabolite in blood and urine, the genetic effects have the same direction in both fluids. In contrast, for the SLC5A11 locus, we found increased levels of myo-inositol in urine whereas mGWAS in blood reported decreased levels for the same genetic variant. This might indicate less effective re-absorption of myo-inositol in the kidneys of carriers. In summary, our study more than doubles the number of known loci that influence urinary phenotypes. It thus allows novel insights into the relationship between blood homeostasis and its regulation through excretion. The newly discovered loci also include variants previously linked to chronic kidney disease (CPS1, SLC6A13), pulmonary hypertension (CPS1), and ischemic stroke (XYLB). By establishing connections from gene to disease via metabolic traits our results provide novel hypotheses about molecular mechanisms involved in the etiology of diseases

DeepWAS: multivariate genotype-phenotype associations by directly integrating regulatory information using deep learning

Genome-wide association studies (GWAS) identify genetic variants associated with traits or diseases. GWAS never directly link variants to regulatory mechanisms. Instead, the functional annotation of variants is typically inferred by post hoc analyses. A specific class of deep learning-based methods allows for the prediction of regulatory effects per variant on several cell type-specific chromatin features. We here describe "DeepWAS", a new approach that integrates these regulatory effect predictions of single variants into a multivariate GWAS setting. Thereby, single variants associated with a trait or disease are directly coupled to their impact on a chromatin feature in a cell type. Up to 61 regulatory SNPs, called dSNPs, were associated with multiple sclerosis (MS, 4,888 cases and 10,395 controls), major depressive disorder (MDD, 1,475 cases and 2,144 controls), and height (5,974 individuals). These variants were mainly non-coding and reached at least nominal significance in classical GWAS. The prediction accuracy was higher for DeepWAS than for classical GWAS models for 91% of the genome-wide significant, MS-specific dSNPs. DSNPs were enriched in public or cohort-matched expression and methylation quantitative trait loci and we demonstrated the potential of DeepWAS to generate testable functional hypotheses based on genotype data alone. DeepWAS is available at https://github.com/cellmapslab/DeepWAS

Gene transcripts associated with muscle strength: a CHARGE meta-analysis of 7,781 persons

This is the author accepted manuscript. The final version is available from the publisher via the DOI in this record.Background: Lower muscle strength in midlife predicts disability and mortality in later life. Bloodborne factors, including growth differentiation factor 11 (GDF11), have been linked to muscle regeneration in animal models. We aimed to identify gene transcripts associated with muscle strength in adults. Methods: Meta-analysis of whole blood gene expression (overall 17,534 unique genes measured by microarray) and hand-grip strength in four independent cohorts (n=7,781, ages: 20-104 years, weighted mean=56), adjusted for age, sex, height, weight, and leukocyte subtypes. Separate analyses were performed in subsets (older/younger than 60, male/female). Results: Expression levels of 221 genes were associated with strength after adjustment for cofactors and for multiple statistical testing, including ALAS2 (rate limiting enzyme in heme synthesis), PRF1 (perforin, a cytotoxic protein associated with inflammation), IGF1R and IGF2BP2 (both insulin like growth factor related). We identified statistical enrichment for hemoglobin biosynthesis, innate immune activation and the stress response. Ten genes were only associated in younger individuals, four in males only and one in females only. For example PIK3R2 (a negative regulator of PI3K/AKT growth pathway) was negatively associated with muscle strength in younger (=60 years). We also show that 115 genes (52%) have not previously been linked to muscle in NCBI PubMed abstracts Conclusions: This first large-scale transcriptome study of muscle strength in human adults confirmed associations with known pathways and provides new evidence for over half of the genes identified. There may be age and sex specific gene expression signatures in blood for muscle strength.Wellcome TrustFHS gene expression profiling was funded through the Division of Intramural Research (Principal Investigator, Daniel Levy), National Heart, Lung, and Blood Institute, National Institutes of Health, Bethesda, MD. Dr. Murabito is supported by NIH grant R01AG029451. Dr. Kiel is supported by NIH R01 AR41398. The Framingham Heart Study is supported by National Heart, Lung, and Blood Institute contract N01-HC-25195.The InCHIANTI study was supported in part by the Intramural Research Program, National Institute on Aging, NIH, Baltimore MD USA. D.M. and L.W.H. were generously supported by a Wellcome Trust Institutional Strategic Support Award (WT097835MF). W.E.H. was funded by the National Institute for Health Research (NIHR) Collaboration for Leadership in Applied Health Research and Care (CLAHRC) for the South West Peninsula. The views expressed in this publication are those of the authors and not necessarily those of the NHS, the NIHR or the Department of Health in EnglandThe infrastructure for the NESDA study (www.nesda.nl) is funded through the Geestkracht program of the Netherlands Organisation for Health Research and Development (Zon-Mw, grant number 10-000-1002) and is supported by participating universities and mental health care organizations (VU University Medical Center, GGZ inGeest, Arkin, Leiden University Medical Center, GGZ Rivierduinen, University Medical Center Groningen, Lentis, GGZ Friesland, GGZ Drenthe, Scientific Institute for Quality of Healthcare (IQ healthcare), Netherlands Institute for Health Services Research (NIVEL) and Netherlands Institute of Mental Health and Addiction (Trimbos Institute).The Rotterdam Study is funded by Erasmus Medical Center and Erasmus University, Rotterdam, Netherlands Organization for the Health Research and Development (ZonMw), the Netherlands Organisation of Scientific Research NWO Investments (nr. 175.010.2005.011, 911-03-012), the Research Institute for Diseases in the Elderly (014-93- 28 015; RIDE2), the Ministry of Education, Culture and Science, the Ministry for Health, Welfare and Sports, the European Commission (DG XII), and the Municipality of Rotterdam. The authors are grateful to the study participants, the staff from the Rotterdam Study and the participating general practitioners and pharmacists. The generation and management of RNA-expression array data for the Rotterdam Study was executed and funded by the Human Genotyping Facility of the Genetic Laboratory of the Department of Internal Medicine, Erasmus MC, the Netherlands. We thank Marjolein Peters, MSc, Ms. Mila Jhamai, Ms. Jeannette M. Vergeer-Drop, Ms. Bernadette van Ast-Copier, Mr. Marijn Verkerk and Jeroen van Rooij, BSc for their help in creating the RNA array expression databaseSHIP is part of the Community Medicine Research net of the University of Greifswald, Germany, which is funded by the Federal Ministry of Education and Research (grants no. 01ZZ9603, 01ZZ0103, and 01ZZ0403), the Ministry of Cultural Affairs as well as the Social Ministry of the Federal State of Mecklenburg-West Pomerania, and the network ‘Greifswald Approach to Individualized Medicine (GANI_MED)’ funded by the Federal Ministry of Education and Research (grant 03IS2061A). The University of Greifswald is a member of the 'Center of Knowledge Interchange' program of the Siemens AG and the Caché Campus program of the InterSystems GmbH

Genome-wide analysis of mitochondrial DNA copy number reveals loci implicated in nucleotide metabolism, platelet activation, and megakaryocyte proliferation

Mitochondrial DNA copy number (mtDNA-CN) measured from blood specimens is a minimally invasive marker of mitochondrial function that exhibits both inter-individual and intercellular variation. To identify genes involved in regulating mitochondrial function, we performed a genome-wide association study (GWAS) in 465,809 White individuals from the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) consortium and the UK Biobank (UKB). We identified 133 SNPs with statistically significant, independent effects associated with mtDNA-CN across 100 loci. A combination of fine-mapping, variant annotation, and co-localization analyses was used to prioritize genes within each of the 133 independent sites. Putative causal genes were enriched for known mitochondrial DNA depletion syndromes (p = 3.09 × 10(–15)) and the gene ontology (GO) terms for mtDNA metabolism (p = 1.43 × 10(–8)) and mtDNA replication (p = 1.2 × 10(–7)). A clustering approach leveraged pleiotropy between mtDNA-CN associated SNPs and 41 mtDNA-CN associated phenotypes to identify functional domains, revealing three distinct groups, including platelet activation, megakaryocyte proliferation, and mtDNA metabolism. Finally, using mitochondrial SNPs, we establish causal relationships between mitochondrial function and a variety of blood cell-related traits, kidney function, liver function and overall (p = 0.044) and non-cancer mortality (p = 6.56 × 10(–4)). SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s00439-021-02394-w

Gene-gene Interaction Analyses for Atrial Fibrillation

Atrial fibrillation (AF) is a heritable disease that affects more than thirty million individuals worldwide. Extensive efforts have been devoted to the study of genetic determinants of AF. The objective of our study is to examine the effect of gene-gene interaction on AF susceptibility. We performed a large-scale association analysis of gene-gene interactions with AF in 8,173 AF cases, and 65,237 AF-free referents collected from 15 studies for discovery. We examined putative interactions between genome-wide SNPs and 17 known AF-related SNPs. The top interactions were then tested for association in a

Meta-analysis of human genome-microbiome association studies: the MiBioGen consortium initiative

Abstract Background In recent years, human microbiota, especially gut microbiota, have emerged as an important yet complex trait influencing human metabolism, immunology, and diseases. Many studies are investigating the forces underlying the observed variation, including the human genetic variants that shape human microbiota. Several preliminary genome-wide association studies (GWAS) have been completed, but more are necessary to achieve a fuller picture. Results Here, we announce the MiBioGen consortium initiative, which has assembled 18 population-level cohorts and some 19,000 participants. Its aim is to generate new knowledge for the rapidly developing field of microbiota research. Each cohort has surveyed the gut microbiome via 16S rRNA sequencing and genotyped their participants with full-genome SNP arrays. We have standardized the analytical pipelines for both the microbiota phenotypes and genotypes, and all the data have been processed using identical approaches. Our analysis of microbiome composition shows that we can reduce the potential artifacts introduced by technical differences in generating microbiota data. We are now in the process of benchmarking the association tests and performing meta-analyses of genome-wide associations. All pipeline and summary statistics results will be shared using public data repositories. Conclusion We present the largest consortium to date devoted to microbiota-GWAS. We have adapted our analytical pipelines to suit multi-cohort analyses and expect to gain insight into host-microbiota cross-talk at the genome-wide level. And, as an open consortium, we invite more cohorts to join us (by contacting one of the corresponding authors) and to follow the analytical pipeline we have developed

The Polygenic and Monogenic Basis of Blood Traits and Diseases

Blood cells play essential roles in human health, underpinning physiological processes such as immunity, oxygen transport, and clotting, which when perturbed cause a significant global health burden. Here we integrate data from UK Biobank and a large-scale international collaborative effort, including data for 563,085 European ancestry participants, and discover 5,106 new genetic variants independently associated with 29 blood cell phenotypes covering a range of variation impacting hematopoiesis. We holistically characterize the genetic architecture of hematopoiesis, assess the relevance of the omnigenic model to blood cell phenotypes, delineate relevant hematopoietic cell states influenced by regulatory genetic variants and gene networks, identify novel splice-altering variants mediating the associations, and assess the polygenic prediction potential for blood traits and clinical disorders at the interface of complex and Mendelian genetics. These results show the power of large-scale blood cell trait GWAS to interrogate clinically meaningful variants across a wide allelic spectrum of human variation. Analysis of blood cell traits in the UK Biobank and other cohorts illuminates the full genetic architecture of hematopoietic phenotypes, with evidence supporting the omnigenic model for complex traits and linking polygenic burden with monogenic blood diseases