A survey of ficolin-3 activity in Systemic Lupus Erythematosus reveals a link to hematological disease manifestations and autoantibody profile.

The complement system plays a central role in the pathogenesis of Systemic Lupus Erythematosus (SLE), but most studies have focused on the classical pathway. Ficolin-3 is the main initiator of the lectin pathway of complement in humans, but its role in systemic autoimmune disease has not been conclusively determined. Here, we combined biochemical and genetic approaches to assess the contribution of ficolin-3 to SLE risk and disease manifestations. Ficolin-3 activity was measured by a functional assay in serum or plasma samples from Swedish SLE patients (n = 786) and controls matched for age and sex (n = 566). Genetic variants in an extended 300 kb genomic region spanning the FCN3 locus were analyzed for their association with ficolin-3 activity and SLE manifestations in a Swedish multicenter cohort (n = 985). Patients with ficolin-3 activity in the highest tertile showed a strong enrichment in an SLE cluster defined by anti-Sm/DNA/nucleosome antibodies (OR 3.0, p < 0.001) and had increased rates of hematological disease (OR 1.4, p = 0.078) and lymphopenia (OR = 1.6, p = 0.039). Genetic variants associated with low ficolin-3 activity mapped to an extended haplotype in high linkage disequilibrium upstream of the FCN3 gene. Patients carrying the lead genetic variant associated with low ficolin-3 activity had a lower frequency of hematological disease (OR 0.67, p = 0.018) and lymphopenia (OR 0.63, p = 0.031) and fewer autoantibodies (p = 0.0019). Loss-of-function variants in the FCN3 gene were not associated with SLE, but four (0.5 %) SLE patients developed acquired ficolin-3 deficiency where ficolin-3 activity in serum was depleted following diagnosis of SLE. Taken together, our results provide genetic and biochemical evidence that implicate the lectin pathway in hematological SLE manifestations. We also identify lectin pathway activation through ficolin-3 as a factor that contributes to the autoantibody response in SLE.

Mitochondrial heteroplasmic shifts reveal a positive selection of breast cancer.

BACKGROUND: Breast cancer is, despite screening, not always detected early enough and is together with other tumor types known to shed genetic information in circulation. Unlike single-copy nuclear DNA, mitochondrial DNA (mtDNA) copies range from 100s to 10,000s per cell, thus providing a potentially alternative to identify potential missing cancer information in circulation at an early stage. METHODS: To characterize mitochondrial mutation landscapes in breast cancer, whole mtDNA sequencing and bioinformatics analyses were performed on 86 breast cancer biopsies and 50 available matched baseline cancer-free whole blood samples from the same individuals, selected from a cohort of middle-aged women in Sweden. To determine whether the mutations can be detected in blood plasma prior to cancer diagnosis, we further designed a nested case-control study (n = 663) and validated the shortlisted mutations using droplet digital PCR. RESULTS: We detected different mutation landscapes between biopsies and matched whole blood samples. Compared to whole blood samples, mtDNA from biopsies had higher heteroplasmic mutations in the D-loop region (P = 0.02), RNR2 (P = 0.005), COX1 (P = 0.037) and CYTB (P = 0.006). Furthermore, the germline mtDNA mutations had higher heteroplasmy level than the lost (P = 0.002) and de novo mutations (P = 0.04). The nonsynonymous to synonymous substitution ratio (dN/dS) was higher for the heteroplasmic mutations (P = 7.25 × 10-12) than that for the homoplasmic mutations, but the de novo (P = 0.06) and lost mutations (P = 0.03) had lower dN/dS than the germline mutations. Interestingly, we found that the critical regions for mitochondrial transcription: MT-HSP1 (odds ratio [OR]: 21.41), MT-TFH (OR: 7.70) and MT-TAS2 (OR: 3.62), had significantly higher heteroplasmic mutations than the rest of the D-loop sub-regions. Finally, we found that the presence of mt.16093T > C mutation increases 67% risk of developing breast cancer. CONCLUSIONS: Our findings show that mitochondrial genetic landscape changes during cancer pathogenesis and positive selection of mtDNA heteroplasmic mutations in breast cancer. Most importantly, the mitochondrial mutations identified in biopsies can be traced back in matched plasma samples and could potentially be used as early breast cancer diagnostic biomarkers.

