Functional Genomic Landscape of Human Breast Cancer Drivers, Vulnerabilities, and Resistance.

Large-scale genomic studies have identified multiple somatic aberrations in breast cancer, including copy number alterations and point mutations. Still, identifying causal variants and emergent vulnerabilities that arise as a consequence of genetic alterations remain major challenges. We performed whole-genome small hairpin RNA (shRNA) "dropout screens" on 77 breast cancer cell lines. Using a hierarchical linear regression algorithm to score our screen results and integrate them with accompanying detailed genetic and proteomic information, we identify vulnerabilities in breast cancer, including candidate "drivers," and reveal general functional genomic properties of cancer cells. Comparisons of gene essentiality with drug sensitivity data suggest potential resistance mechanisms, effects of existing anti-cancer drugs, and opportunities for combination therapy. Finally, we demonstrate the utility of this large dataset by identifying BRD4 as a potential target in luminal breast cancer and PIK3CA mutations as a resistance determinant for BET-inhibitors.

The multimodality cell segmentation challenge: toward universal solutions.

Cell segmentation is a critical step for quantitative single-cell analysis in microscopy images. Existing cell segmentation methods are often tailored to specific modalities or require manual interventions to specify hyper-parameters in different experimental settings. Here, we present a multimodality cell segmentation benchmark, comprising more than 1,500 labeled images derived from more than 50 diverse biological experiments. The top participants developed a Transformer-based deep-learning algorithm that not only exceeds existing methods but can also be applied to diverse microscopy images across imaging platforms and tissue types without manual parameter adjustments. This benchmark and the improved algorithm offer promising avenues for more accurate and versatile cell analysis in microscopy imaging.

IL-7 coordinates proliferation, differentiation and Tcra recombination during thymocyte ß-selection.

Signaling via the pre-T cell antigen receptor (pre-TCR) and the receptor Notch1 induces transient self-renewal (ß-selection) of TCRß(+) CD4(-)CD8(-) double-negative stage 3 (DN3) and DN4 progenitor cells that differentiate into CD4(+)CD8(+) double-positive (DP) thymocytes, which then rearrange the locus encoding the TCR α-chain (Tcra). Interleukin 7 (IL-7) promotes the survival of TCRß(-) DN thymocytes by inducing expression of the pro-survival molecule Bcl-2, but the functions of IL-7 during ß-selection have remained unclear. Here we found that IL-7 signaled TCRß(+) DN3 and DN4 thymocytes to upregulate genes encoding molecules involved in cell growth and repressed the gene encoding the transcriptional repressor Bcl-6. Accordingly, IL-7-deficient DN4 cells lacked trophic receptors and did not proliferate but rearranged Tcra prematurely and differentiated rapidly. Deletion of Bcl6 partially restored the self-renewal of DN4 cells in the absence of IL-7, but overexpression of BCL2 did not. Thus, IL-7 critically acts cooperatively with signaling via the pre-TCR and Notch1 to coordinate proliferation, differentiation and Tcra recombination during ß-selection.

A roadmap for the Human Developmental Cell Atlas.

The Human Developmental Cell Atlas (HDCA) initiative, which is part of the Human Cell Atlas, aims to create a comprehensive reference map of cells during development. This will be critical to understanding normal organogenesis, the effect of mutations, environmental factors and infectious agents on human development, congenital and childhood disorders, and the cellular basis of ageing, cancer and regenerative medicine. Here we outline the HDCA initiative and the challenges of mapping and modelling human development using state-of-the-art technologies to create a reference atlas across gestation. Similar to the Human Genome Project, the HDCA will integrate the output from a growing community of scientists who are mapping human development into a unified atlas. We describe the early milestones that have been achieved and the use of human stem-cell-derived cultures, organoids and animal models to inform the HDCA, especially for prenatal tissues that are hard to acquire. Finally, we provide a roadmap towards a complete atlas of human development.

A clinically applicable integrative molecular classification of meningiomas.

