Integrative analysis of neuroblastoma by single-cell RNA sequencing identifies the NECTIN2-TIGIT axis as a target for immunotherapy.

Pediatric patients with high-risk neuroblastoma have poor survival rates and urgently need more effective treatment options with less side effects. Since novel and improved immunotherapies may fill this need, we dissect the immunoregulatory interactions in neuroblastoma by single-cell RNA-sequencing of 24 tumors (10 pre- and 14 post-chemotherapy, including 5 pairs) to identify strategies for optimizing immunotherapy efficacy. Neuroblastomas are infiltrated by natural killer (NK), T and B cells, and immunosuppressive myeloid populations. NK cells show reduced cytotoxicity and T cells have a dysfunctional profile. Interaction analysis reveals a vast immunoregulatory network and identifies NECTIN2-TIGIT as a crucial immune checkpoint. Combined blockade of TIGIT and PD-L1 significantly reduces neuroblastoma growth, with complete responses (CR) in vivo. Moreover, addition of TIGIT+PD-L1 blockade to standard relapse treatment in a chemotherapy-resistant Th-ALKF1174L/MYCN 129/SvJ syngeneic model induces CR. In conclusion, our integrative analysis provides promising targets and a rationale for immunotherapeutic combination strategies.

An ERK-dependent molecular switch antagonizes fibrosis and promotes regeneration in spiny mice (Acomys).

Although most mammals heal injured tissues and organs with scarring, spiny mice (Acomys) naturally regenerate skin and complex musculoskeletal tissues. Now, the core signaling pathways driving mammalian tissue regeneration are poorly characterized. Here, we show that, while immediate extracellular signal-regulated kinase (ERK) activation is a shared feature of scarring (Mus) and regenerating (Acomys) injuries, ERK activity is only sustained at high levels during complex tissue regeneration. Following ERK inhibition, ear punch regeneration in Acomys shifted toward fibrotic repair. Using single-cell RNA sequencing, we identified ERK-responsive cell types. Loss- and gain-of-function experiments prompted us to uncover fibroblast growth factor and ErbB signaling as upstream ERK regulators of regeneration. The ectopic activation of ERK in scar-prone injuries induced a pro-regenerative response, including cell proliferation, extracellular matrix remodeling, and hair follicle neogenesis. Our data detail an important distinction in ERK activity between regenerating and poorly regenerating adult mammals and open avenues to redirect fibrotic repair toward regenerative healing.

Lysosomal glycogen accumulation in Pompe disease results in disturbed cytoplasmic glycogen metabolism.

Pompe disease is an inherited metabolic myopathy caused by deficiency of acid alpha-glucosidase (GAA), resulting in lysosomal glycogen accumulation. Residual GAA enzyme activity affects disease onset and severity, although other factors, including dysregulation of cytoplasmic glycogen metabolism, are suspected to modulate the disease course. In this study, performed in mice and patient biopsies, we found elevated protein levels of enzymes involved in glucose uptake and cytoplasmic glycogen synthesis in skeletal muscle from mice with Pompe disease, including glycogenin (GYG1), glycogen synthase (GYS1), glucose transporter 4 (GLUT4), glycogen branching enzyme 1 (GBE1), and UDP-glucose pyrophosphorylase (UGP2). Expression levels were elevated before the loss of muscle mass and function. For first time, quantitative mass spectrometry in skeletal muscle biopsies from five adult patients with Pompe disease showed increased expression of GBE1 protein relative to healthy controls at the group level. Paired analysis of individual patients who responded well to treatment with enzyme replacement therapy (ERT) showed reduction of GYS1, GYG1, and GBE1 in all patients after start of ERT compared to baseline. These results indicate that metabolic changes precede muscle wasting in Pompe disease, and imply a positive feedforward loop in Pompe disease, in which lysosomal glycogen accumulation promotes cytoplasmic glycogen synthesis and glucose uptake, resulting in aggravation of the disease phenotype.

