A genome-wide association study of mass spectrometry proteomics using the Seer Proteograph platform.

Genome-wide association studies (GWAS) with proteomics are essential tools for drug discovery. To date, most studies have used affinity proteomics platforms, which have limited discovery to protein panels covered by the available affinity binders. Furthermore, it is not clear to which extent protein epitope changing variants interfere with the detection of protein quantitative trait loci (pQTLs). Mass spectrometry-based (MS) proteomics can overcome some of these limitations. Here we report a GWAS using the MS-based Seer Proteograph™ platform with blood samples from a discovery cohort of 1,260 American participants and a replication in 325 individuals from Asia, with diverse ethnic backgrounds. We analysed 1,980 proteins quantified in at least 80% of the samples, out of 5,753 proteins quantified across the discovery cohort. We identified 252 and replicated 90 pQTLs, where 30 of the replicated pQTLs have not been reported before. We further investigated 200 of the strongest associated cis-pQTLs previously identified using the SOMAscan and the Olink platforms and found that up to one third of the affinity proteomics pQTLs may be affected by epitope effects, while another third were confirmed by MS proteomics to be consistent with the hypothesis that genetic variants induce changes in protein expression. The present study demonstrates the complementarity of the different proteomics approaches and reports pQTLs not accessible to affinity proteomics, suggesting that many more pQTLs remain to be discovered using MS-based platforms.

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.

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.

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.

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.

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.

piTracer - Automatic reconstruction of molecular cascades for the identification of synergistic drug targets.

Cancer cells frequently undergo metabolic reprogramming as a mechanism of resistance against chemotherapeutic drugs. Metabolomic profiling provides a direct readout of metabolic changes and can thus be used to identify these tumor escape mechanisms. Here, we introduce piTracer, a computational tool that uses multi-scale molecular networks to identify potential combination therapies from pre- and post-treatment metabolomics data. We first demonstrate piTracer’s core ability to reconstruct cellular cascades by inspecting well-characterized molecular pathways and previously studied associations between genetic variants and metabolite levels. We then apply a new gene ranking algorithm on differential metabolomic profiles from human breast cancer cells after glutaminase inhibition. Four of the automatically identified gene targets were experimentally tested by simultaneous inhibition of the respective targets and glutaminase. Of these combination treatments, two were be confirmed to induce synthetic lethality in the cell line. In summary, piTracer integrates the molecular monitoring of escape mechanisms into comprehensive pathway networks to accelerate drug target identification. The tool is open source and can be accessed at https://github.com/krumsieklab/pitracer .

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.

Ratios of Acetaminophen Metabolites Identify New Loci of Pharmacogenetic Relevance in a Genome-Wide Association Study.

Genome-wide association studies (GWAS) with non-targeted metabolomics have identified many genetic loci of biomedical interest. However, metabolites with a high degree of missingness, such as drug metabolites and xenobiotics, are often excluded from such studies due to a lack of statistical power and higher uncertainty in their quantification. Here we propose ratios between related drug metabolites as GWAS phenotypes that can drastically increase power to detect genetic associations between pairs of biochemically related molecules. As a proof-of-concept we conducted a GWAS with 520 individuals from the Qatar Biobank for who at least five of the nine available acetaminophen metabolites have been detected. We identified compelling evidence for genetic variance in acetaminophen glucuronidation and methylation by UGT2A15 and COMT, respectively. Based on the metabolite ratio association profiles of these two loci we hypothesized the chemical structure of one of their products or substrates as being 3-methoxyacetaminophen, which we then confirmed experimentally. Taken together, our study suggests a novel approach to analyze metabolites with a high degree of missingness in a GWAS setting with ratios, and it also demonstrates how pharmacological pathways can be mapped out using non-targeted metabolomics measurements in large population-based studies.

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.

Advancing Cancer Treatment by Targeting Glutamine Metabolism-A Roadmap.

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.

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.

Deep sequencing of DNA from urine of kidney allograft recipients to estimate donor/recipient-specific DNA fractions.

