Reconstruction of the mouse otocyst and early neuroblast lineage at single-cell resolution.

The otocyst harbors progenitors for most cell types of the mature inner ear. Developmental lineage analyses and gene expression studies suggest that distinct progenitor populations are compartmentalized to discrete axial domains in the early otocyst. Here, we conducted highly parallel quantitative RT-PCR measurements on 382 individual cells from the developing otocyst and neuroblast lineages to assay 96 genes representing established otic markers, signaling-pathway-associated transcripts, and novel otic-specific genes. By applying multivariate cluster, principal component, and network analyses to the data matrix, we were able to readily distinguish the delaminating neuroblasts and to describe progressive states of gene expression in this population at single-cell resolution. It further established a three-dimensional model of the otocyst in which each individual cell can be precisely mapped into spatial expression domains. Our bioinformatic modeling revealed spatial dynamics of different signaling pathways active during early neuroblast development and prosensory domain specification.

Expanded T lymphocytes in the cerebrospinal fluid of multiple sclerosis patients are specific for Epstein-Barr-virus-infected B cells.

Epstein-Barr virus (EBV) infection has long been associated with multiple sclerosis (MS), but the role of EBV in the pathogenesis of MS is not clear. Our hypothesis is that a major fraction of the expanded clones of T lymphocytes in the cerebrospinal fluid (CSF) are specific for autologous EBV-infected B cells. We obtained blood and CSF samples from eight relapsing-remitting patients in the process of diagnosis. We stimulated cells from the blood with autologous EBV-infected lymphoblastoid cell lines (LCL), EBV, varicella zoster virus, influenza, and candida and sorted the responding cells with flow cytometry after 6 d. We sequenced the RNA for T cell receptors (TCR) from CSF, unselected blood cells, and the antigen-specific cells. We used the TCR Vß CDR3 sequences from the antigen-specific cells to assign antigen specificity to the sequences from the CSF and blood. LCL-specific cells comprised 13.0 ± 4.3% (mean ± SD) of the total reads present in CSF and 13.3 ± 7.5% of the reads present in blood. The next most abundant antigen specificity was flu, which was 4.7 ± 1.7% of the reads in the CSF and 9.3 ± 6.6% in the blood. The prominence of LCL-specific reads was even more marked in the top 1% most abundant CSF clones with statistically significant 47% mean overlap with LCL. We conclude that LCL-specific sequences form a major portion of the TCR repertoire in both CSF and blood and that expanded clones specific for LCL are present in MS CSF. This has important implications for the pathogenesis of MS.

Tissue-specific atlas of trans-models for gene regulation elucidates complex regulation patterns.

BACKGROUND: Deciphering gene regulation is essential for understanding the underlying mechanisms of healthy and disease states. While the regulatory networks formed by transcription factors (TFs) and their target genes has been mostly studied with relation to cis effects such as in TF binding sites, we focused on trans effects of TFs on the expression of their transcribed genes and their potential mechanisms. RESULTS: We provide a comprehensive tissue-specific atlas, spanning 49 tissues of TF variations affecting gene expression through computational models considering two potential mechanisms, including combinatorial regulation by the expression of the TFs, and by genetic variants within the TF. We demonstrate that similarity between tissues based on our discovered genes corresponds to other types of tissue similarity. The genes affected by complex TF regulation, and their modelled TFs, were highly enriched for pharmacogenomic functions, while the TFs themselves were also enriched in several cancer and metabolic pathways. Additionally, genes that appear in multiple clusters are enriched for regulation of immune system while tissue clusters include cluster-specific genes that are enriched for biological functions and diseases previously associated with the tissues forming the cluster. Finally, our atlas exposes multilevel regulation across multiple tissues, where TFs regulate other TFs through the two tested mechanisms. CONCLUSIONS: Our tissue-specific atlas provides hierarchical tissue-specific trans genetic regulations that can be further studied for association with human phenotypes.

Montelukast as a repurposable additive drug for standard-efficacy multiple sclerosis treatment: Emulating clinical trials with retrospective administrative health claims data.

