Genetic associations with ratios between protein levels detect new pQTLs and reveal protein-protein interactions

Protein quantitative trait loci (pQTLs) are an invaluable source of information for drug target development because they provide genetic evidence to support protein function, suggest relationships between cis- and trans-associated proteins, and link proteins to disease endpoints. Using Olink proteomics data for 1,463 proteins measured in over 54,000 samples of the UK Biobank, we identified 4,248 associations with 2,821 ratios between protein levels (rQTLs). rQTLs were 7.6-fold enriched in known protein-protein interactions, suggesting that their ratios reflect biological links between the implicated proteins. Conducting a GWAS on ratios increased the number of discovered genetic signals by 24.7%. The approach can identify novel loci of clinical relevance, support causal gene identification, and reveal complex networks of interacting proteins. Taken together, our study adds significant value to the genetic insights that can be derived from the UKB proteomics data and motivates the wider use of ratios in large-scale GWAS.

Omentin associates with serum metabolite profiles indicating lower diabetes risk: KORA F4 Study

INTRODUCTION

Circulating omentin levels have been positively associated with insulin sensitivity. Although a role for adiponectin in this relationship has been suggested, underlying mechanisms remain elusive. In order to reveal the relationship between omentin and systemic metabolism, this study aimed to investigate associations of serum concentrations of omentin and metabolites.

RESEARCH DESIGN AND METHODS

This study is based on 1124 participants aged 61-82 years from the population-based KORA (Cooperative Health Research in the Region of Augsburg) F4 Study, for whom both serum omentin levels and metabolite concentration profiles were available. Associations were assessed with five multivariable regression models, which were stepwise adjusted for multiple potential confounders, including age, sex, body mass index, waist-to-hip ratio, lifestyle markers (physical activity, smoking behavior and alcohol consumption), serum adiponectin levels, high-density lipoprotein cholesterol, use of lipid-lowering or anti-inflammatory medication, history of myocardial infarction and stroke, homeostasis model assessment 2 of insulin resistance, diabetes status, and use of oral glucose-lowering medication and insulin.

RESULTS

Omentin levels significantly associated with multiple metabolites including amino acids, acylcarnitines, and lipids (eg, sphingomyelins and phosphatidylcholines (PCs)). Positive associations for several PCs, such as diacyl (PC aa C32:1) and alkyl-alkyl (PC ae C32:2), were significant in models 1-4, whereas those with hydroxytetradecenoylcarnitine (C14:1-OH) were significant in all five models. Omentin concentrations were negatively associated with several metabolite ratios, such as the valine-to-PC ae C32:2 and the serine-to-PC ae C32:2 ratios in most models.

CONCLUSIONS

Our results suggest that omentin may influence insulin sensitivity and diabetes risk by changing systemic lipid metabolism, but further mechanistic studies investigating effects of omentin on metabolism of insulin-sensitive tissues are needed.

Nanoparticle enrichment mass-spectrometry proteomics identifies protein-altering variants for precise pQTL mapping

Proteogenomics studies generate hypotheses on protein function and provide genetic evidence for drug target prioritization. Most previous work has been conducted using affinity-based proteomics approaches. These technologies face challenges, such as uncertainty regarding target identity, non-specific binding, and handling of variants that affect epitope affinity binding. Mass spectrometry-based proteomics can overcome some of these challenges. Here we report a pQTL study using the Proteograph™ Product Suite workflow (Seer, Inc.) where we quantify over 18,000 unique peptides from nearly 3000 proteins in more than 320 blood samples from a multi-ethnic cohort in a bottom-up, peptide-centric, mass spectrometry-based proteomics approach. We identify 184 protein-altering variants in 137 genes that are significantly associated with their corresponding variant peptides, confirming target specificity of co-associated affinity binders, identifying putatively causal cis-encoded proteins and providing experimental evidence for their presence in blood, including proteins that may be inaccessible to affinity-based proteomics.

Individual bioenergetic capacity as a potential source of resilience to Alzheimer's disease

Impaired glucose uptake in the brain is one of the earliest presymptomatic manifestations of Alzheimer's disease (AD). The absence of symptoms for extended periods of time suggests that compensatory metabolic mechanisms can provide resilience. Here, we introduce the concept of a systemic 'bioenergetic capacity' as the innate ability to maintain energy homeostasis under pathological conditions, potentially serving as such a compensatory mechanism. We argue that fasting blood acylcarnitine profiles provide an approximate peripheral measure for this capacity that mirrors bioenergetic dysregulation in the brain. Using unsupervised subgroup identification, we show that fasting serum acylcarnitine profiles of participants from the AD Neuroimaging Initiative yields bioenergetically distinct subgroups with significant differences in AD biomarker profiles and cognitive function. To assess the potential clinical relevance of this finding, we examined factors that may offer diagnostic and therapeutic opportunities. First, we identified a genotype affecting the bioenergetic capacity which was linked to succinylcarnitine metabolism and significantly modulated the rate of future cognitive decline. Second, a potentially modifiable influence of beta-oxidation efficiency seemed to decelerate bioenergetic aging and disease progression. Our findings, which are supported by data from more than 9,000 individuals, suggest that interventions tailored to enhance energetic health and to slow bioenergetic aging could mitigate the risk of symptomatic AD, especially in individuals with specific mitochondrial genotypes.

Urinary Cell Gene Signature of Acute Rejection in Kidney Allografts

Introduction

A kidney allograft biopsy may display acute T cell-mediated rejection (TCMR), antibody-mediated rejection (ABMR), or concurrent TCMR + ABMR (MR). Development of noninvasive biomarkers diagnostic of all three types of acute rejection is a useful addition to the diagnostic armamentarium.

Methods

We developed customized RT-qPCR assays and measured urinary cell mRNA copy number in 145 biopsy-matched urine samples from 126 kidney allograft recipients and calculated urinary cell three-gene signature score from log 10 -transformed values for the 18S-normalized CD3E mRNA, 18S-normalized CXCL10 mRNA and 18S rRNA. We determined whether the signature score in biopsy-matched urine specimens discriminates biopsies without rejection (NR, n=50) from biopsies displaying TCMR (n=47), ABMR (n=28) or MR (n=20).

Results

Urinary cell three-gene signature discriminated TCMR, ABMR or MR biopsies from NR biopsies (P <0.0001, One-way ANOVA). Dunnett's multiple comparisons test yielded P<0.0001 for NR vs. TCMR; P <0.001 for NR vs. ABMR; and P <0.0001 for NR vs. MR. By bootstrap resampling, optimism-corrected area under the receiver operating characteristic curve (AUC) was 0.749 (bias-corrected 95% confidence interval [CI], 0.638 to 0.840) for NR vs. TCMR (P<0.0001); 0.780 (95% CI, 0.656 to 0.878) for NR vs. ABMR (P<0.0001); and 0.857 (95% CI, 0.727 to 0.947) for NR vs. MR (P<0.0001). All three rejection biopsy categories were distinguished from NR biopsies with similar accuracy (all AUC comparisons P>0.05).

Conclusion

Urinary cell three-gene signature score may serve as a universal diagnostic signature of acute rejection due to TCMR, ABMR or MR in human kidney allografts with similar performance characteristics.

Genome-wide association of dry (Tamar) date palm fruit color

Date palm (Phoenix dactylifera) fruit (dates) are an economically and culturally significant crop in the Middle East and North Africa. There are hundreds of different commercial cultivars producing dates with distinctive shapes, colors, and sizes. Genetic studies of some date palm traits have been performed, including sex determination, sugar content, and fresh fruit color. In this study, we used genome sequences and image data of 199 dry dates (Tamar) collected from 14 countries to identify genetic loci associated with the color of this fruit stage. Here, we find loci across multiple linkage groups (LG) associated with dry fruit color phenotype. We recover both the previously identified VIRESCENS (VIR) genotype associated with fresh fruit yellow or red color and new associations with the lightness and darkness of dry fruit. This study will add resolution to our understanding of date color phenotype, especially at the most commercially important Tamar stage.

Using trials of caloric restriction and bariatric surgery to explore the effects of body mass index on the circulating proteome

Thousands of proteins circulate in the bloodstream; identifying those which associate with weight and intervention-induced weight loss may help explain mechanisms of diseases associated with adiposity. We aimed to identify consistent protein signatures of weight loss across independent studies capturing changes in body mass index (BMI). We analysed proteomic data from studies implementing caloric restriction (Diabetes Remission Clinical trial) and bariatric surgery (By-Band-Sleeve), using SomaLogic and Olink Explore1536 technologies, respectively. Linear mixed models were used to estimate the effect of the interventions on circulating proteins. Twenty-three proteins were altered in a consistent direction after both bariatric surgery and caloric restriction, suggesting that these proteins are modulated by weight change, independent of intervention type. We also integrated Mendelian randomisation (MR) estimates of the effect of BMI on proteins measured by SomaLogic from a UK blood donor cohort as a third line of causal evidence. These MR estimates provided further corroborative evidence for a role of BMI in regulating the levels of six proteins including alcohol dehydrogenase-4, nogo receptor and interleukin-1 receptor antagonist protein. These results indicate the importance of triangulation in interrogating causal relationships; further study into the role of proteins modulated by weight in disease is now warranted.

