A public resource of single cell transcriptomes and multiscale networks from persons with and without Alzheimer's disease

The emergence of technologies that can support high-throughput profiling of single cell transcriptomes offers to revolutionize the study of brain tissue from persons with and without Alzheimer's disease (AD). Integration of these data with additional complementary multiomics data such as genetics, proteomics and clinical data provides powerful opportunities to link observed cell subpopulations and molecular network features within a broader disease-relevant context. We report here single nucleus RNA sequencing (snRNA-seq) profiles generated from superior frontal gyrus cortical tissue samples from 101 exceptionally well characterized, aged subjects from the Banner Brain and Body Donation Program in combination with whole genome sequences. We report findings that link common AD risk variants with CR1 expression in oligodendrocytes as well as alterations in peripheral hematological lab parameters, with these observations replicated in an independent, prospective cohort study of ageing and dementia. We also observed an AD-associated CD83(+) microglial subtype with unique molecular networks that encompass many known regulators of AD-relevant microglial biology, and which are associated with immunoglobulin IgG4 production in the transverse colon. These findings illustrate the power of multi-tissue molecular profiling to contextualize snRNA-seq brain transcriptomics and reveal novel disease biology. The transcriptomic, genetic, phenotypic, and network data resources described within this study are available for access and utilization by the scientific community.

Validation of a Transcriptome-Based Assay for Classifying Cancers of Unknown Primary Origin

INTRODUCTION

Cancers assume a variety of distinct histologies, and may originate from a myriad of sites including solid organs, hematopoietic cells, and connective tissue. Clinical decision-making based on consensus guidelines such as the National Comprehensive Cancer Network (NCCN) is often predicated on a specific histologic and anatomic diagnosis, supported by clinical features and pathologist interpretation of morphology and immunohistochemical (IHC) staining patterns. However, in patients with nonspecific morphologic and IHC findings-in addition to ambiguous clinical presentations such as recurrence versus new primary-a definitive diagnosis may not be possible, resulting in the patient being categorized as having a cancer of unknown primary (CUP). Therapeutic options and clinical outcomes are poor for patients with CUP, with a median survival of 8-11 months.

METHODS

Here, we describe and validate the Tempus Tumor Origin (Tempus TO) assay, an RNA-sequencing-based machine learning classifier capable of discriminating between 68 clinically relevant cancer subtypes. Model accuracy was assessed using primary and/or metastatic samples with known subtype.

RESULTS

We show that the Tempus TO model is 91% accurate when assessed on both a retrospectively held out cohort and a set of samples sequenced after model freeze that collectively contained 9210 total samples with known diagnoses. When evaluated on a cohort of CUPs, the model recapitulated established associations between genomic alterations and cancer subtype.

DISCUSSION

Combining diagnostic prediction tests (e.g., Tempus TO) with sequencing-based variant reporting (e.g., Tempus xT) may expand therapeutic options for patients with cancers of unknown primary or uncertain histology.

Saracatinib, a Selective Src Kinase Inhibitor, Blocks Fibrotic Responses in Preclinical Models of Pulmonary Fibrosis

Rationale: Idiopathic pulmonary fibrosis (IPF) is a chronic, progressive, and often fatal disorder. Two U.S. Food and Drug Administration-approved antifibrotic drugs, nintedanib and pirfenidone, slow the rate of decline in lung function, but responses are variable and side effects are common. Objectives: Using an in silico data-driven approach, we identified a robust connection between the transcriptomic perturbations in IPF disease and those induced by saracatinib, a selective Src kinase inhibitor originally developed for oncological indications. Based on these observations, we hypothesized that saracatinib would be effective at attenuating pulmonary fibrosis. Methods: We investigated the antifibrotic efficacy of saracatinib relative to nintedanib and pirfenidone in three preclinical models: 1) in vitro in normal human lung fibroblasts; 2) in vivo in bleomycin and recombinant Ad-TGF-β (adenovirus transforming growth factor-β) murine models of pulmonary fibrosis; and 3) ex vivo in mice and human precision-cut lung slices from these two murine models as well as patients with IPF and healthy donors. Measurements and Main Results: In each model, the effectiveness of saracatinib in blocking fibrogenic responses was equal or superior to nintedanib and pirfenidone. Transcriptomic analyses of TGF-β-stimulated normal human lung fibroblasts identified specific gene sets associated with fibrosis, including epithelial-mesenchymal transition, TGF-β, and WNT signaling that was uniquely altered by saracatinib. Transcriptomic analysis of whole-lung extracts from the two animal models of pulmonary fibrosis revealed that saracatinib reverted many fibrogenic pathways, including epithelial-mesenchymal transition, immune responses, and extracellular matrix organization. Amelioration of fibrosis and inflammatory cascades in human precision-cut lung slices confirmed the potential therapeutic efficacy of saracatinib in human lung fibrosis. Conclusions: These studies identify novel Src-dependent fibrogenic pathways and support the study of the therapeutic effectiveness of saracatinib in IPF treatment.

Quantifying the impact of immune history and variant on SARS-CoV-2 viral kinetics and infection rebound: A retrospective cohort study

Background

The combined impact of immunity and SARS-CoV-2 variants on viral kinetics during infections has been unclear.

Methods

We characterized 1,280 infections from the National Basketball Association occupational health cohort identified between June 2020 and January 2022 using serial RT-qPCR testing. Logistic regression and semi-mechanistic viral RNA kinetics models were used to quantify the effect of age, variant, symptom status, infection history, vaccination status and antibody titer to the founder SARS-CoV-2 strain on the duration of potential infectiousness and overall viral kinetics. The frequency of viral rebounds was quantified under multiple cycle threshold (Ct) value-based definitions.

Results

Among individuals detected partway through their infection, 51.0% (95% credible interval [CrI]: 48.3-53.6%) remained potentially infectious (Ct <30) 5 days post detection, with small differences across variants and vaccination status. Only seven viral rebounds (0.7%; N=999) were observed, with rebound defined as 3+days with Ct <30 following an initial clearance of 3+days with Ct ≥30. High antibody titers against the founder SARS-CoV-2 strain predicted lower peak viral loads and shorter durations of infection. Among Omicron BA.1 infections, boosted individuals had lower pre-booster antibody titers and longer clearance times than non-boosted individuals.

Conclusions

SARS-CoV-2 viral kinetics are partly determined by immunity and variant but dominated by individual-level variation. Since booster vaccination protects against infection, longer clearance times for BA.1-infected, boosted individuals may reflect a less effective immune response, more common in older individuals, that increases infection risk and reduces viral RNA clearance rate. The shifting landscape of viral kinetics underscores the need for continued monitoring to optimize isolation policies and to contextualize the health impacts of therapeutics and vaccines.

Funding

Supported in part by CDC contract #200-2016-91779, a sponsored research agreement to Yale University from the National Basketball Association contract #21-003529, and the National Basketball Players Association.

A Wipe-Based Stool Collection and Preservation Kit for Microbiome Community Profiling

While a range of methods for stool collection exist, many require complicated, self-directed protocols and stool transfer. In this study, we introduce and validate a novel, wipe-based approach to fecal sample collection and stabilization for metagenomics analysis. A total of 72 samples were collected across four different preservation types: freezing at -20°C, room temperature storage, a commercial DNA preservation kit, and a dissolvable wipe used with DESS (dimethyl sulfoxide, ethylenediaminetetraacetic acid, sodium chloride) solution. These samples were sequenced and analyzed for taxonomic abundance metrics, bacterial metabolic pathway classification, and diversity analysis. Overall, the DESS wipe results validated the use of a wipe-based capture method to collect stool samples for microbiome analysis, showing an R2 of 0.96 for species across all kingdoms, as well as exhibiting a maintenance of Shannon diversity (3.1-3.3) and species richness (151-159) compared to frozen samples. Moreover, DESS showed comparable performance to the commercially available preservation kit (R2 of 0.98), and samples consistently clustered by subject across each method. These data support that the DESS wipe method can be used for stable, room temperature collection and transport of human stool specimens.