Meningiomas are the most common primary intracranial tumour in adults1. Patients with symptoms are generally treated with surgery as there are no effective medical therapies. The World Health Organization histopathological grade of the tumour and the extent of resection at surgery (Simpson grade) are associated with the recurrence of disease; however, they do not accurately reflect the clinical behaviour of all meningiomas2. Molecular classifications of meningioma that reliably reflect tumour behaviour and inform on therapies are required. Here we introduce four consensus molecular groups of meningioma by combining DNA somatic copy-number aberrations, DNA somatic point mutations, DNA methylation and messenger RNA abundance in a unified analysis. These molecular groups more accurately predicted clinical outcomes compared with existing classification schemes. Each molecular group showed distinctive and prototypical biology (immunogenic, benign NF2 wild-type, hypermetabolic and proliferative) that informed therapeutic options. Proteogenomic characterization reinforced the robustness of the newly defined molecular groups and uncovered highly abundant and group-specific protein targets that we validated using immunohistochemistry. Single-cell RNA sequencing revealed inter-individual variations in meningioma as well as variations in intrinsic expression programs in neoplastic cells that mirrored the biology of the molecular groups identified.

Hypophosphorylated pRb knock-in mice exhibit hallmarks of aging and vitamin C-preventable diabetes.

Despite extensive analysis of pRB phosphorylation in vitro, how this modification influences development and homeostasis in vivo is unclear. Here, we show that homozygous Rb∆K4 and Rb∆K7 knock-in mice, in which either four or all seven phosphorylation sites in the C-terminal region of pRb, respectively, have been abolished by Ser/Thr-to-Ala substitutions, undergo normal embryogenesis and early development, notwithstanding suppressed phosphorylation of additional upstream sites. Whereas Rb∆K4 mice exhibit telomere attrition but no other abnormalities, Rb∆K7 mice are smaller and display additional hallmarks of premature aging including infertility, kyphosis, and diabetes, indicating an accumulative effect of blocking pRb phosphorylation. Diabetes in Rb∆K7 mice is insulin-sensitive and associated with failure of quiescent pancreatic ß-cells to re-enter the cell cycle in response to mitogens, resulting in induction of DNA damage response (DDR), senescence-associated secretory phenotype (SASP), and reduced pancreatic islet mass and circulating insulin level. Pre-treatment with the epigenetic regulator vitamin C reduces DDR, increases cell cycle re-entry, improves islet morphology, and attenuates diabetes. These results have direct implications for cell cycle regulation, CDK-inhibitor therapeutics, diabetes, and longevity.

Functional genomic landscape of cancer-intrinsic evasion of killing by T cells.

The genetic circuits that allow cancer cells to evade destruction by the host immune system remain poorly understood1-3. Here, to identify a phenotypically robust core set of genes and pathways that enable cancer cells to evade killing mediated by cytotoxic T lymphocytes (CTLs), we performed genome-wide CRISPR screens across a panel of genetically diverse mouse cancer cell lines that were cultured in the presence of CTLs. We identify a core set of 182 genes across these mouse cancer models, the individual perturbation of which increases either the sensitivity or the resistance of cancer cells to CTL-mediated toxicity. Systematic exploration of our dataset using genetic co-similarity reveals the hierarchical and coordinated manner in which genes and pathways act in cancer cells to orchestrate their evasion of CTLs, and shows that discrete functional modules that control the interferon response and tumour necrosis factor (TNF)-induced cytotoxicity are dominant sub-phenotypes. Our data establish a central role for genes that were previously identified as negative regulators of the type-II interferon response (for example, Ptpn2, Socs1 and Adar1) in mediating CTL evasion, and show that the lipid-droplet-related gene Fitm2 is required for maintaining cell fitness after exposure to interferon-γ (IFNγ). In addition, we identify the autophagy pathway as a conserved mediator of the evasion of CTLs by cancer cells, and show that this pathway is required to resist cytotoxicity induced by the cytokines IFNγ and TNF. Through the mapping of cytokine- and CTL-based genetic interactions, together with in vivo CRISPR screens, we show how the pleiotropic effects of autophagy control cancer-cell-intrinsic evasion of killing by CTLs and we highlight the importance of these effects within the tumour microenvironment. Collectively, these data expand our knowledge of the genetic circuits that are involved in the evasion of the immune system by cancer cells, and highlight genetic interactions that contribute to phenotypes associated with escape from killing by CTLs.

Drugst.One - a plug-and-play solution for online systems medicine and network-based drug repurposing.