Identification of Endothelial Proteins in Plasma Associated With Cardiovascular Risk Factors.

OBJECTIVE: Endothelial cell (EC) dysfunction is a well-established response to cardiovascular disease risk factors, such as smoking and obesity. Risk factor exposure can modify EC signaling and behavior, leading to arterial and venous disease development. Here, we aimed to identify biomarker panels for the assessment of EC dysfunction, which could be useful for risk stratification or to monitor treatment response. Approach and Results: We used affinity proteomics to identify EC proteins circulating in plasma that were associated with cardiovascular disease risk factor exposure. Two hundred sixteen proteins, which we previously predicted to be EC-enriched across vascular beds, were measured in plasma samples (N=1005) from the population-based SCAPIS (Swedish Cardiopulmonary Bioimage Study) pilot. Thirty-eight of these proteins were associated with body mass index, total cholesterol, low-density lipoprotein, smoking, hypertension, or diabetes. Sex-specific analysis revealed that associations predominantly observed in female- or male-only samples were most frequently with the risk factors body mass index, or total cholesterol and smoking, respectively. We show a relationship between individual cardiovascular disease risk, calculated with the Framingham risk score, and the corresponding biomarker profiles. CONCLUSIONS: EC proteins in plasma could reflect vascular health status.

Circulating proteins associated with allergy development in infants-an exploratory analysis.

BACKGROUND: Protein profiles that can predict allergy development in children are lacking and the ideal sampling age is unknown. By applying an exploratory proteomics approach in the prospective FARMFLORA birth cohort, we sought to identify previously unknown circulating proteins in early life that associate to protection or risk for development of allergy up to 8 years of age. METHODS: We analyzed plasma prepared from umbilical cord blood (n = 38) and blood collected at 1 month (n = 42), 4 months (n = 39), 18 months (n = 42), 36 months (n = 42) and 8 years (n = 44) of age. We profiled 230 proteins with a multiplexed assay and evaluated the global structure of the data with principal component analysis (PCA). Protein profiles informative to allergic disease at 18 months, 36 months and/or 8 years were evaluated using Lasso logistic regression and random forest. RESULTS: Two clusters emerged in the PCA analysis that separated samples obtained at birth and at 1 month of age from samples obtained later. Differences between the clusters were mostly driven by abundant plasma proteins. For the prediction of allergy, both Lasso logistic regression and random forest were most informative with samples collected at 1 month of age. A Lasso model with 27 proteins together with farm environment differentiated children who remained healthy from those developing allergy. This protein panel was primarily composed of antigen-presenting MHC class I molecules, interleukins and chemokines. CONCLUSION: Sampled at one month of age, circulating proteins that reflect processes of the immune system may predict the development of allergic disease later in childhood.

Predicting and elucidating the etiology of fatty liver disease: A machine learning modeling and validation study in the IMI DIRECT cohorts.

