A network of transcriptomic signatures identifies novel comorbidity mechanisms between schizophrenia and somatic disorders.

The clinical burden of mental illness, in particular schizophrenia and bipolar disorder, are driven by frequent chronic courses and increased mortality, as well as the risk for comorbid conditions such as cardiovascular disease and type 2 diabetes. Evidence suggests an overlap of molecular pathways between psychotic disorders and somatic comorbidities. In this study, we developed a computational framework to perform comorbidity modeling via an improved integrative unsupervised machine learning approach based on multi-rank non-negative matrix factorization (mrNMF). Using this procedure, we extracted molecular signatures potentially explaining shared comorbidity mechanisms. For this, 27 case-control microarray transcriptomic datasets across multiple tissues were collected, covering three main categories of conditions including psychotic disorders, cardiovascular diseases and type II diabetes. We addressed the limitation of normal NMF for parameter selection by introducing multi-rank ensembled NMF to identify signatures under various hierarchical levels simultaneously. Analysis of comorbidity signature pairs was performed to identify several potential mechanisms involving activation of inflammatory response auxiliarily interconnecting angiogenesis, oxidative response and GABAergic neuro-action. Overall, we proposed a general cross-cohorts computing workflow for investigating the comorbid pattern across multiple symptoms, applied it to the real-data comorbidity study on schizophrenia, and further discussed the potential for future application of the approach.

Exploring the impact of dexamethasone on gene regulation in myeloma cells.

Among glucocorticoids (GCs), dexamethasone (Dex) is widely used in treatment of multiple myelomas. However, despite a definite benefit, all patients relapse. Moreover, the molecular basis of glucocorticoid efficacy remains elusive. To determine genomic response to Dex in myeloma cells, we generated bulk and single-cell multi-omics data and high-resolution contact maps of active enhancers and target genes. We show that a minority of glucocorticoid receptor-binding sites are associated with enhancer activity gains, increased interaction loops, and transcriptional activity. We identified and characterized a predominant enhancer enriched in cohesin (RAD21) and more accessible upon Dex exposure. Analysis of four gene-specific networks revealed the importance of the CTCF-cohesin couple and the synchronization of regulatory sequence openings for efficient transcription in response to Dex. Notably, these epigenomic changes are associated with cell-to-cell transcriptional heterogeneity, in particular, lineage-specific genes. As consequences, BCL2L11-encoding BIM critical for Dex-induced apoptosis and CXCR4 protective from chemotherapy-induced apoptosis are rather up-regulated in different cells. In summary, our work provides new insights into the molecular mechanisms involved in Dex escape.

Biologically informed variational autoencoders allow predictive modeling of genetic and drug-induced perturbations.

MOTIVATION: Variational autoencoders (VAEs) have rapidly increased in popularity in biological applications and have already successfully been used on many omic datasets. Their latent space provides a low-dimensional representation of input data, and VAEs have been applied, e.g. for clustering of single-cell transcriptomic data. However, due to their non-linear nature, the patterns that VAEs learn in the latent space remain obscure. Hence, the lower-dimensional data embedding cannot directly be related to input features. RESULTS: To shed light on the inner workings of VAE and enable direct interpretability of the model through its structure, we designed a novel VAE, OntoVAE (Ontology guided VAE) that can incorporate any ontology in its latent space and decoder part and, thus, provide pathway or phenotype activities for the ontology terms. In this work, we demonstrate that OntoVAE can be applied in the context of predictive modeling and show its ability to predict the effects of genetic or drug-induced perturbations using different ontologies and both, bulk and single-cell transcriptomic datasets. Finally, we provide a flexible framework, which can be easily adapted to any ontology and dataset. AVAILABILITY AND IMPLEMENTATION: OntoVAE is available as a python package under https://github.com/hdsu-bioquant/onto-vae.

High-Fidelity 3D Stray Magnetic Field Mapping of Smartphones to Address Safety Considerations with Active Implantable Electronic Medical Devices.

