Search

Platinum-based chemotherapy in combination with anti-PD-L1 antibodies has shown promising results in mesothelioma. However, the immunological mechanisms underlying its efficacy are not well understood and there are no predictive biomarkers to guide treatment decisions. Here, we combine time course RNA sequencing (RNA-seq) of peripheral blood mononuclear cells with pre-treatment tumor transcriptome data from the single-arm, phase 2 DREAM trial (N = 54). Single-cell RNA-seq and T cell receptor sequencing (TCR-seq) reveal that CD8 + T effector memory (TEM) cells with stem-like properties are more abundant in peripheral blood of responders and that this population expands upon treatment. These peripheral blood changes are linked to the transcriptional state of the tumor microenvironment. Combining information from both compartments, rather than individually, is most predictive of response. Our study highlights complex interactions between the tumor and immune cells in peripheral blood during objective tumor responses to chemoimmunotherapy. This trial is registered with the Australian New Zealand Clinical Trials Registry, number ACTRN12616001170415., Competing Interests: Declaration of interests W.J.L. declares consultancy for Douglas Pharmaceuticals and research funding from Douglas Pharmaceuticals, AstraZeneca, and ENA therapeutics. W.J.L. is a director of the Cancer Research Trust, which has funded parts of this work, but is not involved in any funding decisions, which are based on independent external peer review. W.J.L. is a founder and director of Setonix Pharmaceuticals, which is not related to the presented work. W.J.L. is an inventor on patent applications unrelated to the presented work. A.K.N. declares consultancy for Douglas Pharmaceuticals and Bristol-Myers Squibb, institutional research funding from AstraZeneca, and consultancy (unrelated to this work) from AstraZeneca., (Copyright © 2024 The Author(s). Published by Elsevier Inc. All rights reserved.)

