Discovering cellular programs of intrinsic and extrinsic drivers of metabolic traits using LipocyteProfiler.

A primary obstacle in translating genetic associations with disease into therapeutic strategies is elucidating the cellular programs affected by genetic risk variants and effector genes. Here, we introduce LipocyteProfiler, a cardiometabolic-disease-oriented high-content image-based profiling tool that enables evaluation of thousands of morphological and cellular profiles that can be systematically linked to genes and genetic variants relevant to cardiometabolic disease. We show that LipocyteProfiler allows surveillance of diverse cellular programs by generating rich context- and process-specific cellular profiles across hepatocyte and adipocyte cell-state transitions. We use LipocyteProfiler to identify known and novel cellular mechanisms altered by polygenic risk of metabolic disease, including insulin resistance, fat distribution, and the polygenic contribution to lipodystrophy. LipocyteProfiler paves the way for large-scale forward and reverse deep phenotypic profiling in lipocytes and provides a framework for the unbiased identification of causal relationships between genetic variants and cellular programs relevant to human disease.

Genome-wide association analysis and Mendelian randomization proteomics identify drug targets for heart failure.

We conduct a large-scale meta-analysis of heart failure genome-wide association studies (GWAS) consisting of over 90,000 heart failure cases and more than 1 million control individuals of European ancestry to uncover novel genetic determinants for heart failure. Using the GWAS results and blood protein quantitative loci, we perform Mendelian randomization and colocalization analyses on human proteins to provide putative causal evidence for the role of druggable proteins in the genesis of heart failure. We identify 39 genome-wide significant heart failure risk variants, of which 18 are previously unreported. Using a combination of Mendelian randomization proteomics and genetic cis-only colocalization analyses, we identify 10 additional putatively causal genes for heart failure. Findings from GWAS and Mendelian randomization-proteomics identify seven (CAMK2D, PRKD1, PRKD3, MAPK3, TNFSF12, APOC3 and NAE1) proteins as potential targets for interventions to be used in primary prevention of heart failure.

Inherited basis of visceral, abdominal subcutaneous and gluteofemoral fat depots.

For any given level of overall adiposity, individuals vary considerably in fat distribution. The inherited basis of fat distribution in the general population is not fully understood. Here, we study up to 38,965 UK Biobank participants with MRI-derived visceral (VAT), abdominal subcutaneous (ASAT), and gluteofemoral (GFAT) adipose tissue volumes. Because these fat depot volumes are highly correlated with BMI, we additionally study six local adiposity traits: VAT adjusted for BMI and height (VATadj), ASATadj, GFATadj, VAT/ASAT, VAT/GFAT, and ASAT/GFAT. We identify 250 independent common variants (39 newly-identified) associated with at least one trait, with many associations more pronounced in female participants. Rare variant association studies extend prior evidence for PDE3B as an important modulator of fat distribution. Local adiposity traits (1) highlight depot-specific genetic architecture and (2) enable construction of depot-specific polygenic scores that have divergent associations with type 2 diabetes and coronary artery disease. These results - using MRI-derived, BMI-independent measures of local adiposity - confirm fat distribution as a highly heritable trait with important implications for cardiometabolic health outcomes.

A single-cell atlas of human and mouse white adipose tissue.

White adipose tissue, once regarded as morphologically and functionally bland, is now recognized to be dynamic, plastic and heterogenous, and is involved in a wide array of biological processes including energy homeostasis, glucose and lipid handling, blood pressure control and host defence1. High-fat feeding and other metabolic stressors cause marked changes in adipose morphology, physiology and cellular composition1, and alterations in adiposity are associated with insulin resistance, dyslipidemia and type 2 diabetes2. Here we provide detailed cellular atlases of human and mouse subcutaneous and visceral white fat at single-cell resolution across a range of body weight. We identify subpopulations of adipocytes, adipose stem and progenitor cells, vascular and immune cells and demonstrate commonalities and differences across species and dietary conditions. We link specific cell types to increased risk of metabolic disease and provide an initial blueprint for a comprehensive set of interactions between individual cell types in the adipose niche in leanness and obesity. These data comprise an extensive resource for the exploration of genes, traits and cell types in the function of white adipose tissue across species, depots and nutritional conditions.

