NINJ1 mediates inflammatory cell death, PANoptosis, and lethality during infection conditions and heat stress

Innate immunity provides the first line of defense through multiple mechanisms, including pyrogen production and cell death. While elevated body temperature during infection is beneficial to clear pathogens, heat stress (HS) can lead to inflammation and pathology. Links between pathogen exposure, HS, cytokine release, and inflammation have been observed, but fundamental innate immune mechanisms driving pathology during pathogen exposure and HS remain unclear. Here, we use multiple genetic approaches to elucidate innate immune pathways in infection or LPS and HS models. Our results show that bacteria and LPS robustly increase inflammatory cell death during HS that is dependent on caspase-1, caspase-11, caspase-8, and RIPK3 through the PANoptosis pathway. Caspase-7 also contributes to PANoptosis in this context. Furthermore, NINJ1 is an important executioner of this cell death to release inflammatory molecules, independent of other pore-forming executioner proteins, gasdermin D, gasdermin E, and MLKL. In an in vivo HS model, mortality is reduced by deleting NINJ1 and fully rescued by deleting key PANoptosis molecules. Our findings suggest that therapeutic strategies blocking NINJ1 or its upstream regulators to prevent PANoptosis may reduce the release of inflammatory mediators and benefit patients.

A community challenge to predict clinical outcomes after immune checkpoint blockade in non-small cell lung cancer

BACKGROUND

Predictive biomarkers of immune checkpoint inhibitor (ICI) efficacy are currently lacking for non-small cell lung cancer (NSCLC). Here, we describe the results from the Anti-PD-1 Response Prediction DREAM Challenge, a crowdsourced initiative that enabled the assessment of predictive models by using data from two randomized controlled clinical trials (RCTs) of ICIs in first-line metastatic NSCLC.

METHODS

Participants developed and trained models using public resources. These were evaluated with data from the CheckMate 026 trial (NCT02041533), according to the model-to-data paradigm to maintain patient confidentiality. The generalizability of the models with the best predictive performance was assessed using data from the CheckMate 227 trial (NCT02477826). Both trials were phase III RCTs with a chemotherapy control arm, which supported the differentiation between predictive and prognostic models. Isolated model containers were evaluated using a bespoke strategy that considered the challenges of handling transcriptome data from clinical trials.

RESULTS

A total of 59 teams participated, with 417 models submitted. Multiple predictive models, as opposed to a prognostic model, were generated for predicting overall survival, progression-free survival, and progressive disease status with ICIs. Variables within the models submitted by participants included tumor mutational burden (TMB), programmed death ligand 1 (PD-L1) expression, and gene-expression-based signatures. The best-performing models showed improved predictive power over reference variables, including TMB or PD-L1.

CONCLUSIONS

This DREAM Challenge is the first successful attempt to use protected phase III clinical data for a crowdsourced effort towards generating predictive models for ICI clinical outcomes and could serve as a blueprint for similar efforts in other tumor types and disease states, setting a benchmark for future studies aiming to identify biomarkers predictive of ICI efficacy.

TRIAL REGISTRATION

CheckMate 026; NCT02041533, registered January 22, 2014. CheckMate 227; NCT02477826, registered June 23, 2015.

Comparative analysis identifies genetic and molecular factors associated with prognostic clusters of PANoptosis in glioma, kidney and melanoma cancer

The importance of inflammatory cell death, PANoptosis, in cancer is increasingly being recognized. PANoptosis can promote or inhibit tumorigenesis in context-dependent manners, and a computational approach leveraging transcriptomic profiling of genes involved in PANoptosis has shown that patients can be stratified into PANoptosis High and PANoptosis Low clusters that have significant differences in overall survival for low grade glioma (LGG), kidney renal cell carcinoma (KIRC) and skin cutaneous melanoma (SKCM). However, the molecular mechanisms that contribute to differential prognosis between PANoptosis clusters require further elucidation. Therefore, we performed a comprehensive comparison of genetic, genomic, tumor microenvironment, and pathway characteristics between the PANoptosis High and PANoptosis Low clusters to determine the relevance of each component in driving the differential associations with prognosis for LGG, KIRC and SKCM. Across these cancer types, we found that activation of the proliferation pathway was significantly different between PANoptosis High and Low clusters. In LGG and SKCM, we also found that aneuploidy and immune cell densities and activations contributed to differences in PANoptosis clusters. In individual cancers, we identified important roles for barrier gene pathway activation (in SKCM) and the somatic mutation profiles of driver oncogenes as well as hedgehog signaling pathway activation (in LGG). By identifying these genetic and molecular factors, we can possibly improve the prognosis for at risk-stratified patient populations based on the PANoptosis phenotype in LGG, KIRC and SKCM. This not only advances our mechanistic understanding of cancer but will allow for the selection of optimal treatment strategies.

Inflammatory cell death, PANoptosis, screen identifies host factors in coronavirus innate immune response as therapeutic targets

The COVID-19 pandemic, caused by the β-coronavirus (β-CoV) severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), continues to cause significant global morbidity and mortality. While vaccines have reduced the overall number of severe infections, there remains an incomplete understanding of viral entry and innate immune activation, which can drive pathology. Innate immune responses characterized by positive feedback between cell death and cytokine release can amplify the inflammatory cytokine storm during β-CoV-mediated infection to drive pathology. Therefore, there remains an unmet need to understand innate immune processes in response to β-CoV infections to identify therapeutic strategies. To address this gap, here we used an MHV model and developed a whole genome CRISPR-Cas9 screening approach to elucidate host molecules required for β-CoV infection and inflammatory cell death, PANoptosis, in macrophages, a sentinel innate immune cell. Our screen was validated through the identification of the known MHV receptor Ceacam1 as the top hit, and its deletion significantly reduced viral replication due to loss of viral entry, resulting in a downstream reduction in MHV-induced cell death. Moreover, this screen identified several other host factors required for MHV infection-induced macrophage cell death. Overall, these findings demonstrate the feasibility and power of using genome-wide PANoptosis screens in macrophage cell lines to accelerate the discovery of key host factors in innate immune processes and suggest new targets for therapeutic development to prevent β-CoV-induced pathology.

Pancancer network analysis reveals key master regulators for cancer invasiveness

BACKGROUND

Tumor invasiveness reflects numerous biological changes, including tumorigenesis, progression, and metastasis. To decipher the role of transcriptional regulators (TR) involved in tumor invasiveness, we performed a systematic network-based pan-cancer assessment of master regulators of cancer invasiveness.

MATERIALS AND METHODS

We stratified patients in The Cancer Genome Atlas (TCGA) into invasiveness high (INV-H) and low (INV-L) groups using consensus clustering based on an established robust 24-gene signature to determine the prognostic association of invasiveness with overall survival (OS) across 32 different cancers. We devise a network-based protocol to identify TRs as master regulators (MRs) unique to INV-H and INV-L phenotypes. We validated the activity of MRs coherently associated with INV-H phenotype and worse OS across cancers in TCGA on a series of additional datasets in the Prediction of Clinical Outcomes from the Genomic Profiles (PRECOG) repository.

RESULTS

Based on the 24-gene signature, we defined the invasiveness score for each patient sample and stratified patients into INV-H and INV-L clusters. We observed that invasiveness was associated with worse survival outcomes in almost all cancers and had a significant association with OS in ten out of 32 cancers. Our network-based framework identified common invasiveness-associated MRs specific to INV-H and INV-L groups across the ten prognostic cancers, including COL1A1, which is also part of the 24-gene signature, thus acting as a positive control. Downstream pathway analysis of MRs specific to INV-H phenotype resulted in the identification of several enriched pathways, including Epithelial into Mesenchymal Transition, TGF-β signaling pathway, regulation of Toll-like receptors, cytokines, and inflammatory response, and selective expression of chemokine receptors during T-cell polarization. Most of these pathways have connotations of inflammatory immune response and feasibility for metastasis.

