Scalable integration of multiomic single-cell data using generative adversarial networks.

MOTIVATION: Single-cell profiling has become a common practice to investigate the complexity of tissues, organs, and organisms. Recent technological advances are expanding our capabilities to profile various molecular layers beyond the transcriptome such as, but not limited to, the genome, the epigenome, and the proteome. Depending on the experimental procedure, these data can be obtained from separate assays or the very same cells. Yet, integration of more than two assays is currently not supported by the majority of the computational frameworks avaiable. RESULTS: We here propose a Multi-Omic data integration framework based on Wasserstein Generative Adversarial Networks suitable for the analysis of paired or unpaired data with a high number of modalities (>2). At the core of our strategy is a single network trained on all modalities together, limiting the computational burden when many molecular layers are evaluated. AVAILABILITY AND IMPLEMENTATION: Source code of our framework is available at https://github.com/vgiansanti/MOWGAN.

Dendritic cell-targeted therapy expands CD8 T cell responses to bona-fide neoantigens in lung tumors.

Cross-presentation by type 1 DCs (cDC1) is critical to induce and sustain antitumoral CD8 T cell responses to model antigens, in various tumor settings. However, the impact of cross-presenting cDC1 and the potential of DC-based therapies in tumors carrying varied levels of bona-fide neoantigens (neoAgs) remain unclear. Here we develop a hypermutated model of non-small cell lung cancer in female mice, encoding genuine MHC-I neoepitopes to study neoAgs-specific CD8 T cell responses in spontaneous settings and upon Flt3L + αCD40 (DC-therapy). We find that cDC1 are required to generate broad CD8 responses against a range of diverse neoAgs. DC-therapy promotes immunogenicity of weaker neoAgs and strongly inhibits the growth of high tumor-mutational burden (TMB) tumors. In contrast, low TMB tumors respond poorly to DC-therapy, generating mild CD8 T cell responses that are not sufficient to block progression. scRNA transcriptional analysis, immune profiling and functional assays unveil the changes induced by DC-therapy in lung tissues, which comprise accumulation of cDC1 with increased immunostimulatory properties and less exhausted effector CD8 T cells. We conclude that boosting cDC1 activity is critical to broaden the diversity of anti-tumoral CD8 T cell responses and to leverage neoAgs content for therapeutic advantage.

Evolutionary fingerprints of EMT in pancreatic cancers.

Mesenchymal plasticity has been extensively described in advanced and metastatic epithelial cancers; however, its functional role in malignant progression, metastatic dissemination and therapy response is controversial. More importantly, the role of epithelial mesenchymal transition (EMT) and cell plasticity in tumor heterogeneity, clonal selection and clonal evolution is poorly understood. Functionally, our work clarifies the contribution of EMT to malignant progression and metastasis in pancreatic cancer. We leveraged ad hoc somatic mosaic genome engineering, lineage tracing and ablation technologies and dynamic genetic reporters to trace and ablate tumor-specific lineages along the phenotypic spectrum of epithelial to mesenchymal plasticity. The experimental evidences clarify the essential contribution of mesenchymal lineages to pancreatic cancer evolution and metastatic dissemination. Spatial genomic analysis combined with single cell transcriptomic and epigenomic profiling of epithelial and mesenchymal lineages reveals that EMT promotes with the emergence of chromosomal instability (CIN). Specifically tumor lineages with mesenchymal features display highly conserved patterns of genomic evolution including complex structural genomic rearrangements and chromotriptic events. Genetic ablation of mesenchymal lineages robustly abolished these mutational processes and evolutionary patterns, as confirmed by cross species analysis of pancreatic and other human epithelial cancers. Mechanistically, we discovered that malignant cells with mesenchymal features display increased chromatin accessibility, particularly in the pericentromeric and centromeric regions, which in turn results in delayed mitosis and catastrophic cell division. Therefore, EMT favors the emergence of high-fitness tumor cells, strongly supporting the concept of a cell-state, lineage-restricted patterns of evolution, where cancer cell sub-clonal speciation is propagated to progenies only through restricted functional compartments. Restraining those evolutionary routes through genetic ablation of clones capable of mesenchymal plasticity and extinction of the derived lineages completely abrogates the malignant potential of one of the most aggressive form of human cancer.