Retinal nerve fibre layer thinning and corneal nerve loss in patients with Bardet-Biedl syndrome

BACKGROUND

Bardet-Biedl syndrome (BBS) is an autosomal recessive, genetically heterogeneous, pleiotropic disorder caused by variants in genes involved in the function of the primary cilium. We have harnessed genomics to identify BBS and ophthalmic technologies to describe novel features of BBS.

CASE PRESENTATION

A patient with an unclear diagnosis of syndromic type 2 diabetes mellitus, another affected sibling and unaffected siblings and parents were sequenced using DNA extracted from saliva samples. Corneal confocal microscopy (CCM) and retinal spectral domain optical coherence tomography (SD-OCT) were used to identify novel ophthalmic features in these patients. The two affected individuals had a homozygous variant in C8orf37 (p.Trp185*). SD-OCT and CCM demonstrated a marked and patchy reduction in the retinal nerve fiber layer thickness and loss of corneal nerve fibers, respectively.

CONCLUSION

This report highlights the use of ophthalmic imaging to identify novel retinal and corneal abnormalities that extend the phenotype of BBS in a patient with syndromic type 2 diabetes.

Metabolomic biomarkers of habitual B vitamin intakes unveil novel differentially methylated positions in the human epigenome

BACKGROUND

B vitamins such as folate (B9), B6, and B12 are key in one carbon metabolism, which generates methyl donors for DNA methylation. Several studies have linked differential methylation to self-reported intakes of folate and B12, but these estimates can be imprecise, while metabolomic biomarkers can offer an objective assessment of dietary intakes. We explored blood metabolomic biomarkers of folate and vitamins B6 and B12, to carry out epigenome-wide analyses across up to three European cohorts. Associations between self-reported habitual daily B vitamin intakes and 756 metabolites (Metabolon Inc.) were assessed in serum samples from 1064 UK participants from the TwinsUK cohort. The identified B vitamin metabolomic biomarkers were then used in epigenome-wide association tests with fasting blood DNA methylation levels at 430,768 sites from the Infinium HumanMethylation450 BeadChip in blood samples from 2182 European participants from the TwinsUK and KORA cohorts. Candidate signals were explored for metabolite associations with gene expression levels in a subset of the TwinsUK sample (n = 297). Metabolomic biomarker epigenetic associations were also compared with epigenetic associations of self-reported habitual B vitamin intakes in samples from 2294 European participants.

RESULTS

Eighteen metabolites were associated with B vitamin intakes after correction for multiple testing (Bonferroni-adj. p < 0.05), of which 7 metabolites were available in both cohorts and tested for epigenome-wide association. Three metabolites - pipecolate (metabolomic biomarker of B6 and folate intakes), pyridoxate (marker of B6 and folate) and docosahexaenoate (DHA, marker of B6) - were associated with 10, 3 and 1 differentially methylated positions (DMPs), respectively. The strongest association was observed between DHA and DMP cg03440556 in the SCD gene (effect = 0.093 ± 0.016, p = 4.07E-09). Pyridoxate, a catabolic product of vitamin B6, was inversely associated with CpG methylation near the SLC1A5 gene promoter region (cg02711608 and cg22304262) and with SLC7A11 (cg06690548), but not with corresponding changes in gene expression levels. The self-reported intake of folate and vitamin B6 had consistent but non-significant associations with the epigenetic signals.

CONCLUSION

Metabolomic biomarkers are a valuable approach to investigate the effects of dietary B vitamin intake on the human epigenome.

Progesterone induces meiosis through two obligate co-receptors with PLA2 activity

The steroid hormone progesterone (P4) regulates multiple aspects of reproductive and metabolic physiology. Classical P4 signaling operates through nuclear receptors that regulate transcription. In addition, P4 signals through membrane P4 receptors (mPRs) in a rapid nongenomic modality. Despite the established physiological importance of P4 nongenomic signaling, its detailed signal transduction remains elusive. Here, using Xenopus oocyte maturation as a well-established physiological readout of nongenomic P4 signaling, we identify the lipid hydrolase ABHD2 (α/β hydrolase domain-containing protein 2) as an essential mPRβ co-receptor to trigger meiosis. We show using functional assays coupled to unbiased and targeted cell-based lipidomics that ABHD2 possesses a phospholipase A2 (PLA2) activity that requires both P4 and mPRβ. This PLA2 activity bifurcates P4 signaling by inducing mPRβ clathrin-dependent endocytosis and producing lipid messengers that are G-protein coupled receptors agonists. Therefore, P4 drives meiosis by inducing the ABHD2 PLA2 activity that requires both mPRβ and ABHD2 as obligate co-receptors.

Plasma proteomic associations with genetics and health in the UK Biobank

The Pharma Proteomics Project is a precompetitive biopharmaceutical consortium characterizing the plasma proteomic profiles of 54,219 UK Biobank participants. Here we provide a detailed summary of this initiative, including technical and biological validations, insights into proteomic disease signatures, and prediction modelling for various demographic and health indicators. We present comprehensive protein quantitative trait locus (pQTL) mapping of 2,923 proteins that identifies 14,287 primary genetic associations, of which 81% are previously undescribed, alongside ancestry-specific pQTL mapping in non-European individuals. The study provides an updated characterization of the genetic architecture of the plasma proteome, contextualized with projected pQTL discovery rates as sample sizes and proteomic assay coverages increase over time. We offer extensive insights into trans pQTLs across multiple biological domains, highlight genetic influences on ligand-receptor interactions and pathway perturbations across a diverse collection of cytokines and complement networks, and illustrate long-range epistatic effects of ABO blood group and FUT2 secretor status on proteins with gastrointestinal tissue-enriched expression. We demonstrate the utility of these data for drug discovery by extending the genetic proxied effects of protein targets, such as PCSK9, on additional endpoints, and disentangle specific genes and proteins perturbed at loci associated with COVID-19 susceptibility. This public-private partnership provides the scientific community with an open-access proteomics resource of considerable breadth and depth to help to elucidate the biological mechanisms underlying proteo-genomic discoveries and accelerate the development of biomarkers, predictive models and therapeutics1.

Genetic determinants of Vitamin D deficiency in the Middle Eastern Qatari population: a genome-wide association study

Introduction

Epidemiological studies have consistently revealed that Vitamin D deficiency is most prevalent in Middle Eastern countries. However, research on the impact of genetic loci and polygenic models related to Vitamin D has primarily focused on European populations.

Methods

We conducted the first genome-wide association study to identify genetic determinants of Vitamin D levels in Middle Easterners using a whole genome sequencing approach in 6,047 subjects from the Qatar Biobank (QBB) project. We performed a GWAS meta-analysis, combining the QBB cohort with recent European GWAS data from the UK Biobank (involving 345,923 individuals). Additionally, we evaluated the performance of European-derived polygenic risk scores using UK Biobank data in the QBB cohort.

Results

Our study identified an association between a variant in a known locus for the group-specific component gene (GC), specifically rs2298850 (p-value = 1.71 × 10-08, Beta = -0.1285), and Vitamin D levels. Furthermore, our GWAS meta-analysis identified two novel variants at a known locus on chromosome 11, rs67609747 and rs1945603, that reached the GWAS significance threshold. Notably, we observed a moderately high heritability of Vitamin D, estimated at 18%, compared to Europeans. Despite the lower predictive performance of Vitamin D levels in Qataris compared to Europeans, the European-derived polygenic risk scores exhibited significant links to Vitamin D deficiency risk within the QBB cohort.

Conclusion

This novel study reveals the genetic architecture contributing to Vitamin D deficiency in the Qatari population, emphasizing the genetic heterogeneity across different populations.

Circulating metabolites modulated by diet are associated with depression

Metabolome reflects the interplay of genome and exposome at molecular level and thus can provide deep insights into the pathogenesis of a complex disease like major depression. To identify metabolites associated with depression we performed a metabolome-wide association analysis in 13,596 participants from five European population-based cohorts characterized for depression, and circulating metabolites using ultra high-performance liquid chromatography/tandem accurate mass spectrometry (UHPLC/MS/MS) based Metabolon platform. We tested 806 metabolites covering a wide range of biochemical processes including those involved in lipid, amino-acid, energy, carbohydrate, xenobiotic and vitamin metabolism for their association with depression. In a conservative model adjusting for life style factors and cardiovascular and antidepressant medication use we identified 8 metabolites, including 6 novel, significantly associated with depression. In individuals with depression, increased levels of retinol (vitamin A), 1-palmitoyl-2-palmitoleoyl-GPC (16:0/16:1) (lecithin) and mannitol/sorbitol and lower levels of hippurate, 4-hydroxycoumarin, 2-aminooctanoate (alpha-aminocaprylic acid), 10-undecenoate (11:1n1) (undecylenic acid), 1-linoleoyl-GPA (18:2) (lysophosphatidic acid; LPA 18:2) are observed. These metabolites are either directly food derived or are products of host and gut microbial metabolism of food-derived products. Our Mendelian randomization analysis suggests that low hippurate levels may be in the causal pathway leading towards depression. Our findings highlight putative actionable targets for depression prevention that are easily modifiable through diet interventions.