Kidney transplantation is the treatment of choice for patients with end-stage kidney failure, but transplanted allograft could be affected by viral and bacterial infections and by immune rejection. The standard test for the diagnosis of acute pathologies in kidney transplants is kidney biopsy. However, noninvasive tests would be desirable. Various methods using different techniques have been developed by the transplantation community. But these methods require improvements. We present here a cost-effective method for kidney rejection diagnosis that estimates donor/recipient-specific DNA fraction in recipient urine by sequencing urinary cell DNA. We hypothesized that in the no-pathology stage, the largest tissue types present in recipient urine are donor kidney cells, and in case of rejection, a larger number of recipient immune cells would be observed. Extensive in-silico simulation was used to tune the sequencing parameters: number of variants and depth of coverage. Sequencing of DNA mixture from 2 healthy individuals showed the method is highly predictive (maximum error < 0.04). We then demonstrated the insignificant impact of familial relationship and ethnicity using an in-house and public database. Lastly, we performed deep DNA sequencing of urinary cell pellets from 32 biopsy-matched samples representing two pathology groups: acute rejection (AR, 11 samples) and acute tubular injury (ATI, 12 samples) and 9 samples with no pathology. We found a significant association between the donor/recipient-specific DNA fraction in the two pathology groups compared to no pathology (P = 0.0064 for AR and P = 0.026 for ATI). We conclude that deep DNA sequencing of urinary cells from kidney allograft recipients offers a noninvasive means of diagnosing acute pathologies in the human kidney allograft.

Metabolic Predictors of Equine Performance in Endurance Racing.

Equine performance in endurance racing depends on the interplay between physiological and metabolic processes. However, there is currently no parameter for estimating the readiness of animals for competition. Our objectives were to provide an in-depth characterization of metabolic consequences of endurance racing and to establish a metabolic performance profile for those animals. We monitored metabolite composition, using a broad non-targeted metabolomics approach, in blood plasma samples from 47 Arabian horses participating in endurance races. The samples were collected before and after the competition and a total of 792 metabolites were measured. We found significant alterations between before and after the race in 417 molecules involved in lipids and amino acid metabolism. Further, even before the race starts, we found metabolic differences between animals who completed the race and those who did not. We identified a set of six metabolite predictors (imidazole propionate, pipecolate, ethylmalonate, 2R-3R-dihydroxybutyrate, ß-hydroxy-isovalerate and X-25455) of animal performance in endurance competition; the resulting model had an area under a receiver operating characteristic (AUC) of 0.92 (95% CI: 0.85-0.98). This study provides an in-depth characterization of metabolic alterations driven by endurance races in equines. Furthermore, we showed the feasibility of identifying potential metabolic signatures as predictors of animal performance in endurance competition.

Association of microRNAs With Embryo Development and Fertilization in Women Undergoing Subfertility Treatments: A Pilot Study.

Objective: Small non-coding RNAs, known as microRNAs (miRNAs), have emerging regulatory functions within the ovary that have been related to fertility. This study was undertaken to determine if circulating miRNAs reflect the changes associated with the parameters of embryo development and fertilization. Methods: In this cross-sectional pilot study. Plasma miRNAs were collected from 48 sequentially presenting women in the follicular phase prior to commencing in vitro fertilization (IVF). Circulating miRNAs were measured using locked nucleic acid (LNA)-based quantitative PCR (qPCR), while an updated miRNA data set was used to determine their level of expression. Results: Body mass index and weight were associated with the miRNAs let7b-3p and miR-375, respectively (p < 0.05), with the same relationship being found between endometrium thickness at oocyte retrieval and miR-885-5p and miR-34a-5p (p < 0.05). In contrast, miR-1260a was found to be inversely associated with anti-Mullerian hormone (AMH; p = 0.007), while miR-365a-3p, miR122-5p, and miR-34a-5p correlated with embryo fertilization rates (p < 0.05). However, when omitting cases of male infertility (n = 15), miR122-5p remained significant (p < 0.05), while miR-365a-3p and miR-34a-5p no longer differed; interestingly, however, miR1260a and mir93.3p became significant (p = 0.0087/0.02, respectively). Furthermore, age was negatively associated with miR-335-3p, miR-28-5p, miR-155-5p, miR-501-3p, and miR-497-5p (p < 0.05). Live birth rate was negatively associated with miR-335-3p, miR-100-5p, miR-497-5p, let-7d, and miR-574-3p (p < 0.05), but these were not significant when age was accounted for.However, with the exclusion of male factor infertility, all those miRNAs were no longer significant, though miR.150.5p emerged as significant (p = 0.042). A beta-regression model identified miR-1260a, miR-486-5p, and miR-132-3p (p < 0.03, p = 0.0003, p < 0.00001, respectively) as the most predictive for fertilization rate. Notably, changes in detectable miRNAs were not linked to cleavage rate, top quality embryos (G3D3), and blastocyst or antral follicle count. An ingenuity pathway analysis showed that miRNAs associated with age were also associated with the variables found in reproductive system diseases. Conclusion: Plasma miRNAs prior to the IVF cycle were associated with differing demographic and IVF parameters, including age, and may be predictive biomarkers of fertilization rate.