Deep learning-based brain transcriptomic signatures associated with the neuropathological and clinical severity of Alzheimer's disease

Brain tissue gene expression from donors with and without Alzheimer's disease has been used to help inform the molecular changes associated with the development and potential treatment of this disorder. Here, we use a deep learning method to analyse RNA-seq data from 1114 brain donors from the Accelerating Medicines Project for Alzheimer's Disease consortium to characterize post-mortem brain transcriptome signatures associated with amyloid-β plaque, tau neurofibrillary tangles and clinical severity in multiple Alzheimer's disease dementia populations. Starting from the cross-sectional data in the Religious Orders Study and Memory and Aging Project cohort (n = 634), a deep learning framework was built to obtain a trajectory that mirrors Alzheimer's disease progression. A severity index was defined to quantitatively measure the progression based on the trajectory. Network analysis was then carried out to identify key gene (index gene) modules present in the model underlying the progression. Within this data set, severity indexes were found to be very closely correlated with all Alzheimer's disease neuropathology biomarkers (R ∼ 0.5, P < 1e-11) and global cognitive function (R = -0.68, P < 2.2e-16). We then applied the model to additional transcriptomic data sets from different brain regions (MAYO, n = 266; Mount Sinai Brain Bank, n = 214), and observed that the model remained significantly predictive (P < 1e-3) of neuropathology and clinical severity. The index genes that significantly contributed to the model were integrated with Alzheimer's disease co-expression regulatory networks, resolving four discrete gene modules that are implicated in vascular and metabolic dysfunction in different cell types, respectively. Our work demonstrates the generalizability of this signature to frontal and temporal cortex measurements and additional brain donors with Alzheimer's disease, other age-related neurological disorders and controls, and revealed that the transcriptomic network modules contribute to neuropathological and clinical disease severity. This study illustrates the promise of using deep learning methods to analyse heterogeneous omics data and discover potentially targetable molecular networks that can inform the development, treatment and prevention of neurodegenerative diseases like Alzheimer's disease.

Predictive Modelling of Susceptibility to Substance Abuse, Mortality and Drug-Drug Interactions in Opioid Patients

Objective: Opioids are a class of drugs that are known for their use as pain relievers. They bind to opioid receptors on nerve cells in the brain and the nervous system to mitigate pain. Addiction is one of the chronic and primary adverse events of prolonged usage of opioids. They may also cause psychological disorders, muscle pain, depression, anxiety attacks etc. In this study, we present a collection of predictive models to identify patients at risk of opioid abuse and mortality by using their prescription histories. Also, we discover particularly threatening drug-drug interactions in the context of opioid usage.

Methods and Materials: Using a publicly available dataset from MIMIC-III, two models were trained, Logistic Regression with L2 regularization (baseline) and Extreme Gradient Boosting (enhanced model), to classify the patients of interest into two categories based on their susceptibility to opioid abuse. We’ve also used K-Means clustering, an unsupervised algorithm, to explore drug-drug interactions that might be of concern.

Results: The baseline model for classifying patients susceptible to opioid abuse has an F1 score of 76.64% (accuracy 77.16%) while the enhanced model has an F1 score of 94.45% (accuracy 94.35%). These models can be used as a preliminary step towards inferring the causal effect of opioid usage and can help monitor the prescription practices to minimize the opioid abuse.

Discussion and Conclusion: Results suggest that the enhanced model provides a promising approach in preemptive identification of patients at risk for opioid abuse. By discovering and correlating the patterns contributing to opioid overdose or abuse among a variety of patients, machine learning models can be used as an efficient tool to help uncover the existing gaps and/or fraudulent practices in prescription writing. To quote an example of one such incidental finding, our study discovered that insulin might possibly be interacting with opioids in an unfavourable way leading to complications in diabetic patients. This indicates that diabetic patients under long term opioid usage might need to take increased amounts of insulin to make it more effective. This observation backs up prior research studies done on a similar aspect. To increase the translational value of our work, the predictive models and the associated software code are made available under the MIT License.

Patient similarity network of newly diagnosed multiple myeloma identifies patient subgroups with distinct genetic features and clinical implications

The remarkable genetic heterogeneity of multiple myeloma poses a substantial challenge for proper prognostication and clinical management of patients. Here, we introduce MM-PSN, the first multiomics patient similarity network of myeloma. MM-PSN enabled accurate dissection of the genetic and molecular landscape of the disease and determined 12 distinct subgroups defined by five data types generated from genomic and transcriptomic profiling of 655 patients. MM-PSN identified patient subgroups not previously described defined by specific patterns of alterations, enriched for specific gene vulnerabilities, and associated with potential therapeutic options. Our analysis revealed that co-occurrence of t(4;14) and 1q gain identified patients at significantly higher risk of relapse and shorter survival as compared to t(4;14) as a single lesion. Furthermore, our results show that 1q gain is the most important single lesion conferring high risk of relapse and that it can improve on the current International Staging Systems (ISS and R-ISS).

Repositioning of a novel GABA-B receptor agonist, AZD3355 (Lesogaberan), for the treatment of non-alcoholic steatohepatitis

Non-alcoholic steatohepatitis (NASH) is a rising health challenge, with no approved drugs. We used a computational drug repositioning strategy to uncover a novel therapy for NASH, identifying a GABA-B receptor agonist, AZD3355 (Lesogaberan) previously evaluated as a therapy for esophageal reflux. AZD3355's potential efficacy in NASH was tested in human stellate cells, human precision cut liver slices (hPCLS), and in vivo in a well-validated murine model of NASH. In human stellate cells AZD3355 significantly downregulated profibrotic gene and protein expression. Transcriptomic analysis of these responses identified key regulatory nodes impacted by AZD3355, including Myc, as well as MAP and ERK kinases. In PCLS, AZD3355 down-regulated collagen1α1, αSMA and TNF-α mRNAs as well as secreted collagen1α1. In vivo, the drug significantly improved histology, profibrogenic gene expression, and tumor development, which was comparable to activity of obeticholic acid in a robust mouse model of NASH, but awaits further testing to determine its relative efficacy in patients. These data identify a well-tolerated clinical stage asset as a novel candidate therapy for human NASH through its hepatoprotective, anti-inflammatory and antifibrotic mechanisms of action. The approach validates computational methods to identify novel therapies in NASH in uncovering new pathways of disease development that can be rapidly translated into clinical trials.

Phe2vec: Automated disease phenotyping based on unsupervised embeddings from electronic health records

Robust phenotyping of patients from electronic health records (EHRs) at scale is a challenge in clinical informatics. Here, we introduce Phe2vec, an automated framework for disease phenotyping from EHRs based on unsupervised learning and assess its effectiveness against standard rule-based algorithms from Phenotype KnowledgeBase (PheKB). Phe2vec is based on pre-computing embeddings of medical concepts and patients' clinical history. Disease phenotypes are then derived from a seed concept and its neighbors in the embedding space. Patients are linked to a disease if their embedded representation is close to the disease phenotype. Comparing Phe2vec and PheKB cohorts head-to-head using chart review, Phe2vec performed on par or better in nine out of ten diseases. Differently from other approaches, it can scale to any condition and was validated against widely adopted expert-based standards. Phe2vec aims to optimize clinical informatics research by augmenting current frameworks to characterize patients by condition and derive reliable disease cohorts.

Integrative chemogenomic analysis identifies small molecules that partially rescue ΔF508-CFTR for cystic fibrosis

Rare diseases affect 10% of the first‐world population, yet over 95% lack even a single pharmaceutical treatment. In the present age of information, we need ways to leverage our vast data and knowledge to streamline therapeutic development and lessen this gap. Here, we develop and implement an innovative informatic approach to identify therapeutic molecules, using the Connectivity Map and LINCS L1000 databases and disease‐associated transcriptional signatures and pathways. We apply this to cystic fibrosis (CF), the most common genetic disease in people of northern European ancestry leading to chronic lung disease and reduced lifespan. We selected and tested 120 small molecules in a CF cell line, finding 8 with activity, and confirmed 3 in primary CF airway epithelia. Although chemically diverse, the transcriptional profiles of the hits suggest a common mechanism associated with the unfolded protein response and/or TNFα signaling. This study highlights the power of informatics to help identify new therapies and reveal mechanistic insights while moving beyond target‐centric drug discovery.