Case reports indicate that magnets in smartphones could be a source of electromagnetic interference (EMI) for active implantable medical devices (AIMD), which could lead to device malfunction, compromising patient safety. Recognizing this challenge, we implemented a high-fidelity 3D magnetic field mapping (spatial resolution 1 mm) setup using a three-axis Hall probe and teslameter, controlled by a robot (COSI Measure). With this setup, we examined the stray magnetic field of an iPhone 13 Pro, iPhone 12, and MagSafe charger to identify sources of magnetic fields for the accurate risk assessment of potential interferences with AIMDs. Our measurements revealed that the stray fields of the annular array of magnets, the wide-angle camera, and the speaker of the smartphones exceeded the 1 mT limit defined by ISO 14117:2019. Our data-driven safety recommendation is that an iPhone 13 Pro should be kept at least 25 mm away from an AIMD to protect it from unwanted EMI interactions. Our study addresses safety concerns due to potential device-device interactions between smartphones and AIMDs and will help to define data-driven safety guidelines. We encourage vendors of electronic consumer products (ECP) to provide information on the magnetic fields of their products and advocate for the inclusion of smartphones in the risk assessment of EMI with AIMDs.

Genomic profiling of HIV-1 integration in microglia cells links viral integration to the topologically associated domains.

HIV-1 encounters the hierarchically organized host chromatin to stably integrate and persist in anatomically distinct latent reservoirs. The contribution of genome organization in HIV-1 infection has been largely understudied across different HIV-1 targets. Here, we determine HIV-1 integration sites (ISs), associate them with chromatin and expression signatures at different genomic scales in a microglia cell model, and profile them together with the primary T cell reservoir. HIV-1 insertions into introns of actively transcribed genes with IS hotspots in genic and super-enhancers, characteristic of blood cells, are maintained in the microglia cell model. Genome organization analysis reveals dynamic CCCTC-binding factor (CTCF) clusters in cells with active and repressed HIV-1 transcription, whereas CTCF removal impairs viral integration. We identify CTCF-enriched topologically associated domain (TAD) boundaries with signatures of transcriptionally active chromatin as HIV-1 integration determinants in microglia and CD4+ T cells, highlighting the importance of host genome organization in HIV-1 infection.

Accelerated Simultaneous T2 and T2* Mapping of Multiple Sclerosis Lesions Using Compressed Sensing Reconstruction of Radial RARE-EPI MRI.

(1) Background: Radial RARE-EPI MRI facilitates simultaneous T2 and T2* mapping (2in1-RARE-EPI). With modest undersampling (R = 2), the speed gain of 2in1-RARE-EPI relative to Multi-Spin-Echo and Multi-Gradient-Recalled-Echo references is limited. Further reduction in scan time is crucial for clinical studies investigating T2 and T2* as imaging biomarkers. We demonstrate the feasibility of further acceleration, utilizing compressed sensing (CS) reconstruction of highly undersampled 2in1-RARE-EPI. (2) Methods: Two-fold radially-undersampled 2in1-RARE-EPI data from phantoms, healthy volunteers (n = 3), and multiple sclerosis patients (n = 4) were used as references, and undersampled (Rextra = 1-12, effective undersampling Reff = 2-24). For each echo time, images were reconstructed using CS-reconstruction. For T2 (RARE module) and T2* mapping (EPI module), a linear least-square fit was applied to the images. T2 and T2* from CS-reconstruction of undersampled data were benchmarked against values from CS-reconstruction of the reference data. (3) Results: We demonstrate accelerated simultaneous T2 and T2* mapping using undersampled 2in1-RARE-EPI with CS-reconstruction is feasible. For Rextra = 6 (TA = 01:39 min), the overall MAPE was ≤8% (T2*) and ≤4% (T2); for Rextra = 12 (TA = 01:06 min), the overall MAPE was <13% (T2*) and <5% (T2). (4) Conclusion: Substantial reductions in scan time are achievable for simultaneous T2 and T2* mapping of the brain using highly undersampled 2in1-RARE-EPI with CS-reconstruction.