In recent decades, the development of new drugs has become increasingly expensive and inefficient, and the molecular mechanisms of most pharmaceuticals remain poorly understood. In response, computational systems and network medicine tools have emerged to identify potential drug repurposing candidates. However, these tools often require complex installation and lack intuitive visual network mining capabilities. To tackle these challenges, we introduce Drugst.One, a platform that assists specialized computational medicine tools in becoming user-friendly, web-based utilities for drug repurposing. With just three lines of code, Drugst.One turns any systems biology software into an interactive web tool for modeling and analyzing complex protein-drug-disease networks. Demonstrating its broad adaptability, Drugst.One has been successfully integrated with 21 computational systems medicine tools. Available at https://drugst.one, Drugst.One has significant potential for streamlining the drug discovery process, allowing researchers to focus on essential aspects of pharmaceutical treatment research.

BIONIC: biological network integration using convolutions.

Biological networks constructed from varied data can be used to map cellular function, but each data type has limitations. Network integration promises to address these limitations by combining and automatically weighting input information to obtain a more accurate and comprehensive representation of the underlying biology. We developed a deep learning-based network integration algorithm that incorporates a graph convolutional network framework. Our method, BIONIC (Biological Network Integration using Convolutions), learns features that contain substantially more functional information compared to existing approaches. BIONIC has unsupervised and semisupervised learning modes, making use of available gene function annotations. BIONIC is scalable in both size and quantity of the input networks, making it feasible to integrate numerous networks on the scale of the human genome. To demonstrate the use of BIONIC in identifying new biology, we predicted and experimentally validated essential gene chemical-genetic interactions from nonessential gene profiles in yeast.

Single-cell, single-nucleus, and spatial transcriptomics characterization of the immunological landscape in the healthy and PSC human liver.

BACKGROUND & AIMS: Primary sclerosing cholangitis (PSC) is an immune-mediated cholestatic liver disease for which there is an unmet need to understand the cellular composition of the affected liver and how it underlies disease pathogenesis. We aimed to generate a comprehensive atlas of the PSC liver using multi-omic modalities and protein-based functional validation. METHODS: We employed single-cell and single-nucleus RNA sequencing (47,156 cells and 23,000 nuclei) and spatial transcriptomics (one sample by 10x Visium and five samples with Nanostring GeoMx DSP) to profile the cellular ecosystem in 10 PSC livers. Transcriptomic profiles were compared to 24 neurologically deceased donor livers (107,542 cells) and spatial transcriptomics controls, as well as 18,240 cells and 20,202 nuclei from three PBC livers. Flow cytometry was performed to validate PSC-specific differences in immune cell phenotype and function. RESULTS: PSC explants with parenchymal cirrhosis and prominent periductal fibrosis contained a population of cholangiocyte-like hepatocytes that were surrounded by diverse immune cell populations. PSC-associated biliary, mesenchymal, and endothelial populations expressed chemokine and cytokine transcripts involved in immune cell recruitment. Additionally, expanded CD4+ T cells and recruited myeloid populations in the PSC liver expressed the corresponding receptors to these chemokines and cytokines, suggesting potential recruitment. Tissue-resident macrophages, by contrast, were reduced in number and exhibited a dysfunctional and downregulated inflammatory response to lipopolysaccharide and interferon-γ stimulation. CONCLUSIONS: We present a comprehensive atlas of the PSC liver and demonstrate an exhaustion-like phenotype of myeloid cells and markers of chronic cytokine expression in late-stage PSC lesions. This atlas expands our understanding of the cellular complexity of PSC and has potential to guide the development of novel treatments. IMPACT AND IMPLICATIONS: Primary sclerosing cholangitis (PSC) is a rare liver disease characterized by chronic inflammation and irreparable damage to the bile ducts, which eventually results in liver failure. Due to a limited understanding of the underlying pathogenesis of disease, treatment options are limited. To address this, we sequenced healthy and diseased livers to compare the activity, interactions, and localization of immune and non-immune cells. This revealed that hepatocytes lining PSC scar regions co-express cholangiocyte markers, whereas immune cells infiltrate the scar lesions. Of these cells, macrophages, which typically contribute to tissue repair, were enriched in immunoregulatory genes and demonstrated a lack of responsiveness to stimulation. These cells may be involved in maintaining hepatic inflammation and could be a target for novel therapies.