The HuMet Repository: Watching human metabolism at work

The human metabolism constantly responds to stimuli such as food intake, fasting, exercise, and stress, triggering adaptive biochemical processes across multiple metabolic pathways. To understand the role of these processes and disruptions thereof in health and disease, detailed documentation of healthy metabolic responses is needed but still scarce on a time-resolved metabolome-wide level. Here, we present the HuMet Repository, a web-based resource for exploring dynamic metabolic responses to six physiological challenges (exercise, 36 h fasting, oral glucose and lipid loads, mixed meal, cold stress) in healthy subjects. For building this resource, we integrated existing and newly derived metabolomics data measured in blood, urine, and breath samples of 15 young healthy men at up to 56 time points during the six highly standardized challenge tests conducted over four days. The data comprise 1.1 million data points acquired on multiple platforms with temporal profiles of 2,656 metabolites from a broad range of biochemical pathways. By embedding the dataset into an interactive web application, we enable users to easily access, search, filter, analyze, and visualize the time-resolved metabolomic readouts and derived results. Users can put metabolites into their larger context by identifying metabolites with similar trajectories or by visualizing metabolites within holistic metabolic networks to pinpoint pathways of interest. In three showcases, we outline the value of the repository for gaining biological insights and generating hypotheses by analyzing the wash-out of dietary markers, the complementarity of metabolomics platforms in dynamic versus cross-sectional data, and similarities and differences in systemic metabolic responses across challenges. With its comprehensive collection of time-resolved metabolomics data, the HuMet Repository, freely accessible at https://humet.org/, is a reference for normal, healthy responses to metabolic challenges in young males. It will enable researchers with and without computational expertise, to flexibly query the data for their own research into the dynamics of human metabolism.

Proteomic analysis of 92 circulating proteins and their effects in cardiometabolic diseases

BACKGROUND

Human plasma contains a wide variety of circulating proteins. These proteins can be important clinical biomarkers in disease and also possible drug targets. Large scale genomics studies of circulating proteins can identify genetic variants that lead to relative protein abundance.

METHODS

We conducted a meta-analysis on genome-wide association studies of autosomal chromosomes in 22,997 individuals of primarily European ancestry across 12 cohorts to identify protein quantitative trait loci (pQTL) for 92 cardiometabolic associated plasma proteins.

RESULTS

We identified 503 (337 cis and 166 trans) conditionally independent pQTLs, including several novel variants not reported in the literature. We conducted a sex-stratified analysis and found that 118 (23.5%) of pQTLs demonstrated heterogeneity between sexes. The direction of effect was preserved but there were differences in effect size and significance. Additionally, we annotate trans-pQTLs with nearest genes and report plausible biological relationships. Using Mendelian randomization, we identified causal associations for 18 proteins across 19 phenotypes, of which 10 have additional genetic colocalization evidence. We highlight proteins associated with a constellation of cardiometabolic traits including angiopoietin-related protein 7 (ANGPTL7) and Semaphorin 3F (SEMA3F).

CONCLUSION

Through large-scale analysis of protein quantitative trait loci, we provide a comprehensive overview of common variants associated with plasma proteins. We highlight possible biological relationships which may serve as a basis for further investigation into possible causal roles in cardiometabolic diseases.

Systemic Markers of Lung Function and Forced Expiratory Volume in 1 Second Decline across Diverse Cohorts

Rationale: Chronic obstructive pulmonary disease (COPD) is a complex disease characterized by airway obstruction and accelerated lung function decline. Our understanding of systemic protein biomarkers associated with COPD remains incomplete. Objectives: To determine what proteins and pathways are associated with impaired pulmonary function in a diverse population. Methods: We studied 6,722 participants across six cohort studies with both aptamer-based proteomic and spirometry data (4,566 predominantly White participants in a discovery analysis and 2,156 African American cohort participants in a validation). In linear regression models, we examined protein associations with baseline forced expiratory volume in 1 second (FEV1) and FEV1/forced vital capacity (FVC). In linear mixed effects models, we investigated the associations of baseline protein levels with rate of FEV1 decline (ml/yr) in 2,777 participants with up to 7 years of follow-up spirometry. Results: We identified 254 proteins associated with FEV1 in our discovery analyses, with 80 proteins validated in the Jackson Heart Study. Novel validated protein associations include kallistatin serine protease inhibitor, growth differentiation factor 2, and tumor necrosis factor-like weak inducer of apoptosis (discovery β = 0.0561, Q = 4.05 × 10-10; β  = 0.0421, Q = 1.12 × 10-3; and β = 0.0358, Q = 1.67 × 10-3, respectively). In longitudinal analyses within cohorts with follow-up spirometry, we identified 15 proteins associated with FEV1 decline (Q < 0.05), including elafin leukocyte elastase inhibitor and mucin-associated TFF2 (trefoil factor 2; β = -4.3 ml/yr, Q = 0.049; β = -6.1 ml/yr, Q = 0.032, respectively). Pathways and processes highlighted by our study include aberrant extracellular matrix remodeling, enhanced innate immune response, dysregulation of angiogenesis, and coagulation. Conclusions: In this study, we identify and validate novel biomarkers and pathways associated with lung function traits in a racially diverse population. In addition, we identify novel protein markers associated with FEV1 decline. Several protein findings are supported by previously reported genetic signals, highlighting the plausibility of certain biologic pathways. These novel proteins might represent markers for risk stratification, as well as novel molecular targets for treatment of COPD.

Auto-immunoproteomics analysis of COVID-19 ICU patients revealed increased levels of autoantibodies related to the male reproductive system

Background: Coronavirus disease (COVID-19) manifests many clinical symptoms, including an exacerbated immune response and cytokine storm. Autoantibodies in COVID-19 may have severe prodromal effects that are poorly understood. The interaction between these autoantibodies and self-antigens can result in systemic inflammation and organ dysfunction. However, the role of autoantibodies in COVID-19 complications has yet to be fully understood. Methods: The current investigation screened two independent cohorts of 97 COVID-19 patients [discovery (Disc) cohort from Qatar (case = 49 vs. control = 48) and replication (Rep) cohort from New York (case = 48 vs. control = 28)] utilizing high-throughput KoRectly Expressed (KREX) Immunome protein-array technology. Total IgG autoantibody responses were evaluated against 1,318 correctly folded and full-length human proteins. Samples were randomly applied on the precoated microarray slides for 2 h. Cy3-labeled secondary antibodies were used to detect IgG autoantibody response. Slides were scanned at a fixed gain setting using the Agilent fluorescence microarray scanner, generating a 16-bit TIFF file. Group comparisons were performed using a linear model and Fisher's exact test. Differentially expressed proteins were used for KEGG and WIKIpathway annotation to determine pathways in which the proteins of interest were significantly over-represented. Results and conclusion: Autoantibody responses to 57 proteins were significantly altered in the COVID-19 Disc cohort compared to healthy controls (p ≤ 0.05). The Rep cohort had altered autoantibody responses against 26 proteins compared to non-COVID-19 ICU patients who served as controls. Both cohorts showed substantial similarities (r 2 = 0.73) and exhibited higher autoantibody responses to numerous transcription factors, immunomodulatory proteins, and human disease markers. Analysis of the combined cohorts revealed elevated autoantibody responses against SPANXN4, STK25, ATF4, PRKD2, and CHMP3 proteins in COVID-19 patients. The sequences for SPANXN4 and STK25 were cross-validated using sequence alignment tools. ELISA and Western blot further verified the autoantigen-autoantibody response of SPANXN4. SPANXN4 is essential for spermiogenesis and male fertility, which may predict a potential role for this protein in COVID-19-associated male reproductive tract complications, and warrants further research.

Identification of candidate metabolite biomarkers for metabolic syndrome and its five components in population-based human cohorts

BACKGROUND

Metabolic Syndrome (MetS) is characterized by risk factors such as abdominal obesity, hypertriglyceridemia, low high-density lipoprotein cholesterol (HDL-C), hypertension, and hyperglycemia, which contribute to the development of cardiovascular disease and type 2 diabetes. Here, we aim to identify candidate metabolite biomarkers of MetS and its associated risk factors to better understand the complex interplay of underlying signaling pathways.

METHODS