Improved Gene Fusion Detection in Childhood Cancer Diagnostics Using RNA Sequencing.

PURPOSE: Gene fusions play a significant role in cancer etiology, making their detection crucial for accurate diagnosis, prognosis, and determining therapeutic targets. Current diagnostic methods largely focus on either targeted or low-resolution genome-wide techniques, which may be unable to capture rare events or both fusion partners. We investigate if RNA sequencing can overcome current limitations with traditional diagnostic techniques to identify gene fusion events. METHODS: We first performed RNA sequencing on a validation cohort of 24 samples with a known gene fusion event, after which a prospective pan-pediatric cancer cohort (n = 244) was tested by RNA sequencing in parallel to existing diagnostic procedures. This cohort included hematologic malignancies, tumors of the CNS, solid tumors, and suspected neoplastic samples. All samples were processed in the routine diagnostic workflow and analyzed for gene fusions using standard-of-care methods and RNA sequencing. RESULTS: We identified a clinically relevant gene fusion in 83 of 244 cases in the prospective cohort. Sixty fusions were detected by both routine diagnostic techniques and RNA sequencing, and one fusion was detected only in routine diagnostics, but an additional 24 fusions were detected solely by RNA sequencing. RNA sequencing, therefore, increased the diagnostic yield by 38%-39%. In addition, RNA sequencing identified both gene partners involved in the gene fusion, in contrast to most routine techniques. For two patients, the newly identified fusion by RNA sequencing resulted in treatment with targeted agents. CONCLUSION: We show that RNA sequencing is sufficiently robust for gene fusion detection in routine diagnostics of childhood cancers and can make a difference in treatment decisions.

LGR6 marks nephron progenitor cells.

BACKGROUND: Nephron progenitor cells (NPCs) undergo a stepwise process to generate all mature nephron structures. Mesenchymal to epithelial transition (MET) is considered a multistep process of NPC differentiation to ensure progressive establishment of new nephrons. However, despite this important role, to date, no marker for NPCs undergoing MET in the nephron exists. RESULTS: Here, we identify LGR6 as a NPC marker, expressed in very early cap mesenchyme, pre-tubular aggregates, renal vesicles, and in segments of S-shaped bodies, following the trajectory of MET. By using a lineage tracing approach in embryonic explants in combination with confocal imaging and single-cell RNA sequencing, we provide evidence for the multiple fates of LGR6+ cells during embryonic nephrogenesis. Moreover, by using long-term in vivo lineage tracing, we show that postnatal LGR6+ cells are capable of generating the multiple lineages of the nephrons. CONCLUSIONS: Given the profound early mesenchymal expression and MET signature of LGR6+ cells, together with the lineage tracing of mesenchymal LGR6+ cells, we conclude that LGR6+ cells contribute to all nephrogenic segments by undergoing MET. LGR6+ cells can therefore be considered an early committed NPC population during embryonic and postnatal nephrogenesis with potential regenerative capability.

Tumor to normal single-cell mRNA comparisons reveal a pan-neuroblastoma cancer cell.

Neuroblastoma is a childhood cancer that resembles developmental stages of the neural crest. It is not established what developmental processes neuroblastoma cancer cells represent. Here, we sought to reveal the phenotype of neuroblastoma cancer cells by comparing cancer (n = 19,723) with normal fetal adrenal single-cell transcriptomes (n = 57,972). Our principal finding was that the neuroblastoma cancer cell resembled fetal sympathoblasts, but no other fetal adrenal cell type. The sympathoblastic state was a universal feature of neuroblastoma cells, transcending cell cluster diversity, individual patients, and clinical phenotypes. We substantiated our findings in 650 neuroblastoma bulk transcriptomes and by integrating canonical features of the neuroblastoma genome with transcriptional signals. Overall, our observations indicate that a pan-neuroblastoma cancer cell state exists, which may be attractive for novel immunotherapeutic and targeted avenues.

Genome-wide off-rates reveal how DNA binding dynamics shape transcription factor function.

