Time-course RNAseq data of murine AB1 mesothelioma and Renca renal cancer following immune checkpoint therapy.

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.

Spatial transcriptomics reveals discrete tumour microenvironments and autocrine loops within ovarian cancer subclones.

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.

Immune checkpoint therapy responders display early clonal expansion of tumor infiltrating lymphocytes.

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.

Therapeutic blockade of ER stress and inflammation prevents NASH and progression to HCC.

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.

Single-nucleus RNA sequencing of pre-malignant liver reveals disease-associated hepatocyte state with HCC prognostic potential.

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.

Characterization of the humanized FRG mouse model and development of an AAV-LK03 variant with improved liver lobular biodistribution.

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.

Human prefrontal cortex gene regulatory dynamics from gestation to adulthood at single-cell resolution.

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.

Correction: MaxHiC: A robust background correction model to identify biologically relevant chromatin interactions in Hi-C and capture Hi-C experiments.

[This corrects the article DOI: 10.1371/journal.pcbi.1010241.].

Temporally restricted activation of IFNß signaling underlies response to immune checkpoint therapy in mice.

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.

Deep sequencing of short capped RNAs reveals novel families of noncoding RNAs.

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.

MaxHiC: A robust background correction model to identify biologically relevant chromatin interactions in Hi-C and capture Hi-C experiments.

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.

Functional validation of variants of unknown significance using CRISPR gene editing and transcriptomics: A Kleefstra syndrome case study.

There are an estimated > 400 million people living with a rare disease globally, with genetic variants the cause of approximately 80% of cases. Next Generation Sequencing (NGS) rapidly identifies genetic variants however they are often of unknown significance. Low throughput functional validation in specialist laboratories is the current ad hoc approach for functional validation of genetic variants, which creating major bottlenecks in patient diagnosis. This study investigates the application of CRISPR gene editing followed by genome wide transcriptomic profiling to facilitate patient diagnosis. As proof-of-concept, we introduced a variant in the Euchromatin histone methyl transferase (EHMT1) gene into HEK293T cells. We identified changes in the regulation of the cell cycle, neural gene expression and suppression of gene expression changes on chromosome 19 and chromosome X, that are in keeping with Kleefstra syndrome clinical phenotype and/or provide insight into disease mechanism. This study demonstrates the utility of genome editing followed by functional readouts to rapidly and systematically validating the function of variants of unknown significance in patients suffering from rare diseases.

Bi-allelic loss-of-function OBSCN variants predispose individuals to severe recurrent rhabdomyolysis.

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.

Improved CRISPR/Cas9 gene editing in primary human myoblasts using low confluency cultures on Matrigel.

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.

Intra- and intertumoral heterogeneity of liver metastases in a patient with uveal melanoma revealed by single-cell RNA sequencing.

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.

Molecular Insights Into Regulatory T-Cell Adaptation to Self, Environment, and Host Tissues: Plasticity or Loss of Function in Autoimmune Disease.

There has been much interest in the ability of regulatory T cells (Treg) to switch function in vivo, either as a result of genetic risk of disease or in response to environmental and metabolic cues. The relationship between levels of FOXP3 and functional fitness plays a significant part in this plasticity. There is an emerging role for Treg in tissue repair that may be less dependent on FOXP3, and the molecular mechanisms underpinning this are not fully understood. As a result of detailed, high-resolution functional genomics, the gene regulatory networks and key functional mediators of Treg phenotype downstream of FOXP3 have been mapped, enabling a mechanistic insight into Treg function. This transcription factor-driven programming of T-cell function to generate Treg requires the switching on and off of key genes that form part of the Treg gene regulatory network and raises the possibility that this is reversible. It is plausible that subtle shifts in expression levels of specific genes, including transcription factors and non-coding RNAs, change the regulation of the Treg gene network. The subtle skewing of gene expression initiates changes in function, with the potential to promote chronic disease and/or to license appropriate inflammatory responses. In the case of autoimmunity, there is an underlying genetic risk, and the interplay of genetic and environmental cues is complex and impacts gene regulation networks frequently involving promoters and enhancers, the regulatory elements that control gene expression levels and responsiveness. These promoter-enhancer interactions can operate over long distances and are highly cell type specific. In autoimmunity, the genetic risk can result in changes in these enhancer/promoter interactions, and this mainly impacts genes which are expressed in T cells and hence impacts Treg/conventional T-cell (Tconv) function. Genetic risk may cause the subtle alterations to the responsiveness of gene regulatory networks which are controlled by or control FOXP3 and its target genes, and the application of assays of the 3D organization of chromatin, enabling the connection of non-coding regulatory regions to the genes they control, is revealing the direct impact of environmental/metabolic/genetic risk on T-cell function and is providing mechanistic insight into susceptibility to inflammatory and autoimmune conditions.

Predicting cell-to-cell communication networks using NATMI.

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.

Brain-Enriched Coding and Long Non-coding RNA Genes Are Overrepresented in Recurrent Neurodevelopmental Disorder CNVs.

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.

Stromal cell diversity associated with immune evasion in human triple-negative breast cancer.

The tumour stroma regulates nearly all stages of carcinogenesis. Stromal heterogeneity in human triple-negative breast cancers (TNBCs) remains poorly understood, limiting the development of stromal-targeted therapies. Single-cell RNA sequencing of five TNBCs revealed two cancer-associated fibroblast (CAF) and two perivascular-like (PVL) subpopulations. CAFs clustered into two states: the first with features of myofibroblasts and the second characterised by high expression of growth factors and immunomodulatory molecules. PVL cells clustered into two states consistent with a differentiated and immature phenotype. We showed that these stromal states have distinct morphologies, spatial relationships and functional properties in regulating the extracellular matrix. Using cell signalling predictions, we provide evidence that stromal-immune crosstalk acts via a diverse array of immunoregulatory molecules. Importantly, the investigation of gene signatures from inflammatory-CAFs and differentiated-PVL cells in independent TNBC patient cohorts revealed strong associations with cytotoxic T-cell dysfunction and exclusion, respectively. Such insights present promising candidates to further investigate for new therapeutic strategies in the treatment of TNBCs.

Comparative transcriptomics of primary cells in vertebrates.

Gene expression profiles in homologous tissues have been observed to be different between species, which may be due to differences between species in the gene expression program in each cell type, but may also reflect differences in cell type composition of each tissue in different species. Here, we compare expression profiles in matching primary cells in human, mouse, rat, dog, and chicken using Cap Analysis Gene Expression (CAGE) and short RNA (sRNA) sequencing data from FANTOM5. While we find that expression profiles of orthologous genes in different species are highly correlated across cell types, in each cell type many genes were differentially expressed between species. Expression of genes with products involved in transcription, RNA processing, and transcriptional regulation was more likely to be conserved, while expression of genes encoding proteins involved in intercellular communication was more likely to have diverged during evolution. Conservation of expression correlated positively with the evolutionary age of genes, suggesting that divergence in expression levels of genes critical for cell function was restricted during evolution. Motif activity analysis showed that both promoters and enhancers are activated by the same transcription factors in different species. An analysis of expression levels of mature miRNAs and of primary miRNAs identified by CAGE revealed that evolutionary old miRNAs are more likely to have conserved expression patterns than young miRNAs. We conclude that key aspects of the regulatory network are conserved, while differential expression of genes involved in cell-to-cell communication may contribute greatly to phenotypic differences between species.