We quantified serum samples of the KORA F4 study participants (N = 2815) and analyzed 121 metabolites. Multiple regression models adjusted for clinical and lifestyle covariates were used to identify metabolites that were Bonferroni significantly associated with MetS. These findings were replicated in the SHIP-TREND-0 study (N = 988) and further analyzed for the association of replicated metabolites with the five components of MetS. Database-driven networks of the identified metabolites and their interacting enzymes were also constructed.

RESULTS

We identified and replicated 56 MetS-specific metabolites: 13 were positively associated (e.g., Val, Leu/Ile, Phe, and Tyr), and 43 were negatively associated (e.g., Gly, Ser, and 40 lipids). Moreover, the majority (89%) and minority (23%) of MetS-specific metabolites were associated with low HDL-C and hypertension, respectively. One lipid, lysoPC a C18:2, was negatively associated with MetS and all of its five components, indicating that individuals with MetS and each of the risk factors had lower concentrations of lysoPC a C18:2 compared to corresponding controls. Our metabolic networks elucidated these observations by revealing impaired catabolism of branched-chain and aromatic amino acids, as well as accelerated Gly catabolism.

CONCLUSION

Our identified candidate metabolite biomarkers are associated with the pathophysiology of MetS and its risk factors. They could facilitate the development of therapeutic strategies to prevent type 2 diabetes and cardiovascular disease. For instance, elevated levels of lysoPC a C18:2 may protect MetS and its five risk components. More in-depth studies are necessary to determine the mechanism of key metabolites in the MetS pathophysiology.

Differences and commonalities in the genetic architecture of protein quantitative trait loci in European and Arab populations

Polygenic scores (PGS) can identify individuals at risk of adverse health events and guide genetics-based personalized medicine. However, it is not clear how well PGS translate between different populations, limiting their application to well-studied ethnicities. Proteins are intermediate traits linking genetic predisposition and environmental factors to disease, with numerous blood circulating protein levels representing functional readouts of disease-related processes. We hypothesized that studying the genetic architecture of a comprehensive set of blood-circulating proteins between a European and an Arab population could shed fresh light on the translatability of PGS to understudied populations. We therefore conducted a genome-wide association study with whole-genome sequencing data using 1301 proteins measured on the SOMAscan aptamer-based affinity proteomics platform in 2935 samples of Qatar Biobank and evaluated the replication of protein quantitative traits (pQTLs) from European studies in an Arab population. Then, we investigated the colocalization of shared pQTL signals between the two populations. Finally, we compared the performance of protein PGS derived from a Caucasian population in a European and an Arab cohort. We found that the majority of shared pQTL signals (81.8%) colocalized between both populations. About one-third of the genetic protein heritability was explained by protein PGS derived from a European cohort, with protein PGS performing ~20% better in Europeans when compared to Arabs. Our results are relevant for the translation of PGS to non-Caucasian populations, as well as for future efforts to extend genetic research to understudied populations.

Multi-ancestry genome-wide association analyses improve resolution of genes and pathways influencing lung function and chronic obstructive pulmonary disease risk

Lung-function impairment underlies chronic obstructive pulmonary disease (COPD) and predicts mortality. In the largest multi-ancestry genome-wide association meta-analysis of lung function to date, comprising 580,869 participants, we identified 1,020 independent association signals implicating 559 genes supported by ≥2 criteria from a systematic variant-to-gene mapping framework. These genes were enriched in 29 pathways. Individual variants showed heterogeneity across ancestries, age and smoking groups, and collectively as a genetic risk score showed strong association with COPD across ancestry groups. We undertook phenome-wide association studies for selected associated variants as well as trait and pathway-specific genetic risk scores to infer possible consequences of intervening in pathways underlying lung function. We highlight new putative causal variants, genes, proteins and pathways, including those targeted by existing drugs. These findings bring us closer to understanding the mechanisms underlying lung function and COPD, and should inform functional genomics experiments and potentially future COPD therapies.

Urine-based multi-omic comparative analysis of COVID-19 and bacterial sepsis-induced ARDS

BACKGROUND

Acute respiratory distress syndrome (ARDS), a life-threatening condition during critical illness, is a common complication of COVID-19. It can originate from various disease etiologies, including severe infections, major injury, or inhalation of irritants. ARDS poses substantial clinical challenges due to a lack of etiology-specific therapies, multisystem involvement, and heterogeneous, poor patient outcomes. A molecular comparison of ARDS groups holds the potential to reveal common and distinct mechanisms underlying ARDS pathogenesis.

METHODS

We performed a comparative analysis of urine-based metabolomics and proteomics profiles from COVID-19 ARDS patients (n = 42) and bacterial sepsis-induced ARDS patients (n = 17). To this end, we used two different approaches, first we compared the molecular omics profiles between ARDS groups, and second, we correlated clinical manifestations within each group with the omics profiles.

RESULTS

The comparison of the two ARDS etiologies identified 150 metabolites and 70 proteins that were differentially abundant between the two groups. Based on these findings, we interrogated the interplay of cell adhesion/extracellular matrix molecules, inflammation, and mitochondrial dysfunction in ARDS pathogenesis through a multi-omic network approach. Moreover, we identified a proteomic signature associated with mortality in COVID-19 ARDS patients, which contained several proteins that had previously been implicated in clinical manifestations frequently linked with ARDS pathogenesis.

CONCLUSION

In summary, our results provide evidence for significant molecular differences in ARDS patients from different etiologies and a potential synergy of extracellular matrix molecules, inflammation, and mitochondrial dysfunction in ARDS pathogenesis. The proteomic mortality signature should be further investigated in future studies to develop prediction models for COVID-19 patient outcomes.

Characterization of exercise-induced hemolysis in endurance horses

Exercise-induced hemolysis occurs as the result of intense physical exercise and is caused by metabolic and mechanical factors including repeated muscle contractions leading to capillary vessels compression, vasoconstriction of internal organs and foot strike among others. We hypothesized that exercise-induced hemolysis occurred in endurance racehorses and its severity was associated with the intensity of exercise. To provide further insight into the hemolysis of endurance horses, the aim of the study was to deployed a strategy for small molecules (metabolites) profiling, beyond standard molecular methods. The study included 47 Arabian endurance horses competing for either 80, 100, or 120 km distances. Blood plasma was collected before and after the competition and analyzed macroscopically, by ELISA and non-targeted metabolomics with liquid chromatography-mass spectrometry. A significant increase in all hemolysis parameters was observed after the race, and an association was found between the measured parameters, average speed, and distance completed. Levels of hemolysis markers were highest in horses eliminated for metabolic reasons in comparison to finishers and horses eliminated for lameness (gait abnormality), which may suggest a connection between the intensity of exercise, metabolic challenges, and hemolysis. Utilization of omics methods alongside conventional methods revealed a broader insight into the exercise-induced hemolysis process by displaying, apart from commonly measured hemoglobin and haptoglobin, levels of hemoglobin degradation metabolites. Obtained results emphasized the importance of respecting horse limitations in regard to speed and distance which, if underestimated, may lead to severe damages.

The Genetic Spectrum of Maturity-Onset Diabetes of the Young (MODY) in Qatar, a Population-Based Study

Maturity-onset diabetes of the young (MODY) is a rare monogenic form of diabetes mellitus. In this study, we estimated the prevalence and genetic spectrum of MODY in the Middle Eastern population of Qatar using whole-genome sequencing (WGS) of 14,364 subjects from the population-based Qatar biobank (QBB) cohort. We focused our investigations on 14 previously identified genes ascribed to the cause of MODY and two potentially novel MODY-causing genes, RFX6 and NKX6-1. Genetic variations within the 16 MODY-related genes were assessed for their pathogenicity to identify disease-causing mutations. Analysis of QBB phenotype data revealed 72 subjects (0.5%) with type 1 diabetes, 2915 subjects (20.3%) with type 2 diabetes and 11,377 (79.2%) without diabetes. We identified 22 mutations in 67 subjects that were previously reported in the Human Genetic Mutation Database (HGMD) as disease-causing (DM) or likely disease causing (DM?) for MODY. We also identified 28 potentially novel MODY-causing mutations, predicted to be among the top 1% most deleterious mutations in the human genome, which showed complete (100%) disease penetrance in 34 subjects. Overall, we estimated that MODY accounts for around 2.2-3.4% of diabetes patients in Qatar. This is the first population-based study to determine the genetic spectrum and estimate the prevalence of MODY in the Middle East. Further research to characterize the newly identified mutations is warranted.

Metabolic and proteomic signatures of type 2 diabetes subtypes in an Arab population