Protein-DNA interactions are dynamic, and these dynamics are an important aspect of chromatin-associated processes such as transcription or replication. Due to a lack of methods to study on- and off-rates across entire genomes, protein-DNA interaction dynamics have not been studied extensively. Here, we determine in vivo off-rates for the Saccharomyces cerevisiae chromatin organizing factor Abf1, at 191 sites simultaneously across the yeast genome. Average Abf1 residence times span a wide range, varying between 4.2 and 33 min. Sites with different off-rates are associated with different functional characteristics. This includes their transcriptional dependency on Abf1, nucleosome positioning and the size of the nucleosome-free region, as well as the ability to roadblock RNA polymerase II for termination. The results show how off-rates contribute to transcription factor function and that DIVORSEQ (Determining In Vivo Off-Rates by SEQuencing) is a meaningful way of investigating protein-DNA binding dynamics genome-wide.

An organoid biobank for childhood kidney cancers that captures disease and tissue heterogeneity.

Kidney tumours are among the most common solid tumours in children, comprising distinct subtypes differing in many aspects, including cell-of-origin, genetics, and pathology. Pre-clinical cell models capturing the disease heterogeneity are currently lacking. Here, we describe the first paediatric cancer organoid biobank. It contains tumour and matching normal kidney organoids from over 50 children with different subtypes of kidney cancer, including Wilms tumours, malignant rhabdoid tumours, renal cell carcinomas, and congenital mesoblastic nephromas. Paediatric kidney tumour organoids retain key properties of native tumours, useful for revealing patient-specific drug sensitivities. Using single cell RNA-sequencing and high resolution 3D imaging, we further demonstrate that organoid cultures derived from Wilms tumours consist of multiple different cell types, including epithelial, stromal and blastemal-like cells. Our organoid biobank captures the heterogeneity of paediatric kidney tumours, providing a representative collection of well-characterised models for basic cancer research, drug-screening and personalised medicine.

CHETAH: a selective, hierarchical cell type identification method for single-cell RNA sequencing.

Cell type identification is essential for single-cell RNA sequencing (scRNA-seq) studies, currently transforming the life sciences. CHETAH (CHaracterization of cEll Types Aided by Hierarchical classification) is an accurate cell type identification algorithm that is rapid and selective, including the possibility of intermediate or unassigned categories. Evidence for assignment is based on a classification tree of previously available scRNA-seq reference data and includes a confidence score based on the variance in gene expression per cell type. For cell types represented in the reference data, CHETAH's accuracy is as good as existing methods. Its specificity is superior when cells of an unknown type are encountered, such as malignant cells in tumor samples which it pinpoints as intermediate or unassigned. Although designed for tumor samples in particular, the use of unassigned and intermediate types is also valuable in other exploratory studies. This is exemplified in pancreas datasets where CHETAH highlights cell populations not well represented in the reference dataset, including cells with profiles that lie on a continuum between that of acinar and ductal cell types. Having the possibility of unassigned and intermediate cell types is pivotal for preventing misclassification and can yield important biological information for previously unexplored tissues.

The ability of transcription factors to differentially regulate gene expression is a crucial component of the mechanism underlying inversion, a frequently observed genetic interaction pattern.