Multiplex qPCR discriminates variants of concern to enhance global surveillance of SARS-CoV-2

With the emergence of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) variants that may increase transmissibility and/or cause escape from immune responses, there is an urgent need for the targeted surveillance of circulating lineages. It was found that the B.1.1.7 (also 501Y.V1) variant, first detected in the United Kingdom, could be serendipitously detected by the Thermo Fisher TaqPath COVID-19 PCR assay because a key deletion in these viruses, spike Δ69-70, would cause a "spike gene target failure" (SGTF) result. However, a SGTF result is not definitive for B.1.1.7, and this assay cannot detect other variants of concern (VOC) that lack spike Δ69-70, such as B.1.351 (also 501Y.V2), detected in South Africa, and P.1 (also 501Y.V3), recently detected in Brazil. We identified a deletion in the ORF1a gene (ORF1a Δ3675-3677) in all 3 variants, which has not yet been widely detected in other SARS-CoV-2 lineages. Using ORF1a Δ3675-3677 as the primary target and spike Δ69-70 to differentiate, we designed and validated an open-source PCR assay to detect SARS-CoV-2 VOC. Our assay can be rapidly deployed in laboratories around the world to enhance surveillance for the local emergence and spread of B.1.1.7, B.1.351, and P.1.

Early introductions and transmission of SARS-CoV-2 variant B.1.1.7 in the United States

The emergence and spread of SARS-CoV-2 lineage B.1.1.7, first detected in the United Kingdom, has become a global public health concern because of its increased transmissibility. Over 2,500 COVID-19 cases associated with this variant have been detected in the United States (US) since December 2020, but the extent of establishment is relatively unknown. Using travel, genomic, and diagnostic data, we highlight that the primary ports of entry for B.1.1.7 in the US were in New York, California, and Florida. Furthermore, we found evidence for many independent B.1.1.7 establishments starting in early December 2020, followed by interstate spread by the end of the month. Finally, we project that B.1.1.7 will be the dominant lineage in many states by mid- to late March. Thus, genomic surveillance for B.1.1.7 and other variants urgently needs to be enhanced to better inform the public health response.

PCR assay to enhance global surveillance for SARS-CoV-2 variants of concern

With the emergence of SARS-CoV-2 variants that may increase transmissibility and/or cause escape from immune responses 1-3 , there is an urgent need for the targeted surveillance of circulating lineages. It was found that the B.1.1.7 (also 501Y.V1) variant first detected in the UK 4,5 could be serendipitously detected by the ThermoFisher TaqPath COVID-19 PCR assay because a key deletion in these viruses, spike Δ69-70, would cause a "spike gene target failure" (SGTF) result. However, a SGTF result is not definitive for B.1.1.7, and this assay cannot detect other variants of concern that lack spike Δ69-70, such as B.1.351 (also 501Y.V2) detected in South Africa 6 and P.1 (also 501Y.V3) recently detected in Brazil 7 . We identified a deletion in the ORF1a gene (ORF1a Δ3675-3677) in all three variants, which has not yet been widely detected in other SARS-CoV-2 lineages. Using ORF1a Δ3675-3677 as the primary target and spike Δ69-70 to differentiate, we designed and validated an open source PCR assay to detect SARS-CoV-2 variants of concern 8 . Our assay can be rapidly deployed in laboratories around the world to enhance surveillance for the local emergence spread of B.1.1.7, B.1.351, and P.1.

Early introductions and community transmission of SARS-CoV-2 variant B.1.1.7 in the United States

The emergence and spread of SARS-CoV-2 lineage B.1.1.7, first detected in the United Kingdom, has become a global public health concern because of its increased transmissibility. Over 2500 COVID-19 cases associated with this variant have been detected in the US since December 2020, but the extent of establishment is relatively unknown. Using travel, genomic, and diagnostic data, we highlight the primary ports of entry for B.1.1.7 in the US and locations of possible underreporting of B.1.1.7 cases. Furthermore, we found evidence for many independent B.1.1.7 establishments starting in early December 2020, followed by interstate spread by the end of the month. Finally, we project that B.1.1.7 will be the dominant lineage in many states by mid to late March. Thus, genomic surveillance for B.1.1.7 and other variants urgently needs to be enhanced to better inform the public health response.

Longitudinal data in peripheral blood confirm that PM20D1 is a quantitative trait locus (QTL) for Alzheimer's disease and implicate its dynamic role in disease progression

BACKGROUND

While Alzheimer's disease (AD) remains one of the most challenging diseases to tackle, genome-wide genetic/epigenetic studies reveal many disease-associated risk loci, which sheds new light onto disease heritability, provides novel insights to understand its underlying mechanism and potentially offers easily measurable biomarkers for early diagnosis and intervention.

METHODS

We analyzed whole-genome DNA methylation data collected from peripheral blood in a cohort (n = 649) from the Alzheimer's Disease Neuroimaging Initiative (ADNI) and compared the DNA methylation level at baseline among participants diagnosed with AD (n = 87), mild cognitive impairment (MCI, n = 175) and normal controls (n = 162), to identify differentially methylated regions (DMRs). We also leveraged up to 4 years of longitudinal DNA methylation data, sampled at approximately 1 year intervals to model alterations in methylation levels at DMRs to delineate methylation changes associated with aging and disease progression, by linear mixed-effects (LME) modeling for the unchanged diagnosis groups (AD, MCI and control, respectively) and U-shape testing for those with changed diagnosis (converters).

RESULTS

When compared with controls, patients with MCI consistently displayed promoter hypomethylation at methylation QTL (mQTL) gene locus PM20D1. This promoter hypomethylation was even more prominent in patients with mild to moderate AD. This is in stark contrast with previously reported hypermethylation in hippocampal and frontal cortex brain tissues in patients with advanced-stage AD at this locus. From longitudinal data, we show that initial promoter hypomethylation of PM20D1 during MCI and early stage AD is reversed to eventual promoter hypermethylation in late stage AD, which helps to complete a fuller picture of methylation dynamics. We also confirm this observation in an independent cohort from the Religious Orders Study and Memory and Aging Project (ROSMAP) Study using DNA methylation and gene expression data from brain tissues as neuropathological staging (Braak score) advances.

CONCLUSIONS

Our results confirm that PM20D1 is an mQTL in AD and demonstrate that it plays a dynamic role at different stages of the disease. Further in-depth study is thus warranted to fully decipher its role in the evolution of AD and potentially explore its utility as a blood-based biomarker for AD.

Integrative genomic meta-analysis reveals novel molecular insights into cystic fibrosis and ΔF508-CFTR rescue

Cystic fibrosis (CF), caused by mutations to CFTR, leads to severe and progressive lung disease. The most common mutant, ΔF508-CFTR, undergoes proteasomal degradation, extinguishing its anion channel function. Numerous in vitro interventions have been identified to partially rescue ΔF508-CFTR function yet remain poorly understood. Improved understanding of both the altered state of CF cells and the mechanisms of existing rescue strategies could reveal novel therapeutic strategies. Toward this aim, we measured transcriptional profiles of established temperature, genetic, and chemical interventions that rescue ΔF508-CFTR and also re-analyzed public datasets characterizing transcription in human CF vs. non-CF samples from airway and whole blood. Meta-analysis yielded a core disease signature and two core rescue signatures. To interpret these through the lens of prior knowledge, we compiled a "CFTR Gene Set Library" from literature. The core disease signature revealed remarkably strong connections to genes with established effects on CFTR trafficking and function and suggested novel roles of EGR1 and SGK1 in the disease state. Our data also revealed an unexpected mechanistic link between several genetic rescue interventions and the unfolded protein response. Finally, we found that C18, an analog of the CFTR corrector compound Lumacaftor, induces almost no transcriptional perturbation despite its rescue activity.

