Pervasive mislocalization of pathogenic coding variants underlying human disorders.

Widespread sequencing has yielded thousands of missense variants predicted or confirmed as disease causing. This creates a new bottleneck: determining the functional impact of each variant-typically a painstaking, customized process undertaken one or a few genes and variants at a time. Here, we established a high-throughput imaging platform to assay the impact of coding variation on protein localization, evaluating 3,448 missense variants of over 1,000 genes and phenotypes. We discovered that mislocalization is a common consequence of coding variation, affecting about one-sixth of all pathogenic missense variants, all cellular compartments, and recessive and dominant disorders alike. Mislocalization is primarily driven by effects on protein stability and membrane insertion rather than disruptions of trafficking signals or specific interactions. Furthermore, mislocalization patterns help explain pleiotropy and disease severity and provide insights on variants of uncertain significance. Our publicly available resource extends our understanding of coding variation in human diseases.

AD-Syn-Net: systematic identification of Alzheimer's disease-associated mutation and co-mutation vulnerabilities via deep learning.

Alzheimer's disease (AD) is one of the most challenging neurodegenerative diseases because of its complicated and progressive mechanisms, and multiple risk factors. Increasing research evidence demonstrates that genetics may be a key factor responsible for the occurrence of the disease. Although previous reports identified quite a few AD-associated genes, they were mostly limited owing to patient sample size and selection bias. There is a lack of comprehensive research aimed to identify AD-associated risk mutations systematically. To address this challenge, we hereby construct a large-scale AD mutation and co-mutation framework ('AD-Syn-Net'), and propose deep learning models named Deep-SMCI and Deep-CMCI configured with fully connected layers that are capable of predicting cognitive impairment of subjects effectively based on genetic mutation and co-mutation profiles. Next, we apply the customized frameworks to data sets to evaluate the importance scores of the mutations and identified mutation effectors and co-mutation combination vulnerabilities contributing to cognitive impairment. Furthermore, we evaluate the influence of mutation pairs on the network architecture to dissect the genetic organization of AD and identify novel co-mutations that could be responsible for dementia, laying a solid foundation for proposing future targeted therapy for AD precision medicine. Our deep learning model codes are available open access here: https://github.com/Pan-Bio/AD-mutation-effectors.

mi-IsoNet: systems-scale microRNA landscape reveals rampant isoform-mediated gain of target interaction diversity and signaling specificity.

MicroRNA (miRNA) is not a single sequence, but a series of multiple variants (also termed isomiRs) with sequence and expression heterogeneity. Whether and how these isoforms contribute to functional variation and complexity at the systems and network levels remain largely unknown. To explore this question systematically, we comprehensively analyzed the expression of small RNAs and their target sites to interrogate functional variations between novel isomiRs and their canonical miRNA sequences. Our analyses of the pan-cancer landscape of miRNA expression indicate that multiple isomiRs generated from the same miRNA locus often exhibit remarkable variation in their sequence, expression and function. We interrogated abundant and differentially expressed 5' isomiRs with novel seed sequences via seed shifting and identified many potential novel targets of these 5' isomiRs that would expand interaction capabilities between small RNAs and mRNAs, rewiring regulatory networks and increasing signaling circuit complexity. Further analyses revealed that some miRNA loci might generate diverse dominant isomiRs that often involved isomiRs with varied seeds and arm-switching, suggesting a selective advantage of multiple isomiRs in regulating gene expression. Finally, experimental validation indicated that isomiRs with shifted seed sequences could regulate novel target mRNAs and therefore contribute to regulatory network rewiring. Our analysis uncovers a widespread expansion of isomiR and mRNA interaction networks compared with those seen in canonical small RNA analysis; this expansion suggests global gene regulation network perturbations by alternative small RNA variants or isoforms. Taken together, the variations in isomiRs that occur during miRNA processing and maturation are likely to play a far more complex and plastic role in gene regulation than previously anticipated.

e-MutPath: computational modeling reveals the functional landscape of genetic mutations rewiring interactome networks.