Genetic interactions, a phenomenon whereby combinations of mutations lead to unexpected effects, reflect how cellular processes are wired and play an important role in complex genetic diseases. Understanding the molecular basis of genetic interactions is crucial for deciphering pathway organization as well as understanding the relationship between genetic variation and disease. Several hypothetical molecular mechanisms have been linked to different genetic interaction types. However, differences in genetic interaction patterns and their underlying mechanisms have not yet been compared systematically between different functional gene classes. Here, differences in the occurrence and types of genetic interactions are compared for two classes, gene-specific transcription factors (GSTFs) and signaling genes (kinases and phosphatases). Genome-wide gene expression data for 63 single and double deletion mutants in baker's yeast reveals that the two most common genetic interaction patterns are buffering and inversion. Buffering is typically associated with redundancy and is well understood. In inversion, genes show opposite behavior in the double mutant compared to the corresponding single mutants. The underlying mechanism is poorly understood. Although both classes show buffering and inversion patterns, the prevalence of inversion is much stronger in GSTFs. To decipher potential mechanisms, a Petri Net modeling approach was employed, where genes are represented as nodes and relationships between genes as edges. This allowed over 9 million possible three and four node models to be exhaustively enumerated. The models show that a quantitative difference in interaction strength is a strict requirement for obtaining inversion. In addition, this difference is frequently accompanied with a second gene that shows buffering. Taken together, these results provide a mechanistic explanation for inversion. Furthermore, the ability of transcription factors to differentially regulate expression of their targets provides a likely explanation why inversion is more prevalent for GSTFs compared to kinases and phosphatases.

Developing a strategy for computational lab skills training through Software and Data Carpentry: Experiences from the ELIXIR Pilot action.

Quality training in computational skills for life scientists is essential to allow them to deliver robust, reproducible and cutting-edge research. A pan-European bioinformatics programme, ELIXIR, has adopted a well-established and progressive programme of computational lab and data skills training from Software and Data Carpentry, aimed at increasing the number of skilled life scientists and building a sustainable training community in this field. This article describes the Pilot action, which introduced the Carpentry training model to the ELIXIR community.

Molecular mechanisms that distinguish TFIID housekeeping from regulatable SAGA promoters.

An important distinction is frequently made between constitutively expressed housekeeping genes versus regulated genes. Although generally characterized by different DNA elements, chromatin architecture and cofactors, it is not known to what degree promoter classes strictly follow regulatability rules and which molecular mechanisms dictate such differences. We show that SAGA-dominated/TATA-box promoters are more responsive to changes in the amount of activator, even compared to TFIID/TATA-like promoters that depend on the same activator Hsf1. Regulatability is therefore an inherent property of promoter class. Further analyses show that SAGA/TATA-box promoters are more dynamic because TATA-binding protein recruitment through SAGA is susceptible to removal by Mot1. In addition, the nucleosome configuration upon activator depletion shifts on SAGA/TATA-box promoters and seems less amenable to preinitiation complex formation. The results explain the fundamental difference between housekeeping and regulatable genes, revealing an additional facet of combinatorial control: an activator can elicit a different response dependent on core promoter class.

A high-resolution gene expression atlas of epistasis between gene-specific transcription factors exposes potential mechanisms for genetic interactions.

BACKGROUND: Genetic interactions, or non-additive effects between genes, play a crucial role in many cellular processes and disease. Which mechanisms underlie these genetic interactions has hardly been characterized. Understanding the molecular basis of genetic interactions is crucial in deciphering pathway organization and understanding the relationship between genotype, phenotype and disease. RESULTS: To investigate the nature of genetic interactions between gene-specific transcription factors (GSTFs) in Saccharomyces cerevisiae, we systematically analyzed 72 GSTF pairs by gene expression profiling double and single deletion mutants. These pairs were selected through previously published growth-based genetic interactions as well as through similarity in DNA binding properties. The result is a high-resolution atlas of gene expression-based genetic interactions that provides systems-level insight into GSTF epistasis. The atlas confirms known genetic interactions and exposes new ones. Importantly, the data can be used to investigate mechanisms that underlie individual genetic interactions. Two molecular mechanisms are proposed, "buffering by induced dependency" and "alleviation by derepression". CONCLUSIONS: These mechanisms indicate how negative genetic interactions can occur between seemingly unrelated parallel pathways and how positive genetic interactions can indirectly expose parallel rather than same-pathway relationships. The focus on GSTFs is important for understanding the transcription regulatory network of yeast as it uncovers details behind many redundancy relationships, some of which are completely new. In addition, the study provides general insight into the complex nature of epistasis and proposes mechanistic models for genetic interactions, the majority of which do not fall into easily recognizable within- or between-pathway relationships.