The Human Cell Atlas (HCA) is a global partnership "to create comprehensive reference maps of all human cells-the fundamental units of life - as a basis for both understanding human health and diagnosing, monitoring, and treating disease." ( https://www.humancellatlas.org/ ) The atlas shall characterize cells from diverse individuals across the globe to better understand human biology. HCA proactively considers the priorities of, and benefits accrued to, contributing communities. Here, we lay out principles and action items that have been adopted to affirm HCA's commitment to equity so that the atlas is beneficial to all of humanity., Competing Interests: Competing interests: A.R. is an employee of Genentech, a member of the Roche group, and has equity in Roche. A.R. was a co-founder and equity holder of Celsius Therapeutics, is an equity holder in Immunitas, and until July 31, 2020 was an S.A.B. member of Thermo Fisher Scientific, Syros Pharmaceuticals, Neogene Therapeutics and Asimov. A.R. is an inventor on multiple patents related to single cell and spatial genomics. All other authors declare no competing interests., (© 2024. The Author(s).)

Breast disseminated cancer cells (DCCs) can remain dormant in the lungs for extended periods, but the mechanisms limiting their expansion are not well understood. Research indicates that tissue-resident alveolar macrophages suppress breast cancer metastasis in lung alveoli by inducing dormancy. Through ligand-receptor mapping and intravital imaging, it was found that alveolar macrophages express transforming growth factor (TGF)-β2. This expression, along with persistent macrophage-cancer cell interactions via the TGF-βRIII receptor, maintains cancer cells in a dormant state. Depleting alveolar macrophages or losing the TGF-β2 receptor in cancer cells triggers metastatic awakening. Aggressive breast cancer cells are either suppressed by alveolar macrophages or evade this suppression by avoiding interaction and downregulating the TGF-β2 receptor. Restoring TGF-βRIII in aggressive cells reinstates TGF-β2-mediated macrophage growth suppression. Thus, alveolar macrophages act as a metastasis immune barrier, and downregulation of TGF-β2 signaling allows cancer cells to overcome macrophage-mediated growth suppression., Competing Interests: Declaration of interests J.A.A.-G. is a scientific co-founder and scientific advisory board member of and an equity owner in HiberCell and receives financial compensation as a consultant for HiberCell, a Mount Sinai spin-off company focused on the research and development of therapeutics that prevent or delay the recurrence of cancer. J.A.A.-G. is also a consultant for Astrin Biosciences and Chief Mission Advisor for the Samuel Waxman Cancer Research Foundation. M.M. serves on the scientific advisory board of and holds stock from Compugen Inc., Myeloid Therapeutics Inc., Morphic Therapeutic Inc., Asher Bio Inc., Dren Bio Inc., Nirogy Inc., Oncoresponse Inc., Owkin Inc., Pionyr Inc., OSE Inc., and Larkspur Inc. M.M. serves on the scientific advisory board of Innate Pharma Inc., DBV Inc., and Genenta Inc. M.M. receives funding for contracted research from Regeneron Inc. and Boerhinger Ingelheim Inc. M.M. is a named co-inventor on an issued patent for multiplex immunohistochemistry to characterize tumors and treatment responses. The technology is filed through Icahn School of Medicine at Mount Sinai (ISMMS) and is currently unlicensed. This technology was used to evaluate tissue in this study, and the results could impact the value of this technology. M.M. has ownership interest (<5%) in Compugen Inc., Morphic Therapeutic Inc., Myeloid Therapeutics Inc., Asher Bio Inc., Dren Bio Inc., Nirogy Inc., Owkin Inc., Pionyr Inc., and Larkspur Inc. M.M. and J.A.A.-G. declare no ownership interest greater than or equal to 5%., (Copyright © 2024 Elsevier Inc. All rights reserved.)

Time-critical transcriptional events in the immune microenvironment are important for response to immune checkpoint blockade (ICB), yet these events are difficult to characterise and remain incompletely understood. Here, we present whole tumor RNA sequencing data in the context of treatment with ICB in murine models of AB1 mesothelioma and Renca renal cell cancer. We sequenced 144 bulk RNAseq samples from these two cancer types across 4 time points prior and after treatment with ICB. We also performed single-cell sequencing on 12 samples of AB1 and Renca tumors an hour before ICB administration. Our samples were equally distributed between responders and non-responders to treatment. Additionally, we sequenced AB1-HA mesothelioma tumors treated with two sample dissociation protocols to assess the impact of these protocols on the quality transcriptional information in our samples. These datasets provide time-course information to transcriptionally characterize the ICB response and provide detailed information at the single-cell level of the early tumor microenvironment prior to ICB therapy., (© 2024. The Author(s).)

Immune checkpoint therapy (ICT) causes durable tumour responses in a subgroup of patients, but it is not well known how T cell receptor beta (TCRβ) repertoire dynamics contribute to the therapeutic response. Using murine models that exclude variation in host genetics, environmental factors and tumour mutation burden, limiting variation between animals to naturally diverse TCRβ repertoires, we applied TCRseq, single cell RNAseq and flow cytometry to study TCRβ repertoire dynamics in ICT responders and non-responders. Increased oligoclonal expansion of TCRβ clonotypes was observed in responding tumours. Machine learning identified TCRβ CDR3 signatures unique to each tumour model, and signatures associated with ICT response at various timepoints before or during ICT. Clonally expanded CD8+ T cells in responding tumours post ICT displayed effector T cell gene signatures and phenotype. An early burst of clonal expansion during ICT is associated with response, and we report unique dynamics in TCRβ signatures associated with ICT response., Competing Interests: LB is employed by the company JJP Biologics. WJL received research funding from Douglas Pharmaceuticals, AstraZeneca, ENA therapeutics, consultancy for Douglas Pharmaceuticals and MSD. AN is on the advisory board of Boehringer Ingelheim, Bayer, Roche, BMS; and received research funding from AstraZeneca. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (© 2024 The Author(s). Published with license by Taylor & Francis Group, LLC.)

High-grade serous ovarian carcinoma (HGSOC) is genetically unstable and characterised by the presence of subclones with distinct genotypes. Intratumoural heterogeneity is linked to recurrence, chemotherapy resistance, and poor prognosis. Here, we use spatial transcriptomics to identify HGSOC subclones and study their association with infiltrating cell populations. Visium spatial transcriptomics reveals multiple tumour subclones with different copy number alterations present within individual tumour sections. These subclones differentially express various ligands and receptors and are predicted to differentially associate with different stromal and immune cell populations. In one sample, CosMx single molecule imaging reveals subclones differentially associating with immune cell populations, fibroblasts, and endothelial cells. Cell-to-cell communication analysis identifies subclone-specific signalling to stromal and immune cells and multiple subclone-specific autocrine loops. Our study highlights the high degree of subclonal heterogeneity in HGSOC and suggests that subclone-specific ligand and receptor expression patterns likely modulate how HGSOC cells interact with their local microenvironment., (© 2024. The Author(s).)

The incidence of hepatocellular carcinoma (HCC) is rapidly rising largely because of increased obesity leading to nonalcoholic steatohepatitis (NASH), a known HCC risk factor. There are no approved treatments to treat NASH. Here, we first used single-nucleus RNA sequencing to characterize a mouse model that mimics human NASH-driven HCC, the MUP-uPA mouse fed a high-fat diet. Activation of endoplasmic reticulum (ER) stress and inflammation was observed in a subset of hepatocytes that was enriched in mice that progress to HCC. We next treated MUP-uPA mice with the ER stress inhibitor BGP-15 and soluble gp130Fc, a drug that blocks inflammation by preventing interleukin-6 trans-signaling. Both drugs have progressed to phase 2/3 human clinical trials for other indications. We show that this combined therapy reversed NASH and reduced NASH-driven HCC. Our data suggest that these drugs could provide a potential therapy for NASH progression to HCC.

Current approaches to staging chronic liver diseases have limited utility for predicting liver cancer risk. Here, we employed single-nucleus RNA sequencing (snRNA-seq) to characterize the cellular microenvironment of healthy and pre-malignant livers using two distinct mouse models. Downstream analyses unraveled a previously uncharacterized disease-associated hepatocyte (daHep) transcriptional state. These cells were absent in healthy livers but increasingly prevalent as chronic liver disease progressed. Copy number variation (CNV) analysis of microdissected tissue demonstrated that daHep-enriched regions are riddled with structural variants, suggesting these cells represent a pre-malignant intermediary. Integrated analysis of three recent human snRNA-seq datasets confirmed the presence of a similar phenotype in human chronic liver disease and further supported its enhanced mutational burden. Importantly, we show that high daHep levels precede carcinogenesis and predict a higher risk of hepatocellular carcinoma development. These findings may change the way chronic liver disease patients are staged, surveilled, and risk stratified., Competing Interests: M.A.F. is the founder and shareholder of Celesta Therapeutics., (© 2023 The Author(s).)

Recent clinical successes have intensified interest in using adeno-associated virus (AAV) vectors for therapeutic gene delivery. The liver is a key clinical target, given its critical physiological functions and involvement in a wide range of genetic diseases. In the present study, we first investigated the validity of a liver xenograft mouse model repopulated with primary hepatocytes using single-nucleus RNA sequencing (sn-RNA-seq) by studying the transcriptomic profile of human hepatocytes pre- and post-engraftment. Complementary immunofluorescence analyses performed in highly engrafted animals confirmed that the human hepatocytes organize and present appropriate patterns of zone-dependent enzyme expression in this model. Next, we tested a set of rationally designed HSPG de-targeted AAV-LK03 variants for relative transduction performance in human hepatocytes. We used immunofluorescence, next-generation sequencing, and single-nucleus transcriptomics data from highly engrafted FRG mice to demonstrate that the optimally HSPG de-targeted AAV-LK03 displayed a significantly improved lobular transduction profile in this model., Competing Interests: L.L. is a cofounder of LogicBio Therapeutics, and L.L. and I.E.A. are co-founders of Exigen Biotherapeutics, companies that utilize similar technologies broadly discussed in this paper., (© 2023 The Authors.)

Rhabdomyolysis is the acute breakdown of skeletal myofibres in response to an initiating factor, most commonly toxins and over exertion. A variety of genetic disorders predispose to rhabdomyolysis through different pathogenic mechanisms, particularly in patients with recurrent episodes. However, most cases remain without a genetic diagnosis. Here we present six patients who presented with severe and recurrent rhabdomyolysis, usually with onset in the teenage years; other features included a history of myalgia and muscle cramps. We identified 10 bi-allelic loss-of-function variants in the gene encoding obscurin (OBSCN) predisposing individuals to recurrent rhabdomyolysis. We show reduced expression of OBSCN and loss of obscurin protein in patient muscle. Obscurin is proposed to be involved in sarcoplasmic reticulum function and Ca2+ handling. Patient cultured myoblasts appear more susceptible to starvation as evidenced by a greater decreased in sarcoplasmic reticulum Ca2+ content compared to control myoblasts. This likely reflects a lower efficiency when pumping Ca2+ back into the sarcoplasmic reticulum and/or a decrease in Ca2+ sarcoplasmic reticulum storage ability when metabolism is diminished. OSBCN variants have previously been associated with cardiomyopathies. None of the patients presented with a cardiomyopathy and cardiac examinations were normal in all cases in which cardiac function was assessed. There was also no history of cardiomyopathy in first degree relatives, in particular in any of the carrier parents. This cohort is relatively young, thus follow-up studies and the identification of additional cases with bi-allelic null OBSCN variants will further delineate OBSCN-related disease and the clinical course of disease., (© The Author(s) 2021. Published by Oxford University Press on behalf of the Guarantors of Brain. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Human brain development is underpinned by cellular and molecular reconfigurations continuing into the third decade of life. To reveal cell dynamics orchestrating neural maturation, we profiled human prefrontal cortex gene expression and chromatin accessibility at single-cell resolution from gestation to adulthood. Integrative analyses define the dynamic trajectories of each cell type, revealing major gene expression reconfiguration at the prenatal-to-postnatal transition in all cell types followed by continuous reconfiguration into adulthood and identifying regulatory networks guiding cellular developmental programs, states, and functions. We uncover links between expression dynamics and developmental milestones, characterize the diverse timing of when cells acquire adult-like states, and identify molecular convergence from distinct developmental origins. We further reveal cellular dynamics and their regulators implicated in neurological disorders. Finally, using this reference, we benchmark cell identities and maturation states in organoid models. Together, this captures the dynamic regulatory landscape of human cortical development., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2022 The Authors. Published by Elsevier Inc. All rights reserved.)

In eukaryotes, capped RNAs include long transcripts such as messenger RNAs and long noncoding RNAs, as well as shorter transcripts such as spliceosomal RNAs, small nucleolar RNAs, and enhancer RNAs. Long capped transcripts can be profiled using cap analysis gene expression (CAGE) sequencing and other methods. Here, we describe a sequencing library preparation protocol for short capped RNAs, apply it to a differentiation time course of the human cell line THP-1, and systematically compare the landscape of short capped RNAs to that of long capped RNAs. Transcription initiation peaks associated with genes in the sense direction have a strong preference to produce either long or short capped RNAs, with one out of six peaks detected in the short capped RNA libraries only. Gene-associated short capped RNAs have highly specific 3' ends, typically overlapping splice sites. Enhancers also preferentially generate either short or long capped RNAs, with 10% of enhancers observed in the short capped RNA libraries only. Enhancers producing either short or long capped RNAs show enrichment for GWAS-associated disease SNPs. We conclude that deep sequencing of short capped RNAs reveals new families of noncoding RNAs and elucidates the diversity of transcripts generated at known and novel promoters and enhancers., (© 2022 de Hoon et al.; Published by Cold Spring Harbor Laboratory Press.)

The biological determinants of the response to immune checkpoint blockade (ICB) in cancer remain incompletely understood. Little is known about dynamic biological events that underpin therapeutic efficacy due to the inability to frequently sample tumours in patients. Here, we map the transcriptional profiles of 144 responding and non-responding tumours within two mouse models at four time points during ICB. We find that responding tumours display on/fast-off kinetics of type-I-interferon (IFN) signaling. Phenocopying of this kinetics using time-dependent sequential dosing of recombinant IFNs and neutralizing antibodies markedly improves ICB efficacy, but only when IFNβ is targeted, not IFNα. We identify Ly6C + /CD11b + inflammatory monocytes as the primary source of IFNβ and find that active type-I-IFN signaling in tumour-infiltrating inflammatory monocytes is associated with T cell expansion in patients treated with ICB. Together, our results suggest that on/fast-off modulation of IFNβ signaling is critical to the therapeutic response to ICB, which can be exploited to drive clinical outcomes towards response., (© 2022. The Author(s).)

Hi-C is a genome-wide chromosome conformation capture technology that detects interactions between pairs of genomic regions and exploits higher order chromatin structures. Conceptually Hi-C data counts interaction frequencies between every position in the genome and every other position. Biologically functional interactions are expected to occur more frequently than transient background and artefactual interactions. To identify biologically relevant interactions, several background models that take biases such as distance, GC content and mappability into account have been proposed. Here we introduce MaxHiC, a background correction tool that deals with these complex biases and robustly identifies statistically significant interactions in both Hi-C and capture Hi-C experiments. MaxHiC uses a negative binomial distribution model and a maximum likelihood technique to correct biases in both Hi-C and capture Hi-C libraries. We systematically benchmark MaxHiC against major Hi-C background correction tools including Hi-C significant interaction callers (SIC) and Hi-C loop callers using published Hi-C, capture Hi-C, and Micro-C datasets. Our results demonstrate that 1) Interacting regions identified by MaxHiC have significantly greater levels of overlap with known regulatory features (e.g. active chromatin histone marks, CTCF binding sites, DNase sensitivity) and also disease-associated genome-wide association SNPs than those identified by currently existing models, 2) the pairs of interacting regions are more likely to be linked by eQTL pairs and 3) more likely to link known regulatory features including known functional enhancer-promoter pairs validated by CRISPRi than any of the existing methods. We also demonstrate that interactions between different genomic region types have distinct distance distributions only revealed by MaxHiC. MaxHiC is publicly available as a python package for the analysis of Hi-C, capture Hi-C and Micro-C data., Competing Interests: The authors have declared that no competing interests exist.