Glycosylation and Cardiovascular Diseases.

Cardiovascular disease (CVD) is the leading cause of death worldwide, accounting for approximately 18 million deaths in 2017. Coronary artery disease is the predominant cause of death from CVD, followed by stroke. Owing to recent technological advancements, glycans and glycosylation patterns of proteins have been investigated in association with CVD risk factors and clinical events. These studies have found significant associations of glycans as biomarkers of systemic inflammation and major CVD risk factors and events. While more limited, studies have also shown that glycans may be useful for monitoring response to anti-inflammatory therapies and may be responsive to changes in lifestyle, particularly in patients with chronic inflammatory diseases. Glycans capture summative risk information related to inflammatory, immune, and signaling pathways and are promising biomarkers for CVD risk prediction and therapeutic monitoring.

SARS2 simplified scores to estimate risk of hospitalization and death among patients with COVID-19.

Although models have been developed for predicting severity of COVID-19 from the medical history of patients, simplified models with good accuracy could be more practical. In this study, we examined utility of simpler models for estimating risk of hospitalization of patients with COVID-19 and mortality of these patients based on demographic characteristics (sex, age, race, median household income based on zip code) and smoking status of 12,347 patients who tested positive at Mass General Brigham centers. The corresponding electronic records were queried (02/26-07/14/2020) to construct derivation and validation cohorts. The derivation cohort was used to fit generalized linear models for estimating risk of hospitalization within 30 days of COVID-19 diagnosis and mortality within approximately 3 months for the hospitalized patients. In the validation cohort, the model resulted in c-statistics of 0.77 [95% CI 0.73-0.80] for hospitalization, and 0.84 [95% CI 0.74-0.94] for mortality among hospitalized patients. Higher risk was associated with older age, male sex, Black ethnicity, lower socioeconomic status, and current/past smoking status. The models can be applied to predict the absolute risks of hospitalization and mortality, and could aid in individualizing the decision making when detailed medical history of patients is not readily available.

NUScon: a community-driven platform for quantitative evaluation of nonuniform sampling in NMR.

Although the concepts of nonuniform sampling (NUS​​​​​​​) and non-Fourier spectral reconstruction in multidimensional NMR began to emerge 4 decades ago , it is only relatively recently that NUS has become more commonplace. Advantages of NUS include the ability to tailor experiments to reduce data collection time and to improve spectral quality, whether through detection of closely spaced peaks (i.e., "resolution") or peaks of weak intensity (i.e., "sensitivity"). Wider adoption of these methods is the result of improvements in computational performance, a growing abundance and flexibility of software, support from NMR spectrometer vendors, and the increased data sampling demands imposed by higher magnetic fields. However, the identification of best practices still remains a significant and unmet challenge. Unlike the discrete Fourier transform, non-Fourier methods used to reconstruct spectra from NUS data are nonlinear, depend on the complexity and nature of the signals, and lack quantitative or formal theory describing their performance. Seemingly subtle algorithmic differences may lead to significant variabilities in spectral qualities and artifacts. A community-based critical assessment of NUS challenge problems has been initiated, called the "Nonuniform Sampling Contest" (NUScon), with the objective of determining best practices for processing and analyzing NUS experiments. We address this objective by constructing challenges from NMR experiments that we inject with synthetic signals, and we process these challenges using workflows submitted by the community. In the initial rounds of NUScon our aim is to establish objective criteria for evaluating the quality of spectral reconstructions. We present here a software package for performing the quantitative analyses, and we present the results from the first two rounds of NUScon. We discuss the challenges that remain and present a roadmap for continued community-driven development with the ultimate aim of providing best practices in this rapidly evolving field. The NUScon software package and all data from evaluating the challenge problems are hosted on the NMRbox platform.

Historical Controls in Randomized Clinical Trials: Opportunities and Challenges.

Randomized control trials (RCTs) with placebo are the gold standard for determining efficacy of novel pharmaceutical treatments. Since their inception, over 75 years ago, researchers have amassed a large body of underutilized data on outcomes in the placebo control arms of these trials. Although rare disease indications have used these historical placebo data as synthetic controls to reduce burden on patients and accelerate drug discovery, broad use of historical controls is in its infancy. Large-scale historical placebo data could be leveraged to benefit both drug developers and patients if warehoused and made more available to guide trial design and analysis. Here, we examine challenges in utilizing historical controls related to heterogeneity in trial design, outcome ascertainment, patient characteristics, and unmeasured pharmacogenomic effects. We then discuss the advantages and disadvantages of current approaches and propose a path forward to broader use of historical controls in RCTs.