In Vitro Selection of Engineered Transcriptional Repressors for Targeted Epigenetic Silencing.

Gene inactivation is instrumental to study gene function and represents a promising strategy for the treatment of a broad range of diseases. Among traditional technologies, RNA interference suffers from partial target abrogation and the requirement for life-long treatments. In contrast, artificial nucleases can impose stable gene inactivation through induction of a DNA double strand break (DSB), but recent studies are questioning the safety of this approach. Targeted epigenetic editing via engineered transcriptional repressors (ETRs) may represent a solution, as a single administration of specific ETR combinations can lead to durable silencing without inducing DNA breaks. ETRs are proteins containing a programmable DNA-binding domain (DBD) and effectors from naturally occurring transcriptional repressors. Specifically, a combination of three ETRs equipped with the KRAB domain of human ZNF10, the catalytic domain of human DNMT3A and human DNMT3L, was shown to induce heritable repressive epigenetic states on the ETR-target gene. The hit-and-run nature of this platform, the lack of impact on the DNA sequence of the target, and the possibility to revert to the repressive state by DNA demethylation on demand, make epigenetic silencing a game-changing tool. A critical step is the identification of the proper ETRs' position on the target gene to maximize on-target and minimize off-target silencing. Performing this step in the final ex vivo or in vivo preclinical setting can be cumbersome. Taking the CRISPR/catalytically dead Cas9 system as a paradigmatic DBD for ETRs, this paper describes a protocol consisting of the in vitro screen of guide RNAs (gRNAs) coupled to the triple-ETR combination for efficient on-target silencing, followed by evaluation of the genome-wide specificity profile of top hits. This allows for reduction of the initial repertoire of candidate gRNAs to a short list of promising ones, whose complexity is suitable for their final evaluation in the therapeutically relevant setting of interest.

Enhancing clinical potential of liquid biopsy through a multi-omic approach: A systematic review.

In the last years, liquid biopsy gained increasing clinical relevance for detecting and monitoring several cancer types, being minimally invasive, highly informative and replicable over time. This revolutionary approach can be complementary and may, in the future, replace tissue biopsy, which is still considered the gold standard for cancer diagnosis. "Classical" tissue biopsy is invasive, often cannot provide sufficient bioptic material for advanced screening, and can provide isolated information about disease evolution and heterogeneity. Recent literature highlighted how liquid biopsy is informative of proteomic, genomic, epigenetic, and metabolic alterations. These biomarkers can be detected and investigated using single-omic and, recently, in combination through multi-omic approaches. This review will provide an overview of the most suitable techniques to thoroughly characterize tumor biomarkers and their potential clinical applications, highlighting the importance of an integrated multi-omic, multi-analyte approach. Personalized medical investigations will soon allow patients to receive predictable prognostic evaluations, early disease diagnosis, and subsequent ad hoc treatments.

Analyzing genomic and epigenetic profiles in single cells by hybrid transposase (scGET-seq).

scGET-seq simultaneously profiles euchromatin and heterochromatin. scGET-seq exploits the concurrent action of transposase Tn5 and its hybrid form TnH, which targets H3K9me3 domains. Here we present a step-by-step protocol to profile single cells by scGET-seq using a 10× Chromium Controller. We describe steps for transposomes preparation and validation. We detail nuclei preparation and transposition, followed by encapsulation, library preparation, sequencing, and data analysis. For complete details on the use and execution of this protocol, please refer to Tedesco et al. (2022)1 and de Pretis and Cittaro (2022).2.

Neural precursor cells tune striatal connectivity through the release of IGFBPL1.