Type 2 diabetes (T2D) has a heterogeneous etiology influencing its progression, treatment, and complications. A data driven cluster analysis in European individuals with T2D previously identified four subtypes: severe insulin deficient (SIDD), severe insulin resistant (SIRD), mild obesity-related (MOD), and mild age-related (MARD) diabetes. Here, the clustering approach was applied to individuals with T2D from the Qatar Biobank and validated in an independent set. Cluster-specific signatures of circulating metabolites and proteins were established, revealing subtype-specific molecular mechanisms, including activation of the complement system with features of autoimmune diabetes and reduced 1,5-anhydroglucitol in SIDD, impaired insulin signaling in SIRD, and elevated leptin and fatty acid binding protein levels in MOD. The MARD cluster was the healthiest with metabolomic and proteomic profiles most similar to the controls. We have translated the T2D subtypes to an Arab population and identified distinct molecular signatures to further our understanding of the etiology of these subtypes.

Identification of PCSK9-like human gene knockouts using metabolomics, proteomics, and whole-genome sequencing in a consanguineous population

Natural human knockouts of genes associated with desirable outcomes, such as PCSK9 with low levels of LDL-cholesterol, can lead to the discovery of new drug targets and treatments. Rare loss-of-function variants are more likely to be found in the homozygous state in consanguineous populations, and deep molecular phenotyping of blood samples from homozygous carriers can help to discriminate between silent and functional variants. Here, we combined whole-genome sequencing with proteomics and metabolomics for 2,935 individuals from the Qatar Biobank (QBB) to evaluate the power of this approach for finding genes of clinical and pharmaceutical interest. As proof-of-concept, we identified a homozygous carrier of a very rare PCSK9 variant with extremely low circulating PCSK9 levels and low LDL. Our study demonstrates that the chances of finding such variants are about 168 times higher in QBB compared with GnomAD and emphasizes the potential of consanguineous populations for drug discovery.

Rare and common genetic determinants of metabolic individuality and their effects on human health

Garrod's concept of 'chemical individuality' has contributed to comprehension of the molecular origins of human diseases. Untargeted high-throughput metabolomic technologies provide an in-depth snapshot of human metabolism at scale. We studied the genetic architecture of the human plasma metabolome using 913 metabolites assayed in 19,994 individuals and identified 2,599 variant-metabolite associations (P < 1.25 × 10^-11) within 330 genomic regions, with rare variants (minor allele frequency ≤ 1%) explaining 9.4% of associations. Jointly modeling metabolites in each region, we identified 423 regional, co-regulated, variant-metabolite clusters called genetically influenced metabotypes. We assigned causal genes for 62.4% of these genetically influenced metabotypes, providing new insights into fundamental metabolite physiology and clinical relevance, including metabolite-guided discovery of potential adverse drug effects (DPYD and SRD5A2). We show strong enrichment of inborn errors of metabolism-causing genes, with examples of metabolite associations and clinical phenotypes of non-pathogenic variant carriers matching characteristics of the inborn errors of metabolism. Systematic, phenotypic follow-up of metabolite-specific genetic scores revealed multiple potential etiological relationships.

Machine learning integration of multimodal data identifies key features of blood pressure regulation

BACKGROUND

Association studies have identified several biomarkers for blood pressure and hypertension, but a thorough understanding of their mutual dependencies is lacking. By integrating two different high-throughput datasets, biochemical and dietary data, we aim to understand the multifactorial contributors of blood pressure (BP).

METHODS

We included 4,863 participants from TwinsUK with concurrent BP, metabolomics, genomics, biochemical measures, and dietary data. We used 5-fold cross-validation with the machine learning XGBoost algorithm to identify features of importance in context of one another in TwinsUK (80% training, 20% test). The features tested in TwinsUK were then probed using the same algorithm in an independent dataset of 2,807 individuals from the Qatari Biobank (QBB).

FINDINGS

Our model explained 39·2% [4·5%, MAE:11·32 mmHg (95%CI, +/- 0·65)] of the variance in systolic BP (SBP) in TwinsUK. Of the top 50 features, the most influential non-demographic variables were dihomo-linolenate, cis-4-decenoyl carnitine, lactate, chloride, urate, and creatinine along with dietary intakes of total, trans and saturated fat. We also highlight the incremental value of each included dimension. Furthermore, we replicated our model in the QBB [SBP variance explained = 45·2% (13·39%)] cohort and 30 of the top 50 features overlapped between cohorts.

INTERPRETATION

We show that an integrated analysis of omics, biochemical and dietary data improves our understanding of their in-between relationships and expands the range of potential biomarkers for blood pressure. Our results point to potentially key biological pathways to be prioritised for mechanistic studies.

FUNDING

Chronic Disease Research Foundation, Medical Research Council, Wellcome Trust, Qatar Foundation.

Dynamic patterns of postprandial metabolic responses to three dietary challenges

Food intake triggers extensive changes in the blood metabolome. The kinetics of these changes depend on meal composition and on intrinsic, health-related characteristics of each individual, making the assessment of changes in the postprandial metabolome an opportunity to assess someone's metabolic status. To enable the usage of dietary challenges as diagnostic tools, profound knowledge about changes that occur in the postprandial period in healthy individuals is needed. In this study, we characterize the time-resolved changes in plasma levels of 634 metabolites in response to an oral glucose tolerance test (OGTT), an oral lipid tolerance test (OLTT), and a mixed meal (SLD) in healthy young males (n = 15). Metabolite levels for samples taken at different time points (20 per individual) during the challenges were available from targeted (132 metabolites) and non-targeted (502 metabolites) metabolomics. Almost half of the profiled metabolites (n = 308) showed a significant change in at least one challenge, thereof 111 metabolites responded exclusively to one particular challenge. Examples include azelate, which is linked to ω-oxidation and increased only in OLTT, and a fibrinogen cleavage peptide that has been linked to a higher risk of cardiovascular events in diabetes patients and increased only in OGTT, making its postprandial dynamics a potential target for risk management. A pool of 89 metabolites changed their plasma levels during all three challenges and represents the core postprandial response to food intake regardless of macronutrient composition. We used fuzzy c-means clustering to group these metabolites into eight clusters based on commonalities of their dynamic response patterns, with each cluster following one of four primary response patterns: (i) “decrease-increase” (valley-like) with fatty acids and acylcarnitines indicating the suppression of lipolysis, (ii) “increase-decrease” (mountain-like) including a cluster of conjugated bile acids and the glucose/insulin cluster, (iii) “steady decrease” with metabolites reflecting a carryover from meals prior to the study, and (iv) “mixed” decreasing after the glucose challenge and increasing otherwise. Despite the small number of subjects, the diversity of the challenges and the wealth of metabolomic data make this study an important step toward the characterization of postprandial responses and the identification of markers of metabolic processes regulated by food intake.

Development and characterization of microsatellite primers for Triops granarius (Branchiopoda: Notostraca) using MiSeq technology

Background

Next-generation sequencing technology has allowed for the rapid development of microsatellites, neutral polymorphic markers that can be used for the analysis of population structure.

Methods and Results

In this study, we performed whole-genome sequencing using the Illumina MiSeq system and de novo assembly to design microsatellite primers for Triops granarius populations in Qatar. The developed microsatellites are suitable for future studies of genetic structuring among geographically isolated freshwater pools. A total of 23 different primer pairs produced typical microsatellite results, with each pair successfully amplified in up to 40 individuals. Only five of the loci produced a significant departure from Hardy-Weinberg equilibrium.

Conclusions

Some of the underlying mechanisms regarding the few loci that deviated from HWE may be further investigated to determine the source of deviation. As T. granarius is the most widely distributed species of the family, the development of these molecular markers would be useful for conducting population genetics and biogeographical studies broadly.

Multi-omic comparative analysis of COVID-19 and bacterial sepsis-induced ARDS

BACKGROUND

Acute respiratory distress syndrome (ARDS), a life-threatening condition characterized by hypoxemia and poor lung compliance, is associated with high mortality. ARDS induced by COVID-19 has similar clinical presentations and pathological manifestations as non-COVID-19 ARDS. However, COVID-19 ARDS is associated with a more protracted inflammatory respiratory failure compared to traditional ARDS. Therefore, a comprehensive molecular comparison of ARDS of different etiologies groups may pave the way for more specific clinical interventions.

METHODS AND FINDINGS

In this study, we compared COVID-19 ARDS (n = 43) and bacterial sepsis-induced (non-COVID-19) ARDS (n = 24) using multi-omic plasma profiles covering 663 metabolites, 1,051 lipids, and 266 proteins. To address both between- and within- ARDS group variabilities we followed two approaches. First, we identified 706 molecules differently abundant between the two ARDS etiologies, revealing more than 40 biological processes differently regulated between the two groups. From these processes, we assembled a cascade of therapeutically relevant pathways downstream of sphingosine metabolism. The analysis suggests a possible overactivation of arginine metabolism involved in long-term sequelae of ARDS and highlights the potential of JAK inhibitors to improve outcomes in bacterial sepsis-induced ARDS. The second part of our study involved the comparison of the two ARDS groups with respect to clinical manifestations. Using a data-driven multi-omic network, we identified signatures of acute kidney injury (AKI) and thrombocytosis within each ARDS group. The AKI-associated network implicated mitochondrial dysregulation which might lead to post-ARDS renal-sequalae. The thrombocytosis-associated network hinted at a synergy between prothrombotic processes, namely IL-17, MAPK, TNF signaling pathways, and cell adhesion molecules. Thus, we speculate that combination therapy targeting two or more of these processes may ameliorate thrombocytosis-mediated hypercoagulation.