Cytoscape.js 2023 update: a graph theory library for visualization and analysis.

SUMMARY: Cytoscape.js is an open-source JavaScript-based graph library. Its most common use case is as a visualization software component, so it can be used to render interactive graphs in a web browser. It also can be used in a headless manner, useful for graph operations on a server, such as Node.js. This update describes new features and enhancements introduced over many new versions from 2015 to 2022. AVAILABILITY AND IMPLEMENTATION: Cytoscape.js is implemented in JavaScript. Documentation, downloads and source code are available at http://js.cytoscape.org. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

PLAG1 dampens protein synthesis to promote human hematopoietic stem cell self-renewal.

Hematopoietic stem cell (HSC) dormancy is understood as supportive of HSC function and its long-term integrity. Although regulation of stress responses incurred as a result of HSC activation is recognized as important in maintaining stem cell function, little is understood of the preventive machinery present in human HSCs that may serve to resist their activation and promote HSC self-renewal. We demonstrate that the transcription factor PLAG1 is essential for long-term HSC function and, when overexpressed, endows a 15.6-fold enhancement in the frequency of functional HSCs in stimulatory conditions. Genome-wide measures of chromatin occupancy and PLAG1-directed gene expression changes combined with functional measures reveal that PLAG1 dampens protein synthesis, restrains cell growth and division, and enhances survival, with the primitive cell advantages it imparts being attenuated by addition of the potent translation activator, c-MYC. We find PLAG1 capitalizes on multiple regulatory factors to ensure protective diminished protein synthesis including 4EBP1 and translation-targeting miR-127 and does so independently of stress response signaling. Overall, our study identifies PLAG1 as an enforcer of human HSC dormancy and self-renewal through its highly context-specific regulation of protein biosynthesis and classifies PLAG1 among a rare set of bona fide regulators of messenger RNA translation in these cells. Our findings showcase the importance of regulated translation control underlying human HSC physiology, its dysregulation under activating demands, and the potential if its targeting for therapeutic benefit.

Integrated transcriptomics uncovers an enhanced association between the prion protein gene expression and vesicle dynamics signatures in glioblastomas.

BACKGROUND: Glioblastoma (GBM) is an aggressive brain tumor that exhibits resistance to current treatment, making the identification of novel therapeutic targets essential. In this context, cellular prion protein (PrPC) stands out as a potential candidate for new therapies. Encoded by the PRNP gene, PrPC can present increased expression levels in GBM, impacting cell proliferation, growth, migration, invasion and stemness. Nevertheless, the exact molecular mechanisms through which PRNP/PrPC modulates key aspects of GBM biology remain elusive. METHODS: To elucidate the implications of PRNP/PrPC in the biology of this cancer, we analyzed publicly available RNA sequencing (RNA-seq) data of patient-derived GBMs from four independent studies. First, we ranked samples profiled by bulk RNA-seq as PRNPhigh and PRNPlow and compared their transcriptomic landscape. Then, we analyzed PRNP+ and PRNP- GBM cells profiled by single-cell RNA-seq to further understand the molecular context within which PRNP/PrPC might function in this tumor. We explored an additional proteomics dataset, applying similar comparative approaches, to corroborate our findings. RESULTS: Functional profiling revealed that vesicular dynamics signatures are strongly correlated with PRNP/PrPC levels in GBM. We found a panel of 73 genes, enriched in vesicle-related pathways, whose expression levels are increased in PRNPhigh/PRNP+ cells across all RNA-seq datasets. Vesicle-associated genes, ANXA1, RAB31, DSTN and SYPL1, were found to be upregulated in vitro in an in-house collection of patient-derived GBM. Moreover, proteome analysis of patient-derived samples reinforces the findings of enhanced vesicle biogenesis, processing and trafficking in PRNPhigh/PRNP+ GBM cells. CONCLUSIONS: Together, our findings shed light on a novel role for PrPC as a potential modulator of vesicle biology in GBM, which is pivotal for intercellular communication and cancer maintenance. We also introduce GBMdiscovery, a novel user-friendly tool that allows the investigation of specific genes in GBM biology.