BACKGROUND: Non-alcoholic fatty liver disease (NAFLD) is highly prevalent and causes serious health complications in individuals with and without type 2 diabetes (T2D). Early diagnosis of NAFLD is important, as this can help prevent irreversible damage to the liver and, ultimately, hepatocellular carcinomas. We sought to expand etiological understanding and develop a diagnostic tool for NAFLD using machine learning. METHODS AND FINDINGS: We utilized the baseline data from IMI DIRECT, a multicenter prospective cohort study of 3,029 European-ancestry adults recently diagnosed with T2D (n = 795) or at high risk of developing the disease (n = 2,234). Multi-omics (genetic, transcriptomic, proteomic, and metabolomic) and clinical (liver enzymes and other serological biomarkers, anthropometry, measures of beta-cell function, insulin sensitivity, and lifestyle) data comprised the key input variables. The models were trained on MRI-image-derived liver fat content (<5% or ≥5%) available for 1,514 participants. We applied LASSO (least absolute shrinkage and selection operator) to select features from the different layers of omics data and random forest analysis to develop the models. The prediction models included clinical and omics variables separately or in combination. A model including all omics and clinical variables yielded a cross-validated receiver operating characteristic area under the curve (ROCAUC) of 0.84 (95% CI 0.82, 0.86; p < 0.001), which compared with a ROCAUC of 0.82 (95% CI 0.81, 0.83; p < 0.001) for a model including 9 clinically accessible variables. The IMI DIRECT prediction models outperformed existing noninvasive NAFLD prediction tools. One limitation is that these analyses were performed in adults of European ancestry residing in northern Europe, and it is unknown how well these findings will translate to people of other ancestries and exposed to environmental risk factors that differ from those of the present cohort. Another key limitation of this study is that the prediction was done on a binary outcome of liver fat quantity (<5% or ≥5%) rather than a continuous one. CONCLUSIONS: In this study, we developed several models with different combinations of clinical and omics data and identified biological features that appear to be associated with liver fat accumulation. In general, the clinical variables showed better prediction ability than the complex omics variables. However, the combination of omics and clinical variables yielded the highest accuracy. We have incorporated the developed clinical models into a web interface (see: https://www.predictliverfat.org/) and made it available to the community. TRIAL REGISTRATION: ClinicalTrials.gov NCT03814915.

Affinity proteomic profiling of plasma for proteins associated to area-based mammographic breast density.

BACKGROUND: Mammographic breast density is one of the strongest risk factors for breast cancer, but molecular understanding of how breast density relates to cancer risk is less complete. Studies of proteins in blood plasma, possibly associated with mammographic density, are well-suited as these allow large-scale analyses and might shed light on the association between breast cancer and breast density. METHODS: Plasma samples from 1329 women in the Swedish KARMA project, without prior history of breast cancer, were profiled with antibody suspension bead array (SBA) assays. Two sample sets comprising 729 and 600 women were screened by two different SBAs targeting a total number of 357 proteins. Protein targets were selected through searching the literature, for either being related to breast cancer or for being linked to the extracellular matrix. Association between proteins and absolute area-based breast density (AD) was assessed by quantile regression, adjusting for age and body mass index (BMI). RESULTS: Plasma profiling revealed linear association between 20 proteins and AD, concordant in the two sets of samples (p < 0.05). Plasma levels of seven proteins were positively associated and 13 proteins negatively associated with AD. For eleven of these proteins evidence for gene expression in breast tissue existed. Among these, ABCC11, TNFRSF10D, F11R and ERRF were positively associated with AD, and SHC1, CFLAR, ACOX2, ITGB6, RASSF1, FANCD2 and IRX5 were negatively associated with AD. CONCLUSIONS: Screening proteins in plasma indicates associations between breast density and processes of tissue homeostasis, DNA repair, cancer development and/or progression in breast cancer. Further validation and follow-up studies of the shortlisted protein candidates in independent cohorts will be needed to infer their role in breast density and its progression in premenopausal and postmenopausal women.

PDGFB, a new candidate plasma biomarker for venous thromboembolism: results from the VEREMA affinity proteomics study.

There is a clear clinical need for high-specificity plasma biomarkers for predicting risk of venous thromboembolism (VTE), but thus far, such markers have remained elusive. Utilizing affinity reagents from the Human Protein Atlas project and multiplexed immuoassays, we extensively analyzed plasma samples from 2 individual studies to identify candidate protein markers associated with VTE risk. We screened plasma samples from 88 VTE cases and 85 matched controls, collected as part of the Swedish "Venous Thromboembolism Biomarker Study," using suspension bead arrays composed of 755 antibodies targeting 408 candidate proteins. We identified significant associations between VTE occurrence and plasma levels of human immunodeficiency virus type I enhancer binding protein 1 (HIVEP1), von Willebrand factor (VWF), glutathione peroxidase 3 (GPX3), and platelet-derived growth factor ß (PDGFB). For replication, we profiled plasma samples of 580 cases and 589 controls from the French FARIVE study. These results confirmed the association of VWF and PDGFB with VTE after correction for multiple testing, whereas only weak trends were observed for HIVEP1 and GPX3. Although plasma levels of VWF and PDGFB correlated modestly (ρ ∼ 0.30) with each other, they were independently associated with VTE risk in a joint model in FARIVE (VWF P < .001; PDGFB P = .002). PDGFΒ was verified as the target of the capture antibody by immunocapture mass spectrometry and sandwich enzyme-linked immunosorbent assay. In conclusion, we demonstrate that high-throughput affinity plasma proteomic profiling is a valuable research strategy to identify potential candidate biomarkers for thrombosis-related disorders, and our study suggests a novel association of PDGFB plasma levels with VTE.

Identification of a Novel Autoimmune Peptide Epitope of Prostein in Prostate Cancer.

There is a demand for novel targets and approaches to diagnose and treat prostate cancer (PCA). In this context, serum and plasma samples from a total of 609 individuals from two independent patient cohorts were screened for IgG reactivity against a sum of 3833 human protein fragments. Starting from planar protein arrays with 3786 protein fragments to screen 80 patients with and without PCA diagnosis, 161 fragments (4%) were chosen for further analysis based on their reactivity profiles. Adding 71 antigens from literature, the selection of antigens was corroborated for their reactivity in a set of 550 samples using suspension bead arrays. The antigens prostein (SLC45A3), TATA-box binding protein (TBP), and insulin-like growth factor 2 mRNA binding protein 2 (IGF2BP2) showed higher reactivity in PCA patients with late disease compared with early disease. Because of its prostate tissue specificity, we focused on prostein and continued with mapping epitopes of the 66-mer protein fragment using patient samples. Using bead-based assays and 15-mer peptides, a minimal peptide epitope was identified and refined by alanine scanning to the KPxAPFP. Further sequence alignment of this motif revealed homology to transmembrane protein 79 (TMEM79) and TGF-beta-induced factor 2 (TGIF2), thus providing a reasoning for cross-reactivity found in females. A comprehensive workflow to discover and validate IgG reactivity against prostein and homologous targets in human serum and plasma was applied. This study provides useful information when searching for novel biomarkers or drug targets that are guided by the reactivity of the immune system against autoantigens.

Multidimensional Normalization to Minimize Plate Effects of Suspension Bead Array Data.

Enhanced by the growing number of biobanks, biomarker studies can now be performed with reasonable statistical power by using large sets of samples. Antibody-based proteomics by means of suspension bead arrays offers one attractive approach to analyze serum, plasma, or CSF samples for such studies in microtiter plates. To expand measurements beyond single batches, with either 96 or 384 samples per plate, suitable normalization methods are required to minimize the variation between plates. Here we propose two normalization approaches utilizing MA coordinates. The multidimensional MA (multi-MA) and MA-loess both consider all samples of a microtiter plate per suspension bead array assay and thus do not require any external reference samples. We demonstrate the performance of the two MA normalization methods with data obtained from the analysis of 384 samples including both serum and plasma. Samples were randomized across 96-well sample plates, processed, and analyzed in assay plates, respectively. Using principal component analysis (PCA), we could show that plate-wise clusters found in the first two components were eliminated by multi-MA normalization as compared with other normalization methods. Furthermore, we studied the correlation profiles between random pairs of antibodies and found that both MA normalization methods substantially reduced the inflated correlation introduced by plate effects. Normalization approaches using multi-MA and MA-loess minimized batch effects arising from the analysis of several assay plates with antibody suspension bead arrays. In a simulated biomarker study, multi-MA restored associations lost due to plate effects. Our normalization approaches, which are available as R package MDimNormn, could also be useful in studies using other types of high-throughput assay data.

Immunocapture strategies in translational proteomics.

Aiming at clinical studies of human diseases, antibody-assisted assays have been applied to biomarker discovery and toward a streamlined translation from patient profiling to assays supporting personalized treatments. In recent years, integrated strategies to couple and combine antibodies with mass spectrometry-based proteomic efforts have emerged, allowing for novel possibilities in basic and clinical research. Described in this review are some of the field's current and emerging immunocapture approaches from an affinity proteomics perspective. Discussed are some of their advantages, pitfalls and opportunities for the next phase in clinical and translational proteomics.

Affinity proteomic profiling of plasma, cerebrospinal fluid, and brain tissue within multiple sclerosis.

The brain is a vital organ and because it is well shielded from the outside environment, possibilities for noninvasive analysis are often limited. Instead, fluids taken from the spinal cord or circulatory system are preferred sources for the discovery of candidate markers within neurological diseases. In the context of multiple sclerosis (MS), we applied an affinity proteomic strategy and screened 22 plasma samples with 4595 antibodies (3450 genes) on bead arrays, then defined 375 antibodies (334 genes) for targeted analysis in a set of 172 samples and finally used 101 antibodies (43 genes) on 443 plasma as well as 573 cerebrospinal spinal fluid (CSF) samples. This revealed alteration of protein profiles in relation to MS subtypes for IRF8, IL7, METTL14, SLC30A7, and GAP43. Respective antibodies were subsequently used for immunofluorescence on human post-mortem brain tissue with MS pathology for expression and association analysis. There, antibodies for IRF8, IL7, and METTL14 stained neurons in proximity of lesions, which highlighted these candidate protein targets for further studies within MS and brain tissue. The affinity proteomic translation of profiles discovered by profiling human body fluids and tissue provides a powerful strategy to suggest additional candidates to studies of neurological disorders.

HIV-2 mediated effects on target and bystander cells induce plasma proteome remodeling.

Despite low or undetectable plasma viral load, people living with HIV-2 (PLWH2) typically progress toward AIDS. The driving forces behind HIV-2 disease progression and the role of viremia are still not known, but low-level replication in tissues is believed to play a role. To investigate the impact of viremic and aviremic HIV-2 infection on target and bystander cell pathology, we used data-independent acquisition mass spectrometry to determine plasma signatures of tissue and cell type engagement. Proteins derived from target and bystander cells in multiple tissues, such as the gastrointestinal tract and brain, were detected at elevated levels in plasma of PLWH2, compared with HIV negative controls. Moreover, viremic HIV-2 infection appeared to induce enhanced release of proteins from a broader range of tissues compared to aviremic HIV-2 infection. This study expands the knowledge on the link between plasma proteome remodeling and the pathological cell engagement in tissues during HIV-2 infection.

A genome-wide assessment of variability in human serum metabolism.

The study of the genetic regulation of metabolism in human serum samples can contribute to a better understanding of the intermediate biological steps that lead from polymorphism to disease. Here, we conducted a genome-wide association study (GWAS) to discover metabolic quantitative trait loci (mQTLs) utilizing samples from a study of prostate cancer in Swedish men, consisting of 402 individuals (214 cases and 188 controls) in a discovery set and 489 case-only samples in a replication set. A global nontargeted metabolite profiling approach was utilized resulting in the detection of 6,138 molecular features followed by targeted identification of associated metabolites. Seven replicating loci were identified (PYROXD2, FADS1, PON1, CYP4F2, UGT1A8, ACADL, and LIPC) with associated sequence variants contributing significantly to trait variance for one or more metabolites (P = 10(-13) -10(-91)). Regional mQTL enrichment analyses implicated two loci that included FADS1 and a novel locus near PDGFC. Biological pathway analysis implicated ACADM, ACADS, ACAD8, ACAD10, ACAD11, and ACOXL, reflecting significant enrichment of genes with acyl-CoA dehydrogenase activity. mQTL SNPs and mQTL-harboring genes were over-represented across GWASs conducted to date, suggesting that these data may have utility in tracing the molecular basis of some complex disease associations.

Characterization of the Mitochondrial Genetic Landscape in Abdominal Aortic Aneurysm.

Background Abdominal aortic aneurysm (AAA) is a vascular disease with a mortality rate of >80% if ruptured. Mitochondrial dysfunction has been previously implicated in AAA pathogenesis. In this study, we aimed to characterize the mitochondrial genetic landscape in AAA. Methods and Results Whole mitochondrial genome sequencing and bioinformatics analysis were performed in comorbidity matched 48 cases without AAA and 48 cases with AAA, objectively diagnosed, and selected from a cohort of 65-year-old men recruited for a screening program. We identified differential mutational landscapes in men with and without AAA, with errors in mitochondrial DNA replication or repair as potential sources. Heteroplasmic insertions and overall heteroplasmy of structural rearrangements were significantly elevated in AAA cases. Three heteroplasmic variants were associated with risk factors of AAA: leukocyte concentration, plasma glucose, and cholesterol levels, respectively. Interestingly, mutations were more prevalent in regulatory part of the mitochondria, the displacement loop region, in AAA as compared with controls (P value <0.05), especially in the conserved and critical mitochondrial extended termination-associated sequence region. Moreover, we report a novel 24 bp mitochondrial DNA duplication present exclusively in cases with AAA (4%) and 75% of the unmatched AAA biopsies. Finally, the haplogroup cluster JTU was overrepresented in AAA and significantly associated with a positive family history of AAA (odds ratio, 2.9 [95% CI, 1.1-8.1]). Conclusions This is the first study investigating the mitochondrial genome in AAA, where important genetic alterations and haplogroups associated with AAA and clinical risk factors were identified. Our findings have the potential to fill in gaps in the missing genetic information on AAA.

Soluble CD14 and Osteoprotegerin Associate with Ankle-Brachial Index as a Measure of Arterial Stiffness in Patients with Mild-to-Moderate Chronic Kidney Disease in a Five-Year Prospective Study.

INTRODUCTION: Vascular lesions and arterial stiffness appear at early stages of chronic kidney disease (CKD) and follow an accelerated course with disease progression, contributing to high cardiovascular mortality. There are limited prospective data on mechanisms contributing to progression of arterial stiffness in mild-to-moderate CKD (stages 2-3). METHODS: We applied an affinity proteomics approach to identify candidates of circulating biomarkers with potential impact on vascular lesions in CKD and selected soluble cluster of differentiation 14 (sCD14), angiogenin (ANG), and osteoprotegerin (OPG) for further analysis. We studied their association with ankle-brachial index (ABI) and carotid intima-media thickness, as measures of arteriosclerosis and atherosclerosis, respectively, in 48 patients with CKD stages 2-3, who were prospectively followed and intensively treated for 5 years, and 44 healthy controls. RESULTS: Concentrations of sCD14 (p < 0.001), ANG (p < 0.001), and OPG (p < 0.05) were higher in patients with CKD 2-3 at baseline, and sCD14 (p < 0.001) and ANG (p < 0.001) remained elevated in CKD patients at follow-up. There were positive correlations between ABI and sCD14 levels (r = 0.36, p = 0.01) and between ABI and OPG (r = 0.31, p = 0.03) at 5 years. The changes in sCD14 during follow-up correlated to changes in ABI from baseline to 5 years (r = 0.41, p = 0.004). CONCLUSION: Elevated levels of circulating sCD14 and OPG in patients with CKD 2-3 were significantly associated with ABI, a measure of arterial stiffness. An increase in sCD14 over time in CKD 2-3 patients was associated with a corresponding increase in ABI. Further studies are needed to examine if early intensive multifactorial medication to align with international treatment targets may influence cardiovascular outcomes.

Screening for Circulating Inflammatory Proteins Does Not Reveal Plasma Biomarkers of Constant Tinnitus.

BACKGROUND AND OBJECTIVE: Tinnitus would benefit from an objective biomarker. The goal of this study is to identify plasma biomarkers of constant and chronic tinnitus among selected circulating inflammatory proteins. METHODS: A case-control retrospective study on 548 cases with constant tinnitus and 548 matched controls from the Swedish Tinnitus Outreach Project (STOP), whose plasma samples were examined using Olink's Inflammatory panel. Replication and meta-analysis were performed using the same method on samples from the TwinsUK cohort. Participants from LifeGene, whose blood was collected in Stockholm and Umeå, were recruited to STOP for a tinnitus subtyping study. An age and sex matching was performed at the individual level. TwinsUK participants (n = 928) were selected based on self-reported tinnitus status over 2 to 10 years. Primary outcomes include normalized levels for 96 circulating proteins, which were used as an index test. No reference standard was available in this study. RESULTS: After adjustment for age, sex, BMI, smoking, hearing loss, and laboratory site, the top proteins identified were FGF-21, MCP4, GDNF, CXCL9, and MCP-1; however, these were no longer statistically significant after correction for multiple testing. Stratification by sex did not yield any significant associations. Similarly, associations with hearing loss or other tinnitus-related comorbidities such as stress, anxiety, depression, hyperacusis, temporomandibular joint disorders, and headache did not yield any significant associations. Analysis in the TwinsUK failed in replicating the top candidates. Meta-analysis of STOP and TwinsUK did not reveal any significant association. Using elastic net regularization, models exhibited poor predictive capacity tinnitus based on inflammatory markers [sensitivity = 0.52 (95% CI 0.47-0.57), specificity = 0.53 (0.48-0.58), positive predictive value = 0.52 (0.47-0.56), negative predictive values = 0.53 (0.49-0.58), and AUC = 0.53 (0.49-0.56)]. DISCUSSION: Our results did not identify significant associations of the selected inflammatory proteins with constant tinnitus. Future studies examining longitudinal relations among those with more severe tinnitus and using more recent expanded proteomics platforms and sampling of cerebrospinal fluid could increase the likelihood of identifying relevant molecular biomarkers.

Pharmacogenomics of GLP-1 receptor agonists: a genome-wide analysis of observational data and large randomised controlled trials.

BACKGROUND: In the treatment of type 2 diabetes, GLP-1 receptor agonists lower blood glucose concentrations, body weight, and have cardiovascular benefits. The efficacy and side effects of GLP-1 receptor agonists vary between people. Human pharmacogenomic studies of this inter-individual variation can provide both biological insight into drug action and provide biomarkers to inform clinical decision making. We therefore aimed to identify genetic variants associated with glycaemic response to GLP-1 receptor agonist treatment. METHODS: In this genome-wide analysis we included adults (aged ≥18 years) with type 2 diabetes treated with GLP-1 receptor agonists with baseline HbA1c of 7% or more (53 mmol/mol) from four prospective observational cohorts (DIRECT, PRIBA, PROMASTER, and GoDARTS) and two randomised clinical trials (HARMONY phase 3 and AWARD). The primary endpoint was HbA1c reduction at 6 months after starting GLP-1 receptor agonists. We evaluated variants in GLP1R, then did a genome-wide association study and gene-based burden tests. FINDINGS: 4571 adults were included in our analysis, of these, 3339 (73%) were White European, 449 (10%) Hispanic, 312 (7%) American Indian or Alaskan Native, and 471 (10%) were other, and around 2140 (47%) of the participants were women. Variation in HbA1c reduction with GLP-1 receptor agonists treatment was associated with rs6923761G→A (Gly168Ser) in the GLP1R (0·08% [95% CI 0·04-0·12] or 0·9 mmol/mol lower reduction in HbA1c per serine, p=6·0 × 10-5) and low frequency variants in ARRB1 (optimal sequence kernel association test p=6·7 × 10-8), largely driven by rs140226575G→A (Thr370Met; 0·25% [SE 0·06] or 2·7 mmol/mol  [SE 0·7] greater HbA1c reduction per methionine, p=5·2 × 10-6). A similar effect size for the ARRB1 Thr370Met was seen in Hispanic and American Indian or Alaska Native populations who have a higher frequency of this variant (6-11%) than in White European populations. Combining these two genes identified 4% of the population who had a 30% greater reduction in HbA1c than the 9% of the population with the worse response. INTERPRETATION: This genome-wide pharmacogenomic study of GLP-1 receptor agonists provides novel biological and clinical insights. Clinically, when genotype is routinely available at the point of prescribing, individuals with ARRB1 variants might benefit from earlier initiation of GLP-1 receptor agonists. FUNDING: Innovative Medicines Initiative and the Wellcome Trust.

Discovery of drug-omics associations in type 2 diabetes with generative deep-learning models.

The application of multiple omics technologies in biomedical cohorts has the potential to reveal patient-level disease characteristics and individualized response to treatment. However, the scale and heterogeneous nature of multi-modal data makes integration and inference a non-trivial task. We developed a deep-learning-based framework, multi-omics variational autoencoders (MOVE), to integrate such data and applied it to a cohort of 789 people with newly diagnosed type 2 diabetes with deep multi-omics phenotyping from the DIRECT consortium. Using in silico perturbations, we identified drug-omics associations across the multi-modal datasets for the 20 most prevalent drugs given to people with type 2 diabetes with substantially higher sensitivity than univariate statistical tests. From these, we among others, identified novel associations between metformin and the gut microbiota as well as opposite molecular responses for the two statins, simvastatin and atorvastatin. We used the associations to quantify drug-drug similarities, assess the degree of polypharmacy and conclude that drug effects are distributed across the multi-omics modalities.

Genetic analysis of blood molecular phenotypes reveals common properties in the regulatory networks affecting complex traits.

We evaluate the shared genetic regulation of mRNA molecules, proteins and metabolites derived from whole blood from 3029 human donors. We find abundant allelic heterogeneity, where multiple variants regulate a particular molecular phenotype, and pleiotropy, where a single variant associates with multiple molecular phenotypes over multiple genomic regions. The highest proportion of share genetic regulation is detected between gene expression and proteins (66.6%), with a further median shared genetic associations across 49 different tissues of 78.3% and 62.4% between plasma proteins and gene expression. We represent the genetic and molecular associations in networks including 2828 known GWAS variants, showing that GWAS variants are more often connected to gene expression in trans than other molecular phenotypes in the network. Our work provides a roadmap to understanding molecular networks and deriving the underlying mechanism of action of GWAS variants using different molecular phenotypes in an accessible tissue.

Elevated plasma complement factor H related 5 protein is associated with venous thromboembolism.

Venous thromboembolism (VTE) is a common, multi-causal disease with potentially serious short- and long-term complications. In clinical practice, there is a need for improved plasma biomarker-based tools for VTE diagnosis and risk prediction. Here we show, using proteomics profiling to screen plasma from patients with suspected acute VTE, and several case-control studies for VTE, how Complement Factor H Related 5 protein (CFHR5), a regulator of the alternative pathway of complement activation, is a VTE-associated plasma biomarker. In plasma, higher CFHR5 levels are associated with increased thrombin generation potential and recombinant CFHR5 enhanced platelet activation in vitro. GWAS analysis of ~52,000 participants identifies six loci associated with CFHR5 plasma levels, but Mendelian randomization do not demonstrate causality between CFHR5 and VTE. Our results indicate an important role for the regulation of the alternative pathway of complement activation in VTE and that CFHR5 represents a potential diagnostic and/or risk predictive plasma biomarker.