CONCLUSION

We present a first comprehensive molecular characterization of differences between two ARDS etiologies-COVID-19 and bacterial sepsis. Further investigation into the identified pathways will lead to a better understanding of the pathophysiological processes, potentially enabling novel therapeutic interventions.

Dynamics and determinants of human plasma bile acid profiles during dietary challenges

In recent years, bile acids (BA) have received great interest due to their pleiotropic biological activity and the presence of plasma membrane-bound and nuclear receptors. Moreover, BA in blood have been identified by metabolite screening approaches as biomarkers that are associated with various diseases and even with a human longevity phenotype. With the growing interest in the microbiota contribution to the health-disease trajectory, BA that undergo deconjugation and other modifications by bacteria in the large intestine have become a prime target as a microbiome diversity modifier. We here profiled BA by a quantitative and a semiquantitative approach in 15 healthy and phenotypically very similar young individuals for over a 36-h fasting period, an oral glucose tolerance test (OGTT), and an oral lipid tolerance test (OLTT). We demonstrate a remarkable heterogeneity of the responses and describe the different dynamics of the plasma changes that likely originate from different routes by which BA enters the peripheral blood, and that may represent a direct secretion from the liver into the blood and a route that reaches the blood as a spill-over after passing from the gallbladder through the intestine and the portal system. We discuss the finding that an individual transport process involved in the passage of BA could be a critical determinant in the kinetics of plasma appearance and the overall phenotypic variability found.

Multi-omic comparative analysis of COVID-19 and bacterial sepsis-induced ARDS

Background

Acute respiratory distress syndrome (ARDS), a life-threatening condition characterized by hypoxemia and poor lung compliance, is associated with high mortality. ARDS induced by COVID-19 has similar clinical presentations and pathological manifestations as non-COVID-19 ARDS. However, COVID-19 ARDS is associated with a more protracted inflammatory respiratory failure compared to traditional ARDS. Therefore, a comprehensive molecular comparison of ARDS of different etiologies groups may pave the way for more specific clinical interventions.

Methods and Findings

In this study, we compared COVID-19 ARDS (n=43) and bacterial sepsis-induced (non-COVID-19) ARDS (n=24) using multi-omic plasma profiles covering 663 metabolites, 1,051 lipids, and 266 proteins. To address both between- and within-ARDS group variabilities we followed two approaches. First, we identified 706 molecules differently abundant between the two ARDS etiologies, revealing more than 40 biological processes differently regulated between the two groups. From these processes, we assembled a cascade of therapeutically relevant pathways downstream of sphingosine metabolism. The analysis suggests a possible overactivation of arginine metabolism involved in long-term sequelae of ARDS and highlights the potential of JAK inhibitors to improve outcomes in bacterial sepsis-induced ARDS. The second part of our study involved the comparison of the two ARDS groups with respect to clinical manifestations. Using a data-driven multi-omic network, we identified signatures of acute kidney injury (AKI) and thrombocytosis within each ARDS group. The AKI-associated network implicated mitochondrial dysregulation which might lead to post-ARDS renal-sequalae. The thrombocytosis-associated network hinted at a synergy between prothrombotic processes, namely IL-17, MAPK, TNF signaling pathways, and cell adhesion molecules. Thus, we speculate that combination therapy targeting two or more of these processes may ameliorate thrombocytosis-mediated hypercoagulation.

Conclusion

We present a first comprehensive molecular characterization of differences between two ARDS etiologies - COVID-19 and bacterial sepsis. Further investigation into the identified pathways will lead to a better understanding of the pathophysiological processes, potentially enabling novel therapeutic interventions.

Urine-based multi-omic comparative analysis of COVID-19 and bacterial sepsis-induced ARDS

Acute respiratory distress syndrome (ARDS), a life-threatening condition during critical illness, is a common complication of COVID-19. It can originate from various disease etiologies, including severe infections, major injury, or inhalation of irritants. ARDS poses substantial clinical challenges due to a lack of etiology-specific therapies, multisystem involvement, and heterogeneous, poor patient outcomes. A molecular comparison of ARDS groups holds the potential to reveal common and distinct mechanisms underlying ARDS pathogenesis. In this study, we performed a comparative analysis of urine-based metabolomics and proteomics profiles from COVID-19 ARDS patients (n = 42) and bacterial sepsis-induced ARDS patients (n = 17). The comparison of these ARDS etiologies identified 150 metabolites and 70 proteins that were differentially abundant between the two groups. Based on these findings, we interrogated the interplay of cell adhesion/extracellular matrix molecules, inflammation, and mitochondrial dysfunction in ARDS pathogenesis through a multi-omic network approach. Moreover, we identified a proteomic signature associated with mortality in COVID-19 ARDS patients, which contained several proteins that had previously been implicated in clinical manifestations frequently linked with ARDS pathogenesis. In summary, our results provide evidence for significant molecular differences in ARDS patients from different etiologies and a potential synergy of extracellular matrix molecules, inflammation, and mitochondrial dysfunction in ARDS pathogenesis. The proteomic mortality signature should be further investigated in future studies to develop prediction models for COVID-19 patient outcomes.

Integrative metabolomic and proteomic signatures define clinical outcomes in severe COVID-19

The coronavirus disease-19 (COVID-19) pandemic has ravaged global healthcare with previously unseen levels of morbidity and mortality. In this study, we performed large-scale integrative multi-omics analyses of serum obtained from COVID-19 patients with the goal of uncovering novel pathogenic complexities of this disease and identifying molecular signatures that predict clinical outcomes. We assembled a network of protein-metabolite interactions through targeted metabolomic and proteomic profiling in 330 COVID-19 patients compared to 97 non-COVID, hospitalized controls. Our network identified distinct protein-metabolite cross talk related to immune modulation, energy and nucleotide metabolism, vascular homeostasis, and collagen catabolism. Additionally, our data linked multiple proteins and metabolites to clinical indices associated with long-term mortality and morbidity. Finally, we developed a novel composite outcome measure for COVID-19 disease severity based on metabolomics data. The model predicts severe disease with a concordance index of around 0.69, and shows high predictive power of 0.83-0.93 in two independent datasets.

Variational autoencoders learn transferrable representations of metabolomics data

Dimensionality reduction approaches are commonly used for the deconvolution of high-dimensional metabolomics datasets into underlying core metabolic processes. However, current state-of-the-art methods are widely incapable of detecting nonlinearities in metabolomics data. Variational Autoencoders (VAEs) are a deep learning method designed to learn nonlinear latent representations which generalize to unseen data. Here, we trained a VAE on a large-scale metabolomics population cohort of human blood samples consisting of over 4500 individuals. We analyzed the pathway composition of the latent space using a global feature importance score, which demonstrated that latent dimensions represent distinct cellular processes. To demonstrate model generalizability, we generated latent representations of unseen metabolomics datasets on type 2 diabetes, acute myeloid leukemia, and schizophrenia and found significant correlations with clinical patient groups. Notably, the VAE representations showed stronger effects than latent dimensions derived by linear and non-linear principal component analysis. Taken together, we demonstrate that the VAE is a powerful method that learns biologically meaningful, nonlinear, and transferrable latent representations of metabolomics data.

Cross-Sectional Blood Metabolite Markers of Hypertension: A Multicohort Analysis of 44,306 Individuals from the COnsortium of METabolomics Studies

Hypertension is the main modifiable risk factor for cardiovascular morbidity and mortality but discovering molecular mechanisms for targeted treatment has been challenging. Here we investigate associations of blood metabolite markers with hypertension by integrating data from nine intercontinental cohorts from the COnsortium of METabolomics Studies. We included 44,306 individuals with circulating metabolites (up to 813). Metabolites were aligned and inverse normalised to allow intra-platform comparison. Logistic models adjusting for covariates were performed in each cohort and results were combined using random-effect inverse-variance meta-analyses adjusting for multiple testing. We further conducted canonical pathway analysis to investigate the pathways underlying the hypertension-associated metabolites. In 12,479 hypertensive cases and 31,827 controls without renal impairment, we identified 38 metabolites, associated with hypertension after adjusting for age, sex, body mass index, ethnicity, and multiple testing. Of these, 32 metabolite associations, predominantly lipid (steroids and fatty acyls) and organic acids (amino-, hydroxy-, and keto-acids) remained after further adjusting for comorbidities and dietary intake. Among the identified metabolites, 5 were novel, including 2 bile acids, 2 glycerophospholipids, and ketoleucine. Pathway analysis further implicates the role of the amino-acids, serine/glycine, and bile acids in hypertension regulation. In the largest cross-sectional hypertension-metabolomics study to date, we identify 32 circulating metabolites (of which 5 novel and 27 confirmed) that are potentially actionable targets for intervention. Further in-vivo studies are needed to identify their specific role in the aetiology or progression of hypertension.