Tumor heterogeneity is a major obstacle to the success of cancer treatment. An accurate understanding and recognition of tumor heterogeneity is critical in the clinical management of cancer patients. Here, we utilized single-cell RNA sequencing (scRNA-seq) to uncover the intra- and intertumoral heterogeneity of liver metastases from a patient with metastatic uveal melanoma. The two metastases analyzed were largely infiltrated by noncancerous cells with significant variability in the proportion of different cell types. Analysis of copy-number variations (CNVs) showed gain of 8q and loss of 6q in both tumors, but loss of Chromosome 3 was only detected in one of the tumors. Single-nucleotide polymorphism (SNP) array revealed a uniparental isodisomy 3 in the tumor with two copies of Chromosome 3, indicating a regain of Chromosome 3 during the development of the metastatic disease. In addition, both tumors harbored subclones with additional CNVs. Pathway enrichment analysis of differentially expressed genes revealed that cancer cells in the metastasis with isodisomy 3 showed up-regulation in epithelial-mesenchymal transition and myogenesis related genes. In contrast, up-regulation in interferon signaling was observed in the metastasis with monosomy 3 and increased T-cell infiltrate. This study highlights the complexity and heterogeneity of different metastases within an individual case of uveal melanoma., (© 2021 Lin et al.; Published by Cold Spring Harbor Laboratory Press.)