Multi-scale systems genomics analysis predicts pathways, cell types, and drug targets involved in normative variation in peri-adolescent human cognition.

An open challenge in human genetics is to better understand the systems-level impact of genotype variation on developmental cognition. To characterize the genetic underpinnings of peri-adolescent cognition, we performed genotype-phenotype and systems analysis for binarized accuracy in nine cognitive tasks from the Philadelphia Neurodevelopmental Cohort (~2,200 individuals of European continental ancestry aged 8-21 years). We report a region of genome-wide significance within the 3' end of the Fibulin-1 gene (P = 4.6 × 10-8), associated with accuracy in nonverbal reasoning, a heritable form of complex reasoning ability. Diffusion tensor imaging data from a subset of these participants identified a significant association of white matter fractional anisotropy with FBLN1 genotypes (P < 0.025); poor performers show an increase in the C and A allele for rs77601382 and rs5765534, respectively, which is associated with increased fractional anisotropy. Integration of published human brain-specific 'omic maps, including single-cell transcriptomes of the developing human brain, shows that FBLN1 demonstrates greatest expression in the fetal brain, as a marker of intermediate progenitor cells, demonstrates negligible expression in the adolescent and adult human brain, and demonstrates increased expression in the brain in schizophrenia. Collectively these findings warrant further study of this gene and genetic locus in cognition, neurodevelopment, and disease. Separately, genotype-pathway analysis identified an enrichment of variants associated with working memory accuracy in pathways related to development and to autonomic nervous system dysfunction. Top-ranking pathway genes include those genetically associated with diseases with working memory deficits, such as schizophrenia and Parkinson's disease. This work advances the "molecules-to-behavior" view of cognition and provides a framework for using systems-level organization of data for other biomedical domains.

The reactome pathway knowledgebase 2022.