The metabolic footprint of compromised insulin sensitivity under fasting and hyperinsulinemic-euglycemic clamp conditions in an Arab population.

Metabolic pathways that are corrupted at early stages of insulin resistance (IR) remain elusive. This study investigates changes in body metabolism in clinically healthy and otherwise asymptomatic subjects that may become apparent already under compromised insulin sensitivity (IS) and prior to IR. 47 clinically healthy Arab male subjects with a broad range of IS, determined by hyperinsulinemic-euglycemic clamp (HIEC), were investigated. Untargeted metabolomics and complex lipidomics were conducted on serum samples collected under fasting and HIEC conditions. Linear models were used to identify associations between metabolites concentrations and IS levels. Among 1896 identified metabolites, 551 showed significant differences between fasting and HIEC, reflecting the metabolic switch in energy utilization. At fasting, 336 metabolites, predominantly di- and tri-acylglycerols, showed significant differences between subjects with low and high levels of IS. Changes in amino acid, carbohydrate and fatty acid metabolism in response to insulin were impaired in subjects with low IS. Association of altered mannose and amino acids with IS was also replicated in an independent cohort of T2D patients. We identified metabolic phenotypes that characterize clinically healthy Arab subjects with low levels of IS at their fasting state. Our study is providing further insights into the metabolic pathways that precede IR.

Identification of genetic variants controlling RNA editing and their effect on RNA structure stabilization.

Post-transcriptional modification of RNA (RNA editing, RNAe) results in differences between the RNA transcript and the genomic DNA sequence (RDD). Enzymatic modification of adenosine to inosine (A2I) by ADAR is the most studied type of RNAe. However, few genetic association studies with A2I RNAe events have been conducted. Some studies have analyzed the inter-population RNAe-QTL diversity in humans, but the sample size of these studies was limited. Other types of RNA and DNA differences have been reported but are largely understudied. Here, we report a comprehensive analysis of all types of RDD, based on two independent datasets. We found that A2I was by far the most observed type of RDD. Moreover, manual curation suggests that A2I is likely the only enzymatically driven RNAe type observed in blood derived DNA, all other non-A2I RDD could either be attributed to sequencing and processing artifacts, or are a result of somatic DNA rearrangements. We then conducted an in-cis genetic association study and identified 472 genetic associations (RNAe-QTL), that were replicated in both datasets. We confirm the potential effect of the RNAe-QTL on RNA structure by showing that allele specific RNAe occurs in heterozygotes. Although the generally assumed function of RNAe is to destabilize double stranded RNA structure, we found clear evidence for the potential additional involvement of RNAe in maintaining RNA hairpin that has been altered by the RNAe-QTL. Our study confirms, in two independent datasets, the potential role of RNAe in maintaining RNA structure in the presence of genetic variation.

Metabolic Signatures of Tumor Responses to Doxorubicin Elucidated by Metabolic Profiling in Ovo.

BACKGROUND: Dysregulated cancer metabolism is associated with acquired resistance to chemotherapeutic treatment and contributes to the activation of cancer survival mechanisms. However, which metabolic pathways are activated following treatment often remains elusive. The combination of chicken embryo tumor models (in ovo) with metabolomics phenotyping could offer a robust platform for drug testing. Here, we assess the potential of this approach in the treatment of an in ovo triple negative breast cancer with doxorubicin. METHODS: MB-MDA-231 cells were grafted in ovo. The resulting tumors were then treated with doxorubicin or dimethyl sulfoxide (DMSO) for six days. Tumors were collected and analyzed using a global untargeted metabolomics and comprehensive lipidomics. RESULTS: We observed a significant suppression of tumor growth in the doxorubicin treated group. The metabolic profiles of doxorubicin and DMSO-treated tumors were clearly separated in a principle component analysis. Inhibition of glycolysis, nucleotide synthesis, and glycerophospholipid metabolism appear to be triggered by doxorubicin treatment, which could explain the observed suppressed tumor growth. In addition, metabolic cancer survival mechanisms could be supported by an acceleration of antioxidative pathways. CONCLUSIONS: Metabolomics in combination with in ovo tumor models provide a robust platform for drug testing to reveal tumor specific treatment targets such as the antioxidative tumor capacity.

Author Correction: Effect of induced hypoglycemia on inflammation and oxidative stress in type 2 diabetes and control subjects.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.