Understanding the functional impact of cancer somatic mutations represents a critical knowledge gap for implementing precision oncology. It has been increasingly appreciated that the interaction profile mediated by a genomic mutation provides a fundamental link between genotype and phenotype. However, specific effects on biological signaling networks for the majority of mutations are largely unknown by experimental approaches. To resolve this challenge, we developed e-MutPath (edgetic Mutation-mediated Pathway perturbations), a network-based computational method to identify candidate 'edgetic' mutations that perturb functional pathways. e-MutPath identifies informative paths that could be used to distinguish disease risk factors from neutral elements and to stratify disease subtypes with clinical relevance. The predicted targets are enriched in cancer vulnerability genes, known drug targets but depleted for proteins associated with side effects, demonstrating the power of network-based strategies to investigate the functional impact and perturbation profiles of genomic mutations. Together, e-MutPath represents a robust computational tool to systematically assign functions to genetic mutations, especially in the context of their specific pathway perturbation effect.

Timing of Resumption of Anticoagulation After Polypectomy and Frequency of Post-procedural Complications: A Post-hoc Analysis.

BACKGROUND: Optimal timing for anticoagulation resumption after polypectomy is unclear. We explored the association between timing of anticoagulation resumption and occurrence of delayed post-polypectomy bleeding (PPB) and thromboembolic (TE) events. METHODS: We performed a post-hoc analysis of patients in an earlier study whose anticoagulants were interrupted for polypectomy. We compared rates of clinically important delayed PPB and TE events in relationship to timing of anticoagulant resumption. Late resumption was defined as > 2 days after polypectomy. RESULTS: Among 437 patients, 351 had early and 86 late resumption. Compared to early resumers, late resumers had greater polypectomy complexity. PPB rate was higher (but not significantly) in the late versus early resumers (2.3% vs. 0.9%, 1.47% greater, 95% CI [- 2.58 to 5.52], p = 0.26). TE events were more frequent in late versus early resumers [0% vs. 1.2% at 30 days, 0% vs. 2.3%, 95% CI 0.3-8, (p = 0.04) at 90 days]. On multivariate analysis, timing of restarting anticoagulation was not a significant predictor of PPB (OR 0.97, 95% CI 0.61-1.44, p = 0.897). Significant predictors were number of polyps ≥ 1 cm (OR 4.14, 95% CI 1.27-13.66, p = 0.014) and use of fulguration (OR 11.43, 95% CI 1.35-80.80, p = 0.014). CONCLUSIONS: Physicians delayed anticoagulation resumption more commonly after complex polypectomies. The timing of restarting anticoagulation was not a significant risk factor for PPB and late resumers had significantly higher rates of TE events within 90 days. Considering the potentially catastrophic consequences of TE events and the generally benign outcome of PPBs, clinicians should be cautious about delaying resumption of anticoagulation after polypectomy.

Exploration of the intracellular chiral metabolome in pediatric BCP-ALL: a pilot study investigating the metabolic phenotype of IgH locus aberrations.

Background and aims: Aberrations in the immunoglobulin heavy chain (IgH) locus are associated with poor prognosis in pediatric precursor B-cell acute lymphoblastic leukemia (BCP-ALL) patients. The primary objective of this pilot study is to enhance our understanding of the IgH phenotype by exploring the intracellular chiral metabolome. Materials and methods: Leukemia cells were isolated from the bone marrow of BCP-ALL pediatric patients at diagnosis. The samples' metabolome and transcriptome were characterized using untargeted chiral metabolomic and next-generation sequencing transcriptomic analyses. Results: For the first time D- amino acids were identified in the leukemic cells' intracellular metabolome from the bone marrow niche. Chiral metabolic signatures at diagnosis was indicative of a resistant phenotype. Through integrated network analysis and Pearson correlation, confirmation was obtained regarding the association of the IgH phenotype with several genes linked to poor prognosis. Conclusion: The findings of this study have contributed to the understanding that the chiral metabolome plays a role in the poor prognosis observed in an exceptionally rare patient cohort. The findings include elevated D-amino acid incorporation in the IgH group, the emergence of several unknown, potentially enantiomeric, metabolites, and insights into metabolic pathways that all warrant further exploration.

Spatial imaging of glycoRNA in single cells with ARPLA.