CONCLUSION

Our pan-cancer study provides a comprehensive master regulator analysis of tumor invasiveness and can suggest more precise therapeutic strategies by targeting the identified MRs and downstream enriched pathways for patients across multiple cancers.

NLRP12-PANoptosome activates PANoptosis and pathology in response to heme and PAMPs

Cytosolic innate immune sensors are critical for host defense and form complexes, such as inflammasomes and PANoptosomes, that induce inflammatory cell death. The sensor NLRP12 is associated with infectious and inflammatory diseases, but its activating triggers and roles in cell death and inflammation remain unclear. Here, we discovered that NLRP12 drives inflammasome and PANoptosome activation, cell death, and inflammation in response to heme plus PAMPs or TNF. TLR2/4-mediated signaling through IRF1 induced Nlrp12 expression, which led to inflammasome formation to induce maturation of IL-1β and IL-18. The inflammasome also served as an integral component of a larger NLRP12-PANoptosome that drove inflammatory cell death through caspase-8/RIPK3. Deletion of Nlrp12 protected mice from acute kidney injury and lethality in a hemolytic model. Overall, we identified NLRP12 as an essential cytosolic sensor for heme plus PAMPs-mediated PANoptosis, inflammation, and pathology, suggesting that NLRP12 and molecules in this pathway are potential drug targets for hemolytic and inflammatory diseases.

Whole-genome CRISPR screen identifies RAVER1 as a key regulator of RIPK1-mediated inflammatory cell death, PANoptosis

Transforming growth factor-β-activated kinase 1 (TAK1) is a central regulator of innate immunity, cell death, inflammation, and cellular homeostasis. Therefore, many pathogens carry TAK1 inhibitors (TAK1i). As a host strategy to counteract this, inhibition or deletion of TAK1 induces spontaneous inflammatory cell death, PANoptosis, through the RIPK1-PANoptosome complex, containing the NLRP3 inflammasome and caspase-8/FADD/RIPK3 as integral components; however, PANoptosis also promotes pathological inflammation. Therefore, understanding molecular mechanisms that regulate TAK1i-induced cell death is essential. Here, we report a genome-wide CRISPR screen in macrophages that identified TAK1i-induced cell death regulators, including polypyrimidine tract-binding (PTB) protein 1 (PTBP1), a known regulator of RIPK1, and a previously unknown regulator RAVER1. RAVER1 blocked alternative splicing of Ripk1, and its genetic depletion inhibited TAK1i-induced, RIPK1-mediated inflammasome activation and PANoptosis. Overall, our CRISPR screen identified several positive regulators of PANoptosis. Moreover, our study highlights the utility of genome-wide CRISPR-Cas9 screens in myeloid cells for comprehensive characterization of complex cell death pathways to discover therapeutic targets.

An integrated tumor, immune and microbiome atlas of colon cancer

The lack of multi-omics cancer datasets with extensive follow-up information hinders the identification of accurate biomarkers of clinical outcome. In this cohort study, we performed comprehensive genomic analyses on fresh-frozen samples from 348 patients affected by primary colon cancer, encompassing RNA, whole-exome, deep T cell receptor and 16S bacterial rRNA gene sequencing on tumor and matched healthy colon tissue, complemented with tumor whole-genome sequencing for further microbiome characterization. A type 1 helper T cell, cytotoxic, gene expression signature, called Immunologic Constant of Rejection, captured the presence of clonally expanded, tumor-enriched T cell clones and outperformed conventional prognostic molecular biomarkers, such as the consensus molecular subtype and the microsatellite instability classifications. Quantification of genetic immunoediting, defined as a lower number of neoantigens than expected, further refined its prognostic value. We identified a microbiome signature, driven by Ruminococcus bromii, associated with a favorable outcome. By combining microbiome signature and Immunologic Constant of Rejection, we developed and validated a composite score (mICRoScore), which identifies a group of patients with excellent survival probability. The publicly available multi-omics dataset provides a resource for better understanding colon cancer biology that could facilitate the discovery of personalized therapeutic approaches.

An integrated multi-omic approach demonstrates distinct molecular signatures between human obesity with and without metabolic complications: a case-control study

OBJECTIVES

To examine the hypothesis that obesity complicated by the metabolic syndrome, compared to uncomplicated obesity, has distinct molecular signatures and metabolic pathways.

METHODS

We analyzed a cohort of 39 participants with obesity that included 21 with metabolic syndrome, age-matched to 18 without metabolic complications. We measured in whole blood samples 754 human microRNAs (miRNAs), 704 metabolites using unbiased mass spectrometry metabolomics, and 25,682 transcripts, which include both protein coding genes (PCGs) as well as non-coding transcripts. We then identified differentially expressed miRNAs, PCGs, and metabolites and integrated them using databases such as mirDIP (mapping between miRNA-PCG network), Human Metabolome Database (mapping between metabolite-PCG network) and tools like MetaboAnalyst (mapping between metabolite-metabolic pathway network) to determine dysregulated metabolic pathways in obesity with metabolic complications.

RESULTS

We identified 8 significantly enriched metabolic pathways comprising 8 metabolites, 25 protein coding genes and 9 microRNAs which are each differentially expressed between the subjects with obesity and those with obesity and metabolic syndrome. By performing unsupervised hierarchical clustering on the enrichment matrix of the 8 metabolic pathways, we could approximately segregate the uncomplicated obesity strata from that of obesity with metabolic syndrome.

CONCLUSIONS

The data suggest that at least 8 metabolic pathways, along with their various dysregulated elements, identified via our integrative bioinformatics pipeline, can potentially differentiate those with obesity from those with obesity and metabolic complications.

AI-driven drug repurposing and binding pose meta dynamics identifies novel targets for Monkeypox virus

Monkeypox virus (MPXV) was confirmed in May 2022 and designated a global health emergency by WHO in July 2022. MPX virions are big, enclosed, brick-shaped, and contain a linear, double-stranded DNA genome as well as enzymes. MPXV particles bind to the host cell membrane via a variety of viral-host protein interactions. As a result, the wrapped structure is a potential therapeutic target. DeepRepurpose, an artificial intelligence-based compound-viral proteins interaction framework, was used via a transfer learning setting to prioritize a set of FDA approved and investigational drugs which can potentially inhibit MPXV viral proteins. To filter and narrow down the lead compounds from curated collections of pharmaceutical compounds, we used a rigorous computational framework that included homology modeling, molecular docking, dynamic simulations, binding free energy calculations, and binding pose metadynamics. We identified Elvitegravir as a potential inhibitor of MPXV virus using our comprehensive pipeline.

Pancancer transcriptomic profiling identifies key PANoptosis markers as therapeutic targets for oncology

Resistance to programmed cell death (PCD) is a hallmark of cancer. While some PCD components are prognostic in cancer, the roles of many molecules can be masked by redundancies and crosstalks between PCD pathways, impeding the development of targeted therapeutics. Recent studies characterizing these redundancies have identified PANoptosis, a unique innate immune-mediated inflammatory PCD pathway that integrates components from other PCD pathways. Here, we designed a systematic computational framework to determine the pancancer clinical significance of PANoptosis and identify targetable biomarkers. We found that high expression of PANoptosis genes was detrimental in low grade glioma (LGG) and kidney renal cell carcinoma (KIRC). ZBP1, ADAR, CASP2, CASP3, CASP4, CASP8 and GSDMD expression consistently had negative effects on prognosis in LGG across multiple survival models, while AIM2, CASP3, CASP4 and TNFRSF10 expression had negative effects for KIRC. Conversely, high expression of PANoptosis genes was beneficial in skin cutaneous melanoma (SKCM), with ZBP1, NLRP1, CASP8 and GSDMD expression consistently having positive prognostic effects. As a therapeutic proof-of-concept, we treated melanoma cells with combination therapy that activates ZBP1 and showed that this treatment induced PANoptosis. Overall, through our systematic framework, we identified and validated key innate immune biomarkers from PANoptosis which can be targeted to improve patient outcomes in cancers.