Global hypomethylation in childhood asthma identified by genome-wide DNA-methylation sequencing preferentially affects enhancer regions.

BACKGROUND: Childhood asthma is a result of a complex interaction of genetic and environmental components causing epigenetic and immune dysregulation, airway inflammation and impaired lung function. Although different microarray based EWAS studies have been conducted, the impact of epigenetic regulation in asthma development is still widely unknown. We have therefore applied unbiased whole genome bisulfite sequencing (WGBS) to characterize global DNA-methylation profiles of asthmatic children compared to healthy controls. METHODS: Peripheral blood samples of 40 asthmatic and 42 control children aged 5-15 years from three birth cohorts were sequenced together with paired cord blood samples. Identified differentially methylated regions (DMRs) were categorized in genotype-associated, cell-type-dependent, or prenatally primed. Network analysis and subsequent natural language processing of DMR-associated genes was complemented by targeted analysis of functional translation of epigenetic regulation on the transcriptional and protein level. RESULTS: In total, 158 DMRs were identified in asthmatic children compared to controls of which 37% were related to the eosinophil content. A global hypomethylation was identified affecting predominantly enhancer regions and regulating key immune genes such as IL4, IL5RA, and EPX. These DMRs were confirmed in n = 267 samples and could be linked to aberrant gene expression. Out of the 158 DMRs identified in the established phenotype, 56 were perturbed already at birth and linked, at least in part, to prenatal influences such as tobacco smoke exposure or phthalate exposure. CONCLUSION: This is the first epigenetic study based on whole genome sequencing to identify marked dysregulation of enhancer regions as a hallmark of childhood asthma.

dsMTL: a computational framework for privacy-preserving, distributed multi-task machine learning.

MOTIVATION: In multi-cohort machine learning studies, it is critical to differentiate between effects that are reproducible across cohorts and those that are cohort-specific. Multi-task learning (MTL) is a machine learning approach that facilitates this differentiation through the simultaneous learning of prediction tasks across cohorts. Since multi-cohort data can often not be combined into a single storage solution, there would be the substantial utility of an MTL application for geographically distributed data sources. RESULTS: Here, we describe the development of 'dsMTL', a computational framework for privacy-preserving, distributed multi-task machine learning that includes three supervised and one unsupervised algorithms. First, we derive the theoretical properties of these methods and the relevant machine learning workflows to ensure the validity of the software implementation. Second, we implement dsMTL as a library for the R programming language, building on the DataSHIELD platform that supports the federated analysis of sensitive individual-level data. Third, we demonstrate the applicability of dsMTL for comorbidity modeling in distributed data. We show that comorbidity modeling using dsMTL outperformed conventional, federated machine learning, as well as the aggregation of multiple models built on the distributed datasets individually. The application of dsMTL was computationally efficient and highly scalable when applied to moderate-size (n < 500), real expression data given the actual network latency. AVAILABILITY AND IMPLEMENTATION: dsMTL is freely available at https://github.com/transbioZI/dsMTLBase (server-side package) and https://github.com/transbioZI/dsMTLClient (client-side package). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

MYCN mediates cysteine addiction and sensitizes neuroblastoma to ferroptosis.

Aberrant expression of MYC transcription factor family members predicts poor clinical outcome in many human cancers. Oncogenic MYC profoundly alters metabolism and mediates an antioxidant response to maintain redox balance. Here we show that MYCN induces massive lipid peroxidation on depletion of cysteine, the rate-limiting amino acid for glutathione (GSH) biosynthesis, and sensitizes cells to ferroptosis, an oxidative, non-apoptotic and iron-dependent type of cell death. The high cysteine demand of MYCN-amplified childhood neuroblastoma is met by uptake and transsulfuration. When uptake is limited, cysteine usage for protein synthesis is maintained at the expense of GSH triggering ferroptosis and potentially contributing to spontaneous tumor regression in low-risk neuroblastomas. Pharmacological inhibition of both cystine uptake and transsulfuration combined with GPX4 inactivation resulted in tumor remission in an orthotopic MYCN-amplified neuroblastoma model. These findings provide a proof of concept of combining multiple ferroptosis targets as a promising therapeutic strategy for aggressive MYCN-amplified tumors.