Little is known about the biological roles of glycosylated RNAs (glycoRNAs), a recently discovered class of glycosylated molecules, because of a lack of visualization methods. We report sialic acid aptamer and RNA in situ hybridization-mediated proximity ligation assay (ARPLA) to visualize glycoRNAs in single cells with high sensitivity and selectivity. The signal output of ARPLA occurs only when dual recognition of a glycan and an RNA triggers in situ ligation, followed by rolling circle amplification of a complementary DNA, which generates a fluorescent signal by binding fluorophore-labeled oligonucleotides. Using ARPLA, we detect spatial distributions of glycoRNAs on the cell surface and their colocalization with lipid rafts as well as the intracellular trafficking of glycoRNAs through SNARE protein-mediated secretory exocytosis. Studies in breast cell lines suggest that surface glycoRNA is inversely associated with tumor malignancy and metastasis. Investigation of the relationship between glycoRNAs and monocyte-endothelial cell interactions suggests that glycoRNAs may mediate cell-cell interactions during the immune response.

Population-level comparisons of gene regulatory networks modeled on high-throughput single-cell transcriptomics data.

Single-cell technologies enable high-resolution studies of phenotype-defining molecular mechanisms. However, data sparsity and cellular heterogeneity make modeling biological variability across single-cell samples difficult. Here we present SCORPION, a tool that uses a message-passing algorithm to reconstruct comparable gene regulatory networks from single-cell/nuclei RNA-sequencing data that are suitable for population-level comparisons by leveraging the same baseline priors. Using synthetic data, we found that SCORPION outperformed 12 existing gene regulatory network reconstruction techniques. Using supervised experiments, we show that SCORPION can accurately identify differences in regulatory networks between wild-type and transcription factor-perturbed cells. We demonstrate SCORPION's scalability to population-level analyses using a single-cell RNA-sequencing atlas containing 200,436 cells from colorectal cancer and adjacent healthy tissues. The differences between tumor regions detected by SCORPION are consistent across multiple cohorts as well as with our understanding of disease progression, and elucidate phenotypic regulators that may impact patient survival.

Integrated multi-omics analyses identify anti-viral host factors and pathways controlling SARS-CoV-2 infection.

Host anti-viral factors are essential for controlling SARS-CoV-2 infection but remain largely unknown due to the biases of previous large-scale studies toward pro-viral host factors. To fill in this knowledge gap, we perform a genome-wide CRISPR dropout screen and integrate analyses of the multi-omics data of the CRISPR screen, genome-wide association studies, single-cell RNA-Seq, and host-virus proteins or protein/RNA interactome. This study uncovers many host factors that are currently underappreciated, including the components of V-ATPases, ESCRT, and N-glycosylation pathways that modulate viral entry and/or replication. The cohesin complex is also identified as an anti-viral pathway, suggesting an important role of three-dimensional chromatin organization in mediating host-viral interaction. Furthermore, we discover another anti-viral regulator KLF5, a transcriptional factor involved in sphingolipid metabolism, which is up-regulated, and harbors genetic variations linked to COVID-19 patients with severe symptoms. Anti-viral effects of three identified candidates (DAZAP2/VTA1/KLF5) are confirmed individually. Molecular characterization of DAZAP2/VTA1/KLF5-knockout cells highlights the involvement of genes related to the coagulation system in determining the severity of COVID-19. Together, our results provide further resources for understanding the host anti-viral network during SARS-CoV-2 infection and may help develop new countermeasure strategies.

Rapid iPSC inclusionopathy models shed light on formation, consequence, and molecular subtype of α-synuclein inclusions.

The heterogeneity of protein-rich inclusions and its significance in neurodegeneration is poorly understood. Standard patient-derived iPSC models develop inclusions neither reproducibly nor in a reasonable time frame. Here, we developed screenable iPSC "inclusionopathy" models utilizing piggyBac or targeted transgenes to rapidly induce CNS cells that express aggregation-prone proteins at brain-like levels. Inclusions and their effects on cell survival were trackable at single-inclusion resolution. Exemplar cortical neuron α-synuclein inclusionopathy models were engineered through transgenic expression of α-synuclein mutant forms or exogenous seeding with fibrils. We identified multiple inclusion classes, including neuroprotective p62-positive inclusions versus dynamic and neurotoxic lipid-rich inclusions, both identified in patient brains. Fusion events between these inclusion subtypes altered neuronal survival. Proteome-scale α-synuclein genetic- and physical-interaction screens pinpointed candidate RNA-processing and actin-cytoskeleton-modulator proteins like RhoA whose sequestration into inclusions could enhance toxicity. These tractable CNS models should prove useful in functional genomic analysis and drug development for proteinopathies.