MicroRNA-195 controls MICU1 expression and tumor growth in ovarian cancer

MICU1 is a mitochondrial inner membrane protein that inhibits mitochondrial calcium entry; elevated MICU1 expression is characteristic of many cancers, including ovarian cancer. MICU1 induces both glycolysis and chemoresistance and is associated with poor clinical outcomes. However, there are currently no available interventions to normalize aberrant MICU1 expression. Here, we demonstrate that microRNA-195-5p (miR-195) directly targets the 3' UTR of the MICU1 mRNA and represses MICU1 expression. Additionally, miR-195 is under-expressed in ovarian cancer cell lines, and restoring miR-195 expression reestablishes native MICU1 levels and the associated phenotypes. Stable expression of miR-195 in a human xenograft model of ovarian cancer significantly reduces tumor growth, increases tumor doubling times, and enhances overall survival. In conclusion, miR-195 controls MICU1 levels in ovarian cancer and could be exploited to normalize aberrant MICU1 expression, thus reversing both glycolysis and chemoresistance and consequently improving patient outcomes.

miR155 regulation of behavior, neuropathology, and cortical transcriptomics in Alzheimer's disease

MicroRNAs are recognized as important regulators of many facets of physiological brain function while also being implicated in the pathogenesis of several neurological disorders. Dysregulation of miR155 is widely reported across a variety of neurodegenerative conditions, including Alzheimer's disease (AD), Parkinson's disease, amyotrophic lateral sclerosis, and traumatic brain injury. In previous work, we observed that experimentally validated miR155 gene targets were consistently enriched among genes identified as differentially expressed across multiple brain tissue and disease contexts. In particular, we found that human herpesvirus-6A (HHV-6A) suppressed miR155, recapitulating reports of miR155 inhibition by HHV-6A in infected T-cells, thyrocytes, and natural killer cells. In earlier studies, we also reported the effects of constitutive deletion of miR155 on accelerating the accumulation of Aβ deposits in 4-month-old APP/PSEN1 mice. Herein, we complete the cumulative characterization of transcriptomic, electrophysiological, neuropathological, and learning behavior profiles from 4-, 8- and 10-month-old WT and APP/PSEN1 mice in the absence or presence of miR155. We also integrated human post-mortem brain RNA-sequences from four independent AD consortium studies, together comprising 928 samples collected from six brain regions. We report that gene expression perturbations associated with miR155 deletion in mouse cortex are in aggregate observed to be concordant with AD-associated changes across these independent human late-onset AD (LOAD) data sets, supporting the relevance of our findings to human disease. LOAD has recently been formulated as the clinicopathological manifestation of a multiplex of genetic underpinnings and pathophysiological mechanisms. Our accumulated data are consistent with such a formulation, indicating that miR155 may be uniquely positioned at the intersection of at least four components of this LOAD "multiplex": (1) innate immune response pathways; (2) viral response gene networks; (3) synaptic pathology; and (4) proamyloidogenic pathways involving the amyloid β peptide (Aβ).

A Novel Approach to Safer Glucocorticoid Receptor-Targeted Anti-lymphoma Therapy via REDD1 (Regulated in Development and DNA Damage 1) Inhibition

Glucocorticoids are widely used for therapy of hematologic malignancies. Unfortunately, chronic treatment with glucocorticoids commonly leads to adverse effects including skin and muscle atrophy and osteoporosis. We found recently that REDD1 (regulated in development and DNA damage 1) plays central role in steroid atrophy. Here, we tested whether REDD1 suppression makes glucocorticoid-based therapy of blood cancer safer. Unexpectedly, approximately 50% of top putative REDD1 inhibitors selected by bioinformatics screening of Library of Integrated Network-Based Cellular Signatures database (LINCS) were PI3K/Akt/mTOR inhibitors. We selected Wortmannin, LY294002, and AZD8055 for our studies and showed that they blocked basal and glucocorticoid-induced REDD1 expression. Moreover, all PI3K/mTOR/Akt inhibitors modified glucocorticoid receptor function shifting it toward therapeutically important transrepression. PI3K/Akt/mTOR inhibitors enhanced anti-lymphoma effects of Dexamethasone in vitro and in vivo, in lymphoma xenograft model. The therapeutic effects of PI3K inhibitor+Dexamethasone combinations ranged from cooperative to synergistic, especially in case of LY294002 and Rapamycin, used as a previously characterized reference REDD1 inhibitor. We found that coadministration of LY294002 or Rapamycin with Dexamethasone protected skin against Dexamethasone-induced atrophy, and normalized RANKL/OPG ratio indicating a reduction of Dexamethasone-induced osteoporosis. Together, our results provide foundation for further development of safer and more effective glucocorticoid-based combination therapy of hematologic malignancies using PI3K/Akt/mTOR inhibitors.

Deep representation learning of electronic health records to unlock patient stratification at scale

Deriving disease subtypes from electronic health records (EHRs) can guide next-generation personalized medicine. However, challenges in summarizing and representing patient data prevent widespread practice of scalable EHR-based stratification analysis. Here we present an unsupervised framework based on deep learning to process heterogeneous EHRs and derive patient representations that can efficiently and effectively enable patient stratification at scale. We considered EHRs of 1,608,741 patients from a diverse hospital cohort comprising a total of 57,464 clinical concepts. We introduce a representation learning model based on word embeddings, convolutional neural networks, and autoencoders (i.e., ConvAE) to transform patient trajectories into low-dimensional latent vectors. We evaluated these representations as broadly enabling patient stratification by applying hierarchical clustering to different multi-disease and disease-specific patient cohorts. ConvAE significantly outperformed several baselines in a clustering task to identify patients with different complex conditions, with 2.61 entropy and 0.31 purity average scores. When applied to stratify patients within a certain condition, ConvAE led to various clinically relevant subtypes for different disorders, including type 2 diabetes, Parkinson's disease, and Alzheimer's disease, largely related to comorbidities, disease progression, and symptom severity. With these results, we demonstrate that ConvAE can generate patient representations that lead to clinically meaningful insights. This scalable framework can help better understand varying etiologies in heterogeneous sub-populations and unlock patterns for EHR-based research in the realm of personalized medicine.

Sepsis in the era of data-driven medicine: personalizing risks, diagnoses, treatments and prognoses

Sepsis is a series of clinical syndromes caused by the immunological response to infection. The clinical evidence for sepsis could typically attribute to bacterial infection or bacterial endotoxins, but infections due to viruses, fungi or parasites could also lead to sepsis. Regardless of the etiology, rapid clinical deterioration, prolonged stay in intensive care units and high risk for mortality correlate with the incidence of sepsis. Despite its prevalence and morbidity, improvement in sepsis outcomes has remained limited. In this comprehensive review, we summarize the current landscape of risk estimation, diagnosis, treatment and prognosis strategies in the setting of sepsis and discuss future challenges. We argue that the advent of modern technologies such as in-depth molecular profiling, biomedical big data and machine intelligence methods will augment the treatment and prevention of sepsis. The volume, variety, veracity and velocity of heterogeneous data generated as part of healthcare delivery and recent advances in biotechnology-driven therapeutics and companion diagnostics may provide a new wave of approaches to identify the most at-risk sepsis patients and reduce the symptom burden in patients within shorter turnaround times. Developing novel therapies by leveraging modern drug discovery strategies including computational drug repositioning, cell and gene-therapy, clustered regularly interspaced short palindromic repeats -based genetic editing systems, immunotherapy, microbiome restoration, nanomaterial-based therapy and phage therapy may help to develop treatments to target sepsis. We also provide empirical evidence for potential new sepsis targets including FER and STARD3NL. Implementing data-driven methods that use real-time collection and analysis of clinical variables to trace, track and treat sepsis-related adverse outcomes will be key. Understanding the root and route of sepsis and its comorbid conditions that complicate treatment outcomes and lead to organ dysfunction may help to facilitate identification of most at-risk patients and prevent further deterioration. To conclude, leveraging the advances in precision medicine, biomedical data science and translational bioinformatics approaches may help to develop better strategies to diagnose and treat sepsis in the next decade.