Investigation of the Effect of Exendin-4 on Oleic Acid-Induced Steatosis in HepG2 Cells Using Fourier Transform Infrared Spectroscopy

Non-alcoholic fatty liver disease (NAFLD) is a common liver lesion that is untreatable with medications. Glucagon-like peptide-1 receptor (GLP-1R) agonists have recently emerged as a potential NAFLD pharmacotherapy. However, the molecular mechanisms underlying these drugs’ beneficial effects are not fully understood. Using Fourier transform infrared (FTIR) spectroscopy, we sought to investigate the biochemical changes in a steatosis cell model treated or not with the GLP-1R agonist Exendin-4 (Ex-4). HepG2 cells were made steatotic with 400 µM of oleic acid and then treated with 200 nM Ex-4 in order to reduce lipid accumulation. We quantified steatosis using the Oil Red O staining method. We investigated the biochemical alterations induced by steatosis and Ex-4 treatment using Fourier transform infrared (FTIR) spectroscopy and chemometric analyses. Analysis of the Oil Red O staining showed that Ex-4 significantly reduces steatosis. This reduction was confirmed by FTIR analysis, as the phospholipid band (C=O) at 1740 cm⁻¹ in Ex-4 treated cells is significantly decreased compared to steatotic cells. The principal component analysis score plots for both the lipid and protein regions showed that the untreated and Ex-4-treated samples, while still separated, are clustered close to each other, far from the steatotic cells. The biochemical and structural changes induced by OA-induced lipotoxicity are at least partially reversed upon Ex-4 treatment. FTIR and chemometric analyses revealed that Ex-4 significantly reduces OA-induced lipid accumulation, and Ex-4 also restored the lipid and protein biochemical alterations caused by lipotoxicity-induced oxidative stress. In combination with chemometric analyses, FTIR spectroscopy may offer new approaches for investigating the mechanisms underpinning NAFLD.

Molecular mechanism of RIPK1 and caspase-8 in homeostatic type I interferon production and regulation

Type I interferons (IFNs) are essential innate immune proteins that maintain tissue homeostasis through tonic expression and can be upregulated to drive antiviral resistance and inflammation upon stimulation. However, the mechanisms that inhibit aberrant IFN upregulation in homeostasis and the impacts of tonic IFN production on health and disease remain enigmatic. Here, we report that caspase-8 negatively regulates type I IFN production by inhibiting the RIPK1-TBK1 axis during homeostasis across multiple cell types and tissues. When caspase-8 is deleted or inhibited, RIPK1 interacts with TBK1 to drive elevated IFN production, leading to heightened resistance to norovirus infection in macrophages but also early onset lymphadenopathy in mice. Combined deletion of caspase-8 and RIPK1 reduces the type I IFN signaling and lymphadenopathy, highlighting the critical role of RIPK1 in this process. Overall, our study identifies a mechanism to constrain tonic type I IFN during homeostasis which could be targeted for infectious and inflammatory diseases.

Wang et al. report the mechanistic regulation of homeostatic type I IFN production by caspase-8 through the RIPK1-TBK1 axis. Hyper-activation of this pathway due to loss of caspase-8 has profound physiological impacts on natural resistance to viral infection and the pathogenesis of lymphadenopathy.

Dysregulated Metabolic Pathways in Subjects with Obesity and Metabolic Syndrome

Background: Obesity coexists with variable features of metabolic syndrome, which is associated with dysregulated metabolic pathways. We assessed potential associations between serum metabolites and features of metabolic syndrome in Arabic subjects with obesity. Methods: We analyzed a dataset of 39 subjects with obesity only (OBO, n = 18) age-matched to subjects with obesity and metabolic syndrome (OBM, n = 21). We measured 1069 serum metabolites and correlated them to clinical features. Results: A total of 83 metabolites, mostly lipids, were significantly different (p < 0.05) between the two groups. Among lipids, 22 sphingomyelins were decreased in OBM compared to OBO. Among non-lipids, quinolinate, kynurenine, and tryptophan were also decreased in OBM compared to OBO. Sphingomyelin is negatively correlated with glucose, HbA1C, insulin, and triglycerides but positively correlated with HDL, LDL, and cholesterol. Differentially enriched pathways include lysine degradation, amino sugar and nucleotide sugar metabolism, arginine and proline metabolism, fructose and mannose metabolism, and galactose metabolism. Conclusions: Metabolites and pathways associated with chronic inflammation are differentially expressed in subjects with obesity and metabolic syndrome compared to subjects with obesity but without the clinical features of metabolic syndrome.

ZBP1-dependent inflammatory cell death, PANoptosis, and cytokine storm disrupt IFN therapeutic efficacy during coronavirus infection

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the virus responsible for coronavirus disease 2019 (COVID-19), continues to cause substantial morbidity and mortality in the ongoing global pandemic. Understanding the fundamental mechanisms that govern innate immune and inflammatory responses during SARS-CoV-2 infection is critical for developing effective therapeutic strategies. Whereas interferon (IFN)-based therapies are generally expected to be beneficial during viral infection, clinical trials in COVID-19 have shown limited efficacy and potential detrimental effects of IFN treatment during SARS-CoV-2 infection. However, the underlying mechanisms responsible for this failure remain unknown. In this study, we found that IFN induced Z-DNA-binding protein 1 (ZBP1)-mediated inflammatory cell death, PANoptosis, in human and murine macrophages and in the lungs of mice infected with β-coronaviruses, including SARS-CoV-2 and mouse hepatitis virus (MHV). In patients with COVID-19, expression of the innate immune sensor ZBP1 was increased in immune cells from those who succumbed to the disease compared with those who recovered, further suggesting a link between ZBP1 and pathology. In mice, IFN-β treatment after β-coronavirus infection increased lethality, and genetic deletion of Zbp1 or its Zα domain suppressed cell death and protected the mice from IFN-mediated lethality during β-coronavirus infection. Overall, our results identify that ZBP1 induced during coronavirus infection limits the efficacy of IFN therapy by driving inflammatory cell death and lethality. Therefore, inhibiting ZBP1 activity may improve the efficacy of IFN therapy, paving the way for the development of new and critically needed therapeutics for COVID-19 as well as other infections and inflammatory conditions where IFN-mediated cell death and pathology occur.

Graphene oxide activates B cells with upregulation of granzyme B expression: evidence at the single-cell level for its immune-modulatory properties and anticancer activity