The DNA methylation landscape of multiple myeloma shows extensive inter- and intrapatient heterogeneity that fuels transcriptomic variability.

BACKGROUND: Cancer evolution depends on epigenetic and genetic diversity. Historically, in multiple myeloma (MM), subclonal diversity and tumor evolution have been investigated mostly from a genetic perspective. METHODS: Here, we performed an analysis of 42 MM samples from 21 patients by using enhanced reduced representation bisulfite sequencing (eRRBS). We combined several metrics of epigenetic heterogeneity to analyze DNA methylation heterogeneity in MM patients. RESULTS: We show that MM is characterized by the continuous accumulation of stochastic methylation at the promoters of development-related genes. High combinatorial entropy change is associated with poor outcomes in our pilot study and depends predominantly on partially methylated domains (PMDs). These PMDs, which represent the major source of inter- and intrapatient DNA methylation heterogeneity in MM, are linked to other key epigenetic aberrations, such as CpG island (CGI)/transcription start site (TSS) hypermethylation and H3K27me3 redistribution as well as 3D organization alterations. In addition, transcriptome analysis revealed that intratumor methylation heterogeneity was associated with low-level expression and high variability. CONCLUSIONS: We propose that disrupted DNA methylation in MM is responsible for high epigenetic and transcriptomic instability allowing tumor cells to adapt to environmental changes by tapping into a pool of evolutionary trajectories.

The endogenous cellular protease inhibitor SPINT2 controls SARS-CoV-2 viral infection and is associated to disease severity.

COVID-19 outbreak is the biggest threat to human health in recent history. Currently, there are over 1.5 million related deaths and 75 million people infected around the world (as of 22/12/2020). The identification of virulence factors which determine disease susceptibility and severity in different cell types remains an essential challenge. The serine protease TMPRSS2 has been shown to be important for S protein priming and viral entry, however, little is known about its regulation. SPINT2 is a member of the family of Kunitz type serine protease inhibitors and has been shown to inhibit TMPRSS2. Here, we explored the existence of a co-regulation between SPINT2/TMPRSS2 and found a tightly regulated protease/inhibitor expression balance across tissues. We found that SPINT2 negatively correlates with SARS-CoV-2 expression in Calu-3 and Caco-2 cell lines and was down-regulated in secretory cells from COVID-19 patients. We validated our findings using Calu-3 cell lines and observed a strong increase in viral load after SPINT2 knockdown, while overexpression lead to a drastic reduction of the viral load. Additionally, we evaluated the expression of SPINT2 in datasets from comorbid diseases using bulk and scRNA-seq data. We observed its down-regulation in colon, kidney and liver tumors as well as in alpha pancreatic islets cells from diabetes Type 2 patients, which could have implications for the observed comorbidities in COVID-19 patients suffering from chronic diseases.

Integrative Ranking of Enhancer Networks Facilitates the Discovery of Epigenetic Markers in Cancer.

Regulation of gene expression through multiple epigenetic components is a highly combinatorial process. Alterations in any of these layers, as is commonly found in cancer diseases, can lead to a cascade of downstream effects on tumor suppressor or oncogenes. Hence, deciphering the effects of epigenetic alterations on regulatory elements requires innovative computational approaches that can benefit from the huge amounts of epigenomic datasets that are available from multiple consortia, such as Roadmap or BluePrint. We developed a software tool named IRENE (Integrative Ranking of Epigenetic Network of Enhancers), which performs quantitative analyses on differential epigenetic modifications through an integrated, network-based approach. The method takes into account the additive effect of alterations on multiple regulatory elements of a gene. Applying this tool to well-characterized test cases, it successfully found many known cancer genes from publicly available cancer epigenome datasets.

Simultaneous T2 and T2∗ mapping of multiple sclerosis lesions with radial RARE-EPI.