Background: CRISPR/Cas9 is an invaluable tool for studying cell biology and the development of molecular therapies. However, delivery of CRISPR/Cas9 components into some cell types remains a major hurdle. Primary human myoblasts are a valuable cell model for muscle studies, but are notoriously difficult to transfect. There are currently no commercial lipofection protocols tailored for primary myoblasts, and most generic guidelines simply recommend transfecting healthy cells at high confluency. This study aimed to maximize CRISPR/Cas9 transfection and editing in primary human myoblasts., Methods: Since increased cell proliferation is associated with increased transfection efficiency, we investigated two factors known to influence myoblast proliferation: cell confluency, and a basement membrane matrix, Matrigel. CRISPR/Cas9 editing was performed by delivering Cas9 ribonucleoprotein complexes via lipofection into primary human myoblasts, cultured in wells with or without a Matrigel coating, at low (~ 40%) or high (~ 80%) confluency., Results: Cells transfected at low confluency on Matrigel-coated wells had the highest levels of transfection, and were most effectively edited across three different target loci, achieving a maximum editing efficiency of 93.8%. On average, editing under these conditions was >4-fold higher compared to commercial recommendations (high confluency, uncoated wells)., Conclusion: This study presents a simple, effective and economical method of maximizing CRISPR/Cas9-mediated gene editing in primary human myoblasts. This protocol could be a valuable tool for improving the genetic manipulation of cultured human skeletal muscle cells, and potentially be adapted for use in other cell types., (© 2021. The Author(s).)

Autism spectrum disorder (ASD) is a neurodevelopmental condition with substantial phenotypic and etiological heterogeneity. Although 10%-20% of ASD cases are attributable to copy number variation (CNV), causative genomic loci and constituent genes remain unclarified. We have developed SNATCNV, a tool that outperforms existing tools, to identify 47 recurrent ASD CNV regions from 19,663 cases and 6,479 controls documented in the AutDB database. Analysis of ASD CNV gene content using FANTOM5 shows that constituent coding genes and long non-coding RNAs have brain-enriched patterns of expression. Notably, such enrichment is not observed for regions identified by using other tools. We also find evidence of sexual dimorphism, one locus uniquely comprising a single lncRNA gene, and correlation of CNVs to distinct clinical and behavioral traits. Finally, we analyze a large dataset for schizophrenia to further demonstrate that SNATCNV is an effective, publicly available tool to define genomic loci and causative genes for multiple CNV-associated conditions., Competing Interests: Declaration of Interests The authors declare no competing interests., (Copyright © 2020 The Author(s). Published by Elsevier Inc. All rights reserved.)