BACKGROUND: Effective and safe treatment options for multiple sclerosis (MS) are still needed. Montelukast, a leukotriene receptor antagonist (LTRA) currently indicated for asthma or allergic rhinitis, may provide an additional therapeutic approach. OBJECTIVE: The study aimed to evaluate the effects of montelukast on the relapses of people with MS (pwMS). METHODS: In this retrospective case-control study, two independent longitudinal claims datasets were used to emulate randomized clinical trials (RCTs). We identified pwMS aged 18-65 years, on MS disease-modifying therapies concomitantly, in de-identified claims from Optum's Clinformatics® Data Mart (CDM) and IQVIA PharMetrics® Plus for Academics. Cases included 483 pwMS on montelukast and with medication adherence in CDM and 208 in PharMetrics Plus for Academics. We randomly sampled controls from 35,330 pwMS without montelukast prescriptions in CDM and 10,128 in PharMetrics Plus for Academics. Relapses were measured over a 2-year period through inpatient hospitalization and corticosteroid claims. A doubly robust causal inference model estimated the effects of montelukast, adjusting for confounders and censored patients. RESULTS: pwMS treated with montelukast demonstrated a statistically significant 23.6% reduction in relapses compared to non-users in 67.3% of emulated RCTs. CONCLUSION: Real-world evidence suggested that montelukast reduces MS relapses, warranting future clinical trials and further research on LTRAs' potential mechanism in MS.

Comparing diagnostic accuracy of current practice guidelines in predicting choledocholithiasis: outcomes from a large healthcare system comprising both academic and community settings.

BACKGROUND AND AIMS: The American Society for Gastrointestinal Endoscopy (ASGE) 2010 guidelines for suspected choledocholithiasis were recently updated by proposing more specific criteria for selection of high-risk patients to undergo direct ERCP while advocating the use of additional imaging studies for intermediate- and low-risk individuals. We aim to compare the performance and diagnostic accuracy of 2019 versus 2010 ASGE criteria for suspected choledocholithiasis. METHODS: We performed a retrospective chart review of a prospectively maintained database (2013-2019) of over 10,000 ERCPs performed by 70 gastroenterologists in our 14-hospital system. We randomly selected 744 ERCPs in which the primary indication was suspected choledocholithiasis. Patients with a history of cholecystectomy or prior sphincterotomy were excluded. The same patient cohort was assigned as low, intermediate, or high risk according to the 2010 and 2019 guideline criteria. Overall accuracy of both guidelines was compared against the presence of stones and/or sludge on ERCP. RESULTS: Of 744 patients who underwent ERCP, 544 patients (73.1%) had definite stones during ERCP and 696 patients (93.5%) had stones and/or sludge during ERCP. When classified according to the 2019 guidelines, fewer patients were high risk (274/744, 36.8%) compared with 2010 guidelines (449/744, 60.4%; P < .001). Within the high-risk group per both guidelines, definitive stone was found during ERCP more frequently in the 2019 guideline cohort (226/274, 82.5%) compared with the 2010 guideline cohort (342/449, 76.2%; P < .001). In our patient cohort, overall specificity of the 2010 guideline was 46.5%, which improved to 76.0% as per 2019 guideline criteria (P < .001). However, no significant change was noted for either positive predictive value or negative predictive value between 2019 and 2010 guidelines. CONCLUSIONS: The 2019 ASGE guidelines are more specific for detection of choledocholithiasis during ERCP when compared with the 2010 guidelines. However, a large number of patients are categorized as intermediate risk per 2019 guidelines and will require an additional confirmatory imaging study.

Population stratification enables modeling effects of reopening policies on mortality and hospitalization rates.