PURPOSE: The characteristic MRI features of multiple sclerosis (MS) lesions make it conceptually appealing to pursue parametric mapping techniques that support simultaneous generation of quantitative maps of 2 or more MR contrast mechanisms. We present a modular rapid acquisition with relaxation enhancement (RARE)-EPI hybrid that facilitates simultaneous T2 and T2∗ mapping (2in1-RARE-EPI). METHODS: In 2in1-RARE-EPI the first echoes in the echo train are acquired with a RARE module, later echoes are acquired with an EPI module. To define the fraction of echoes covered by the RARE and EPI module, an error analysis of T2 and T2∗ was conducted with Monte Carlo simulations. Radial k-space (under)sampling was implemented for acceleration (R = 2). The feasibility of 2in1-RARE-EPI for simultaneous T2 and T2∗ mapping was examined in a phantom study mimicking T2 and T2∗ relaxation times of the brain. For validation, 2in1-RARE-EPI was benchmarked versus multi spin-echo (MSE) and multi gradient-echo (MGRE) techniques. The clinical applicability of 2in1-RARE-EPI was demonstrated in healthy subjects and MS patients. RESULTS: There was a good agreement between T2 / T2∗ values derived from 2in1-RARE-EPI and T2 / T2∗ reference values obtained from MSE and MGRE in both phantoms and healthy subjects. In patients, MS lesions in T2 and T2∗ maps deduced from 2in1-RARE-EPI could be just as clearly delineated as in reference maps calculated from MSE/MGRE. CONCLUSION: This work demonstrates the feasibility of radially (under)sampled 2in1-RARE-EPI for simultaneous T2 and T2∗ mapping in MS patients.

Single-cell analyses reveal SARS-CoV-2 interference with intrinsic immune response in the human gut.

Exacerbated pro-inflammatory immune response contributes to COVID-19 pathology. However, despite the mounting evidence about SARS-CoV-2 infecting the human gut, little is known about the antiviral programs triggered in this organ. To address this gap, we performed single-cell transcriptomics of SARS-CoV-2-infected intestinal organoids. We identified a subpopulation of enterocytes as the prime target of SARS-CoV-2 and, interestingly, found the lack of positive correlation between susceptibility to infection and the expression of ACE2. Infected cells activated strong pro-inflammatory programs and produced interferon, while expression of interferon-stimulated genes was limited to bystander cells due to SARS-CoV-2 suppressing the autocrine action of interferon. These findings reveal that SARS-CoV-2 curtails the immune response and highlights the gut as a pro-inflammatory reservoir that should be considered to fully understand SARS-CoV-2 pathogenesis.

Identification of therapeutic targets of the hijacked super-enhancer complex in EVI1-rearranged leukemia.

Deregulation of the EVI1 proto-oncogene by the GATA2 distal hematopoietic enhancer (G2DHE) is a key event in high-risk acute myeloid leukemia carrying 3q21q26 aberrations (3q-AML). Upon chromosomal rearrangement, G2DHE acquires characteristics of a super-enhancer and causes overexpression of EVI1 at 3q26.2. However, the transcription factor (TF) complex of G2DHE remains poorly characterized. The aim of this study was to unravel key components of G2DHE-bound TFs involved in the deregulation of EVI1. We have identified several CEBPA and RUNX1 binding sites to be enriched and critical for G2DHE function in 3q-AML cells. Using ChIP-SICAP (ChIP followed by selective isolation of chromatin-associated proteins), a panel of chromatin interactors of RUNX1 and CEBPA were detected in 3q-AML, including PARP1 and IKZF1. PARP1 inhibition (PARPi) caused a reduction of EVI1 expression and a decrease in EVI1-G2DHE interaction frequency, highlighting the involvement of PARP1 in oncogenic super-enhancer formation. Furthermore, 3q-AML cells were highly sensitive to PARPi and displayed morphological changes with higher rates of differentiation and apoptosis as well as depletion of CD34 + cells. In summary, integrative analysis of the 3q-AML super-enhancer complex identified CEBPA and RUNX1 associated proteins and nominated PARP1 as a potential new therapeutic target in EVI1 + 3q-AML.