Transcriptomic analyses of patient peripheral blood with hemoglobin depletion reveal glioblastoma biomarkers.

Peripheral blood is gaining prominence as a noninvasive alternative to tissue biopsy to develop biomarkers for glioblastoma (GBM); however, widely utilized blood-based biomarkers in clinical settings have not yet been identified due to the lack of a robust detection approach. Here, we describe the application of globin reduction in RNA sequencing of whole blood (i.e., WBGR) and perform transcriptomic analysis to identify GBM-associated transcriptomic changes. By using WBGR, we improved the detection sensitivity of informatic reads and identified differential gene expression in GBM blood. By analyzing tumor tissues, we identified transcriptomic traits of GBM blood. Further functional enrichment analyses retained the most changed genes in GBM. Subsequent validation elicited a 10-gene panel covering mRNA, long noncoding RNA, and microRNA (i.e., GBM-Dx panel) that has translational potential to aid in the early detection or clinical management of GBM. Here, we report an integrated approach, WBGR, with comprehensive analytic capacity for blood-based marker identification.

The HAPSTR2 retrogene buffers stress signaling and resilience in mammals.

We recently identified HAPSTR1 (C16orf72) as a key component in a novel pathway which regulates the cellular response to molecular stressors, such as DNA damage, nutrient scarcity, and protein misfolding. Here, we identify a functional paralog to HAPSTR1: HAPSTR2. HAPSTR2 formed early in mammalian evolution, via genomic integration of a reverse transcribed HAPSTR1 transcript, and has since been preserved under purifying selection. HAPSTR2, expressed primarily in neural and germline tissues and a subset of cancers, retains established biochemical features of HAPSTR1 to achieve two functions. In normal physiology, HAPSTR2 directly interacts with HAPSTR1, markedly augmenting HAPSTR1 protein stability in a manner independent from HAPSTR1's canonical E3 ligase, HUWE1. Alternatively, in the context of HAPSTR1 loss, HAPSTR2 expression is sufficient to buffer stress signaling and resilience. Thus, we discover a mammalian retrogene which safeguards fitness.

AI-DrugNet: A network-based deep learning model for drug repurposing and combination therapy in neurological disorders.

Discovering effective therapies is difficult for neurological and developmental disorders in that disease progression is often associated with a complex and interactive mechanism. Over the past few decades, few drugs have been identified for treating Alzheimer's disease (AD), especially for impacting the causes of cell death in AD. Although drug repurposing is gaining more success in developing therapeutic efficacy for complex diseases such as common cancer, the complications behind AD require further study. Here, we developed a novel prediction framework based on deep learning to identify potential repurposed drug therapies for AD, and more importantly, our framework is broadly applicable and may generalize to identifying potential drug combinations in other diseases. Our prediction framework is as follows: we first built a drug-target pair (DTP) network based on multiple drug features and target features, as well as the associations between DTP nodes where drug-target pairs are the DTP nodes and the associations between DTP nodes are represented as the edges in the AD disease network; furthermore, we incorporated the drug-target feature from the DTP network and the relationship information between drug-drug, target-target, drug-target within and outside of drug-target pairs, representing each drug-combination as a quartet to generate corresponding integrated features; finally, we developed an AI-based Drug discovery Network (AI-DrugNet), which exhibits robust predictive performance. The implementation of our network model help identify potential repurposed and combination drug options that may serve to treat AD and other diseases.

Gain-of-Function Variomics and Multi-omics Network Biology for Precision Medicine.