Development of high throughput single-cell sequencing technologies has made it cost-effective to profile thousands of cells from diverse samples containing multiple cell types. To study how these different cell types work together, here we develop NATMI (Network Analysis Toolkit for Multicellular Interactions). NATMI uses connectomeDB2020 (a database of 2293 manually curated ligand-receptor pairs with literature support) to predict and visualise cell-to-cell communication networks from single-cell (or bulk) expression data. Using multiple published single-cell datasets we demonstrate how NATMI can be used to identify (i) the cell-type pairs that are communicating the most (or most specifically) within a network, (ii) the most active (or specific) ligand-receptor pairs active within a network, (iii) putative highly-communicating cellular communities and (iv) differences in intercellular communication when profiling given cell types under different conditions. Furthermore, analysis of the Tabula Muris (organism-wide) atlas confirms our previous prediction that autocrine signalling is a major feature of cell-to-cell communication networks, while also revealing that hundreds of ligands and their cognate receptors are co-expressed in individual cells suggesting a substantial potential for self-signalling.

The recent coronavirus disease 2019 (COVID-19) pandemic has caused worldwide disruption which also extends to the arena of scientific meetings around the world. Here, we explore the lessons learned from moving two human genetics and genomics meetings quickly to an online format in early 2020. The tips presented herein may be useful not only for future virtual meetings but may also enrich future physical if not hybrid meetings once they resume.

Long noncoding RNAs (lncRNAs) constitute the majority of transcripts in the mammalian genomes, and yet, their functions remain largely unknown. As part of the FANTOM6 project, we systematically knocked down the expression of 285 lncRNAs in human dermal fibroblasts and quantified cellular growth, morphological changes, and transcriptomic responses using Capped Analysis of Gene Expression (CAGE). Antisense oligonucleotides targeting the same lncRNAs exhibited global concordance, and the molecular phenotype, measured by CAGE, recapitulated the observed cellular phenotypes while providing additional insights on the affected genes and pathways. Here, we disseminate the largest-to-date lncRNA knockdown data set with molecular phenotyping (over 1000 CAGE deep-sequencing libraries) for further exploration and highlight functional roles for ZNF213-AS1 and lnc-KHDC3L-2 ., (© 2020 Ramilowski et al.; Published by Cold Spring Harbor Laboratory Press.)

Background: Single-cell RNA sequencing has been widely adopted to estimate the cellular composition of heterogeneous tissues and obtain transcriptional profiles of individual cells. Multiple approaches for optimal sample dissociation and storage of single cells have been proposed as have single-nuclei profiling methods. What has been lacking is a systematic comparison of their relative biases and benefits., Results: Here, we compare gene expression and cellular composition of single-cell suspensions prepared from adult mouse kidney using two tissue dissociation protocols. For each sample, we also compare fresh cells to cryopreserved and methanol-fixed cells. Lastly, we compare this single-cell data to that generated using three single-nucleus RNA sequencing workflows. Our data confirms prior reports that digestion on ice avoids the stress response observed with 37 °C dissociation. It also reveals cell types more abundant either in the cold or warm dissociations that may represent populations that require gentler or harsher conditions to be released intact. For cell storage, cryopreservation of dissociated cells results in a major loss of epithelial cell types; in contrast, methanol fixation maintains the cellular composition but suffers from ambient RNA leakage. Finally, cell type composition differences are observed between single-cell and single-nucleus RNA sequencing libraries. In particular, we note an underrepresentation of T, B, and NK lymphocytes in the single-nucleus libraries., Conclusions: Systematic comparison of recovered cell types and their transcriptional profiles across the workflows has highlighted protocol-specific biases and thus enables researchers starting single-cell experiments to make an informed choice.

Macrophages in the lung detect and respond to influenza A virus (IAV), determining the nature of the immune response. Using terminal-depth cap analysis of gene expression (CAGE), we quantified transcriptional activity of both host and pathogen over a 24-h time course of IAV infection in primary human monocyte-derived macrophages (MDMs). This method allowed us to observe heterogenous host sequences incorporated into IAV mRNA, "snatched" 5' RNA caps, and corresponding RNA sequences from host RNAs. In order to determine whether cap-snatching is random or exhibits a bias, we systematically compared host sequences incorporated into viral mRNA ("snatched") against a complete survey of all background host RNA in the same cells, at the same time. Using a computational strategy designed to eliminate sources of bias due to read length, sequencing depth, and multimapping, we were able to quantify overrepresentation of host RNA features among the sequences that were snatched by IAV. We demonstrate biased snatching of numerous host RNAs, particularly small nuclear RNAs (snRNAs), and avoidance of host transcripts encoding host ribosomal proteins, which are required by IAV for replication. We then used a systems approach to describe the transcriptional landscape of the host response to IAV, observing many new features, including a failure of IAV-treated MDMs to induce feedback inhibitors of inflammation, seen in response to other treatments. IMPORTANCE Infection with influenza A virus (IAV) infection is responsible for an estimated 500,000 deaths and up to 5 million cases of severe respiratory illness each year. In this study, we looked at human primary immune cells (macrophages) infected with IAV. Our method allows us to look at both the host and the virus in parallel. We used these data to explore a process known as "cap-snatching," where IAV snatches a short nucleotide sequence from capped host RNA. This process was believed to be random. We demonstrate biased snatching of numerous host RNAs, including those associated with snRNA transcription, and avoidance of host transcripts encoding host ribosomal proteins, which are required by IAV for replication. We then describe the transcriptional landscape of the host response to IAV, observing new features, including a failure of IAV-treated MDMs to induce feedback inhibitors of inflammation, seen in response to other treatments., (Copyright © 2020 Clohisey et al.)