We recently found by single-cell mass cytometry that ex vivo human B cells internalize graphene oxide (GO). The functional impact of such uptake on B cells remains unexplored. Here, we disclosed the effects of GO and amino-functionalized GO (GONH₂) interacting with human B cells in vitro and ex vivo at the protein and gene expression levels. Moreover, our study considered three different subpopulations of B cells and their functionality in terms of: (i) cytokine production, (ii) activation markers, (iii) killing activity towards cancer cells. Single-cell mass cytometry screening revealed the higher impact of GO on cell viability towards naïve, memory, and plasma B cell subsets. Different cytokines such as granzyme B (GrB) and activation markers, like CD69, CD80, CD138, and CD38, were differently regulated by GONH₂ compared to GO, supporting possible diverse B cell activation paths. Moreover, co-culture experiments also suggest the functional ability of both GOs to activate B cells and therefore enhance the toxicity towards HeLa cancer cell line. Complete transcriptomic analysis on a B cell line highlighted the distinctive GO and GONH₂ elicited responses, inducing pathways such as B cell receptor and CD40 signaling pathways, key players for GrB secretion. B cells were regularly left behind the scenes in graphene biological studies; our results may open new horizons in the development of GO-based immune-modulatory strategies having B cell as main actors.

Fasting-Mimicking Diet Is Safe and Reshapes Metabolism and Antitumor Immunity in Patients with Cancer

In tumor-bearing mice, cyclic fasting or fasting-mimicking diets (FMD) enhance the activity of antineoplastic treatments by modulating systemic metabolism and boosting antitumor immunity. Here we conducted a clinical trial to investigate the safety and biological effects of cyclic, five-day FMD in combination with standard antitumor therapies. In 101 patients, the FMD was safe, feasible, and resulted in a consistent decrease of blood glucose and growth factor concentration, thus recapitulating metabolic changes that mediate fasting/FMD anticancer effects in preclinical experiments. Integrated transcriptomic and deep-phenotyping analyses revealed that FMD profoundly reshapes anticancer immunity by inducing the contraction of peripheral blood immunosuppressive myeloid and regulatory T-cell compartments, paralleled by enhanced intratumor Th1/cytotoxic responses and an enrichment of IFNγ and other immune signatures associated with better clinical outcomes in patients with cancer. Our findings lay the foundations for phase II/III clinical trials aimed at investigating FMD antitumor efficacy in combination with standard antineoplastic treatments. SIGNIFICANCE: Cyclic FMD is well tolerated and causes remarkable systemic metabolic changes in patients with different tumor types and treated with concomitant antitumor therapies. In addition, the FMD reshapes systemic and intratumor immunity, finally activating several antitumor immune programs. Phase II/III clinical trials are needed to investigate FMD antitumor activity/efficacy.This article is highlighted in the In This Issue feature, p. 1.

Characteristic MicroRNAs Linked to Dysregulated Metabolic Pathways in Qatari Adult Subjects With Obesity and Metabolic Syndrome

BACKGROUND

Obesity-associated dysglycemia is associated with metabolic disorders. MicroRNAs (miRNAs) are known regulators of metabolic homeostasis. We aimed to assess the relationship of circulating miRNAs with clinical features in obese Qatari individuals.

METHODS

We analyzed a dataset of 39 age-matched patients that includes 18 subjects with obesity only (OBO) and 21 subjects with obesity and metabolic syndrome (OBM). We measured 754 well-characterized human microRNAs (miRNAs) and identified differentially expressed miRNAs along with their significant associations with clinical markers in these patients.

RESULTS

A total of 64 miRNAs were differentially expressed between metabolically healthy obese (OBO) versus metabolically unhealthy obese (OBM) patients. Thirteen out of 64 miRNAs significantly correlated with at least one clinical trait of the metabolic syndrome. Six out of the thirteen demonstrated significant association with HbA1c levels; miR-331-3p, miR-452-3p, and miR-485-5p were over-expressed, whereas miR-153-3p, miR-182-5p, and miR-433-3p were under-expressed in the OBM patients with elevated HbA1c levels. We also identified, miR-106b-3p, miR-652-3p, and miR-93-5p that showed a significant association with creatinine; miR-130b-5p, miR-363-3p, and miR-636 were significantly associated with cholesterol, whereas miR-130a-3p was significantly associated with LDL. Additionally, miR-652-3p's differential expression correlated significantly with HDL and creatinine.

CONCLUSIONS

MicroRNAs associated with metabolic syndrome in obese subjects may have a pathophysiologic role and can serve as markers for obese individuals predisposed to various metabolic diseases like diabetes.

Simple risk score to screen for prediabetes: A cross-sectional study from the Qatar Biobank cohort

AIMS/INTRODUCTION

The progression from prediabetes to type 2 diabetes is preventable by lifestyle intervention and/or pharmacotherapy in a large fraction of individuals with prediabetes. Our objective was to develop a risk score to screen for prediabetes in the Middle East, where diabetes prevalence is one of the highest in the world.

MATERIALS AND METHODS

In this cross-sectional, case-control study, we used data of 4,895 controls and 2,373 prediabetic adults obtained from the Qatar Biobank cohort. Significant risk factors were identified by logistic regression and other machine learning methods. The receiver operating characteristic was used to calculate the area under curve, cut-off point, sensitivity, specificity, positive and negative predictive values. The prediabetes risk score was developed from data of Qatari citizens, as well as long-term (≥15 years) residents.

RESULTS

The significant risk factors for the Prediabetes Risk Score in Qatar were age, sex, body mass index, waist circumference and blood pressure. The risk score ranges from 0 to 45. The area under the curve of the score was 80% (95% confidence interval 78-83%), and the cut-off point of 16 yielded sensitivity and specificity of 86.2% (95% confidence interval 82.7-89.2%) and 57.9% (95% confidence interval 65.5-71.4%), respectively. Prediabetes Risk Score in Qatar performed equally in Qatari nationals and long-term residents.

CONCLUSIONS

Prediabetes Risk Score in Qatar is the first prediabetes screening score developed in a Middle Eastern population. It only uses risk factors measured non-invasively, is simple, cost-effective, and can be easily understood by the general public and health providers. Prediabetes Risk Score in Qatar is an important tool for early detection of prediabetes, and can help tremendously in curbing the diabetes epidemic in the region.

Network-based identification of key master regulators associated with an immune-silent cancer phenotype

A cancer immune phenotype characterized by an active T-helper 1 (Th1)/cytotoxic response is associated with responsiveness to immunotherapy and favorable prognosis across different tumors. However, in some cancers, such an intratumoral immune activation does not confer protection from progression or relapse. Defining mechanisms associated with immune evasion is imperative to refine stratification algorithms, to guide treatment decisions and to identify candidates for immune-targeted therapy. Molecular alterations governing mechanisms for immune exclusion are still largely unknown. The availability of large genomic datasets offers an opportunity to ascertain key determinants of differential intratumoral immune response. We follow a network-based protocol to identify transcription regulators (TRs) associated with poor immunologic antitumor activity. We use a consensus of four different pipelines consisting of two state-of-the-art gene regulatory network inference techniques, regularized gradient boosting machines and ARACNE to determine TR regulons, and three separate enrichment techniques, including fast gene set enrichment analysis, gene set variation analysis and virtual inference of protein activity by enriched regulon analysis to identify the most important TRs affecting immunologic antitumor activity. These TRs, referred to as master regulators (MRs), are unique to immune-silent and immune-active tumors, respectively. We validated the MRs coherently associated with the immune-silent phenotype across cancers in The Cancer Genome Atlas and a series of additional datasets in the Prediction of Clinical Outcomes from Genomic Profiles repository. A downstream analysis of MRs specific to the immune-silent phenotype resulted in the identification of several enriched candidate pathways, including NOTCH1, TGF-$\beta $, Interleukin-1 and TNF-$\alpha $ signaling pathways. TGFB1I1 emerged as one of the main negative immune modulators preventing the favorable effects of a Th1/cytotoxic response.

Data-driven enhancement of cubic phase stability in mixed-cation perovskites