Quantum approximate Bayesian computation for NMR model inference.

Recent technological advances may lead to the development of small scale quantum computers capable of solving problems that cannot be tackled with classical computers. A limited number of algorithms has been proposed and their relevance to real world problems is a subject of active investigation. Analysis of many-body quantum system is particularly challenging for classical computers due to the exponential scaling of Hilbert space dimension with the number of particles. Hence, solving problems relevant to chemistry and condensed matter physics are expected to be the first successful applications of quantum computers. In this paper, we propose another class of problems from the quantum realm that can be solved efficiently on quantum computers: model inference for nuclear magnetic resonance (NMR) spectroscopy, which is important for biological and medical research. Our results are based on three interconnected studies. Firstly, we use methods from classical machine learning to analyze a dataset of NMR spectra of small molecules. We perform a stochastic neighborhood embedding and identify clusters of spectra, and demonstrate that these clusters are correlated with the covalent structure of the molecules. Secondly, we propose a simple and efficient method, aided by a quantum simulator, to extract the NMR spectrum of any hypothetical molecule described by a parametric Heisenberg model. Thirdly, we propose a simple variational Bayesian inference procedure for extracting Hamiltonian parameters of experimentally relevant NMR spectra.

One-Year Effects of Omega-3 Treatment on Fatty Acids, Oxylipins, and Related Bioactive Lipids and Their Associations with Clinical Lipid and Inflammatory Biomarkers: Findings from a Substudy of the Vitamin D and Omega-3 Trial (VITAL).

Omega-3 (n-3) treatment may lower cardiovascular risk, yet its effects on the circulating lipidome and relation to cardiovascular risk biomarkers are unclear. We hypothesized that n-3 treatment is associated with favorable changes in downstream fatty acids (FAs), oxylipins, bioactive lipids, clinical lipid and inflammatory biomarkers. We examined these VITAL200, a nested substudy of 200 subjects balanced on demographics and treatment and randomly selected from the Vitamin D and Omega-3 Trial (VITAL). VITAL is a randomized double-blind trial of 840 mg/d eicosapentaenoic acid (EPA) + docosahexaenoic acid (DHA) vs. placebo among 25,871 individuals. Small polar bioactive lipid features, oxylipins and FAs from plasma and red blood cells were measured using three independent assaying techniques at baseline and one year. The Women's Health Study (WHS) was used for replication with dietary n-3 intake. Randomized n-3 treatment led to changes in 143 FAs, oxylipins and bioactive lipids (False Discovery Rate (FDR) < 0.05 in VITAL200, validated (p-values < 0.05)) in WHS with increases in 95 including EPA, DHA, n-3 docosapentaenoic acid (DPA-n3), and decreases in 48 including DPA-n6, dihomo gamma linolenic (DGLA), adrenic and arachidonic acids. N-3 related changes in the bioactive lipidome were heterogeneously associated with changes in clinical lipid and inflammatory biomarkers. N-3 treatment significantly modulates the bioactive lipidome, which may contribute to its clinical benefits.

SARS2 simplified scores to estimate risk of hospitalization and death among patients with COVID-19.

Although models have been developed for predicting severity of COVID-19 based on the medical history of patients, simplified risk prediction models with good accuracy could be more practical. In this study, we examined utility of simpler models for estimating risk of hospitalization of patients with COVID-19 and mortality of these patients based on demographic characteristics (sex, age, race, median household income based on zip code) and smoking status of 12,347 patients who tested positive at Mass General Brigham centers. The corresponding electronic health records were queried from 02/26/2020 to 07/14/2020 to construct derivation and validation cohorts. The derivation cohort was used to fit a generalized linear model for estimating risk of hospitalization within 30 days of COVID-19 diagnosis and mortality within approximately 3 months for the hospitalized patients. On the validation cohort, the model resulted in c-statistics of 0.77 [95% CI: 0.73-0.80] for hospitalization outcome, and 0.72 [95% CI: 0.69-0.74] for mortality among hospitalized patients. Higher risk was associated with older age, male sex, black ethnicity, lower socioeconomic status, and current/past smoking status. The model can be applied to predict risk of hospitalization and mortality, and could aid decision making when detailed medical history of patients is not easily available.