The adult brain retains over life endogenous neural stem/precursor cells (eNPCs) within the subventricular zone (SVZ). Whether or not these cells exert physiological functions is still unclear. In the present work, we provide evidence that SVZ-eNPCs tune structural, electrophysiological, and behavioural aspects of striatal function via secretion of insulin-like growth factor binding protein-like 1 (IGFBPL1). In mice, selective ablation of SVZ-eNPCs or selective abrogation of IGFBPL1 determined an impairment of striatal medium spiny neuron morphology, a higher failure rate in GABAergic transmission mediated by fast-spiking interneurons, and striatum-related behavioural dysfunctions. We also found IGFBPL1 expression in the human SVZ, foetal and induced-pluripotent stem cell-derived NPCs. Finally, we found a significant correlation between SVZ damage, reduction of striatum volume, and impairment of information processing speed in neurological patients. Our results highlight the physiological role of adult SVZ-eNPCs in supporting cognitive functions by regulating striatal neuronal activity.

Into the multiverse: advances in single-cell multiomic profiling.

Single-cell transcriptomic approaches have revolutionised the study of complex biological systems, with the routine measurement of gene expression in thousands of cells enabling construction of whole-organism cell atlases. However, the transcriptome is just one layer amongst many that coordinate to define cell type and state and, ultimately, function. In parallel with the widespread uptake of single-cell RNA-seq (scRNA-seq), there has been a rapid emergence of methods that enable multiomic profiling of individual cells, enabling parallel measurement of intercellular heterogeneity in the genome, epigenome, transcriptome, and proteomes. Linking measurements from each of these layers has the potential to reveal regulatory and functional mechanisms underlying cell behaviour in healthy development and disease.

Chromatin Velocity reveals epigenetic dynamics by single-cell profiling of heterochromatin and euchromatin.

Recent efforts have succeeded in surveying open chromatin at the single-cell level, but high-throughput, single-cell assessment of heterochromatin and its underlying genomic determinants remains challenging. We engineered a hybrid transposase including the chromodomain (CD) of the heterochromatin protein-1α (HP-1α), which is involved in heterochromatin assembly and maintenance through its binding to trimethylation of the lysine 9 on histone 3 (H3K9me3), and developed a single-cell method, single-cell genome and epigenome by transposases sequencing (scGET-seq), that, unlike single-cell assay for transposase-accessible chromatin with sequencing (scATAC-seq), comprehensively probes both open and closed chromatin and concomitantly records the underlying genomic sequences. We tested scGET-seq in cancer-derived organoids and human-derived xenograft (PDX) models and identified genetic events and plasticity-driven mechanisms contributing to cancer drug resistance. Next, building upon the differential enrichment of closed and open chromatin, we devised a method, Chromatin Velocity, that identifies the trajectories of epigenetic modifications at the single-cell level. Chromatin Velocity uncovered paths of epigenetic reorganization during stem cell reprogramming and identified key transcription factors driving these developmental processes. scGET-seq reveals the dynamics of genomic and epigenetic landscapes underlying any cellular processes.

Aging, inflammation and DNA damage in the somatic testicular niche with idiopathic germ cell aplasia.

Molecular mechanisms associated with human germ cell aplasia in infertile men remain undefined. Here we perform single-cell transcriptome profiling to highlight differentially expressed genes and pathways in each somatic cell type in testes of men with idiopathic germ cell aplasia. We identify immaturity of Leydig cells, chronic tissue inflammation, fibrosis, and senescence phenotype of the somatic cells, as well markers of chronic inflammation in the blood. We find that deregulated expression of parentally imprinted genes in myoid and immature Leydig cells, with relevant changes in the ratio of Lamin A/C transcripts and an active DNA damage response in Leydig and peritubular myoid cells are also indicative of senescence of the testicular niche. This study offers molecular insights into the pathogenesis of idiopathic germ cell aplasia.

Transcriptional dynamics of induced pluripotent stem cell differentiation into ß cells reveals full endodermal commitment and homology with human islets.