Detection of infiltrating fibroblasts by single-cell transcriptomics in human kidney allografts

We tested the hypothesis that single-cell RNA-sequencing (scRNA-seq) analysis of human kidney allograft biopsies will reveal distinct cell types and states and yield insights to decipher the complex heterogeneity of alloimmune injury. We selected 3 biopsies of kidney cortex from 3 individuals for scRNA-seq and processed them fresh using an identical protocol on the 10x Chromium platform; (i) HK: native kidney biopsy from a living donor, (ii) AK1: allograft kidney with transplant glomerulopathy, tubulointerstitial fibrosis, and worsening graft function, and (iii) AK2: allograft kidney after successful treatment of active antibody-mediated rejection. We did not study T-cell-mediated rejections. We generated 7217 high-quality single cell transcriptomes. Taking advantage of the recipient-donor sex mismatches revealed by X and Y chromosome autosomal gene expression, we determined that in AK1 with fibrosis, 42 months after transplantation, more than half of the kidney allograft fibroblasts were recipient-derived and therefore likely migratory and graft infiltrative, whereas in AK2 without fibrosis, 84 months after transplantation, most fibroblasts were donor-organ-derived. Furthermore, AK1 was enriched for tubular progenitor cells overexpressing profibrotic extracellular matrix genes. AK2, eight months after successful treatment of rejection, contained plasmablast cells with high expression of immunoglobulins, endothelial cell elaboration of T cell chemoattractant cytokines, and persistent presence of cytotoxic T cells. In addition to these key findings, our analysis revealed unique cell types and states in the kidney. Altogether, single-cell transcriptomics yielded novel mechanistic insights, which could pave the way for individualizing the care of transplant recipients.

Matching Drug Metabolites from Non-Targeted Metabolomics to Self-Reported Medication in the Qatar Biobank Study

Modern metabolomics platforms are able to identify many drug-related metabolites in blood samples. Applied to population-based biobank studies, the detection of drug metabolites can then be used as a proxy for medication use or serve as a validation tool for questionnaire-based health assessments. However, it is not clear how well detection of drug metabolites in blood samples matches information on self-reported medication provided by study participants. Here, we curate free-text responses to a drug-usage questionnaire from 6000 participants of the Qatar Biobank (QBB) using standardized WHO Anatomical Therapeutic Chemical (ATC) Classification System codes and compare the occurrence of these ATC terms to the detection of drug-related metabolites in matching blood plasma samples from 2807 QBB participants for which we collected non-targeted metabolomics data. We found that the detection of 22 drug-related metabolites significantly associated with the self-reported use of the corresponding medication. Good agreement of self-reported medication with non-targeted metabolomics was observed, with self-reported drugs and their metabolites being detected in a same blood sample in 79.4% of the cases. On the other hand, only 29.5% of detected drug metabolites matched to self-reported medication. Possible explanations for differences include under-reporting of over-the-counter medications from the study participants, such as paracetamol, misannotation of low abundance metabolites, such as metformin, and inability of the current methods to detect them. Taken together, our study provides a broad real-world view of what to expect from large non-targeted metabolomics measurements in population-based biobank studies and indicates areas where further improvements can be made.

Qatar Genome: Insights on Genomics from the Middle East

Despite recent biomedical breakthroughs and large genomic studies growing momentum, the Middle Eastern population, home to over 400 million people, is under-represented in the human genome variation databases. Here we describe insights from phase 1 of the Qatar Genome Program with whole genome sequenced 6,047 individuals from Qatar. We identified more than 88 million variants of which 24 million are novel and 23 million are singletons. Consistent with the high consanguinity and founder effects in the region, we found that several rare deleterious variants were more common in the Qatari population while others seem to provide protection against diseases and have shaped the genetic architecture of adaptive phenotypes. These results highlight the value of our data as a resource to advance genetic studies in the Arab and neighbouring Middle Eastern populations and will significantly boost the current efforts to improve our understanding of global patterns of human variations, human history and genetic contributions to health and diseases in diverse populations. This article is protected by copyright. All rights reserved.

Metabolic and Metabo-Clinical Signatures of Type 2 Diabetes, Obesity, Retinopathy, and Dyslipidemia

Macro- and microvascular complications of type 2 diabetes (T2D), obesity, and dyslipidemia share common metabolic pathways. In this study, using a total of 1,300 metabolites from 996 Qatari adults (57% with T2D) and 1,159 metabolites from an independent cohort of 2,618 individuals from the Qatar BioBank (11% with T2D), we identified 373 metabolites associated with T2D, obesity, retinopathy, dyslipidemia, and lipoprotein levels, 161 of which were novel. Novel metabolites included phospholipids, sphingolipids, lysolipids, fatty acids, dipeptides, and metabolites of the urea cycle and xanthine, steroid, and glutathione metabolism. The identified metabolites enrich pathways of oxidative stress, lipotoxicity, glucotoxicity, and proteolysis. Second, we identified 15 patterns we defined as "metabo-clinical signatures." These are clusters of patients with T2D who group together based on metabolite levels and reveal the same clustering in two or more clinical variables (obesity, LDL, HDL, triglycerides, and retinopathy). These signatures revealed metabolic pathways associated with different clinical patterns and identified patients with extreme (very high/low) clinical variables associated with extreme metabolite levels in specific pathways. Among our novel findings are the role of N-acetylmethionine in retinopathy in conjunction with dyslipidemia and the possible roles of N-acetylvaline and pyroglutamine in association with high cholesterol levels and kidney function.

Identification of Robust Protein Associations With COVID-19 Disease Based on Five Clinical Studies

Multiple studies have investigated the role of blood circulating proteins in COVID-19 disease using the Olink affinity proteomics platform. However, study inclusion criteria and sample collection conditions varied between studies, leading to sometimes incongruent associations. To identify the most robust protein markers of the disease and the underlying pathways that are relevant under all conditions, it is essential to identify proteins that replicate most widely. Here we combined the Olink proteomics profiles of two newly recruited COVID-19 studies (N=68 and N=98) with those of three previously published COVID-19 studies (N=383, N=83, N=57). For these studies, three Olink panels (Inflammation and Cardiovascular II & III) with 253 unique proteins were compared. Case/control analysis revealed thirteen proteins (CCL16, CCL7, CXCL10, CCL8, LGALS9, CXCL11, IL1RN, CCL2, CD274, IL6, IL18, MERTK, IFNγ, and IL18R1) that were differentially expressed in COVID-19 patients in all five studies. Except CCL16, which was higher in controls, all proteins were overexpressed in COVID-19 patients. Pathway analysis revealed concordant trends across all studies with pathways related to cytokine-cytokine interaction, IL18 signaling, fluid shear stress and rheumatoid arthritis. Our results reaffirm previous findings related to a COVID-19 cytokine storm syndrome. Cross-study robustness of COVID-19 specific protein expression profiles support the utility of affinity proteomics as a tool and for the identification of potential therapeutic targets.

Advancing Cancer Treatment by Targeting Glutamine Metabolism—A Roadmap

Simple Summary

Dysregulated glutamine metabolism is one of the metabolic features evident in cancer cells when compared to normal cells. Cancer cells utilize glutamine for energy generation as well as the synthesis of other molecules that are critical for cancer growth and progression. Therefore, drugs targeting glutamine metabolism have been extensively investigated. However, inhibition of glutamine metabolism in cancer cells results in the activation of other metabolic pathways enabling cancer cells to survive. In this review, we summarize and discuss the targets in glutamine metabolism, which has been probed in the development of anticancer drugs in preclinical and clinical studies. We further discuss pathways activated in response to glutamine metabolism inhibition, enabling cancer cells to survive the challenge. Finally, we put into perspective combined treatment strategies targeting glutamine metabolism along with other pathways as potential treatment options.

Abstract

Tumor growth and metastasis strongly depend on adapted cell metabolism. Cancer cells adjust their metabolic program to their specific energy needs and in response to an often challenging tumor microenvironment. Glutamine metabolism is one of the metabolic pathways that can be successfully targeted in cancer treatment. The dependence of many hematological and solid tumors on glutamine is associated with mitochondrial glutaminase (GLS) activity that enables channeling of glutamine into the tricarboxylic acid (TCA) cycle, generation of ATP and NADPH, and regulation of glutathione homeostasis and reactive oxygen species (ROS). Small molecules that target glutamine metabolism through inhibition of GLS therefore simultaneously limit energy availability and increase oxidative stress. However, some cancers can reprogram their metabolism to evade this metabolic trap. Therefore, the effectiveness of treatment strategies that rely solely on glutamine inhibition is limited. In this review, we discuss the metabolic and molecular pathways that are linked to dysregulated glutamine metabolism in multiple cancer types. We further summarize and review current clinical trials of glutaminolysis inhibition in cancer patients. Finally, we put into perspective strategies that deploy a combined treatment targeting glutamine metabolism along with other molecular or metabolic pathways and discuss their potential for clinical applications.