FOXO target gene CTDSP2 regulates cell cycle progression through Ras and p21(Cip1/Waf1).

Activity of FOXO (forkhead box O) transcription factors is inhibited by growth factor-PI3K (phosphoinositide 3-kinase)-PKB (protein kinase B)/Akt signalling to control a variety of cellular processes including cell cycle progression. Through comparative analysis of a number of microarray datasets we identified a set of genes commonly regulated by FOXO proteins and PI3K-PKB/Akt, which includes CTDSP2 (C-terminal domain small phosphatase 2). We validated CTDSP2 as a genuine FOXO target gene and show that ectopic CTDSP2 can induce cell cycle arrest. We analysed transcriptional regulation after CTDSP2 expression and identified extensive regulation of genes involved in cell cycle progression, which depends on the phosphatase activity of CTDSP2. The most notably regulated gene is the CDK (cyclin-dependent kinase) inhibitor p21(Cip1/Waf1) and in the present study we show that p21(Cip1/Waf1) is partially responsible for the cell cycle arrest through decreasing cyclin-CDK activity. Our data suggest that CTDSP2 induces p21(Cip1/Waf1) through increasing the activity of Ras. As has been described previously, Ras induces p21(Cip1/Waf1) through p53-dependent and p53-independent pathways and indeed both p53 and MEK inhibition can mitigate the CTDSP2-induced p21(Cip1/Waf1) mRNA up-regulation. In support of Ras activation by CTDSP2, depletion of endogenous CTDSP2 results in reduced Ras activity and thus CTDSP2 seems to be part of a larger set of genes regulated by FOXO proteins, which increase growth factor signalling upon FOXO activation.

Cell cycle population effects in perturbation studies.

Growth condition perturbation or gene function disruption are commonly used strategies to study cellular systems. Although it is widely appreciated that such experiments may involve indirect effects, these frequently remain uncharacterized. Here, analysis of functionally unrelated Saccharyomyces cerevisiae deletion strains reveals a common gene expression signature. One property shared by these strains is slower growth, with increased presence of the signature in more slowly growing strains. The slow growth signature is highly similar to the environmental stress response (ESR), an expression response common to diverse environmental perturbations. Both environmental and genetic perturbations result in growth rate changes. These are accompanied by a change in the distribution of cells over different cell cycle phases. Rather than representing a direct expression response in single cells, both the slow growth signature and ESR mainly reflect a redistribution of cells over different cell cycle phases, primarily characterized by an increase in the G1 population. The findings have implications for any study of perturbation that is accompanied by growth rate changes. Strategies to counter these effects are presented and discussed.

Large-scale genetic perturbations reveal regulatory networks and an abundance of gene-specific repressors.

To understand regulatory systems, it would be useful to uniformly determine how different components contribute to the expression of all other genes. We therefore monitored mRNA expression genome-wide, for individual deletions of one-quarter of yeast genes, focusing on (putative) regulators. The resulting genetic perturbation signatures reflect many different properties. These include the architecture of protein complexes and pathways, identification of expression changes compatible with viability, and the varying responsiveness to genetic perturbation. The data are assembled into a genetic perturbation network that shows different connectivities for different classes of regulators. Four feed-forward loop (FFL) types are overrepresented, including incoherent type 2 FFLs that likely represent feedback. Systematic transcription factor classification shows a surprisingly high abundance of gene-specific repressors, suggesting that yeast chromatin is not as generally restrictive to transcription as is often assumed. The data set is useful for studying individual genes and for discovering properties of an entire regulatory system.

A novel microscopy-based high-throughput screening method to identify proteins that regulate global histone modification levels.