BACKGROUND AIMS: Induced pluripotent stem cells (iPSCs) have the capacity to generate ß cells in vitro, but the differentiation is incomplete and generates a variable percentage of off-target cells. Single-cell RNA sequencing offers the possibility of characterizing the transcriptional dynamics throughout differentiation and determining the identity of the final differentiation product. METHODS: Single-cell transcriptomics data were obtained from four stages across differentiation of iPSCs into ß cells and from human donor islets. RESULTS: Clustering analysis revealed that iPSCs undertake a full endoderm commitment, and the obtained endocrine pancreatic cells have high homology with mature islets. The iPSC-derived ß cells were devoid of pluripotent residual cells, and the differentiation was pancreas-specific, as it did not generate ectodermal or mesodermal cells. Pseudotime trajectory identified a dichotomic endocrine/non-endocrine cell fate and distinct subgroups in the endocrine branch. CONCLUSIONS: Future efforts to produce ß cells from iPSCs must aim not only to improve the resulting endocrine cell but also to avoid differentiation into non-pancreatic endoderm cells.

MSH6 gene pathogenic variant identified in familial pancreatic cancer in the absence of colon cancer.

OBJECTIVES: Lynch syndrome is characterized by pathogenetic variants in the mismatch repair genes and autosomal dominant inheritance with incomplete penetrance. Lynch syndrome is characterized by colorectal and, with lesser and variable extent, extracolonic cancers. We describe a family with MSH6-dependent Lynch syndrome and familial pancreatic cancer and other tumours (gastric and endometrial), in the absence of colorectal neoplasia. METHODS: Patients were analysed by sequencing, Next Generation or Sanger, to identify germinal pathogenic variants in hereditary cancer genes. RESULTS: We identified the MSH6 gene pathogenic variant c.2194C>T, p.(Arg732Ter) in a family with hereditary pancreatic cancer without diagnosed cases of colorectal adenocarcinoma. Seven family members were affected by the MSH6 pathogenic variant. Three had pancreatic adenocarcinoma at 65, 57 and 44 years; one had endometrial cancer at 36 years. None of the remaining three subjects (75, 45 and 17 years old) had developed any cancer yet. CONCLUSIONS: Lynch syndrome should be suspected in families with familial pancreatic cancer, even in the absence of colon cancers. Specifically, our observation supports the association between the MSH6 c.2194C>T pathogenic variant and extracolonic tumours and it suggests that MSH6 pathogenic variants are associated with familial pancreatic cancer more frequently than assumed.

Conserved DNA Methylation Signatures in Early Maternal Separation and in Twins Discordant for CO2 Sensitivity.

Respiratory and emotional responses to blood-acidifying inhalation of CO2 are markers of some human anxiety disorders, and can be enhanced by repeatedly cross-fostering (RCF) mouse pups from their biological mother to unrelated lactating females. Yet, these dynamics remain poorly understood. We show RCF-associated intergenerational transmission of CO2 sensitivity in normally-reared mice descending from RCF-exposed females, and describe the accompanying alterations in brain DNA methylation patterns. These epigenetic signatures were compared to DNA methylation profiles of monozygotic twins discordant for emotional reactivity to a CO2 challenge. Altered methylation was consistently associated with repeated elements and transcriptional regulatory regions among RCF-exposed animals, their normally-reared offspring, and humans with CO2 hypersensitivity. In both species, regions bearing differential methylation were associated with neurodevelopment, circulation, and response to pH acidification processes, and notably included the ASIC2 gene. Our data show that CO2 hypersensitivity is associated with specific methylation clusters and genes that subserve chemoreception and anxiety. The methylation status of genes implicated in acid-sensing functions can inform etiological and therapeutic research in this field.

Structures of purine nucleosidase from Trypanosoma brucei bound to isozyme-specific trypanocidals and a novel metalorganic inhibitor.