OBJECTIVE: Study the impact of local policies on near-future hospitalization and mortality rates. MATERIALS AND METHODS: We introduce a novel risk-stratified SIR-HCD model that introduces new variables to model the dynamics of low-contact (e.g., work from home) and high-contact (e.g., work on-site) subpopulations while sharing parameters to control their respective R0(t) over time. We test our model on data of daily reported hospitalizations and cumulative mortality of COVID-19 in Harris County, Texas, from May 1, 2020, until October 4, 2020, collected from multiple sources (USA FACTS, U.S. Bureau of Labor Statistics, Southeast Texas Regional Advisory Council COVID-19 report, TMC daily news, and Johns Hopkins University county-level mortality reporting). RESULTS: We evaluated our model's forecasting accuracy in Harris County, TX (the most populated county in the Greater Houston area) during Phase-I and Phase-II reopening. Not only does our model outperform other competing models, but it also supports counterfactual analysis to simulate the impact of future policies in a local setting, which is unique among existing approaches. DISCUSSION: Mortality and hospitalization rates are significantly impacted by local quarantine and reopening policies. Existing models do not directly account for the effect of these policies on infection, hospitalization, and death rates in an explicit and explainable manner. Our work is an attempt to improve prediction of these trends by incorporating this information into the model, thus supporting decision-making. CONCLUSION: Our work is a timely effort to attempt to model the dynamics of pandemics under the influence of local policies.

Environmental stresses disrupt telomere length homeostasis.

Telomeres protect the chromosome ends from degradation and play crucial roles in cellular aging and disease. Recent studies have additionally found a correlation between psychological stress, telomere length, and health outcome in humans. However, studies have not yet explored the causal relationship between stress and telomere length, or the molecular mechanisms underlying that relationship. Using yeast as a model organism, we show that stresses may have very different outcomes: alcohol and acetic acid elongate telomeres, whereas caffeine and high temperatures shorten telomeres. Additional treatments, such as oxidative stress, show no effect. By combining genome-wide expression measurements with a systematic genetic screen, we identify the Rap1/Rif1 pathway as the central mediator of the telomeric response to environmental signals. These results demonstrate that telomere length can be manipulated, and that a carefully regulated homeostasis may become markedly deregulated in opposing directions in response to different environmental cues.

Ranking adverse drug reactions with crowdsourcing.

BACKGROUND: There is no publicly available resource that provides the relative severity of adverse drug reactions (ADRs). Such a resource would be useful for several applications, including assessment of the risks and benefits of drugs and improvement of patient-centered care. It could also be used to triage predictions of drug adverse events. OBJECTIVE: The intent of the study was to rank ADRs according to severity. METHODS: We used Internet-based crowdsourcing to rank ADRs according to severity. We assigned 126,512 pairwise comparisons of ADRs to 2589 Amazon Mechanical Turk workers and used these comparisons to rank order 2929 ADRs. RESULTS: There is good correlation (rho=.53) between the mortality rates associated with ADRs and their rank. Our ranking highlights severe drug-ADR predictions, such as cardiovascular ADRs for raloxifene and celecoxib. It also triages genes associated with severe ADRs such as epidermal growth-factor receptor (EGFR), associated with glioblastoma multiforme, and SCN1A, associated with epilepsy. CONCLUSIONS: ADR ranking lays a first stepping stone in personalized drug risk assessment. Ranking of ADRs using crowdsourcing may have useful clinical and financial implications, and should be further investigated in the context of health care decision making.

Generalizable pipeline for constructing HIV risk prediction models across electronic health record systems.

OBJECTIVE: The HIV epidemic remains a significant public health issue in the United States. HIV risk prediction models could be beneficial for reducing HIV transmission by helping clinicians identify patients at high risk for infection and refer them for testing. This would facilitate initiation on treatment for those unaware of their status and pre-exposure prophylaxis for those uninfected but at high risk. Existing HIV risk prediction algorithms rely on manual construction of features and are limited in their application across diverse electronic health record systems. Furthermore, the accuracy of these models in predicting HIV in females has thus far been limited. MATERIALS AND METHODS: We devised a pipeline for automatic construction of prediction models based on automatic feature engineering to predict HIV risk and tested our pipeline on a local electronic health records system and a national claims data. We also compared the performance of general models to female-specific models. RESULTS: Our models obtain similarly good performance on both health record datasets despite difference in represented populations and data availability (AUC = 0.87). Furthermore, our general models obtain good performance on females but are also improved by constructing female-specific models (AUC between 0.81 and 0.86 across datasets). DISCUSSION AND CONCLUSIONS: We demonstrated that flexible construction of prediction models performs well on HIV risk prediction across diverse health records systems and perform as well in predicting HIV risk in females, making deployment of such models into existing health care systems tangible.