Mixing cations has been a successful strategy in perovskite synthesis by solution-processing, delivering improvements in the thermodynamic stability as well as in the lattice parameter control. Unfortunately, the relation between a given cation mixture and the associated structural deformation is not well-established, a fact that hinders an adequate identification of the optimum chemical compositions. Such difficulty arises since local distortion and microscopic disorder influence structural stability and also determine phase segregation. Hence, the search for an optimum composition is currently based on experimental trial and error, a tedious and high-cost process. Here, we report on a machine-learning-reinforced cubic-phase-perovskite stability predictor that has been constructed over an extensive dataset of first-principles calculations. Such a predictor allows us to determine the cubic phase stability at a given cation mixture regardless of the various cations’ pair and concentration, even assessing very dilute concentrations, a notoriously challenging task for first-principles calculations. In particular, we construct machine learning models, predicting multiple target quantities such as the enthalpy of mixing and various octahedral distortions. It is then the combination of these targets that guide the laboratory synthesis. Our theoretical analysis is also validated by the experimental synthesis and characterization of methylammonium–dimethylammonium-mixed perovskite thin films, demonstrating the ability of the stability predictor to drive the chemical design of this class of materials.

A Modelling Framework for Embedding-based Predictions for Compound-Viral Protein Activity

Motivation

A global effort is underway to identify compounds for the treatment of COVID-19. Since de novo compound design is an extremely long, time-consuming, and expensive process, efforts are underway to discover existing compounds that can be repurposed for COVID-19 and new viral diseases.

Model

We propose a machine learning representation framework that uses deep learning induced vector embeddings of compounds and viral proteins as features to predict compound-viral protein activity. The prediction model in-turn uses a consensus framework to rank approved compounds against viral proteins of interest.

Results

Our consensus framework achieves a highmean Pearson correlation of 0.916, mean R2 of 0.840 and a low mean squared error of 0.313 for the task of compound-viral protein activity prediction on an independent test set. As a use case, we identify a ranked list of 47 compounds common to three main proteins of SARS-COV-2 virus (PL-PRO, 3CL-PRO and Spike protein) as potential targets including 21 antivirals, 15 anticancer, 5 antibiotics and 6 other investigationalhuman compounds.We performadditional molecular docking simulations to demonstrate thatmajority of these compounds have low binding energies and thus high binding affinity with the potential to be effective against the SARS-COV-2 virus.

Availability

All the source code and data is available at: https://github.com/raghvendra5688/Drug-Repurposing and https://dx.doi.org/10.17632/8rrwnbcgmx.3. We also implemented a web-server at: https://machinelearning-protein.qcri.org/index.html.

Supplementary information

Supplementary data are available at Bioinformatics online.

Ancestry-associated transcriptomic profiles of breast cancer in patients of African, Arab, and European ancestry

Breast cancer largely dominates the global cancer burden statistics; however, there are striking disparities in mortality rates across countries. While socioeconomic factors contribute to population-based differences in mortality, they do not fully explain disparity among women of African ancestry (AA) and Arab ancestry (ArA) compared to women of European ancestry (EA). In this study, we sought to identify molecular differences that could provide insight into the biology of ancestry-associated disparities in clinical outcomes. We applied a unique approach that combines the use of curated survival data from The Cancer Genome Atlas (TCGA) Pan-Cancer clinical data resource, improved single-nucleotide polymorphism-based inferred ancestry assignment, and a novel breast cancer subtype classification to interrogate the TCGA and a local Arab breast cancer dataset. We observed an enrichment of BasalMyo tumors in AA patients (38 vs 16.5% in EA, p = 1.30E - 10), associated with a significant worse overall (hazard ratio (HR) = 2.39, p = 0.02) and disease-specific survival (HR = 2.57, p = 0.03). Gene set enrichment analysis of BasalMyo AA and EA samples revealed differences in the abundance of T-regulatory and T-helper type 2 cells, and enrichment of cancer-related pathways with prognostic implications (AA: PI3K-Akt-mTOR and ErbB signaling; EA: EGF, estrogen-dependent and DNA repair signaling). Strikingly, AMPK signaling was associated with opposing prognostic connotation (AA: 10-year HR = 2.79, EA: 10-year HR = 0.34). Analysis of ArA patients suggests enrichment of BasalMyo tumors with a trend for differential enrichment of T-regulatory cells and AMPK signaling. Together, our findings suggest that the disparity in the clinical outcome of AA breast cancer patients is likely related to differences in cancer-related and microenvironmental features.

Artificial Intelligence (AI) based machine learning models predict glucose variability and hypoglycaemia risk in patients with type 2 diabetes on a multiple drug regimen who fast during ramadan (The PROFAST - IT Ramadan study)

OBJECTIVE

To develop a machine-based algorithm from clinical and demographic data, physical activity and glucose variability to predict hyperglycaemic and hypoglycaemic excursions in patients with type 2 diabetes on multiple glucose lowering therapies who fast during Ramadan.

PATIENTS AND METHODS

Thirteen patients (10 males and three females) with type 2 diabetes on 3 or more anti-diabetic medications were studied with a Fitbit-2 pedometer device and Freestyle Libre (Abbott Diagnostics) 2 weeks before and 2 weeks during Ramadan. Several machine learning techniques were trained to predict blood glucose levels in a regression framework utilising physical activity and contemporaneous blood glucose levels, comparing Ramadan to non-Ramadan days.

RESULTS

The median age of participants was 51 years (IQR 49-52); median BMI was 33.2 kg/m² (IQR 33.0-35.9) and median HbA1c was 7.3% (IQR 6.7-7.8). The optimal model using physical activity achieved an R² of 0.548 and a mean absolute error (MAE) of 30.30. The addition of electronic health record (ehr) information increased R² to 0.636 and reduced MAE to 26.89 and the time of the day feature further increased R² to 0.768 and reduced MAE to 20.55. Combining all the features together resulted in an optimal XGBoost model with an R² of 0.836 and MAE of 17.47. This model accurately estimated normal glucose levels in 2584/2715 (95.2%) readings and hyperglycaemic events in 852/1031 (82.6%) readings, but fewer hypoglycaemic events (48/172 (27.9%)). The optimal XGBoost model prioritized age, gender, BMI and HbA1c followed by glucose levels and physical activity. Interestingly, the blood glucose level prediction by our model was influenced by use of SGLT2i.

CONCLUSION

XGBoost, a machine learning AI algorithm achieves high predictive performance for normal and hyperglycaemic excursions, but has limited predictive value for hypoglycaemia in patients on multiple therapies who fast during Ramadan.

BCrystal: an interpretable sequence-based protein crystallization predictor

MOTIVATION

X-ray crystallography has facilitated the majority of protein structures determined to date. Sequence-based predictors that can accurately estimate protein crystallization propensities would be highly beneficial to overcome the high expenditure, large attrition rate, and to reduce the trial-and-error settings required for crystallization.

RESULTS

In this study, we present a novel model, BCrystal, which uses an optimized gradient boosting machine (XGBoost) on sequence, structural and physio-chemical features extracted from the proteins of interest. BCrystal also provides explanations, highlighting the most important features for the predicted crystallization propensity of an individual protein using the SHAP algorithm. On three independent test sets, BCrystal outperforms state-of-the-art sequence-based methods by more than 12.5% in accuracy, 18% in recall and 0.253 in Matthew's correlation coefficient, with an average accuracy of 93.7%, recall of 96.63% and Matthew's correlation coefficient of 0.868. For relative solvent accessibility of exposed residues, we observed higher values to associate positively with protein crystallizability and the number of disordered regions, fraction of coils and tripeptide stretches that contain multiple histidines associate negatively with crystallizability. The higher accuracy of BCrystal enables it to accurately screen for sequence variants with enhanced crystallizability.