Epigenetic scores for the circulating proteome as tools for disease prediction

Protein biomarkers have been identified across many age-related morbidities. However, characterising epigenetic influences could further inform disease predictions. Here, we leverage epigenome-wide data to study links between the DNA methylation (DNAm) signatures of the circulating proteome and incident diseases. Using data from four cohorts, we trained and tested epigenetic scores (EpiScores) for 953 plasma proteins, identifying 109 scores that explained between 1% and 58% of the variance in protein levels after adjusting for known protein quantitative trait loci (pQTL) genetic effects. By projecting these EpiScores into an independent sample (Generation Scotland; n = 9537) and relating them to incident morbidities over a follow-up of 14 years, we uncovered 137 EpiScore-disease associations. These associations were largely independent of immune cell proportions, common lifestyle and health factors, and biological aging. Notably, we found that our diabetes-associated EpiScores highlighted previous top biomarker associations from proteome-wide assessments of diabetes. These EpiScores for protein levels can therefore be a valuable resource for disease prediction and risk stratification.

TIGER: technical variation elimination for metabolomics data using ensemble learning architecture

Large metabolomics datasets inevitably contain unwanted technical variations which can obscure meaningful biological signals and affect how this information is applied to personalized healthcare. Many methods have been developed to handle unwanted variations. However, the underlying assumptions of many existing methods only hold for a few specific scenarios. Some tools remove technical variations with models trained on quality control (QC) samples which may not generalize well on subject samples. Additionally, almost none of the existing methods supports datasets with multiple types of QC samples, which greatly limits their performance and flexibility. To address these issues, a non-parametric method TIGER (Technical variation elImination with ensemble learninG architEctuRe) is developed in this study and released as an R package (https://CRAN.R-project.org/package=TIGERr). TIGER integrates the random forest algorithm into an adaptable ensemble learning architecture. Evaluation results show that TIGER outperforms four popular methods with respect to robustness and reliability on three human cohort datasets constructed with targeted or untargeted metabolomics data. Additionally, a case study aiming to identify age-associated metabolites is performed to illustrate how TIGER can be used for cross-kit adjustment in a longitudinal analysis with experimental data of three time-points generated by different analytical kits. A dynamic website is developed to help evaluate the performance of TIGER and examine the patterns revealed in our longitudinal analysis (https://han-siyu.github.io/TIGER_web/). Overall, TIGER is expected to be a powerful tool for metabolomics data analysis.

Proteomic Analysis of Plasma Markers in Patients Maintained on Antipsychotics: Comparison to Patients Off Antipsychotics and Normal Controls

BACKGROUND

Schizophrenia (SZ) and bipolar disorder (BD) share many features: overlap in mood and psychotic symptoms, common genetic predisposition, treatment with antipsychotics (APs), and similar metabolic comorbidities. The pathophysiology of both is still not well defined, and no biomarkers can be used clinically for diagnosis and management. This study aimed to assess the plasma proteomics profile of patients with SZ and BD maintained on APs compared to those who had been off APs for 6 months and to healthy controls (HCs).

METHODS

We analyzed the data using functional enrichment, random forest modeling to identify potential biomarkers, and multivariate regression for the associations with metabolic abnormalities.

RESULTS

We identified several proteins known to play roles in the differentiation of the nervous system like NTRK2, CNTN1, ROBO2, and PLXNC1, which were downregulated in AP-free SZ and BD patients but were "normalized" in those on APs. Other proteins (like NCAM1 and TNFRSF17) were "normal" in AP-free patients but downregulated in patients on APs, suggesting that these changes are related to medication's effects. We found significant enrichment of proteins involved in neuronal plasticity, mainly in SZ patients on APs. Most of the proteins associated with metabolic abnormalities were more related to APs use than having SZ or BD. The biomarkers identification showed specific and sensitive results for schizophrenia, where two proteins (PRL and MRC2) produced adequate results.

CONCLUSIONS

Our results confirmed the utility of blood samples to identify protein signatures and mechanisms involved in the pathophysiology and treatment of SZ and BD.

Proteome-wide associations with short- and long-term weight loss and regain after Roux-en-Y gastric bypass surgery

OBJECTIVE

Gastric bypass surgery results in long-term weight loss. Small studies have examined protein changes during rapid weight loss (up to 1 or 2 years post surgery). This study tested whether short-term changes were maintained after 12 years.

METHODS

A 12-year follow-up, protein-wide association study of 1,297 SomaLogic aptamer-based plasma proteins compared short- (2-year) and long-term (12-year) protein changes in 234 individuals who had gastric bypass surgery with 144 nonintervened individuals with severe obesity.

RESULTS

There were 51 replicated 12-year protein changes that differed between the surgery and nonsurgery groups. Adjusting for change in BMI, only 12 proteins remained significant, suggesting that BMI change was the primary reason for most protein changes and not non-BMI-related surgical effects. Protein changes were related to BMI changes during both weight-loss and weight-regain periods. The significant proteins were associated primarily with lipid, uric acid, or resting energy expenditure clinical variables and metabolic pathways. Eight protein changes were associated with 12-year diabetes remission, including apolipoprotein M, sex hormone binding globulin, and adiponectin (p < 3.5 × 10^-5 ).

CONCLUSIONS

This study showed that most short-term postsurgical changes in proteins were maintained at 12 years. Systemic protection pathways, including inflammation, complement, lipid, and adipocyte pathways, were related to the long-term benefits of gastric bypass surgery.

Defining the landscape of metabolic dysregulations in cancer metastasis

Metastasis is the primary cause of cancer related deaths due to the limited number of efficient druggable targets. Signatures of dysregulated cancer metabolism could serve as a roadmap for the determination of new treatment strategies. However, the metabolic signatures of metastatic cells remain vastly elusive. Our aim was to determine metabolic dysregulations associated with high metastatic potential in breast cancer cell lines. We have selected 5 triple negative breast cancer (TNBC) cell lines including three with high metastatic potential (HMP) (MDA-MB-231, MDA-MB-436, MDA-MB-468) and two with low metastatic potential (LMP) (BT549, HCC1143). The normal epithelial breast cell line (hTERT-HME1) was also investigated. The untargeted metabolic profiling of cells and growth media was conducted and total of 479 metabolites were quantified. First we characterized metabolic features differentiating TNBC cell lines from normal cells as well as identified cell line specific metabolic fingerprints. Next, we determined 92 metabolites in cells and 22 in growth medium that display significant differences between LMP and HMP. The HMP cell lines had elevated level of molecules involved in glycolysis, TCA cycle and lipid metabolism. We identified metabolic advantages of cell lines with HMP beyond enhanced glycolysis by pinpointing the role of branched chain amino acids (BCAA) catabolism as well as molecules supporting coagulation and platelet activation as important contributors to the metastatic cascade. The landscape of metabolic dysregulations, characterized in our study, could serve as a roadmap for the identification of treatment strategies targeting cancer cells with enhanced metastatic potential.

The Proteomic Signature of Recombinant Growth Hormone in Recreational Athletes

OBJECTIVE

Administration of human growth hormone (hGH) is prohibited in competitive sport and its detection in an athlete's sample triggers an adverse analytical finding. However, the biological processes that are modulated by recombinant hGH are not well characterized and associated blood serum proteins may constitute new biomarkers for hGH misuse.

METHODS

Thirty-five recreational athletes were enrolled in a study to investigate the time- and dose-dependent response of serum protein levels to recombinant hGH administration. Participants were randomly assigned to 4 groups, receiving 1 of 3 different doses of recombinant hGH or a placebo. Bio samples were collected at 22 time points over a period of 13 weeks, starting 4 weeks before treatment, during 3 weeks of treatment, and at 6 weeks' follow-up. A total of 749 serum samples were analyzed for 1305 protein markers using the SOMAscan proteomics platform.

RESULTS

We identified 66 proteins that significantly associated with recombinant hGH administration and dosage, including well known hGH targets, such as IGF1, but also previously unknown hGH-related proteins (eg, protease inhibitors, WFIKKN1, and chemokines, CCL2). Network analysis revealed changes in specific biological pathways, mainly related to the immune system and glucose metabolism.

CONCLUSION

Our analysis suggests that hGH administration affects biological processes more strongly than previously acknowledged. Some of the proteins were dysregulated even after hGH treatment and could potentially be developed into biomarkers for hGH misuse. Moreover, our findings suggest new roles for hGH-associated proteins in the etiology of hGH-related diseases and may indicate new risks that may be associated with hGH misuse.