The Reactome Knowledgebase (https://reactome.org), an Elixir core resource, provides manually curated molecular details across a broad range of physiological and pathological biological processes in humans, including both hereditary and acquired disease processes. The processes are annotated as an ordered network of molecular transformations in a single consistent data model. Reactome thus functions both as a digital archive of manually curated human biological processes and as a tool for discovering functional relationships in data such as gene expression profiles or somatic mutation catalogs from tumor cells. Recent curation work has expanded our annotations of normal and disease-associated signaling processes and of the drugs that target them, in particular infections caused by the SARS-CoV-1 and SARS-CoV-2 coronaviruses and the host response to infection. New tools support better simultaneous analysis of high-throughput data from multiple sources and the placement of understudied ('dark') proteins from analyzed datasets in the context of Reactome's manually curated pathways.

Genotyping SARS-CoV-2 Variants Using Ratiometric Nucleic Acid Barcode Panels.

Designing diagnostic assays to genotype rapidly mutating viruses remains a challenge despite the overall improvements in nucleic acid detection technologies. RT-PCR and next-generation sequencing are unsuitable for genotyping during outbreaks or in point-of-care detection due to their infrastructure requirements and longer turnaround times. We developed a quantum dot barcode multiplexing system to genotype mutated viruses. We designed multiple quantum dot barcodes to target conserved, wildtype, and mutated regions of SARS-CoV-2. We calculated ratios of the signal output from different barcodes that enabled SARS-CoV-2 detection and identified SARS-CoV-2 variant strains from a sample. We detected different sequence types, including conserved genes, nucleotide deletions, and single nucleotide substitutions. Our system detected SARS-CoV-2 patient specimens with 98% sensitivity and 94% specificity across 91 patient samples. Further, we leveraged our barcoding and ratio system to track the emergence of the N501Y SARS-CoV-2 mutation from December 2020 to May 2021 and demonstrated that the more transmissible N501Y mutation started to dominate infections by April 2021. Our barcoding and signal ratio approach can genotype viruses and track the emergence of viral mutations in a single diagnostic test. This technology can be extended to tracking other viruses. Combined with smartphone detection technologies, this assay can be adapted for point-of-care tracking of viral mutations in real time.

SBGN Bricks Ontology as a tool to describe recurring concepts in molecular networks.

A comprehensible representation of a molecular network is key to communicating and understanding scientific results in systems biology. The Systems Biology Graphical Notation (SBGN) has emerged as the main standard to represent such networks graphically. It has been implemented by different software tools, and is now largely used to communicate maps in scientific publications. However, learning the standard, and using it to build large maps, can be tedious. Moreover, SBGN maps are not grounded on a formal semantic layer and therefore do not enable formal analysis. Here, we introduce a new set of patterns representing recurring concepts encountered in molecular networks, called SBGN bricks. The bricks are structured in a new ontology, the Bricks Ontology (BKO), to define clear semantics for each of the biological concepts they represent. We show the usefulness of the bricks and BKO for both the template-based construction and the semantic annotation of molecular networks. The SBGN bricks and BKO can be freely explored and downloaded at sbgnbricks.org.

Tracing the origins of relapse in acute myeloid leukaemia to stem cells.

In acute myeloid leukaemia, long-term survival is poor as most patients relapse despite achieving remission. Historically, the failure of therapy has been thought to be due to mutations that produce drug resistance, possibly arising as a consequence of the mutagenic properties of chemotherapy drugs. However, other lines of evidence have pointed to the pre-existence of drug-resistant cells. For example, deep sequencing of paired diagnosis and relapse acute myeloid leukaemia samples has provided direct evidence that relapse in some cases is generated from minor genetic subclones present at diagnosis that survive chemotherapy, suggesting that resistant cells are generated by evolutionary processes before treatment and are selected by therapy. Nevertheless, the mechanisms of therapy failure and capacity for leukaemic regeneration remain obscure, as sequence analysis alone does not provide insight into the cell types that are fated to drive relapse. Although leukaemia stem cells have been linked to relapse owing to their dormancy and self-renewal properties, and leukaemia stem cell gene expression signatures are highly predictive of therapy failure, experimental studies have been primarily correlative and a role for leukaemia stem cells in acute myeloid leukaemia relapse has not been directly proved. Here, through combined genetic and functional analysis of purified subpopulations and xenografts from paired diagnosis/relapse samples, we identify therapy-resistant cells already present at diagnosis and two major patterns of relapse. In some cases, relapse originated from rare leukaemia stem cells with a haematopoietic stem/progenitor cell phenotype, while in other instances relapse developed from larger subclones of immunophenotypically committed leukaemia cells that retained strong stemness transcriptional signatures. The identification of distinct patterns of relapse should lead to improved methods for disease management and monitoring in acute myeloid leukaemia. Moreover, the shared functional and transcriptional stemness properties that underlie both cellular origins of relapse emphasize the importance of developing new therapeutic approaches that target stemness to prevent relapse.

Unraveling the Complexity of Liver Disease One Cell at a Time.

The human liver is a complex organ made up of multiple specialized cell types that carry out key physiological functions. An incomplete understanding of liver biology limits our ability to develop therapeutics to prevent chronic liver diseases, liver cancers, and death as a result of organ failure. Recently, single-cell modalities have expanded our understanding of the cellular phenotypic heterogeneity and intercellular cross-talk in liver health and disease. This review summarizes these findings and looks forward to highlighting new avenues for the application of single-cell genomics to unravel unknown pathogenic pathways and disease mechanisms for the development of new therapeutics targeting liver pathology. As these technologies mature, their integration into clinical data analysis will aid in patient stratification and in developing treatment plans for patients suffering from liver disease.

EPIC: software toolkit for elution profile-based inference of protein complexes.

Protein complexes are key macromolecular machines of the cell, but their description remains incomplete. We and others previously reported an experimental strategy for global characterization of native protein assemblies based on chromatographic fractionation of biological extracts coupled to precision mass spectrometry analysis (chromatographic fractionation-mass spectrometry, CF-MS), but the resulting data are challenging to process and interpret. Here, we describe EPIC (elution profile-based inference of complexes), a software toolkit for automated scoring of large-scale CF-MS data to define high-confidence multi-component macromolecules from diverse biological specimens. As a case study, we used EPIC to map the global interactome of Caenorhabditis elegans, defining 612 putative worm protein complexes linked to diverse biological processes. These included novel subunits and assemblies unique to nematodes that we validated using orthogonal methods. The open source EPIC software is freely available as a Jupyter notebook packaged in a Docker container (https://hub.docker.com/r/baderlab/bio-epic/).