Analysis of cancer mutational signatures have been instrumental in identification of responsible endogenous and exogenous molecular processes in cancer. The quantitative approach used to deconvolute mutational signatures is becoming an integral part of cancer research. Therefore, development of a stand-alone tool with a user-friendly interface for analysis of cancer mutational signatures is necessary. In this manuscript we introduce CANCERSIGN, which enables users to identify 3-mer and 5-mer mutational signatures within whole genome, whole exome or pooled samples. Additionally, this tool enables users to perform clustering on tumor samples based on the proportion of mutational signatures in each sample. Using CANCERSIGN, we analysed all the whole genome somatic mutation datasets profiled by the International Cancer Genome Consortium (ICGC) and identified a number of novel signatures. By examining signatures found in exonic and non-exonic regions of the genome using WGS and comparing this to signatures found in WES data we observe that WGS can identify additional non-exonic signatures that are enriched in the non-coding regions of the genome while the deeper sequencing of WES may help identify weak signatures that are otherwise missed in shallower WGS data.

Cancer precision medicine aims to predict the drug likely to yield the best response for a patient. Genomic sequencing of tumors is currently being used to better inform treatment options; however, this approach has had a limited clinical impact due to the paucity of actionable mutations. An alternative to mutation status is the use of gene expression signatures to predict response. Using data from two large-scale studies, The Genomics of Drug Sensitivity of Cancer (GDSC) and The Cancer Therapeutics Response Portal (CTRP), we investigated the relationship between the sensitivity of hundreds of cell lines to hundreds of drugs, and the relative expression levels of the targets these drugs are directed against. For approximately one third of the drugs considered (73/222 in GDSC and 131/360 in CTRP), sensitivity was significantly correlated with the expression of at least one of the known targets. Surprisingly, for 8% of the annotated targets, there was a significant anticorrelation between target expression and sensitivity. For several cases, this corresponded to drugs targeting multiple genes in the same family, with the expression of one target significantly correlated with sensitivity and another significantly anticorrelated suggesting a possible role in resistance. Furthermore, we identified nontarget genes that are significantly correlated or anticorrelated with drug sensitivity, and find literature linking several to sensitization and resistance. Our analyses provide novel and important insights into both potential mechanisms of resistance and relative efficacy of drugs against the same target., (©2019 American Association for Cancer Research.)

Motivation: Single-cell RNA sequencing (scRNA-seq) measures gene expression at the resolution of individual cells. Massively multiplexed single-cell profiling has enabled large-scale transcriptional analyses of thousands of cells in complex tissues. In most cases, the true identity of individual cells is unknown and needs to be inferred from the transcriptomic data. Existing methods typically cluster (group) cells based on similarities of their gene expression profiles and assign the same identity to all cells within each cluster using the averaged expression levels. However, scRNA-seq experiments typically produce low-coverage sequencing data for each cell, which hinders the clustering process., Results: We introduce scMatch, which directly annotates single cells by identifying their closest match in large reference datasets. We used this strategy to annotate various single-cell datasets and evaluated the impacts of sequencing depth, similarity metric and reference datasets. We found that scMatch can rapidly and robustly annotate single cells with comparable accuracy to another recent cell annotation tool (SingleR), but that it is quicker and can handle larger reference datasets. We demonstrate how scMatch can handle large customized reference gene expression profiles that combine data from multiple sources, thus empowering researchers to identify cell populations in any complex tissue with the desired precision., Availability and Implementation: scMatch (Python code) and the FANTOM5 reference dataset are freely available to the research community here https://github.com/forrest-lab/scMatch., Supplementary Information: Supplementary data are available at Bioinformatics online., (© The Author(s) 2019. Published by Oxford University Press.)

Background: The work of the FANTOM5 Consortium has brought forth a new level of understanding of the regulation of gene transcription and the cellular processes involved in creating diversity of cell types. In this study, we extended the analysis of the FANTOM5 Cap Analysis of Gene Expression (CAGE) transcriptome data to focus on understanding the genetic regulators involved in mouse cerebellar development., Results: We used the HeliScopeCAGE library sequencing on cerebellar samples over 8 embryonic and 4 early postnatal times. This study showcases temporal expression pattern changes during cerebellar development. Through a bioinformatics analysis that focused on transcription factors, their promoters and binding sites, we identified genes that appear as strong candidates for involvement in cerebellar development. We selected several candidate transcriptional regulators for validation experiments including qRT-PCR and shRNA transcript knockdown. We observed marked and reproducible developmental defects in Atf4, Rfx3, and Scrt2 knockdown embryos, which support the role of these genes in cerebellar development., Conclusions: The successful identification of these novel gene regulators in cerebellar development demonstrates that the FANTOM5 cerebellum time series is a high-quality transcriptome database for functional investigation of gene regulatory networks in cerebellar development.