AVAILABILITY AND IMPLEMENTATION

Our BCrystal webserver is at https://machinelearning-protein.qcri.org/ and source code is available at https://github.com/raghvendra5688/BCrystal.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Oncogenic states dictate the prognostic and predictive connotations of intratumoral immune response

BACKGROUND

An immune active cancer phenotype typified by a T helper 1 (Th-1) immune response has been associated with increased responsiveness to immunotherapy and favorable prognosis in some but not all cancer types. The reason of this differential prognostic connotation remains unknown.

METHODS

To explore the contextual prognostic value of cancer immune phenotypes, we applied a multimodal pan-cancer analysis among 31 different histologies (9282 patients), encompassing immune and oncogenic transcriptomic analysis, mutational and neoantigen load and copy number variations.

RESULTS

We demonstrated that the favorable prognostic connotation conferred by the presence of a Th-1 immune response was abolished in tumors displaying specific tumor-cell intrinsic attributes such as high transforming growth factor-beta (TGF-β) signaling and low proliferation capacity. This observation was independent of mutation rate. We validated this observation in the context of immune checkpoint inhibition. WNT-β catenin, barrier molecules, Notch, hedgehog, mismatch repair, telomerase activity and AMPK signaling were the pathways most coherently associated with an immune silent phenotype together with mutations of driver genes including IDH1/2, FOXA2, HDAC3, PSIP1, MAP3K1, KRAS, NRAS, EGFR, FGFR3, WNT5A and IRF7.

CONCLUSIONS

This is the first systematic study demonstrating that the prognostic and predictive role of a bona fide favorable intratumoral immune response is dependent on the disposition of specific oncogenic pathways. This information could be used to refine stratification algorithms and prioritize hierarchically relevant targets for combination therapies.

The future of sleep health: a data-driven revolution in sleep science and medicine

In recent years, there has been a significant expansion in the development and use of multi-modal sensors and technologies to monitor physical activity, sleep and circadian rhythms. These developments make accurate sleep monitoring at scale a possibility for the first time. Vast amounts of multi-sensor data are being generated with potential applications ranging from large-scale epidemiological research linking sleep patterns to disease, to wellness applications, including the sleep coaching of individuals with chronic conditions. However, in order to realise the full potential of these technologies for individuals, medicine and research, several significant challenges must be overcome. There are important outstanding questions regarding performance evaluation, as well as data storage, curation, processing, integration, modelling and interpretation. Here, we leverage expertise across neuroscience, clinical medicine, bioengineering, electrical engineering, epidemiology, computer science, mHealth and human-computer interaction to discuss the digitisation of sleep from a inter-disciplinary perspective. We introduce the state-of-the-art in sleep-monitoring technologies, and discuss the opportunities and challenges from data acquisition to the eventual application of insights in clinical and consumer settings. Further, we explore the strengths and limitations of current and emerging sensing methods with a particular focus on novel data-driven technologies, such as Artificial Intelligence.

Probing the fibrillation of lysozyme by nanoscale-infrared spectroscopy

Amyloid fibrillation is the root cause of several neuro as well as non-neurological disorders. Understanding the molecular basis of amyloid aggregate formation is crucial for deciphering various neurodegenerative diseases. In our study, we have examined the lysozyme fibrillation process using nano-infrared spectroscopy (nanoIR). NanoIR enabled us to investigate both structural and chemical characteristics of lysozyme fibrillar species concurrently. The spectroscopic results indicate that lysozyme transformed into a fibrillar structure having mainly parallel β-sheets, with almost no antiparallel β-sheets. Features such as protein stiffness have a good correlation with obtained secondary structural information showing the state of the protein within the fibrillation state. The structural and chemical details were compared with transmission electron microscopy (TEM) and circular dichroism (CD). We have utilized nanoIR and measured infrared spectra to characterize lysozyme amyloid fibril structures in terms of morphology, molecular structure, secondary structure content, stability, and size of the cross-β core. We have shown that the use of nanoIR can complement other biophysical studies to analyze the aggregation process and is particularly useful for studying proteins involved in aggregation to help in designing molecules against amyloid aggregation. Specifically, the nanoIR spectra afford higher resolution information and a characteristic fingerprint for determining states of aggregation.Communicated by Ramaswamy H. Sarma.

Investigation of the Effect of PD-L1 Blockade on Triple Negative Breast Cancer Cells Using Fourier Transform Infrared Spectroscopy

Interactions between programmed death-1 (PD-1) with its ligand PD-L1 on tumor cells can antagonize T cell responses. Inhibiting these interactions using immune checkpoint inhibitors has shown promise in cancer immunotherapy. MDA-MB-231 is a triple negative breast cancer cell line that expresses PD-L1. In this study, we investigated the biochemical changes in MDA-MB-231 cells following treatment with atezolizumab, a specific PD-L1 blocker. Our readouts were Fourier Transform Infrared (FTIR) spectroscopy and flow cytometric analyses. Chemometrical analysis, such as principal component analysis (PCA), was applied to delineate the spectral differences. We were able to identify the chemical alterations in both protein and lipid structure of the treated cells. We found that there was a shift from random coil and α-helical structure to β-sheet conformation of PD-L1 on tumor cells due to atezolizumab treatment, which could hinder binding with its receptors on immune cells, ensuring sustained T cell activation for potent immune responses. This work provides novel information about the effects of atezolizumab at molecular and cellular levels. FTIR bio-spectroscopy, in combination with chemometric analyses, may expedite research and offer new approaches for cancer immunology.

DeepSol: a deep learning framework for sequence-based protein solubility prediction

Motivation

Protein solubility plays a vital role in pharmaceutical research and production yield. For a given protein, the extent of its solubility can represent the quality of its function, and is ultimately defined by its sequence. Thus, it is imperative to develop novel, highly accurate in silico sequence-based protein solubility predictors. In this work we propose, DeepSol, a novel Deep Learning-based protein solubility predictor. The backbone of our framework is a convolutional neural network that exploits k-mer structure and additional sequence and structural features extracted from the protein sequence.

Results

DeepSol outperformed all known sequence-based state-of-the-art solubility prediction methods and attained an accuracy of 0.77 and Matthew's correlation coefficient of 0.55. The superior prediction accuracy of DeepSol allows to screen for sequences with enhanced production capacity and can more reliably predict solubility of novel proteins.

Availability and implementation

DeepSol's best performing models and results are publicly deposited at https://doi.org/10.5281/zenodo.1162886 (Khurana and Mall, 2018).

Supplementary information

Supplementary data are available at Bioinformatics online.

Benchmark on a large cohort for sleep-wake classification with machine learning techniques

Accurately measuring sleep and its quality with polysomnography (PSG) is an expensive task. Actigraphy, an alternative, has been proven cheap and relatively accurate. However, the largest experiments conducted to date, have had only hundreds of participants. In this work, we processed the data of the recently published Multi-Ethnic Study of Atherosclerosis (MESA) Sleep study to have both PSG and actigraphy data synchronized. We propose the adoption of this publicly available large dataset, which is at least one order of magnitude larger than any other dataset, to systematically compare existing methods for the detection of sleep-wake stages, thus fostering the creation of new algorithms. We also implemented and compared state-of-the-art methods to score sleep-wake stages, which range from the widely used traditional algorithms to recent machine learning approaches. We identified among the traditional algorithms, two approaches that perform better than the algorithm implemented by the actigraphy device used in the MESA Sleep experiments. The performance, in regards to accuracy and F 1 score of the machine learning algorithms, was also superior to the device's native algorithm and comparable to human annotation. Future research in developing new sleep-wake scoring algorithms, in particular, machine learning approaches, will be highly facilitated by the cohort used here. We exemplify this potential by showing that two particular deep-learning architectures, CNN and LSTM, among the many recently created, can achieve accuracy scores significantly higher than other methods for the same tasks.