Prediction of Mild Cognitive Impairment Status: Pilot Study of Machine Learning Models Based on Longitudinal Data From Fitness Trackers.

BACKGROUND: Early signs of Alzheimer disease (AD) are difficult to detect, causing diagnoses to be significantly delayed to time points when brain damage has already occurred and current experimental treatments have little effect on slowing disease progression. Tracking cognitive decline at early stages is critical for patients to make lifestyle changes and consider new and experimental therapies. Frequently studied biomarkers are invasive and costly and are limited for predicting conversion from normal to mild cognitive impairment (MCI). OBJECTIVE: This study aimed to use data collected from fitness trackers to predict MCI status. METHODS: In this pilot study, fitness trackers were worn by 20 participants: 12 patients with MCI and 8 age-matched controls. We collected physical activity, heart rate, and sleep data from each participant for up to 1 month and further developed a machine learning model to predict MCI status. RESULTS: Our machine learning model was able to perfectly separate between MCI and controls (area under the curve=1.0). The top predictive features from the model included peak, cardio, and fat burn heart rate zones; resting heart rate; average deep sleep time; and total light activity time. CONCLUSIONS: Our results suggest that a longitudinal digital biomarker differentiates between controls and patients with MCI in a very cost-effective and noninvasive way and hence may be very useful for identifying patients with very early AD who can benefit from clinical trials and new, disease-modifying therapies.

Unsupervised deep representation learning enables phenotype discovery for genetic association studies of brain imaging.

Understanding the genetic architecture of brain structure is challenging, partly due to difficulties in designing robust, non-biased descriptors of brain morphology. Until recently, brain measures for genome-wide association studies (GWAS) consisted of traditionally expert-defined or software-derived image-derived phenotypes (IDPs) that are often based on theoretical preconceptions or computed from limited amounts of data. Here, we present an approach to derive brain imaging phenotypes using unsupervised deep representation learning. We train a 3-D convolutional autoencoder model with reconstruction loss on 6130 UK Biobank (UKBB) participants' T1 or T2-FLAIR (T2) brain MRIs to create a 128-dimensional representation known as Unsupervised Deep learning derived Imaging Phenotypes (UDIPs). GWAS of these UDIPs in held-out UKBB subjects (n = 22,880 discovery and n = 12,359/11,265 replication cohorts for T1/T2) identified 9457 significant SNPs organized into 97 independent genetic loci of which 60 loci were replicated. Twenty-six loci were not reported in earlier T1 and T2 IDP-based UK Biobank GWAS. We developed a perturbation-based decoder interpretation approach to show that these loci are associated with UDIPs mapped to multiple relevant brain regions. Our results established unsupervised deep learning can derive robust, unbiased, heritable, and interpretable brain imaging phenotypes.

A method for inferring medical diagnoses from patient similarities.

BACKGROUND: Clinical decision support systems assist physicians in interpreting complex patient data. However, they typically operate on a per-patient basis and do not exploit the extensive latent medical knowledge in electronic health records (EHRs). The emergence of large EHR systems offers the opportunity to integrate population information actively into these tools. METHODS: Here, we assess the ability of a large corpus of electronic records to predict individual discharge diagnoses. We present a method that exploits similarities between patients along multiple dimensions to predict the eventual discharge diagnoses. RESULTS: Using demographic, initial blood and electrocardiography measurements, as well as medical history of hospitalized patients from two independent hospitals, we obtained high performance in cross-validation (area under the curve >0.88) and correctly predicted at least one diagnosis among the top ten predictions for more than 84% of the patients tested. Importantly, our method provides accurate predictions (>0.86 precision in cross validation) for major disease categories, including infectious and parasitic diseases, endocrine and metabolic diseases and diseases of the circulatory systems. Our performance applies to both chronic and acute diagnoses. CONCLUSIONS: Our results suggest that one can harness the wealth of population-based information embedded in electronic health records for patient-specific predictive tasks.

INDI: a computational framework for inferring drug interactions and their associated recommendations.