Posttranslational modifications of histones play an important role in the regulation of gene expression and chromatin structure in eukaryotes. The balance between chromatin factors depositing (writers) and removing (erasers) histone marks regulates the steady-state levels of chromatin modifications. Here we describe a novel microscopy-based screening method to identify proteins that regulate histone modification levels in a high-throughput fashion. We named our method CROSS, for Chromatin Regulation Ontology SiRNA Screening. CROSS is based on an siRNA library targeting the expression of 529 proteins involved in chromatin regulation. As a proof of principle, we used CROSS to identify chromatin factors involved in histone H3 methylation on either lysine-4 or lysine-27. Furthermore, we show that CROSS can be used to identify chromatin factors that affect growth in cancer cell lines. Taken together, CROSS is a powerful method to identify the writers and erasers of novel and known chromatin marks and facilitates the identification of drugs targeting epigenetic modifications.

Centromere binding and a conserved role in chromosome stability for SUMO-dependent ubiquitin ligases.

The Saccharomyces cerevisiae Slx5/8 complex is the founding member of a recently defined class of SUMO-targeted ubiquitin ligases (STUbLs). Slx5/8 has been implicated in genome stability and transcription, but the precise contribution is unclear. To characterise Slx5/8 function, we determined genome-wide changes in gene expression upon loss of either subunit. The majority of mRNA changes are part of a general stress response, also exhibited by mutants of other genome integrity pathways and therefore indicative of an indirect effect on transcription. Genome-wide binding analysis reveals a uniquely centromeric location for Slx5. Detailed phenotype analyses of slx5Δ and slx8Δ mutants show severe mitotic defects that include aneuploidy, spindle mispositioning, fish hooks and aberrant spindle kinetics. This is associated with accumulation of the PP2A regulatory subunit Rts1 at centromeres prior to entry into anaphase. Knockdown of the human STUbL orthologue RNF4 also results in chromosome segregation errors due to chromosome bridges. The study shows that STUbLs have a conserved role in maintenance of chromosome stability and links SUMO-dependent ubiquitination to a centromere-specific function during mitosis.

The specificity and topology of chromatin interaction pathways in yeast.

Packaging of DNA into chromatin has a profound impact on gene expression. To understand how changes in chromatin influence transcription, we analyzed 165 mutants of chromatin machinery components in Saccharomyces cerevisiae. mRNA expression patterns change in 80% of mutants, always with specific effects, even for loss of widespread histone marks. The data are assembled into a network of chromatin interaction pathways. The network is function based, has a branched, interconnected topology, and lacks strict one-to-one relationships between complexes. Chromatin pathways are not separate entities for different gene sets, but share many components. The study evaluates which interactions are important for which genes and predicts additional interactions, for example between Paf1C and Set3C, as well as a role for Mediator in subtelomeric silencing. The results indicate the presence of gene-dependent effects that go beyond context-dependent binding of chromatin factors and provide a framework for understanding how specificity is achieved through regulating chromatin.

Gcn4 misregulation reveals a direct role for the evolutionary conserved EKC/KEOPS in the t6A modification of tRNAs.

The EKC/KEOPS complex is universally conserved in Archaea and Eukarya and has been implicated in several cellular processes, including transcription, telomere homeostasis and genomic instability. However, the molecular function of the complex has remained elusive so far. We analyzed the transcriptome of EKC/KEOPS mutants and observed a specific profile that is highly enriched in targets of the Gcn4p transcriptional activator. GCN4 expression was found to be activated at the translational level in mutants via the defective recognition of the inhibitory upstream ORFs (uORFs) present in its leader. We show that EKC/KEOPS mutants are defective for the N6-threonylcarbamoyl adenosine modification at position 37 (t(6)A(37)) of tRNAs decoding ANN codons, which affects initiation at the inhibitory uORFs and provokes Gcn4 de-repression. Structural modeling reveals similarities between Kae1 and bacterial enzymes involved in carbamoylation reactions analogous to t(6)A(37) formation, supporting a direct role for the EKC in tRNA modification. These findings are further supported by strong genetic interactions of EKC mutants with a translation initiation factor and with threonine biosynthesis genes. Overall, our data provide a novel twist to understanding the primary function of the EKC/KEOPS and its impact on several essential cellular functions like transcription and telomere homeostasis.