Traditionally, disease causal mutations were thought to disrupt gene function. However, it becomes more clear that many deleterious mutations could exhibit a "gain-of-function" (GOF) behavior. Systematic investigation of such mutations has been lacking and largely overlooked. Advances in next-generation sequencing have identified thousands of genomic variants that perturb the normal functions of proteins, further contributing to diverse phenotypic consequences in disease. Elucidating the functional pathways rewired by GOF mutations will be crucial for prioritizing disease-causing variants and their resultant therapeutic liabilities. In distinct cell types (with varying genotypes), precise signal transduction controls cell decision, including gene regulation and phenotypic output. When signal transduction goes awry due to GOF mutations, it would give rise to various disease types. Quantitative and molecular understanding of network perturbations by GOF mutations may provide explanations for 'missing heritability" in previous genome-wide association studies. We envision that it will be instrumental to push current paradigm toward a thorough functional and quantitative modeling of all GOF mutations and their mechanistic molecular events involved in disease development and progression. Many fundamental questions pertaining to genotype-phenotype relationships remain unresolved. For example, which GOF mutations are key for gene regulation and cellular decisions? What are the GOF mechanisms at various regulation levels? How do interaction networks undergo rewiring upon GOF mutations? Is it possible to leverage GOF mutations to reprogram signal transduction in cells, aiming to cure disease? To begin to address these questions, we will cover a wide range of topics regarding GOF disease mutations and their characterization by multi-omic networks. We highlight the fundamental function of GOF mutations and discuss the potential mechanistic effects in the context of signaling networks. We also discuss advances in bioinformatic and computational resources, which will dramatically help with studies on the functional and phenotypic consequences of GOF mutations.

Alternative splicing of GSDMB modulates killer lymphocyte-triggered pyroptosis.

Granzyme A from killer lymphocytes cleaves gasdermin B (GSDMB) and triggers pyroptosis in targeted human tumor cells, eliciting antitumor immunity. However, GSDMB has a controversial role in pyroptosis and has been linked to both anti- and protumor functions. Here, we found that GSDMB splicing variants are functionally distinct. Cleaved N-terminal (NT) fragments of GSDMB isoforms 3 and 4 caused pyroptosis, but isoforms 1, 2, and 5 did not. The nonfunctional isoforms have a deleted or modified exon 6 and therefore lack a stable belt motif. The belt likely contributes to the insertion of oligomeric GSDMB-NTs into the membrane. Consistently, noncytotoxic GSDMB-NTs blocked pyroptosis caused by cytotoxic GSDMB-NTs in a dominant-negative manner. Upon natural killer (NK) cell attack, GSDMB3-expressing cells died by pyroptosis, whereas GSDMB4-expressing cells died by mixed pyroptosis and apoptosis, and GSDMB1/2-expressing cells died only by apoptosis. GSDMB4 partially resisted NK cell-triggered cleavage, suggesting that only GSDMB3 is fully functional. GSDMB1-3 were the most abundant isoforms in the tested tumor cell lines and were similarly induced by interferon-γ and the chemotherapy drug methotrexate. Expression of cytotoxic GSDMB3/4 isoforms, but not GSDMB1/2 isoforms that are frequently up-regulated in tumors, was associated with better outcomes in bladder and cervical cancers, suggesting that GSDMB3/4-mediated pyroptosis was protective in those tumors. Our study indicates that tumors may block and evade killer cell-triggered pyroptosis by generating noncytotoxic GSDMB isoforms. Therefore, therapeutics that favor the production of cytotoxic GSDMB isoforms by alternative splicing may improve antitumor immunity.

Distinct astrocytic modulatory roles in sensory transmission during sleep, wakefulness, and arousal states in freely moving mice.

Despite extensive research on astrocytic Ca2+ in synaptic transmission, its contribution to the modulation of sensory transmission during different brain states remains largely unknown. Here, by using two-photon microscopy and whole-cell recordings, we show two distinct astrocytic Ca2+ signals in the murine barrel cortex: a small, long-lasting Ca2+ increase during sleep and a large, widespread but short-lasting Ca2+ spike when aroused. The large Ca2+ wave in aroused mice was inositol trisphosphate (IP3)-dependent, evoked by the locus coeruleus-norepinephrine system, and enhanced sensory input, contributing to reliable sensory transmission. However, the small Ca2+ transient was IP3-independent and contributed to decreased extracellular K+, hyperpolarization of the neurons, and suppression of sensory transmission. These events respond to different pharmacological inputs and contribute to distinct sleep and arousal functions by modulating the efficacy of sensory transmission. Together, our data demonstrate an important function for astrocytes in sleep and arousal states via astrocytic Ca2+ waves.

Widespread BRCA1/2-independent homologous recombination defects are caused by alterations in RNA-binding proteins.