Correction: Exploring new approaches towards the formability of mixed-ion perovskites by DFT and machine learning

Correction for 'Exploring new approaches towards the formability of mixed-ion perovskites by DFT and machine learning' by Heesoo Park et al., Phys. Chem. Chem. Phys., 2019, DOI: 10.1039/c8cp06528d.

PaRSnIP: sequence-based protein solubility prediction using gradient boosting machine

Motivation

Protein solubility can be a decisive factor in both research and production efficiency, and in silico sequence-based predictors that can accurately estimate solubility outcomes are highly sought.

Results

In this study, we present a novel approach termed PRotein SolubIlity Predictor (PaRSnIP), which uses a gradient boosting machine algorithm as well as an approximation of sequence and structural features of the protein of interest. Based on an independent test set, PaRSnIP outperformed other state-of-the-art sequence-based methods by more than 9% in accuracy and 0.17 in Matthew's correlation coefficient, with an overall accuracy of 74% and Matthew's correlation coefficient of 0.48. Additionally, PaRSnIP provides importance scores for all features used in training. We observed higher fractions of exposed residues to associate positively with protein solubility and tripeptide stretches with multiple histidines to associate negatively with solubility. The improved prediction accuracy of PaRSnIP should enable it to predict protein solubility with greater reliability and to screen for sequence variants with enhanced manufacturability.

Availability and implementation

PaRSnIP software is available for download under GitHub (https://github.com/RedaRawi/PaRSnIP).

Contact

gwo-yu.chuang@nih.gov.

Supplementary information

Supplementary data are available at Bioinformatics online.

Exploring new approaches towards the formability of mixed-ion perovskites by DFT and machine learning

Recent years have witnessed a growing effort in engineering and tuning the properties of hybrid halide perovskites as light absorbers.

Comparison and assessment of family- and population-based genotype imputation methods in large pedigrees

Genotype imputation is widely used in genome-wide association studies to boost variant density, allowing increased power in association testing. Many studies currently include pedigree data due to increasing interest in rare variants coupled with the availability of appropriate analysis tools. The performance of population-based (subjects are unrelated) imputation methods is well established. However, the performance of family- and population-based imputation methods on family data has been subject to much less scrutiny. Here, we extensively compare several family- and population-based imputation methods on family data of large pedigrees with both European and African ancestry. Our comparison includes many widely used family- and population-based tools and another method, Ped_Pop, which combines family- and population-based imputation results. We also compare four subject selection strategies for full sequencing to serve as the reference panel for imputation: GIGI-Pick, ExomePicks, PRIMUS, and random selection. Moreover, we compare two imputation accuracy metrics: the Imputation Quality Score and Pearson's correlation R ² for predicting power of association analysis using imputation results. Our results show that (1) GIGI outperforms Merlin; (2) family-based imputation outperforms population-based imputation for rare variants but not for common ones; (3) combining family- and population-based imputation outperforms all imputation approaches for all minor allele frequencies; (4) GIGI-Pick gives the best selection strategy based on the R ² criterion; and (5) R ² is the best measure of imputation accuracy. Our study is the first to extensively evaluate the imputation performance of many available family- and population-based tools on the same family data and provides guidelines for future studies.

DeepCrystal: a deep learning framework for sequence-based protein crystallization prediction

Motivation

Protein structure determination has primarily been performed using X-ray crystallography. To overcome the expensive cost, high attrition rate and series of trial-and-error settings, many in-silico methods have been developed to predict crystallization propensities of proteins based on their sequences. However, the majority of these methods build their predictors by extracting features from protein sequences, which is computationally expensive and can explode the feature space. We propose DeepCrystal, a deep learning framework for sequence-based protein crystallization prediction. It uses deep learning to identify proteins which can produce diffraction-quality crystals without the need to manually engineer additional biochemical and structural features from sequence. Our model is based on convolutional neural networks, which can exploit frequently occurring k-mers and sets of k-mers from the protein sequences to distinguish proteins that will result in diffraction-quality crystals from those that will not.

Results

Our model surpasses previous sequence-based protein crystallization predictors in terms of recall, F-score, accuracy and Matthew’s correlation coefficient (MCC) on three independent test sets. DeepCrystal achieves an average improvement of 1.4, 12.1% in recall, when compared to its closest competitors, Crysalis II and Crysf, respectively. In addition, DeepCrystal attains an average improvement of 2.1, 6.0% for F-score, 1.9, 3.9% for accuracy and 3.8, 7.0% for MCC w.r.t. Crysalis II and Crysf on independent test sets.

Availability and implementation

The standalone source code and models are available at https://github.com/elbasir/DeepCrystal and a web-server is also available at https://deeplearning-protein.qcri.org.

Supplementary information

Supplementary data are available at Bioinformatics online.

Fourier-Transform Infrared Imaging Spectroscopy and Laser Ablation -ICPMS New Vistas for Biochemical Analyses of Ischemic Stroke in Rat Brain

Objective: Stroke is the main cause of adult disability in the world, leaving more than half of the patients dependent on daily assistance. Understanding the post-stroke biochemical and molecular changes are critical for patient survival and stroke management. The aim of this work was to investigate the photo-thrombotic ischemic stroke in male rats with particular focus on biochemical and elemental changes in the primary stroke lesion in the somatosensory cortex and surrounding areas, including the corpus callosum.

Materials and Methods: FT-IR imaging spectroscopy and LA-ICPMS techniques examined stroke brain samples, which were compared with standard immunohistochemistry studies.

Results: The FTIR results revealed that in the lesioned gray matter the relative distribution of lipid, lipid acyl and protein contents decreased significantly. Also at this locus, there was a significant increase in aggregated protein as detected by high-levels Aβ_1-42. Areas close to the stroke focus experienced decrease in the lipid and lipid acyl contents associated with an increase in lipid ester, olefin, and methyl bio-contents with a novel finding of Aβ_1-42 in the PL-GM and L-WM. Elemental analyses realized major changes in the different brain structures that may underscore functionality.

Conclusion: In conclusion, FTIR bio-spectroscopy is a non-destructive, rapid, and a refined technique to characterize oxidative stress markers associated with lipid degradation and protein denaturation not characterized by routine approaches. This technique may expedite research into stroke and offer new approaches for neurodegenerative disorders. The results suggest that a good therapeutic strategy should include a mechanism that provides protective effect from brain swelling (edema) and neurotoxicity by scavenging the lipid peroxidation end products.

RGBM: regularized gradient boosting machines for identification of the transcriptional regulators of discrete glioma subtypes

We propose a generic framework for gene regulatory network (GRN) inference approached as a feature selection problem. GRNs obtained using Machine Learning techniques are often dense, whereas real GRNs are rather sparse. We use a Tikonov regularization inspired optimal L-curve criterion that utilizes the edge weight distribution for a given target gene to determine the optimal set of TFs associated with it. Our proposed framework allows to incorporate a mechanistic active biding network based on cis-regulatory motif analysis. We evaluate our regularization framework in conjunction with two non-linear ML techniques, namely gradient boosting machines (GBM) and random-forests (GENIE), resulting in a regularized feature selection based method specifically called RGBM and RGENIE respectively. RGBM has been used to identify the main transcription factors that are causally involved as master regulators of the gene expression signature activated in the FGFR3-TACC3-positive glioblastoma. Here, we illustrate that RGBM identifies the main regulators of the molecular subtypes of brain tumors. Our analysis reveals the identity and corresponding biological activities of the master regulators characterizing the difference between G-CIMP-high and G-CIMP-low subtypes and between PA-like and LGm6-GBM, thus providing a clue to the yet undetermined nature of the transcriptional events among these subtypes.

Harnessing Qatar Biobank to understand type 2 diabetes and obesity in adult Qataris from the First Qatar Biobank Project

BACKGROUND

Human tissues are invaluable resources for researchers worldwide. Biobanks are repositories of such human tissues and can have a strategic importance for genetic research, clinical care, and future discoveries and treatments. One of the aims of Qatar Biobank is to improve the understanding and treatment of common diseases afflicting Qatari population such as obesity and diabetes.

METHODS

In this study we apply a panorama of state-of-the-art statistical methods and machine learning algorithms to investigate associations and risk factors for diabetes and obesity on a sample of 1000 Qatari population.

RESULTS

Regarding diabetes, we identified pronounced associations and risk factors in Qatari population including magnesium, chloride, c-peptide of insulin, insulin, and uric acid. Similarly, for obesity, significant associations and risk factors include insulin, c-peptide of insulin, albumin, and uric acid. Moreover, our study has revealed interactions of hypomagnesemia with HDL-C, triglycerides, and free thyroxine.

CONCLUSIONS

Our study strongly confirms known associations and risk factors associated with diabetes and obesity in Qatari population as previously found in other population studies in different parts of the world. Moreover, interactions of hypomagnesemia with other associations and risk factors merit further investigations.

Detection of statistically significant network changes in complex biological networks

Background

Biological networks contribute effectively to unveil the complex structure of molecular interactions and to discover driver genes especially in cancer context. It can happen that due to gene mutations, as for example when cancer progresses, the gene expression network undergoes some amount of localized re-wiring. The ability to detect statistical relevant changes in the interaction patterns induced by the progression of the disease can lead to the discovery of novel relevant signatures. Several procedures have been recently proposed to detect sub-network differences in pairwise labeled weighted networks.

Methods

In this paper, we propose an improvement over the state-of-the-art based on the Generalized Hamming Distance adopted for evaluating the topological difference between two networks and estimating its statistical significance. The proposed procedure exploits a more effective model selection criteria to generate p-values for statistical significance and is more efficient in terms of computational time and prediction accuracy than literature methods. Moreover, the structure of the proposed algorithm allows for a faster parallelized implementation.

Results

In the case of dense random geometric networks the proposed approach is 10-15x faster and achieves 5-10% higher AUC, Precision/Recall, and Kappa value than the state-of-the-art. We also report the application of the method to dissect the difference between the regulatory networks of IDH-mutant versus IDH-wild-type glioma cancer. In such a case our method is able to identify some recently reported master regulators as well as novel important candidates.

Conclusions

We show that our network differencing procedure can effectively and efficiently detect statistical significant network re-wirings in different conditions. When applied to detect the main differences between the networks of IDH-mutant and IDH-wild-type glioma tumors, it correctly selects sub-networks centered on important key regulators of these two different subtypes. In addition, its application highlights several novel candidates that cannot be detected by standard single network-based approaches.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-017-0412-6) contains supplementary material, which is available to authorized users.