Natural antisense transcripts are common features of mammalian genes providing additional regulatory layers of gene expression. A comprehensive description of antisense transcription in loci associated to familial neurodegenerative diseases may identify key players in gene regulation and provide tools for manipulating gene expression. We take advantage of the FANTOM5 sequencing datasets that represent the largest collection to date of genome-wide promoter usage in almost 2000 human samples. Transcription start sites (TSSs) are mapped at high resolution by the use of a modified protocol of cap analysis of gene expression (CAGE) for high-throughput single molecule next-generation sequencing with Helicos (hCAGE). Here we present the analysis of antisense transcription at 17 loci associated to hereditary Alzheimer's disease, Frontotemporal Dementia, Parkinson's disease, Amyotrophic Lateral Sclerosis, and Huntington's disease. We focused our analysis on libraries derived from brain tissues and primary cells. We also screened libraries from total blood and blood cell populations in the quest for peripheral biomarkers of neurodegenerative diseases. We identified 63 robust promoters in antisense orientation to genes associated to familial neurodegeneration. When applying a less stringent cutoff, this number increases to over 400. A subset of these promoters represents alternative TSSs for 24 FANTOM5 annotated long noncoding RNA (lncRNA) genes, in antisense orientation to 13 of the loci analyzed here, while the remaining contribute to the expression of additional transcript variants. Intersection with GWAS studies, sample ontology, and dynamic expression reveals association to specific genetic traits as well as cell and tissue types, not limited to neurodegenerative diseases. Antisense transcription was validated for a subset of genes, including those encoding for Microtubule-Associated Protein Tau, α-synuclein, Parkinsonism-associated deglycase DJ-1, and Leucin-Rich Repeat Kinase 2. This work provides evidence for the existence of additional regulatory mechanisms of the expression of neurodegenerative disease-causing genes by previously not-annotated and/or not-validated antisense long noncoding RNAs.

Cancer immunotherapy using antibodies that target immune checkpoints has delivered outstanding results. However, responses only occur in a subset of patients, and it is not fully understood what biological processes determine an effective outcome. This lack of understanding hinders the development of rational combination treatments. We set out to define the pretreatment microenvironment associated with an effective outcome by using the fact that inbred mouse strains bearing monoclonal cancer cell line-derived tumors respond in a dichotomous manner to immune checkpoint blockade (ICB). We compared the cellular composition and gene expression profiles of responsive and nonresponsive tumors from mice before ICB and validated the findings in cohorts of patients with cancer treated with ICB antibodies. We found that responsive tumors were characterized by an inflammatory gene expression signature consistent with up-regulation of signal transducer and activator of transcription 1 (STAT1) and Toll-like receptor 3 (TLR3) signaling and down-regulation of interleukin-10 (IL-10) signaling. In addition, responsive tumors had more infiltrating-activated natural killer (NK) cells, which were necessary for response. Pretreatment of mice with large established tumors using the STAT1-activating cytokine interferon-γ (IFNγ), the TLR3 ligand poly(I:C), and an anti-IL-10 antibody sensitized tumors to ICB by attracting IFNγ-producing NK cells into the tumor, resulting in increased cure rates. Our results identify a pretreatment tumor microenvironment that predicts response to ICB, which can be therapeutically attained. These data suggest a biomarker-driven approach to patient management to establish whether a patient would benefit from treatment with sensitizing therapeutics before ICB., (Copyright © 2019 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works.)

Myoglobin, encoded by MB, is a small cytoplasmic globular hemoprotein highly expressed in cardiac myocytes and oxidative skeletal myofibers. Myoglobin binds O 2, facilitates its intracellular transport and serves as a controller of nitric oxide and reactive oxygen species. Here, we identify a recurrent c.292C>T (p.His98Tyr) substitution in MB in fourteen members of six European families suffering from an autosomal dominant progressive myopathy with highly characteristic sarcoplasmic inclusions in skeletal and cardiac muscle. Myoglobinopathy manifests in adulthood with proximal and axial weakness that progresses to involve distal muscles and causes respiratory and cardiac failure. Biochemical characterization reveals that the mutant myoglobin has altered O 2 binding, exhibits a faster heme dissociation rate and has a lower reduction potential compared to wild-type myoglobin. Preliminary studies show that mutant myoglobin may result in elevated superoxide levels at the cellular level. These data define a recognizable muscle disease associated with MB mutation.