Inferring drug-drug interactions (DDIs) is an essential step in drug development and drug administration. Most computational inference methods focus on modeling drug pharmacokinetics, aiming at interactions that result from a common metabolizing enzyme (CYP). Here, we introduce a novel prediction method, INDI (INferring Drug Interactions), allowing the inference of both pharmacokinetic, CYP-related DDIs (along with their associated CYPs) and pharmacodynamic, non-CYP associated ones. On cross validation, it obtains high specificity and sensitivity levels (AUC (area under the receiver-operating characteristic curve) ≥0.93). In application to the FDA adverse event reporting system, 53% of the drug events could potentially be connected to known (41%) or predicted (12%) DDIs. Additionally, INDI predicts the severity level of each DDI upon co-administration of the involved drugs, suggesting that severe interactions are abundant in the clinical practice. Examining regularly taken medications by hospitalized patients, 18% of the patients receive known or predicted severely interacting drugs and are hospitalized more frequently. Access to INDI and its predictions is provided via a web tool at http://www.cs.tau.ac.il/~bnet/software/INDI, facilitating the inference and exploration of drug interactions and providing important leads for physicians and pharmaceutical companies alike.

PRINCIPLE: a tool for associating genes with diseases via network propagation.

SUMMARY: PRINCIPLE is a Java application implemented as a Cytoscape plug-in, based on a previously published algorithm, PRINCE. Given a query disease, it prioritizes disease-related genes based on their closeness in a protein-protein interaction network to genes causing phenotypically similar disorders to the query disease. AVAILABILITY: Implemented in Java, PRINCIPLE runs over Cytoscape 2.7 or newer versions. Binaries, default input files and documentation are freely available at http://www.cs.tau.ac.il/~bnet/software/PrincePlugin/. CONTACT: roded@tau.ac.il; assafgot@tau.ac.il.

PREDICT: a method for inferring novel drug indications with application to personalized medicine.

Inferring potential drug indications, for either novel or approved drugs, is a key step in drug development. Previous computational methods in this domain have focused on either drug repositioning or matching drug and disease gene expression profiles. Here, we present a novel method for the large-scale prediction of drug indications (PREDICT) that can handle both approved drugs and novel molecules. Our method is based on the observation that similar drugs are indicated for similar diseases, and utilizes multiple drug-drug and disease-disease similarity measures for the prediction task. On cross-validation, it obtains high specificity and sensitivity (AUC=0.9) in predicting drug indications, surpassing existing methods. We validate our predictions by their overlap with drug indications that are currently under clinical trials, and by their agreement with tissue-specific expression information on the drug targets. We further show that disease-specific genetic signatures can be used to accurately predict drug indications for new diseases (AUC=0.92). This lays the computational foundation for future personalized drug treatments, where gene expression signatures from individual patients would replace the disease-specific signatures.

SynPathy: Predicting Drug Synergy through Drug-Associated Pathways Using Deep Learning.

Drug combination therapy has become a promising therapeutic strategy for cancer treatment. While high-throughput drug combination screening is effective for identifying synergistic drug combinations, measuring all possible combinations is impractical due to the vast space of therapeutic agents and cell lines. In this study, we propose a biologically-motivated deep learning approach to identify pathway-level features from drug and cell lines' molecular data for predicting drug synergy and quantifying the interactions in synergistic drug pairs. This method obtained an MSE of 70.6 ± 6.4, significantly surpassing previous approaches while providing potential candidate pathways to explain the prediction. We further demonstrate that drug combinations tend to be more synergistic when their top contributing pathways are closer to each other on a protein interaction network, suggesting a potential strategy for combination therapy with topologically interacting pathways. Our computational approach can thus be utilized both for prescreening of potential drug combinations and for designing new combinations based on proximity of pathways associated with drug targets and cell lines. IMPLICATIONS: Our computational framework may be translated in the future to clinical scenarios where synergistic drugs are tailored to the patient and additionally, drug development could benefit from designing drugs that target topologically close pathways.

Tissue-Specific Variations in Transcription Factors Elucidate Complex Immune System Regulation.