Probabilistic identification of saccharide moieties in biomolecules and their protein complexes.

The chemical composition of saccharide complexes underlies their biomedical activities as biomarkers for cardiometabolic disease, various types of cancer, and other conditions. However, because these molecules may undergo major structural modifications, distinguishing between compounds of saccharide and non-saccharide origin becomes a challenging computational problem that hinders the aggregation of information about their bioactive moieties. We have developed an algorithm and software package called "Cheminformatics Tool for Probabilistic Identification of Carbohydrates" (CTPIC) that analyzes the covalent structure of a compound to yield a probabilistic measure for distinguishing saccharides and saccharide-derivatives from non-saccharides. CTPIC analysis of the RCSB Ligand Expo (database of small molecules found to bind proteins in the Protein Data Bank) led to a substantial increase in the number of ligands characterized as saccharides. CTPIC analysis of Protein Data Bank identified 7.7% of the proteins as saccharide-binding. CTPIC is freely available as a webservice at (http://ctpic.nmrfam.wisc.edu).

Fragment screening targeting Ebola virus nucleoprotein C-terminal domain identifies lead candidates.

The Ebola Virus is a causative agent of viral hemorrhagic fever outbreaks and a potential global health risk. The outbreak in West Africa (2013-2016) led to 11,000+ deaths and 30,000+ Ebola infected individuals. The current outbreak in the Democratic Republic of Congo (DRC) with 3000+ confirmed cases and 2000+ deaths attributed to Ebola virus infections provides a reminder that innovative countermeasures are still needed. Ebola virus encodes 7 open reading frames (ORFs). Of these, the nucleocapsid protein (eNP) encoded by the first ORF plays many significant roles, including a role in viral RNA synthesis. Here we describe efforts to target the C-terminal domain of eNP (eNP-CTD) that contains highly conserved residues 641-739 as a pan-Ebola antiviral target. Interactions of eNP-CTD with Ebola Viral Protein 30 (eVP30) and Viral Protein 40 (eVP40) have been shown to be crucial for viral RNA synthesis, virion formation, and virion transport. We used nuclear magnetic response (NMR)-based methods to screened the eNP-CTD against a fragment library. Perturbations of 1D 1H NMR spectra identified of 48 of the 439 compounds screened as potential eNP CTD interactors. Subsequent analysis of these compounds to measure chemical shift perturbations in 2D 1H,15N NMR spectra of 15N-labeled protein identified six with low millimolar affinities. All six perturbed an area consisting mainly of residues at or near the extreme C-terminus that we named "Site 1" while three other sites were perturbed by other compounds. Our findings here demonstrate the potential utility of eNP as a target, several fragment hits, and provide an experimental pipeline to validate viral-viral interactions as potential panfiloviral inhibitor targets.

Tools for Enhanced NMR-Based Metabolomics Analysis.

Metabolomics is the study of profiles of small molecules in biological fluids, cells, or organs. These profiles can be thought of as the "fingerprints" left behind from chemical processes occurring in biological systems. Because of its potential for groundbreaking applications in disease diagnostics, biomarker discovery, and systems biology, metabolomics has emerged as a rapidly growing area of research. Metabolomics investigations often, but not always, involve the identification and quantification of endogenous and exogenous metabolites in biological samples. Software tools and databases play a crucial role in advancing the rigor, robustness, reproducibility, and validation of these studies. Specifically, the establishment of a robust library of spectral signatures with unique compound descriptors and atom identities plays a key role in profiling studies based on data from nuclear magnetic resonance (NMR) spectroscopy. Here, we discuss developments leading to a rigorous basis for unique identification of compounds, reproducible numbering of atoms, the compact representation of NMR spectra of metabolites and small molecules, tools for improved compound identification, quantification and visualization, and approaches toward the goal of rigorous analysis of metabolomics data.

I-PINE web server: an integrative probabilistic NMR assignment system for proteins.

Various methods for understanding the structural and dynamic properties of proteins rely on the analysis of their NMR chemical shifts. These methods require the initial assignment of NMR signals to particular atoms in the sequence of the protein, a step that can be very time-consuming. The probabilistic interaction network of evidence (PINE) algorithm for automated assignment of backbone and side chain chemical shifts utilizes a Bayesian probabilistic network model that analyzes sequence data and peak lists from multiple NMR experiments. PINE, which is one of the most popular and reliable automated chemical shift assignment algorithms, has been available to the protein NMR community for longer than a decade. We announce here a new web server version of PINE, called Integrative PINE (I-PINE), which supports more types of NMR experiments than PINE (including three-dimensional nuclear Overhauser enhancement and four-dimensional J-coupling experiments) along with more comprehensive visualization of chemical shift based analysis of protein structure and dynamics. The I-PINE server is freely accessible at http://i-pine.nmrfam.wisc.edu . Help pages and tutorial including browser capability are available at: http://i-pine.nmrfam.wisc.edu/instruction.html . Sample data that can be used for testing the web server are available at: http://i-pine.nmrfam.wisc.edu/examples.html .

Automated evaluation of consistency within the PubChem Compound database.

Identification of discrepant data in aggregated databases is a key step in data curation and remediation. We have applied the ALATIS approach, which is based on the international chemical shift identifier (InChI) model, to the full PubChem Compound database to generate unique and reproducible compound and atom identifiers for all entries for which three-dimensional structures were available. This exercise also served to identify entries with discrepancies between structures and chemical formulas or InChI strings. The use of unique compound identifiers and atom nomenclature should support more rigorous links between small-molecule databases including those containing atom-specific information of the type available from crystallography and spectroscopy. The comprehensive results from this analysis are publicly available through our webserver [http://alatis.nmrfam.wisc.edu/].

NMR-STAR: comprehensive ontology for representing, archiving and exchanging data from nuclear magnetic resonance spectroscopic experiments.

The growth of the biological nuclear magnetic resonance (NMR) field and the development of new experimental technology have mandated the revision and enlargement of the NMR-STAR ontology used to represent experiments, spectral and derived data, and supporting metadata. We present here a brief description of the NMR-STAR ontology and software tools for manipulating NMR-STAR data files, editing the files, extracting selected data, and creating data visualizations. Detailed information on these is accessible from the links provided.

Applications of Parametrized NMR Spin Systems of Small Molecules.

We have developed technology for producing accurate spectral fingerprints of small molecules through modeling of NMR spin system matrices to encapsulate their chemical shifts and scalar couplings. We describe here how libraries of these spin systems utilizing unique and reproducible atom numbering can be used to improve NMR-based ligand screening and metabolomics studies. We introduce new Web services that facilitate the analysis of NMR spectra of mixtures of small molecules to yield their identification and quantification. The library of parametrized compounds has been expanded to cover simulations of 1H NMR spectra at a variety of magnetic fields of more than 1100 compounds, included are many common metabolites and a library of drug-like molecular fragments used in ligand screening. The compound library and related Web services are freely available from http://gissmo.nmrfam.wisc.edu/ .

Architectural Features of Human Mitochondrial Cysteine Desulfurase Complexes from Crosslinking Mass Spectrometry and Small-Angle X-Ray Scattering.

Cysteine desulfurase plays a central role in mitochondrial iron-sulfur cluster biogenesis by generating sulfur through the conversion of L-cysteine to L-alanine and by serving as the platform for assembling other components of the biosynthetic machinery, including ISCU, frataxin, and ferredoxin. The human mitochondrial cysteine desulfurase complex consists of two copies each of NFS1, ISD11, and acyl carrier protein. We describe results from chemical crosslinking coupled with tandem mass spectrometry and small-angle X-ray scattering studies that are consistent with a closed NFS1 dimer rather than an open one for both the cysteine desulfurase-ISCU and cysteine desulfurase-ISCU-frataxin complexes. We present a structural model for the cysteine desulfurase-ISCU-frataxin complex derived from chemical crosslinking restraints in conjunction with the recent crystal structure of the cysteine desulfurase-ISCU-zinc complex and distance constraints from nuclear magnetic resonance.