Congenital myopathies are typically characterised by early onset hypotonia, weakness and hallmark features on biopsy. Despite the rapid pace of gene discovery, ∼50% of patients with a congenital myopathy remain without a genetic diagnosis following screening of known disease genes. We performed exome sequencing on two consanguineous probands diagnosed with a congenital myopathy and muscle biopsy showing selective atrophy/hypotrophy or absence of type II myofibres. We identified variants in the gene (MYL1) encoding the skeletal muscle fast-twitch specific myosin essential light chain (ELC) in both probands. A homozygous essential splice acceptor variant (c.479-2A > G, predicted to result in skipping of exon 5 was identified in Proband 1, and a homozygous missense substitution (c.488T>G, p.(Met163Arg)) was identified in Proband 2. Protein modelling of the p.(Met163Arg) substitution predicted it might impede intermolecular interactions that facilitate binding to the IQ domain of myosin heavy chain, thus likely impacting on the structure and functioning of the myosin motor. MYL1 was markedly reduced in skeletal muscle from both probands, suggesting that the missense substitution likely results in an unstable protein. Knock down of myl1 in zebrafish resulted in abnormal morphology, disrupted muscle structure and impaired touch-evoked escape responses, thus confirming that skeletal muscle fast-twitch specific myosin ELC is critical for myofibre development and function. Our data implicate MYL1 as a crucial protein for adequate skeletal muscle function and that MYL1 deficiency is associated with severe congenital myopathy.

The authors regret that Luba M. Pardo was omitted in error from the author list of the original version of this Data Descriptor. This omission has now been corrected in the HTML and PDF versions. The authors also regret that Anemieke Rozemuller was omitted in error from the Acknowledgements of the original version of this Data Descriptor. This omission has now been corrected in the HTML and PDF versions.

The promoters of immediate early genes (IEGs) are rapidly activated in response to an external stimulus. These genes, also known as primary response genes, have been identified in a range of cell types, under diverse extracellular signals and using varying experimental protocols. Whereas genomic dissection on a case-by-case basis has not resulted in a comprehensive catalogue of IEGs, a rigorous meta-analysis of eight genome-wide FANTOM5 CAGE (cap analysis of gene expression) time course datasets reveals successive waves of promoter activation in IEGs, recapitulating known relationships between cell types and stimuli: we obtain a set of 57 (42 protein-coding) candidate IEGs possessing promoters that consistently drive a rapid but transient increase in expression over time. These genes show significant enrichment for known IEGs reported previously, pathways associated with the immediate early response, and include a number of non-coding RNAs with roles in proliferation and differentiation. Surprisingly, we also find strong conservation of the ordering of activation for these genes, such that 77 pairwise promoter activation orderings are conserved. Using the leverage of comprehensive CAGE time series data across cell types, we also document the extensive alternative promoter usage by such genes, which is likely to have been a barrier to their discovery until now. The common activation ordering of the core set of early-responding genes we identify may indicate conserved underlying regulatory mechanisms. By contrast, the considerably larger number of transiently activated genes that are specific to each cell type and stimulus illustrates the breadth of the primary response., (© 2018 The Authors.)

Laser-capture microdissection was used to isolate external germinal layer tissue from three developmental periods of mouse cerebellar development: embryonic days 13, 15, and 18. The cerebellar granule cell-enriched mRNA library was generated with next-generation sequencing using the Helicos technology. Our objective was to discover transcriptional regulators that could be important for the development of cerebellar granule cells-the most numerous neuron in the central nervous system. Through differential expression analysis, we have identified 82 differentially expressed transcription factors (TFs) from a total of 1311 differentially expressed genes. In addition, with TF-binding sequence analysis, we have identified 46 TF candidates that could be key regulators responsible for the variation in the granule cell transcriptome between developmental stages. Altogether, we identified 125 potential TFs (82 from differential expression analysis, 46 from motif analysis with 3 overlaps in the two sets). From this gene set, 37 TFs are considered novel due to the lack of previous knowledge about their roles in cerebellar development. The results from transcriptome-wide analyses were validated with existing online databases, qRT-PCR, and in situ hybridization. This study provides an initial insight into the TFs of cerebellar granule cells that might be important for development and provide valuable information for further functional studies on these transcriptional regulators.

Mycobacterium tuberculosis (Mtb) infection reveals complex and dynamic host-pathogen interactions, leading to host protection or pathogenesis. Using a unique transcriptome technology (CAGE), we investigated the promoter-based transcriptional landscape of IFNγ (M1) or IL-4/IL-13 (M2) stimulated macrophages during Mtb infection in a time-kinetic manner. Mtb infection widely and drastically altered macrophage-specific gene expression, which is far larger than that of M1 or M2 activations. Gene Ontology enrichment analysis for Mtb-induced differentially expressed genes revealed various terms, related to host-protection and inflammation, enriched in up-regulated genes. On the other hand, terms related to dis-regulation of cellular functions were enriched in down-regulated genes. Differential expression analysis revealed known as well as novel transcription factor genes in Mtb infection, many of them significantly down-regulated. IFNγ or IL-4/IL-13 pre-stimulation induce additional differentially expressed genes in Mtb-infected macrophages. Cluster analysis uncovered significant numbers, prolonging their expressional changes. Furthermore, Mtb infection augmented cytokine-mediated M1 and M2 pre-activations. In addition, we identified unique transcriptional features of Mtb-mediated differentially expressed lncRNAs. In summary we provide a comprehensive in depth gene expression/regulation profile in Mtb-infected macrophages, an important step forward for a better understanding of host-pathogen interaction dynamics in Mtb infection.

Genetic variants underlying complex traits, including disease susceptibility, are enriched within the transcriptional regulatory elements, promoters and enhancers. There is emerging evidence that regulatory elements associated with particular traits or diseases share similar patterns of transcriptional activity. Accordingly, shared transcriptional activity (coexpression) may help prioritise loci associated with a given trait, and help to identify underlying biological processes. Using cap analysis of gene expression (CAGE) profiles of promoter- and enhancer-derived RNAs across 1824 human samples, we have analysed coexpression of RNAs originating from trait-associated regulatory regions using a novel quantitative method (network density analysis; NDA). For most traits studied, phenotype-associated variants in regulatory regions were linked to tightly-coexpressed networks that are likely to share important functional characteristics. Coexpression provides a new signal, independent of phenotype association, to enable fine mapping of causative variants. The NDA coexpression approach identifies new genetic variants associated with specific traits, including an association between the regulation of the OCT1 cation transporter and genetic variants underlying circulating cholesterol levels. NDA strongly implicates particular cell types and tissues in disease pathogenesis. For example, distinct groupings of disease-associated regulatory regions implicate two distinct biological processes in the pathogenesis of ulcerative colitis; a further two separate processes are implicated in Crohn's disease. Thus, our functional analysis of genetic predisposition to disease defines new distinct disease endotypes. We predict that patients with a preponderance of susceptibility variants in each group are likely to respond differently to pharmacological therapy. Together, these findings enable a deeper biological understanding of the causal basis of complex traits.