Sleeping sickness is a deadly disease that primarily affects sub-Saharan Africa and is caused by protozoan parasites of the Trypanosoma genus. Trypanosomes are purine auxotrophs and their uptake pathway has long been appreciated as an attractive target for drug design. Recently, one tight-binding competitive inhibitor of the trypanosomal purine-specific nucleoside hydrolase (IAGNH) showed remarkable trypanocidal activity in a murine model of infection. Here, the enzymatic characterization of T. brucei brucei IAGNH is presented, together with its high-resolution structures in the unliganded form and in complexes with different inhibitors, including the trypanocidal compound UAMC-00363. A description of the crucial contacts that account for the high-affinity inhibition of IAGNH by iminoribitol-based compounds is provided and the molecular mechanism underlying the conformational change necessary for enzymatic catalysis is identified. It is demonstrated for the first time that metalorganic complexes can compete for binding at the active site of nucleoside hydrolase enzymes, mimicking the positively charged transition state of the enzymatic reaction. Moreover, we show that divalent metal ions can act as noncompetitive IAGNH inhibitors, stabilizing a nonproductive conformation of the catalytic loop. These results open a path for rational improvement of the potency and the selectivity of existing compounds and suggest new scaffolds that may be used as blueprints for the design of novel antitrypanosomal compounds.

Evaluation of nucleoside hydrolase inhibitors for treatment of African trypanosomiasis.

In this paper, we present the biochemical and biological evaluation of N-arylmethyl-substituted iminoribitol derivatives as potential chemotherapeutic agents against trypanosomiasis. Previously, a library of 52 compounds was designed and synthesized as potent and selective inhibitors of Trypanosoma vivax inosine-adenosine-guanosine nucleoside hydrolase (IAG-NH). However, when the compounds were tested against bloodstream-form Trypanosoma brucei brucei, only one inhibitor, N-(9-deaza-adenin-9-yl)methyl-1,4-dideoxy-1,4-imino-d-ribitol (UAMC-00363), displayed significant activity (mean 50% inhibitory concentration [IC(50)] +/- standard error, 0.49 +/- 0.31 microM). Validation in an in vivo model of African trypanosomiasis showed promising results for this compound. Several experiments were performed to investigate why only UAMC-00363 showed antiparasitic activity. First, the compound library was screened against T. b. brucei IAG-NH and inosine-guanosine nucleoside hydrolase (IG-NH) to confirm the previously demonstrated inhibitory effects of the compounds on T. vivax IAG-NH. Second, to verify the uptake of these compounds by T. b. brucei, their affinities for the nucleoside P1 and nucleoside/nucleobase P2 transporters of T. b. brucei were tested. Only UAMC-00363 displayed significant affinity for the P2 transporter. It was also shown that UAMC-00363 is concentrated in the cell via at least one additional transporter, since P2 knockout mutants of T. b. brucei displayed no resistance to the compound. Consequently, no cross-resistance to the diamidine or the melaminophenyl arsenical classes of trypanocides is expected. Third, three enzymes of the purine salvage pathway of procyclic T. b. brucei (IAG-NH, IG-NH, and methylthioadenosine phosphorylase [MTAP]) were investigated using RNA interference. The findings from all these studies showed that it is probably not sufficient to target only the nucleoside hydrolase activity to block the purine salvage pathway of T. b. brucei and that, therefore, it is possible that UAMC-00363 acts on an additional target.

Structural basis for substrate specificity in group I nucleoside hydrolases.

Enzymes with nucleoside hydrolase activity (NHs) belonging to homology group I either are markedly specific for pyrimidine nucleoside substrates or hydrolyze with comparable efficiencies the N-glycosidic bond in all common nucleosides. The biochemical and structural basis for these differences in substrate specificity is still unknown. Here we characterize the binding interactions between the slowly hydrolyzed substrate inosine and the Escherichia coli pyrimidine-specific NH YeiK using cryotrapping and X-ray crystallography. Guided by the structural features of the Michaelis complex, we show the synergic effect of two specific point mutations in YeiK that increase the catalytic efficiency toward purine nucleosides to values comparable to those of natural nonspecific NHs. We demonstrate that the integrity of an active-site catalytic triad comprised of two hydroxylated amino acids and one histidine residue is a requirement for the highly efficient hydrolysis of inosine by group I NHs. Instead, cleavage of the YeiK-preferred substrate uridine is not affected by mutations at the same locations, suggesting a different fine chemical mechanism for the hydrolysis of the two nucleoside substrates. Our study provides for the first time direct evidence that distinct subsets of amino acid residues are involved in the hydrolysis of purine or pyrimidine nucleosides in group I NHs.