COUSCOus: improved protein contact prediction using an empirical Bayes covariance estimator

Background

The post-genomic era with its wealth of sequences gave rise to a broad range of protein residue-residue contact detecting methods. Although various coevolution methods such as PSICOV, DCA and plmDCA provide correct contact predictions, they do not completely overlap. Hence, new approaches and improvements of existing methods are needed to motivate further development and progress in the field. We present a new contact detecting method, COUSCOus, by combining the best shrinkage approach, the empirical Bayes covariance estimator and GLasso.

Results

Using the original PSICOV benchmark dataset, COUSCOus achieves mean accuracies of 0.74, 0.62 and 0.55 for the top L/10 predicted long, medium and short range contacts, respectively. In addition, COUSCOus attains mean areas under the precision-recall curves of 0.25, 0.29 and 0.30 for long, medium and short contacts and outperforms PSICOV. We also observed that COUSCOus outperforms PSICOV w.r.t. Matthew’s correlation coefficient criterion on full list of residue contacts. Furthermore, COUSCOus achieves on average 10% more gain in prediction accuracy compared to PSICOV on an independent test set composed of CASP11 protein targets. Finally, we showed that when using a simple random forest meta-classifier, by combining contact detecting techniques and sequence derived features, PSICOV predictions should be replaced by the more accurate COUSCOus predictions.

Conclusion

We conclude that the consideration of superior covariance shrinkage approaches will boost several research fields that apply the GLasso procedure, amongst the presented one of residue-residue contact prediction as well as fields such as gene network reconstruction.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-016-1400-3) contains supplementary material, which is available to authorized users.

Very sparse LSSVM reductions for large-scale data

Least squares support vector machines (LSSVMs) have been widely applied for classification and regression with comparable performance with SVMs. The LSSVM model lacks sparsity and is unable to handle large-scale data due to computational and memory constraints. A primal fixed-size LSSVM (PFS-LSSVM) introduce sparsity using Nyström approximation with a set of prototype vectors (PVs). The PFS-LSSVM model solves an overdetermined system of linear equations in the primal. However, this solution is not the sparsest. We investigate the sparsity-error tradeoff by introducing a second level of sparsity. This is done by means of L0 -norm-based reductions by iteratively sparsifying LSSVM and PFS-LSSVM models. The exact choice of the cardinality for the initial PV set is not important then as the final model is highly sparse. The proposed method overcomes the problem of memory constraints and high computational costs resulting in highly sparse reductions to LSSVM models. The approximations of the two models allow to scale the models to large-scale datasets. Experiments on real-world classification and regression data sets from the UCI repository illustrate that these approaches achieve sparse models without a significant tradeoff in errors.

Multiclass semisupervised learning based upon kernel spectral clustering

This paper proposes a multiclass semisupervised learning algorithm by using kernel spectral clustering (KSC) as a core model. A regularized KSC is formulated to estimate the class memberships of data points in a semisupervised setting using the one-versus-all strategy while both labeled and unlabeled data points are present in the learning process. The propagation of the labels to a large amount of unlabeled data points is achieved by adding the regularization terms to the cost function of the KSC formulation. In other words, imposing the regularization term enforces certain desired memberships. The model is then obtained by solving a linear system in the dual. Furthermore, the optimal embedding dimension is designed for semisupervised clustering. This plays a key role when one deals with a large number of clusters.

Multilevel hierarchical kernel spectral clustering for real-life large scale complex networks

Kernel spectral clustering corresponds to a weighted kernel principal component analysis problem in a constrained optimization framework. The primal formulation leads to an eigen-decomposition of a centered Laplacian matrix at the dual level. The dual formulation allows to build a model on a representative subgraph of the large scale network in the training phase and the model parameters are estimated in the validation stage. The KSC model has a powerful out-of-sample extension property which allows cluster affiliation for the unseen nodes of the big data network. In this paper we exploit the structure of the projections in the eigenspace during the validation stage to automatically determine a set of increasing distance thresholds. We use these distance thresholds in the test phase to obtain multiple levels of hierarchy for the large scale network. The hierarchical structure in the network is determined in a bottom-up fashion. We empirically showcase that real-world networks have multilevel hierarchical organization which cannot be detected efficiently by several state-of-the-art large scale hierarchical community detection techniques like the Louvain, OSLOM and Infomap methods. We show that a major advantage of our proposed approach is the ability to locate good quality clusters at both the finer and coarser levels of hierarchy using internal cluster quality metrics on 7 real-life networks.