Gene expression plays a key role in health and disease. Estimating the genetic components underlying gene expression can thus help understand disease etiology. Polygenic models termed "transcriptome imputation" are used to estimate the genetic component of gene expression, but these models typically consider only the cis regions of the gene. However, these cis-based models miss large variability in expression for multiple genes. Transcription factors (TFs) that regulate gene expression are natural candidates for looking for additional sources of the missing variability. We developed a hypothesis-driven approach to identify second-tier regulation by variability in TFs. Our approach tested two models representing possible mechanisms by which variations in TFs can affect gene expression: variability in the expression of the TF and genetic variants within the TF that may affect the binding affinity of the TF to the TF-binding site. We tested our TF models in whole blood and skeletal muscle tissues and identified TF variability that can partially explain missing gene expression for 1035 genes, 76% of which explains more than the cis-based models. While the discovered regulation patterns were tissue-specific, they were both enriched for immune system functionality, elucidating complex regulation patterns. Our hypothesis-driven approach is useful for identifying tissue-specific genetic regulation patterns involving variations in TF expression or binding.

Brain Antigens Stimulate Proliferation of T Lymphocytes With a Pathogenic Phenotype in Multiple Sclerosis Patients.

A method to stimulate T lymphocytes with a broad range of brain antigens would facilitate identification of the autoantigens for multiple sclerosis and enable definition of the pathogenic mechanisms important for multiple sclerosis. In a previous work, we found that the obvious approach of culturing leukocytes with homogenized brain tissue does not work because the brain homogenate suppresses antigen-specific lymphocyte proliferation. We now report a method that substantially reduces the suppressive activity. We used this non-suppressive brain homogenate to stimulate leukocytes from multiple sclerosis patients and controls. We also stimulated with common viruses for comparison. We measured proliferation, selected the responding CD3+ cells with flow cytometry, and sequenced their transcriptomes for mRNA and T-cell receptor sequences. The mRNA expression suggested that the brain-responding cells from MS patients are potentially pathogenic. The T-cell receptor repertoire of the brain-responding cells was clonal with minimal overlap with virus antigens.

Molecular pathways enhance drug response prediction using transfer learning from cell lines to tumors and patient-derived xenografts.

Computational models have been successful in predicting drug sensitivity in cancer cell line data, creating an opportunity to guide precision medicine. However, translating these models to tumors remains challenging. We propose a new transfer learning workflow that transfers drug sensitivity predicting models from large-scale cancer cell lines to both tumors and patient derived xenografts based on molecular pathways derived from genomic features. We further compute feature importance to identify pathways most important to drug response prediction. We obtained good performance on tumors (AUROC = 0.77) and patient derived xenografts from triple negative breast cancers (RMSE = 0.11). Using feature importance, we highlight the association between ER-Golgi trafficking pathway in everolimus sensitivity within breast cancer patients and the role of class II histone deacetylases and interlukine-12 in response to drugs for triple-negative breast cancer. Pathway information support transfer of drug response prediction models from cell lines to tumors and can provide biological interpretation underlying the predictions, serving as a steppingstone towards usage in clinical setting.

A Method for Bridging Population-Specific Genotypes to Detect Gene Modules Associated with Alzheimer's Disease.

BACKGROUND: Genome-wide association studies have successfully identified variants associated with multiple conditions. However, generalizing discoveries across diverse populations remains challenging due to large variations in genetic composition. Methods that perform gene expression imputation have attempted to address the transferability of gene discoveries across populations, but with limited success. METHODS: Here, we introduce a pipeline that combines gene expression imputation with gene module discovery, including a dense gene module search and a gene set variation analysis, to address the transferability issue. Our method feeds association probabilities of imputed gene expression with a selected phenotype into tissue-specific gene-module discovery over protein interaction networks to create higher-level gene modules. RESULTS: We demonstrate our method's utility in three case-control studies of Alzheimer's disease (AD) for three different race/ethnic populations (Whites, African descent and Hispanics). We discovered 182 AD-associated genes from gene modules shared between these populations, highlighting new gene modules associated with AD. CONCLUSIONS: Our innovative framework has the potential to identify robust discoveries across populations based on gene modules, as demonstrated in AD.