CANDI: an R package and Shiny app for annotating radiographs and evaluating computer-aided diagnosis

Motivation

Radiologists have used algorithms for Computer-Aided Diagnosis (CAD) for decades. These algorithms use machine learning with engineered features, and there have been mixed findings on whether they improve radiologists’ interpretations. Deep learning offers superior performance but requires more training data and has not been evaluated in joint algorithm-radiologist decision systems.

Results

We developed the Computer-Aided Note and Diagnosis Interface (CANDI) for collaboratively annotating radiographs and evaluating how algorithms alter human interpretation. The annotation app collects classification, segmentation, and image captioning training data, and the evaluation app randomizes the availability of CAD tools to facilitate clinical trials on radiologist enhancement.

Availability and implementation

Demonstrations and source code are hosted at (https://candi.nextgenhealthcare.org), and (https://github.com/mbadge/candi), respectively, under GPL-3 license.

Supplementary information

Supplementary material is available at Bioinformatics online.

Identifying Acute Low Back Pain Episodes in Primary Care Practice From Clinical Notes: Observational Study

BACKGROUND

Acute and chronic low back pain (LBP) are different conditions with different treatments. However, they are coded in electronic health records with the same International Classification of Diseases, 10th revision (ICD-10) code (M54.5) and can be differentiated only by retrospective chart reviews. This prevents an efficient definition of data-driven guidelines for billing and therapy recommendations, such as return-to-work options.

OBJECTIVE

The objective of this study was to evaluate the feasibility of automatically distinguishing acute LBP episodes by analyzing free-text clinical notes.

METHODS

We used a dataset of 17,409 clinical notes from different primary care practices; of these, 891 documents were manually annotated as acute LBP and 2973 were generally associated with LBP via the recorded ICD-10 code. We compared different supervised and unsupervised strategies for automated identification: keyword search, topic modeling, logistic regression with bag of n-grams and manual features, and deep learning (a convolutional neural network-based architecture [ConvNet]). We trained the supervised models using either manual annotations or ICD-10 codes as positive labels.

RESULTS

ConvNet trained using manual annotations obtained the best results with an area under the receiver operating characteristic curve of 0.98 and an F score of 0.70. ConvNet's results were also robust to reduction of the number of manually annotated documents. In the absence of manual annotations, topic models performed better than methods trained using ICD-10 codes, which were unsatisfactory for identifying LBP acuity.

CONCLUSIONS

This study uses clinical notes to delineate a potential path toward systematic learning of therapeutic strategies, billing guidelines, and management options for acute LBP at the point of care.

Alterations in innate immunity and epithelial cell differentiation are the molecular pillars of hidradenitis suppurativa

BACKGROUND

The large unmet need of hidradenitis suppurativa/acne inversa (HS) therapy requires the elucidation of disease-driving mechanisms and tissue targeting.

OBJECTIVE

Robust characterization of the underlying HS mechanisms and detection of the involved skin compartments.

METHODS

Hidradenitis suppurativa/acne inversa molecular taxonomy and key signalling pathways were studied by whole transcriptome profiling. Dysregulated genes were detected by comparing lesional and non-lesional skin obtained from female HS patients and matched healthy controls using the Agilent array platform. The differential gene expression was confirmed by quantitative real-time PCR and targeted protein characterization via immunohistochemistry in another set of female patients. HS-involved skin compartments were also recognized by immunohistochemistry.

RESULTS

Alterations to key regulatory pathways involving glucocorticoid receptor, atherosclerosis, HIF1α and IL17A signalling as well as inhibition of matrix metalloproteases were detected. From a functional standpoint, cellular assembly, maintenance and movement, haematological system development and function, immune cell trafficking and antimicrobial response were key processes probably being affected in HS. Sixteen genes were found to characterize HS from a molecular standpoint (DEFB4, MMP1, GJB2, PI3, KRT16, MMP9, SERPINB4, SERPINB3, SPRR3, S100A8, S100A9, S100A12, S100A7A (15), KRT6A, TCN1, TMPRSS11D). Among the proteins strongly expressed in HS, calgranulin-A, calgranulin-B and serpin-B4 were detected in the hair root sheath, koebnerisin and connexin-32 in stratum granulosum, transcobalamin-1 in stratum spinosum/hair root sheath, small prolin-rich protein-3 in apocrine sweat gland ducts/sebaceous glands-ducts and matrix metallopeptidase-9 in resident monocytes.

CONCLUSION

Our findings highlight a panel of immune-related drivers in HS, which influence innate immunity and cell differentiation in follicular and epidermal keratinocytes as well as skin glands.

MEWS++: Enhancing the Prediction of Clinical Deterioration in Admitted Patients through a Machine Learning Model

Early detection of patients at risk for clinical deterioration is crucial for timely intervention. Traditional detection systems rely on a limited set of variables and are unable to predict the time of decline. We describe a machine learning model called MEWS++ that enables the identification of patients at risk of escalation of care or death six hours prior to the event. A retrospective single-center cohort study was conducted from July 2011 to July 2017 of adult (age > 18) inpatients excluding psychiatric, parturient, and hospice patients. Three machine learning models were trained and tested: random forest (RF), linear support vector machine, and logistic regression. We compared the models’ performance to the traditional Modified Early Warning Score (MEWS) using sensitivity, specificity, and Area Under the Curve for Receiver Operating Characteristic (AUC-ROC) and Precision-Recall curves (AUC-PR). The primary outcome was escalation of care from a floor bed to an intensive care or step-down unit, or death, within 6 h. A total of 96,645 patients with 157,984 hospital encounters and 244,343 bed movements were included. Overall rate of escalation or death was 3.4%. The RF model had the best performance with sensitivity 81.6%, specificity 75.5%, AUC-ROC of 0.85, and AUC-PR of 0.37. Compared to traditional MEWS, sensitivity increased 37%, specificity increased 11%, and AUC-ROC increased 14%. This study found that using machine learning and readily available clinical data, clinical deterioration or death can be predicted 6 h prior to the event. The model we developed can warn of patient deterioration hours before the event, thus helping make timely clinical decisions.

Mutation-derived Neoantigen-specific T-cell Responses in Multiple Myeloma

PURPOSE

Somatic mutations in cancer cells can give rise to novel protein sequences that can be presented by antigen-presenting cells as neoantigens to the host immune system. Tumor neoantigens represent excellent targets for immunotherapy, due to their specific expression in cancer tissue. Despite the widespread use of immunomodulatory drugs and immunotherapies that recharge T and NK cells, there has been no direct evidence that neoantigen-specific T-cell responses are elicited in multiple myeloma.

EXPERIMENTAL DESIGN

Using next-generation sequencing data we describe the landscape of neo-antigens in 184 patients with multiple myeloma and successfully validate neoantigen-specific T cells in patients with multiple myeloma and support the feasibility of neoantigen-based therapeutic vaccines for use in cancers with intermediate mutational loads such as multiple myeloma.

RESULTS

In this study, we demonstrate an increase in neoantigen load in relapsed patients with multiple myeloma as compared with newly diagnosed patients with multiple myeloma. Moreover, we identify shared neoantigens across multiple patients in three multiple myeloma oncogenic driver genes (KRAS, NRAS, and IRF4). Next, we validate neoantigen T-cell response and clonal expansion in correlation with clinical response in relapsed patients with multiple myeloma. This is the first study to experimentally validate the immunogenicity of predicted neoantigens from next-generation sequencing in relapsed patients with multiple myeloma.

CONCLUSIONS

Our findings demonstrate that somatic mutations in multiple myeloma can be immunogenic and induce neoantigen-specific T-cell activation that is associated with antitumor activity in vitro and clinical response in vivo. Our results provide the foundation for using neoantigen targeting strategies such as peptide vaccines in future trials for patients with multiple myeloma.

Scaling structural learning with NO-BEARS to infer causal transcriptome networks

Constructing gene regulatory networks is a critical step in revealing disease mechanisms from transcriptomic data. In this work, we present NO-BEARS, a novel algorithm for estimating gene regulatory networks. The NO-BEARS algorithm is built on the basis of the NO-TEARS algorithm with two improvements. First, we propose a new constraint and its fast approximation to reduce the computational cost of the NO-TEARS algorithm. Next, we introduce a polynomial regression loss to handle non-linearity in gene expressions. Our implementation utilizes modern GPU computation that can decrease the time of hours-long CPU computation to seconds. Using synthetic data, we demonstrate improved performance, both in processing time and accuracy, on inferring gene regulatory networks from gene expression data.

Automated phenotyping of patients with non-alcoholic fatty liver disease reveals clinically relevant disease subtypes

Non-alcoholic fatty liver disease (NAFLD) is a complex heterogeneous disease which affects more than 20% of the population worldwide. Some subtypes of NAFLD have been clinically identified using hypothesis-driven methods. In this study, we used data mining techniques to search for subtypes in an unbiased fashion. Using electronic signatures of the disease, we identified a cohort of 13,290 patients with NAFLD from a hospital database. We gathered clinical data from multiple sources and applied unsupervised clustering to identify five subtypes among this cohort. Descriptive statistics and survival analysis showed that the subtypes were clinically distinct and were associated with different rates of death, cirrhosis, hepatocellular carcinoma, chronic kidney disease, cardiovascular disease, and myocardial infarction. Novel disease subtypes identified in this manner could be used to risk-stratify patients and guide management.

KRCC1: A potential therapeutic target in ovarian cancer

Using a systems biology approach to prioritize potential points of intervention in ovarian cancer, we identified the lysine rich coiled-coil 1 (KRCC1), as a potential target. High-grade serous ovarian cancer patient tumors and cells express significantly higher levels of KRCC1 which correlates with poor overall survival and chemoresistance. We demonstrate that KRCC1 is predominantly present in the chromatin-bound nuclear fraction, interacts with HDAC1, HDAC2, and with the serine-threonine phosphatase PP1CC. Silencing KRCC1 inhibits cellular plasticity, invasive properties, and potentiates apoptosis resulting in reduced tumor growth. These phenotypes are associated with increased acetylation of histones and with increased phosphorylation of H2AX and CHK1, suggesting the modulation of transcription and DNA damage that may be mediated by the action of HDAC and PP1CC, respectively. Hence, we address an urgent need to develop new targets in cancer.

MetaCyto: A Tool for Automated Meta-analysis of Mass and Flow Cytometry Data

While meta-analysis has demonstrated increased statistical power and more robust estimations in studies, the application of this commonly accepted methodology to cytometry data has been challenging. Different cytometry studies often involve diverse sets of markers. Moreover, the detected values of the same marker are inconsistent between studies due to different experimental designs and cytometer configurations. As a result, the cell subsets identified by existing auto-gating methods cannot be directly compared across studies. We developed MetaCyto for automated meta-analysis of both flow and mass cytometry (CyTOF) data. By combining clustering methods with a silhouette scanning method, MetaCyto is able to identify commonly labeled cell subsets across studies, thus enabling meta-analysis. Applying MetaCyto across a set of ten heterogeneous cytometry studies totaling 2,926 samples enabled us to identify multiple cell populations exhibiting differences in abundance between demographic groups. Software is released to the public through Bioconductor (http://bioconductor.org/packages/release/bioc/html/MetaCyto.html).

PatientExploreR: an extensible application for dynamic visualization of patient clinical history from electronic health records in the OMOP common data model

MOTIVATION

Electronic health records (EHRs) are quickly becoming omnipresent in healthcare, but interoperability issues and technical demands limit their use for biomedical and clinical research. Interactive and flexible software that interfaces directly with EHR data structured around a common data model (CDM) could accelerate more EHR-based research by making the data more accessible to researchers who lack computational expertise and/or domain knowledge.

RESULTS

We present PatientExploreR, an extensible application built on the R/Shiny framework that interfaces with a relational database of EHR data in the Observational Medical Outcomes Partnership CDM format. PatientExploreR produces patient-level interactive and dynamic reports and facilitates visualization of clinical data without any programming required. It allows researchers to easily construct and export patient cohorts from the EHR for analysis with other software. This application could enable easier exploration of patient-level data for physicians and researchers. PatientExploreR can incorporate EHR data from any institution that employs the CDM for users with approved access. The software code is free and open source under the MIT license, enabling institutions to install and users to expand and modify the application for their own purposes.

AVAILABILITY AND IMPLEMENTATION

PatientExploreR can be freely obtained from GitHub: https://github.com/BenGlicksberg/PatientExploreR. We provide instructions for how researchers with approved access to their institutional EHR can use this package. We also release an open sandbox server of synthesized patient data for users without EHR access to explore: http://patientexplorer.ucsf.edu.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Building a Secure Biomedical Data Sharing Decentralized App (DApp): Tutorial

Decentralized apps (DApps) are computer programs that run on a distributed computing system, such as a blockchain network. Unlike the client-server architecture that powers most internet apps, DApps that are integrated with a blockchain network can execute app logic that is guaranteed to be transparent, verifiable, and immutable. This new paradigm has a number of unique properties that are attractive to the biomedical and health care communities. However, instructional resources are scarcely available for biomedical software developers to begin building DApps on a blockchain. Such apps require new ways of thinking about how to build, maintain, and deploy software. This tutorial serves as a complete working prototype of a DApp, motivated by a real use case in biomedical research requiring data privacy. We describe the architecture of a DApp, the implementation details of a smart contract, a sample iPhone operating system (iOS) DApp that interacts with the smart contract, and the development tools and libraries necessary to get started. The code necessary to recreate the app is publicly available.

Privacy-Preserving Methods for Feature Engineering Using Blockchain: Review, Evaluation, and Proof of Concept

Background

The protection of private data is a key responsibility for research studies that collect identifiable information from study participants. Limiting the scope of data collection and preventing secondary use of the data are effective strategies for managing these risks. An ideal framework for data collection would incorporate feature engineering, a process where secondary features are derived from sensitive raw data in a secure environment without a trusted third party.

Objective

This study aimed to compare current approaches based on how they maintain data privacy and the practicality of their implementations. These approaches include traditional approaches that rely on trusted third parties, and cryptographic, secure hardware, and blockchain-based techniques.

Methods

A set of properties were defined for evaluating each approach. A qualitative comparison was presented based on these properties. The evaluation of each approach was framed with a use case of sharing geolocation data for biomedical research.

Results

We found that approaches that rely on a trusted third party for preserving participant privacy do not provide sufficiently strong guarantees that sensitive data will not be exposed in modern data ecosystems. Cryptographic techniques incorporate strong privacy-preserving paradigms but are appropriate only for select use cases or are currently limited because of computational complexity. Blockchain smart contracts alone are insufficient to provide data privacy because transactional data are public. Trusted execution environments (TEEs) may have hardware vulnerabilities and lack visibility into how data are processed. Hybrid approaches combining blockchain and cryptographic techniques or blockchain and TEEs provide promising frameworks for privacy preservation. For reference, we provide a software implementation where users can privately share features of their geolocation data using the hybrid approach combining blockchain with TEEs as a supplement.

Conclusions

Blockchain technology and smart contracts enable the development of new privacy-preserving feature engineering methods by obviating dependence on trusted parties and providing immutable, auditable data processing workflows. The overlap between blockchain and cryptographic techniques or blockchain and secure hardware technologies are promising fields for addressing important data privacy needs. Hybrid blockchain and TEE frameworks currently provide practical tools for implementing experimental privacy-preserving applications.

Artificial Intelligence in Cardiology

Artificial intelligence and machine learning are poised to influence nearly every aspect of the human condition, and cardiology is not an exception to this trend. This paper provides a guide for clinicians on relevant aspects of artificial intelligence and machine learning, reviews selected applications of these methods in cardiology to date, and identifies how cardiovascular medicine could incorporate artificial intelligence in the future. In particular, the paper first reviews predictive modeling concepts relevant to cardiology such as feature selection and frequent pitfalls such as improper dichotomization. Second, it discusses common algorithms used in supervised learning and reviews selected applications in cardiology and related disciplines. Third, it describes the advent of deep learning and related methods collectively called unsupervised learning, provides contextual examples both in general medicine and in cardiovascular medicine, and then explains how these methods could be applied to enable precision cardiology and improve patient outcomes.

Integrative analysis of loss-of-function variants in clinical and genomic data reveals novel genes associated with cardiovascular traits

Background

Genetic loss-of-function variants (LoFs) associated with disease traits are increasingly recognized as critical evidence for the selection of therapeutic targets. We integrated the analysis of genetic and clinical data from 10,511 individuals in the Mount Sinai BioMe Biobank to identify genes with loss-of-function variants (LoFs) significantly associated with cardiovascular disease (CVD) traits, and used RNA-sequence data of seven metabolic and vascular tissues isolated from 600 CVD patients in the Stockholm-Tartu Atherosclerosis Reverse Network Engineering Task (STARNET) study for validation. We also carried out in vitro functional studies of several candidate genes, and in vivo studies of one gene.

Results

We identified LoFs in 433 genes significantly associated with at least one of 10 major CVD traits. Next, we used RNA-sequence data from the STARNET study to validate 115 of the 433 LoF harboring-genes in that their expression levels were concordantly associated with corresponding CVD traits. Together with the documented hepatic lipid-lowering gene, APOC3, the expression levels of six additional liver LoF-genes were positively associated with levels of plasma lipids in STARNET. Candidate LoF-genes were subjected to gene silencing in HepG2 cells with marked overall effects on cellular LDLR, levels of triglycerides and on secreted APOB100 and PCSK9. In addition, we identified novel LoFs in DGAT2 associated with lower plasma cholesterol and glucose levels in BioMe that were also confirmed in STARNET, and showed a selective DGAT2-inhibitor in C57BL/6 mice not only significantly lowered fasting glucose levels but also affected body weight.

Conclusion

In sum, by integrating genetic and electronic medical record data, and leveraging one of the world’s largest human RNA-sequence datasets (STARNET), we identified known and novel CVD-trait related genes that may serve as targets for CVD therapeutics and as such merit further investigation.

Electronic supplementary material

The online version of this article (10.1186/s12920-019-0542-3) contains supplementary material, which is available to authorized users.

A transcriptomic model to predict increase in fibrous cap thickness in response to high-dose statin treatment: Validation by serial intracoronary OCT imaging

Background

Fibrous cap thickness (FCT), best measured by intravascular optical coherence tomography (OCT), is the most important determinant of plaque rupture in the coronary arteries. Statin treatment increases FCT and thus reduces the likelihood of acute coronary events. However, substantial statin-related FCT increase occurs in only a subset of patients. Currently, there are no methods to predict which patients will benefit. We use transcriptomic data from a clinical trial of rosuvastatin to predict if a patient's FCT will increase in response to statin therapy.

Methods

FCT was measured using OCT in 69 patients at (1) baseline and (2) after 8–10 weeks of 40  mg rosuvastatin. Peripheral blood mononuclear cells were assayed via microarray. We constructed machine learning models with baseline gene expression data to predict change in FCT. Finally, we ascertained the biological functions of the most predictive transcriptomic markers.

Findings

Machine learning models were able to predict FCT responders using baseline gene expression with high fidelity (Classification AUC = 0.969 and 0.972). The first model (elastic net) using 73 genes had an accuracy of 92.8%, sensitivity of 94.1%, and specificity of 91.4%. The second model (KTSP) using 18 genes has an accuracy of 95.7%, sensitivity of 94.3%, and specificity of 97.1%. We found 58 enriched gene ontology terms, including many involved with immune cell function and cholesterol biometabolism.

Interpretation

In this pilot study, transcriptomic models could predict if FCT increased following 8–10 weeks of rosuvastatin. These findings may have significance for therapy selection and could supplement invasive imaging modalities.

The next generation of precision medicine: observational studies, electronic health records, biobanks and continuous monitoring

Precision medicine can utilize new techniques in order to more effectively translate research findings into clinical practice. In this article, we first explore the limitations of traditional study designs, which stem from (to name a few): massive cost for the assembly of large patient cohorts; non-representative patient data; and the astounding complexity of human biology. Second, we propose that harnessing electronic health records and mobile device biometrics coupled to longitudinal data may prove to be a solution to many of these problems by capturing a 'real world' phenotype. We envision that future biomedical research utilizing more precise approaches to patient care will utilize continuous and longitudinal data sources.

Systems Pharmacology Identifies an Arterial Wall Regulatory Gene Network Mediating Coronary Artery Disease Side Effects of Antiretroviral Therapy

BACKGROUND

Antiretroviral therapy (ART) for HIV infection increases risk for coronary artery disease (CAD), presumably by causing dyslipidemia and increased atherosclerosis. We applied systems pharmacology to identify and validate specific regulatory gene networks through which ART drugs may promote CAD.

METHODS

Transcriptional responses of human cell lines to 15 ART drugs retrieved from the Library of Integrated Cellular Signatures (overall 1127 experiments) were used to establish consensus ART gene/transcriptional signatures. Next, enrichments of differentially expressed genes and gene-gene connectivity within these ART-consensus signatures were sought in 30 regulatory gene networks associated with CAD and CAD-related phenotypes in the Stockholm Atherosclerosis Gene Expression study.

RESULTS

Ten of 15 ART signatures were significantly enriched both for differential expression and connectivity in a specific atherosclerotic arterial wall regulatory gene network (AR-RGN) causal for CAD involving RNA processing genes. An atherosclerosis in vitro model of cholestryl ester-loaded foam cells was then used for experimental validation. Treatments of these foam cells with ritonavir, nelfinavir, and saquinavir at least doubled cholestryl ester accumulation ( P=0.02, 0.0009, and 0.02, respectively), whereas RNA silencing of the AR-RGN top key driver, PQBP1 (polyglutamine binding protein 1), significantly curbed cholestryl ester accumulation following treatment with any of these ART drugs by >37% ( P<0.05).

CONCLUSIONS

By applying a novel systems pharmacology data analysis framework, 3 commonly used ARTs (ritonavir, nelfinavir, and saquinavir) were found altering the activity of AR-RGN, a regulatory gene network promoting foam cell formation and risk of CAD. Targeting AR-RGN or its top key driver PQBP1 may help reduce CAD side effects of these ART drugs.

Critical period plasticity-related transcriptional aberrations in schizophrenia and bipolar disorder

Childhood critical periods of experience-dependent plasticity are essential for the development of environmentally appropriate behavior and cognition. Disruption of critical periods can alter development of normal function and confer risk for neurodevelopmental disorders. While genes and their expression relevant to neurodevelopment are associated with schizophrenia, the molecular relationship between schizophrenia and critical periods has not been assessed systematically. Here, we apply a transcriptome-based bioinformatics approach to assess whether genes associated with the human critical period for visual cortex plasticity, a well-studied model of cortical critical periods, are aberrantly expressed in schizophrenia and bipolar disorder. Across two dozen datasets encompassing 522 cases and 374 controls, we find that the majority show aberrations in expression of genes associated with the critical period. We observed both hyper- and hypo-critical period plasticity phenotypes at the transcriptome level, which partially mapped to drug candidates that reverse the disorder signatures in silico. Our findings indicate plasticity aberrations in schizophrenia and their treatment may need to be considered in the context of subpopulations with elevated and others reduced plasticity. Future work should leverage ongoing consortia RNA-sequencing efforts to tease out the sources of plasticity-related transcriptional aberrations seen in schizophrenia, including true biological heterogeneity, interaction between normal development/aging and the disorder, and medication history. Our study also urges innovation towards direct assessment of visual cortex plasticity in humans with schizophrenia to precisely deconstruct the role of plasticity in this disorder.

Deep learning predicts hip fracture using confounding patient and healthcare variables

Hip fractures are a leading cause of death and disability among older adults. Hip fractures are also the most commonly missed diagnosis on pelvic radiographs, and delayed diagnosis leads to higher cost and worse outcomes. Computer-aided diagnosis (CAD) algorithms have shown promise for helping radiologists detect fractures, but the image features underpinning their predictions are notoriously difficult to understand. In this study, we trained deep-learning models on 17,587 radiographs to classify fracture, 5 patient traits, and 14 hospital process variables. All 20 variables could be individually predicted from a radiograph, with the best performances on scanner model (AUC = 1.00), scanner brand (AUC = 0.98), and whether the order was marked "priority" (AUC = 0.79). Fracture was predicted moderately well from the image (AUC = 0.78) and better when combining image features with patient data (AUC = 0.86, DeLong paired AUC comparison, p = 2e-9) or patient data plus hospital process features (AUC = 0.91, p = 1e-21). Fracture prediction on a test set that balanced fracture risk across patient variables was significantly lower than a random test set (AUC = 0.67, DeLong unpaired AUC comparison, p = 0.003); and on a test set with fracture risk balanced across patient and hospital process variables, the model performed randomly (AUC = 0.52, 95% CI 0.46-0.58), indicating that these variables were the main source of the model's fracture predictions. A single model that directly combines image features, patient, and hospital process data outperforms a Naive Bayes ensemble of an image-only model prediction, patient, and hospital process data. If CAD algorithms are inexplicably leveraging patient and process variables in their predictions, it is unclear how radiologists should interpret their predictions in the context of other known patient data. Further research is needed to illuminate deep-learning decision processes so that computers and clinicians can effectively cooperate.

Natural Language Processing of Clinical Notes on Chronic Diseases: Systematic Review

Background

Novel approaches that complement and go beyond evidence-based medicine are required in the domain of chronic diseases, given the growing incidence of such conditions on the worldwide population. A promising avenue is the secondary use of electronic health records (EHRs), where patient data are analyzed to conduct clinical and translational research. Methods based on machine learning to process EHRs are resulting in improved understanding of patient clinical trajectories and chronic disease risk prediction, creating a unique opportunity to derive previously unknown clinical insights. However, a wealth of clinical histories remains locked behind clinical narratives in free-form text. Consequently, unlocking the full potential of EHR data is contingent on the development of natural language processing (NLP) methods to automatically transform clinical text into structured clinical data that can guide clinical decisions and potentially delay or prevent disease onset.

Objective

The goal of the research was to provide a comprehensive overview of the development and uptake of NLP methods applied to free-text clinical notes related to chronic diseases, including the investigation of challenges faced by NLP methodologies in understanding clinical narratives.

Methods

Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) guidelines were followed and searches were conducted in 5 databases using “clinical notes,” “natural language processing,” and “chronic disease” and their variations as keywords to maximize coverage of the articles.

Results

Of the 2652 articles considered, 106 met the inclusion criteria. Review of the included papers resulted in identification of 43 chronic diseases, which were then further classified into 10 disease categories using the International Classification of Diseases, 10th Revision. The majority of studies focused on diseases of the circulatory system (n=38) while endocrine and metabolic diseases were fewest (n=14). This was due to the structure of clinical records related to metabolic diseases, which typically contain much more structured data, compared with medical records for diseases of the circulatory system, which focus more on unstructured data and consequently have seen a stronger focus of NLP. The review has shown that there is a significant increase in the use of machine learning methods compared to rule-based approaches; however, deep learning methods remain emergent (n=3). Consequently, the majority of works focus on classification of disease phenotype with only a handful of papers addressing extraction of comorbidities from the free text or integration of clinical notes with structured data. There is a notable use of relatively simple methods, such as shallow classifiers (or combination with rule-based methods), due to the interpretability of predictions, which still represents a significant issue for more complex methods. Finally, scarcity of publicly available data may also have contributed to insufficient development of more advanced methods, such as extraction of word embeddings from clinical notes.

Conclusions

Efforts are still required to improve (1) progression of clinical NLP methods from extraction toward understanding; (2) recognition of relations among entities rather than entities in isolation; (3) temporal extraction to understand past, current, and future clinical events; (4) exploitation of alternative sources of clinical knowledge; and (5) availability of large-scale, de-identified clinical corpora.

Designing Robust N-of-1 Studies for Precision Medicine: Simulation Study and Design Recommendations

Background

Recent advances in molecular biology, sensors, and digital medicine have led to an explosion of products and services for high-resolution monitoring of individual health. The N-of-1 study has emerged as an important methodological tool for harnessing these new data sources, enabling researchers to compare the effectiveness of health interventions at the level of a single individual.

Objective

N-of-1 studies are susceptible to several design flaws. We developed a model that generates realistic data for N-of-1 studies to enable researchers to optimize study designs in advance.

Methods

Our stochastic time-series model simulates an N-of-1 study, incorporating all study-relevant effects, such as carryover and wash-in effects, as well as various sources of noise. The model can be used to produce realistic simulated data for a near-infinite number of N-of-1 study designs, treatment profiles, and patient characteristics.

Results

Using simulation, we demonstrate how the number of treatment blocks, ordering of treatments within blocks, duration of each treatment, and sampling frequency affect our ability to detect true differences in treatment efficacy. We provide a set of recommendations for study designs on the basis of treatment, outcomes, and instrument parameters, and make our simulation software publicly available for use by the precision medicine community.

Conclusions

Simulation can facilitate rapid optimization of N-of-1 study designs and increase the likelihood of study success while minimizing participant burden.

Lyme Disease Patient Trajectories Learned from Electronic Medical Data for Stratification of Disease Risk and Therapeutic Response

Lyme disease (LD) is the most common tick-borne illness in the United States. Although appropriate antibiotic treatment is effective for most cases, up to 20% of patients develop post-treatment Lyme disease syndrome (PTLDS). There is an urgent need to improve clinical management of LD using precise understanding of disease and patient stratification. We applied machine-learning to electronic medical records to better characterize the heterogeneity of LD and developed predictive models for identifying medications that are associated with risks of subsequent comorbidities. For broad disease categories, we identified 3, 16, and 17 comorbidities within 2, 5, and 10 years of diagnosis, respectively. At a higher resolution of ICD-9 codes, we identified known associations with LD including chronic pain and cognitive disorders, as well as particular comorbidities on a timescale that matched PTLDS symptomology. We identified 7, 30, and 35 medications associated with risks of these comorbidities within 2, 5, and 10 years, respectively. For instance, the first-line antibiotic doxycycline exhibited a consistently protective association for typical symptoms of LD, including backache. Our approach and findings may suggest new hypotheses for more personalized treatments regimens for LD patients.

Predicting Future Cardiovascular Events in Patients With Peripheral Artery Disease Using Electronic Health Record Data

BACKGROUND

Patients with peripheral artery disease (PAD) are at risk of major adverse cardiac and cerebrovascular events. There are no readily available risk scores that can accurately identify which patients are most likely to sustain an event, making it difficult to identify those who might benefit from more aggressive intervention. Thus, we aimed to develop a novel predictive model-using machine learning methods on electronic health record data-to identify which PAD patients are most likely to develop major adverse cardiac and cerebrovascular events.

METHODS AND RESULTS

Data were derived from patients diagnosed with PAD at 2 tertiary care institutions. Predictive models were built using a common data model that allowed for utilization of both structured (coded) and unstructured (text) data. Only data from time of entry into the health system up to PAD diagnosis were used for modeling. Models were developed and tested using nested cross-validation. A total of 7686 patients were included in learning our predictive models. Utilizing almost 1000 variables, our best predictive model accurately determined which PAD patients would go on to develop major adverse cardiac and cerebrovascular events with an area under the curve of 0.81 (95% CI, 0.80-0.83).

CONCLUSIONS

Machine learning algorithms applied to data in the electronic health record can learn models that accurately identify PAD patients at risk of future major adverse cardiac and cerebrovascular events, highlighting the great potential of electronic health records to provide automated risk stratification for cardiovascular diseases. Common data models that can enable cross-institution research and technology development could potentially be an important aspect of widespread adoption of newer risk-stratification models.