Defects in homologous recombination DNA repair (HRD) both predispose to cancer development and produce therapeutic vulnerabilities, making it critical to define the spectrum of genetic events that cause HRD. However, we found that mutations in BRCA1/2 and other canonical HR genes only identified 10%-20% of tumors that display genomic evidence of HRD. Using a networks-based approach, we discovered that over half of putative genes causing HRD originated outside of canonical DNA damage response genes, with a particular enrichment for RNA-binding protein (RBP)-encoding genes. These putative drivers of HRD were experimentally validated, cross-validated in an independent cohort, and enriched in cancer-associated genome-wide association study loci. Mechanistic studies indicate that some RBPs are recruited to sites of DNA damage to facilitate repair, whereas others control the expression of canonical HR genes. Overall, this study greatly expands the repertoire of known drivers of HRD, with implications for basic biology, genetic screening, and therapy stratification.

Pervasive mislocalization of pathogenic coding variants underlying human disorders.

Widespread sequencing has yielded thousands of missense variants predicted or confirmed as disease-causing. This creates a new bottleneck: determining the functional impact of each variant - largely a painstaking, customized process undertaken one or a few genes or variants at a time. Here, we established a high-throughput imaging platform to assay the impact of coding variation on protein localization, evaluating 3,547 missense variants of over 1,000 genes and phenotypes. We discovered that mislocalization is a common consequence of coding variation, affecting about one-sixth of all pathogenic missense variants, all cellular compartments, and recessive and dominant disorders alike. Mislocalization is primarily driven by effects on protein stability and membrane insertion rather than disruptions of trafficking signals or specific interactions. Furthermore, mislocalization patterns help explain pleiotropy and disease severity and provide insights on variants of unknown significance. Our publicly available resource will likely accelerate the understanding of coding variation in human diseases.

Landscape of MicroRNA Regulatory Network Architecture and Functional Rerouting in Cancer.

Somatic mutations are a major source of cancer development, and many driver mutations have been identified in protein coding regions. However, the function of mutations located in miRNA and their target binding sites throughout the human genome remains largely unknown. Here, we built detailed cancer-specific miRNA regulatory networks across 30 cancer types to systematically analyze the effect of mutations in miRNAs and their target sites in 3' untranslated region (3' UTR), coding sequence (CDS), and 5' UTR regions. A total of 3,518,261 mutations from 9,819 samples were mapped to miRNA-gene interactions (mGI). Mutations in miRNAs showed a mutually exclusive pattern with mutations in their target genes in almost all cancer types. A linear regression method identified 148 candidate driver mutations that can significantly perturb miRNA regulatory networks. Driver mutations in 3'UTRs played their roles by altering RNA binding energy and the expression of target genes. Finally, mutated driver gene targets in 3' UTRs were significantly downregulated in cancer and functioned as tumor suppressors during cancer progression, suggesting potential miRNA candidates with significant clinical implications. A user-friendly, open-access web portal (mGI-map) was developed to facilitate further use of this data resource. Together, these results will facilitate novel noncoding biomarker identification and therapeutic drug design targeting the miRNA regulatory networks. SIGNIFICANCE: A detailed miRNA-gene interaction map reveals extensive miRNA-mediated gene regulatory networks with mutation-induced perturbations across multiple cancers, serving as a resource for noncoding biomarker discovery and drug development.

PRMT blockade induces defective DNA replication stress response and synergizes with PARP inhibition.

Multiple cancers exhibit aberrant protein arginine methylation by both type I arginine methyltransferases, predominately protein arginine methyltransferase 1 (PRMT1) and to a lesser extent PRMT4, and by type II PRMTs, predominately PRMT5. Here, we perform targeted proteomics following inhibition of PRMT1, PRMT4, and PRMT5 across 12 cancer cell lines. We find that inhibition of type I and II PRMTs suppresses phosphorylated and total ATR in cancer cells. Loss of ATR from PRMT inhibition results in defective DNA replication stress response activation, including from PARP inhibitors. Inhibition of type I and II PRMTs is synergistic with PARP inhibition regardless of homologous recombination function, but type I PRMT inhibition is more toxic to non-malignant cells. Finally, we demonstrate that the combination of PARP and PRMT5 inhibition improves survival in both BRCA-mutant and wild-type patient-derived xenografts without toxicity. Taken together, these results demonstrate that PRMT5 inhibition may be a well-tolerated approach to sensitize tumors to PARP inhibition.