Single-cell transcriptomic analyses provide insights into the developmental origins of neuroblastoma.

Neuroblastoma is a pediatric tumor of the developing sympathetic nervous system. However, the cellular origin of neuroblastoma has yet to be defined. Here we studied the single-cell transcriptomes of neuroblastomas and normal human developing adrenal glands at various stages of embryonic and fetal development. We defined normal differentiation trajectories from Schwann cell precursors over intermediate states to neuroblasts or chromaffin cells and showed that neuroblastomas transcriptionally resemble normal fetal adrenal neuroblasts. Importantly, neuroblastomas with varying clinical phenotypes matched different temporal states along normal neuroblast differentiation trajectories, with the degree of differentiation corresponding to clinical prognosis. Our work highlights the roles of oncogenic MYCN and loss of TFAP2B in blocking differentiation and may provide the basis for designing therapeutic interventions to overcome differentiation blocks.

Kalirin-RAC controls nucleokinetic migration in ADRN-type neuroblastoma.

The migrational propensity of neuroblastoma is affected by cell identity, but the mechanisms behind the divergence remain unknown. Using RNAi and time-lapse imaging, we show that ADRN-type NB cells exhibit RAC1- and kalirin-dependent nucleokinetic (NUC) migration that relies on several integral components of neuronal migration. Inhibition of NUC migration by RAC1 and kalirin-GEF1 inhibitors occurs without hampering cell proliferation and ADRN identity. Using three clinically relevant expression dichotomies, we reveal that most of up-regulated mRNAs in RAC1- and kalirin-GEF1-suppressed ADRN-type NB cells are associated with low-risk characteristics. The computational analysis shows that, in a context of overall gene set poverty, the upregulomes in RAC1- and kalirin-GEF1-suppressed ADRN-type cells are a batch of AU-rich element-containing mRNAs, which suggests a link between NUC migration and mRNA stability. Gene set enrichment analysis-based search for vulnerabilities reveals prospective weak points in RAC1- and kalirin-GEF1-suppressed ADRN-type NB cells, including activities of H3K27- and DNA methyltransferases. Altogether, these data support the introduction of NUC inhibitors into cancer treatment research.

Magnetic Resonance Imaging of Multiple Sclerosis at 7.0 Tesla.

The overall goal of this article is to demonstrate a state-of-the-art ultrahigh field (UHF) magnetic resonance (MR) protocol of the brain at 7.0 Tesla in multiple sclerosis (MS) patients. MS is a chronic inflammatory, demyelinating, neurodegenerative disease that is characterized by white and gray matter lesions. Detection of spatially and temporally disseminated T2-hyperintense lesions by the use of MRI at 1.5 T and 3 T represents a crucial diagnostic tool in clinical practice to establish accurate diagnosis of MS based on the current version of the 2017 McDonald criteria. However, the differentiation of MS lesions from brain white matter lesions of other origins can sometimes be challenging due to their resembling morphology at lower magnetic field strengths (typically 3 T). Ultrahigh field MR (UHF-MR) benefits from increased signal-to-noise ratio and enhanced spatial resolution, both key to superior imaging for more accurate and definitive diagnoses of subtle lesions. Hence, MRI at 7.0 T has shown encouraging results to overcome the challenges of MS differential diagnosis by providing MS-specific neuroimaging markers (e.g., central vein sign, hypointense rim structures and differentiation of MS grey matter lesions). These markers and others can be identified by other MR contrasts other than T1 and T2 (T2*, phase, diffusion) and substantially improve the differentiation of MS lesions from those occurring in other neuroinflammatory conditions such as neuromyelitis optica and Susac syndrome. In this article, we describe our current technical approach to study cerebral white and grey matter lesions in MS patients at 7.0 T using different MR acquisition methods. The up-to-date protocol includes the preparation of the MR setup including the radio-frequency coils customized for UHF-MR, standardized screening, safety and interview procedures with MS patients, patient positioning in the MR scanner and acquisition of dedicated brain scans tailored for examining MS.