Molecular maps of synovial cells in inflammatory arthritis using an optimized synovial tissue dissociation protocol.

In this study, we optimized the dissociation of synovial tissue biopsies for single-cell omics studies and created a single-cell atlas of human synovium in inflammatory arthritis. The optimized protocol allowed consistent isolation of highly viable cells from tiny fresh synovial biopsies, minimizing the synovial biopsy drop-out rate. The synovium scRNA-seq atlas contained over 100,000 unsorted synovial cells from 25 synovial tissues affected by inflammatory arthritis, including 16 structural, 11 lymphoid, and 15 myeloid cell clusters. This synovial cell map expanded the diversity of synovial cell types/states, detected synovial neutrophils, and broadened synovial endothelial cell classification. We revealed tissue-resident macrophage subsets with proposed matrix-sensing (FOLR2+COLEC12high) and iron-recycling (LYVE1+SLC40A1+) activities and identified fibroblast subsets with proposed functions in cartilage breakdown (SOD2highSAA1+SAA2+SDC4+) and extracellular matrix remodeling (SERPINE1+COL5A3+LOXL2+). Our study offers an efficient synovium dissociation method and a reference scRNA-seq resource, that advances the current understanding of synovial cell heterogeneity in inflammatory arthritis.

On the identification of differentially-active transcription factors from ATAC-seq data.

ATAC-seq has emerged as a rich epigenome profiling technique, and is commonly used to identify Transcription Factors (TFs) underlying given phenomena. A number of methods can be used to identify differentially-active TFs through the accessibility of their DNA-binding motif, however little is known on the best approaches for doing so. Here we benchmark several such methods using a combination of curated datasets with various forms of short-term perturbations on known TFs, as well as semi-simulations. We include both methods specifically designed for this type of data as well as some that can be repurposed for it. We also investigate variations to these methods, and identify three particularly promising approaches (chromVAR-limma with critical adjustments, monaLisa and a combination of GC smooth quantile normalization and multivariate modeling). We further investigate the specific use of nucleosome-free fragments, the combination of top methods, and the impact of technical variation. Finally, we illustrate the use of the top methods on a novel dataset to characterize the impact on DNA accessibility of TRAnscription Factor TArgeting Chimeras (TRAFTAC), which can deplete TFs - in our case NFkB - at the protein level.

A model of human neural networks reveals NPTX2 pathology in ALS and FTLD.

Human cellular models of neurodegeneration require reproducibility and longevity, which is necessary for simulating age-dependent diseases. Such systems are particularly needed for TDP-43 proteinopathies1, which involve human-specific mechanisms2-5 that cannot be directly studied in animal models. Here, to explore the emergence and consequences of TDP-43 pathologies, we generated induced pluripotent stem cell-derived, colony morphology neural stem cells (iCoMoNSCs) via manual selection of neural precursors6. Single-cell transcriptomics and comparison to independent neural stem cells7 showed that iCoMoNSCs are uniquely homogenous and self-renewing. Differentiated iCoMoNSCs formed a self-organized multicellular system consisting of synaptically connected and electrophysiologically active neurons, which matured into long-lived functional networks (which we designate iNets). Neuronal and glial maturation in iNets was similar to that of cortical organoids8. Overexpression of wild-type TDP-43 in a minority of neurons within iNets led to progressive fragmentation and aggregation of the protein, resulting in a partial loss of function and neurotoxicity. Single-cell transcriptomics revealed a novel set of misregulated RNA targets in TDP-43-overexpressing neurons and in patients with TDP-43 proteinopathies exhibiting a loss of nuclear TDP-43. The strongest misregulated target encoded the synaptic protein NPTX2, the levels of which are controlled by TDP-43 binding on its 3' untranslated region. When NPTX2 was overexpressed in iNets, it exhibited neurotoxicity, whereas correcting NPTX2 misregulation partially rescued neurons from TDP-43-induced neurodegeneration. Notably, NPTX2 was consistently misaccumulated in neurons from patients with amyotrophic lateral sclerosis and frontotemporal lobar degeneration with TDP-43 pathology. Our work directly links TDP-43 misregulation and NPTX2 accumulation, thereby revealing a TDP-43-dependent pathway of neurotoxicity.

DESpace: spatially variable gene detection via differential expression testing of spatial clusters.

MOTIVATION: Spatially resolved transcriptomics (SRT) enables scientists to investigate spatial context of mRNA abundance, including identifying spatially variable genes (SVGs), i.e. genes whose expression varies across the tissue. Although several methods have been proposed for this task, native SVG tools cannot jointly model biological replicates, or identify the key areas of the tissue affected by spatial variability. RESULTS: Here, we introduce DESpace, a framework, based on an original application of existing methods, to discover SVGs. In particular, our approach inputs all types of SRT data, summarizes spatial information via spatial clusters, and identifies spatially variable genes by performing differential gene expression testing between clusters. Furthermore, our framework can identify (and test) the main cluster of the tissue affected by spatial variability; this allows scientists to investigate spatial expression changes in specific areas of interest. Additionally, DESpace enables joint modeling of multiple samples (i.e. biological replicates); compared to inference based on individual samples, this approach increases statistical power, and targets SVGs with consistent spatial patterns across replicates. Overall, in our benchmarks, DESpace displays good true positive rates, controls for false positive and false discovery rates, and is computationally efficient. AVAILABILITY AND IMPLEMENTATION: DESpace is freely distributed as a Bioconductor R package at https://bioconductor.org/packages/DESpace.

Rbm8a deficiency causes hematopoietic defects by modulating Wnt/PCP signaling.

Defects in blood development frequently occur among syndromic congenital anomalies. Thrombocytopenia-Absent Radius (TAR) syndrome is a rare congenital condition with reduced platelets (hypomegakaryocytic thrombocytopenia) and forelimb anomalies, concurrent with more variable heart and kidney defects. TAR syndrome associates with hypomorphic gene function for RBM8A/Y14 that encodes a component of the exon junction complex involved in mRNA splicing, transport, and nonsense-mediated decay. How perturbing a general mRNA-processing factor causes the selective TAR Syndrome phenotypes remains unknown. Here, we connect zebrafish rbm8a perturbation to early hematopoietic defects via attenuated non-canonical Wnt/Planar Cell Polarity (PCP) signaling that controls developmental cell re-arrangements. In hypomorphic rbm8a zebrafish, we observe a significant reduction of cd41-positive thrombocytes. rbm8a-mutant zebrafish embryos accumulate mRNAs with individual retained introns, a hallmark of defective nonsense-mediated decay; affected mRNAs include transcripts for non-canonical Wnt/PCP pathway components. We establish that rbm8a-mutant embryos show convergent extension defects and that reduced rbm8a function interacts with perturbations in non-canonical Wnt/PCP pathway genes wnt5b, wnt11f2, fzd7a, and vangl2. Using live-imaging, we found reduced rbm8a function impairs the architecture of the lateral plate mesoderm (LPM) that forms hematopoietic, cardiovascular, kidney, and forelimb skeleton progenitors as affected in TAR Syndrome. Both mutants for rbm8a and for the PCP gene vangl2 feature impaired expression of early hematopoietic/endothelial genes including runx1 and the megakaryocyte regulator gfi1aa. Together, our data propose aberrant LPM patterning and hematopoietic defects as consequence of attenuated non-canonical Wnt/PCP signaling upon reduced rbm8a function. These results also link TAR Syndrome to a potential LPM origin and a developmental mechanism.

RNA recognition by Npl3p reveals U2 snRNA-binding compatible with a chaperone role during splicing.

The conserved SR-like protein Npl3 promotes splicing of diverse pre-mRNAs. However, the RNA sequence(s) recognized by the RNA Recognition Motifs (RRM1 & RRM2) of Npl3 during the splicing reaction remain elusive. Here, we developed a split-iCRAC approach in yeast to uncover the consensus sequence bound to each RRM. High-resolution NMR structures show that RRM2 recognizes a 5´-GNGG-3´ motif leading to an unusual mille-feuille topology. These structures also reveal how RRM1 preferentially interacts with a CC-dinucleotide upstream of this motif, and how the inter-RRM linker and the region C-terminal to RRM2 contribute to cooperative RNA-binding. Structure-guided functional studies show that Npl3 genetically interacts with U2 snRNP specific factors and we provide evidence that Npl3 melts U2 snRNA stem-loop I, a prerequisite for U2/U6 duplex formation within the catalytic center of the Bact spliceosomal complex. Thus, our findings suggest an unanticipated RNA chaperoning role for Npl3 during spliceosome active site formation.

Validation of hypermethylated DNA regions found in colorectal cancers as potential aging-independent biomarkers of precancerous colorectal lesions.

BACKGROUND: We previously identified 16,772 colorectal cancer-associated hypermethylated DNA regions that were also detectable in precancerous colorectal lesions (preCRCs) and unrelated to normal mucosal aging. We have now conducted a study to validate 990 of these differentially methylated DNA regions (DMRs) in a new series of preCRCs. METHODS: We used targeted bisulfite sequencing to validate these 990 potential biomarkers in 59 preCRC tissue samples (41 conventional adenomas, 18 sessile serrated lesions), each with a patient-matched normal mucosal sample. Based on differential DNA methylation tests, a panel of candidate DMRs was chosen on a subset of our cohort and then validated on the remaining part of our cohort and two publicly available datasets with respect to their stratifying potential between preCRCs and normal mucosa. RESULTS: Strong statistical significance for the difference in methylation levels was observed across the full set of 990 investigated DMRs. From these, a selected candidate panel of 30 DMRs correctly identified 58/59 tumors (area under the receiver operating curve: 0.998). CONCLUSIONS: These validated DNA hypermethylation markers can be exploited to develop more accurate noninvasive colorectal tumor screening assays.

DifferentialRegulation: a Bayesian hierarchical approach to identify differentially regulated genes.

Motivation: Although transcriptomics data is typically used to analyse mature spliced mRNA, recent attention has focused on jointly investigating spliced and unspliced (or precursor-) mRNA, which can be used to study gene regulation and changes in gene expression production. Nonetheless, most methods for spliced/unspliced inference (such as RNA velocity tools) focus on individual samples, and rarely allow comparisons between groups of samples (e.g., healthy vs. diseased). Furthermore, this kind of inference is challenging, because spliced and unspliced mRNA abundance is characterized by a high degree of quantification uncertainty, due to the prevalence of multi-mapping reads, i.e., reads compatible with multiple transcripts (or genes), and/or with both their spliced and unspliced versions. Results: Here, we present DifferentialRegulation, a Bayesian hierarchical method to discover changes between experimental conditions with respect to the relative abundance of unspliced mRNA (over the total mRNA). We model the quantification uncertainty via a latent variable approach, where reads are allocated to their gene/transcript of origin, and to the respective splice version. We designed several benchmarks where our approach shows good performance, in terms of sensitivity and error control, versus state-of-the-art competitors. Importantly, our tool is flexible, and works with both bulk and single-cell RNA-sequencing data. Availability and implementation: DifferentialRegulation is distributed as a Bioconductor R package.

Meta-analysis of (single-cell method) benchmarks reveals the need for extensibility and interoperability.

Computational methods represent the lifeblood of modern molecular biology. Benchmarking is important for all methods, but with a focus here on computational methods, benchmarking is critical to dissect important steps of analysis pipelines, formally assess performance across common situations as well as edge cases, and ultimately guide users on what tools to use. Benchmarking can also be important for community building and advancing methods in a principled way. We conducted a meta-analysis of recent single-cell benchmarks to summarize the scope, extensibility, and neutrality, as well as technical features and whether best practices in open data and reproducible research were followed. The results highlight that while benchmarks often make code available and are in principle reproducible, they remain difficult to extend, for example, as new methods and new ways to assess methods emerge. In addition, embracing containerization and workflow systems would enhance reusability of intermediate benchmarking results, thus also driving wider adoption.

Conserved stromal-immune cell circuits secure B cell homeostasis and function.

B cell zone reticular cells (BRCs) form stable microenvironments that direct efficient humoral immunity with B cell priming and memory maintenance being orchestrated across lymphoid organs. However, a comprehensive understanding of systemic humoral immunity is hampered by the lack of knowledge of global BRC sustenance, function and major pathways controlling BRC-immune cell interactions. Here we dissected the BRC landscape and immune cell interactome in human and murine lymphoid organs. In addition to the major BRC subsets underpinning the follicle, including follicular dendritic cells, PI16+ RCs were present across organs and species. As well as BRC-produced niche factors, immune cell-driven BRC differentiation and activation programs governed the convergence of shared BRC subsets, overwriting tissue-specific gene signatures. Our data reveal that a canonical set of immune cell-provided cues enforce bidirectional signaling programs that sustain functional BRC niches across lymphoid organs and species, thereby securing efficient humoral immunity.

Gapless provides combined scaffolding, gap filling, and assembly correction with long reads.

Continuity, correctness, and completeness of genome assemblies are important for many biological projects. Long reads represent a major driver towards delivering high-quality genomes, but not everybody can achieve the necessary coverage for good long read-only assemblies. Therefore, improving existing assemblies with low-coverage long reads is a promising alternative. The improvements include correction, scaffolding, and gap filling. However, most tools perform only one of these tasks and the useful information of reads that supported the scaffolding is lost when running separate programs successively. Therefore, we propose a new tool for combined execution of all three tasks using PacBio or Oxford Nanopore reads. gapless is available at: https://github.com/schmeing/gapless.

Transformation of alignment files improves performance of variant callers for long-read RNA sequencing data.

Long-read RNA sequencing (lrRNA-seq) produces detailed information about full-length transcripts, including novel and sample-specific isoforms. Furthermore, there is an opportunity to call variants directly from lrRNA-seq data. However, most state-of-the-art variant callers have been developed for genomic DNA. Here, there are two objectives: first, we perform a mini-benchmark on GATK, DeepVariant, Clair3, and NanoCaller primarily on PacBio Iso-Seq, data, but also on Nanopore and Illumina RNA-seq data; second, we propose a pipeline to process spliced-alignment files, making them suitable for variant calling with DNA-based callers. With such manipulations, high calling performance can be achieved using DeepVariant on Iso-seq data.

The shaky foundations of simulating single-cell RNA sequencing data.

BACKGROUND: With the emergence of hundreds of single-cell RNA-sequencing (scRNA-seq) datasets, the number of computational tools to analyze aspects of the generated data has grown rapidly. As a result, there is a recurring need to demonstrate whether newly developed methods are truly performant-on their own as well as in comparison to existing tools. Benchmark studies aim to consolidate the space of available methods for a given task and often use simulated data that provide a ground truth for evaluations, thus demanding a high quality standard results credible and transferable to real data. RESULTS: Here, we evaluated methods for synthetic scRNA-seq data generation in their ability to mimic experimental data. Besides comparing gene- and cell-level quality control summaries in both one- and two-dimensional settings, we further quantified these at the batch- and cluster-level. Secondly, we investigate the effect of simulators on clustering and batch correction method comparisons, and, thirdly, which and to what extent quality control summaries can capture reference-simulation similarity. CONCLUSIONS: Our results suggest that most simulators are unable to accommodate complex designs without introducing artificial effects, they yield over-optimistic performance of integration and potentially unreliable ranking of clustering methods, and it is generally unknown which summaries are important to ensure effective simulation-based method comparisons.

Complement downregulation promotes an inflammatory signature that renders colorectal cancer susceptible to immunotherapy.

BACKGROUND AND AIMS: The role of inflammatory immune responses in colorectal cancer (CRC) development and response to therapy is a matter of intense debate. While inflammation is a known driver of CRC, inflammatory immune infiltrates are a positive prognostic factor in CRC and predispose to response to immune checkpoint blockade (ICB) therapy. Unfortunately, over 85% of CRC cases are primarily unresponsive to ICB due to the absence of an immune infiltrate, and even the cases that show an initial immune infiltration can become refractory to ICB. The identification of therapy supportive immune responses in the field has been partially hindered by the sparsity of suitable mouse models to recapitulate the human disease. In this study, we aimed to understand how the dysregulation of the complement anaphylatoxin C3a receptor (C3aR), observed in subsets of patients with CRC, affects the immune responses, the development of CRC, and response to ICB therapy. METHODS: We use a comprehensive approach encompassing analysis of publicly available human CRC datasets, inflammation-driven and newly generated spontaneous mouse models of CRC, and multiplatform high-dimensional analysis of immune responses using microbiota sequencing, RNA sequencing, and mass cytometry. RESULTS: We found that patients' regulation of the complement C3aR is associated with epigenetic modifications. Specifically, downregulation of C3ar1 in human CRC promotes a tumor microenvironment characterized by the accumulation of innate and adaptive immune cells that support antitumor immunity. In addition, in vivo studies in our newly generated mouse model revealed that the lack of C3a in the colon activates a microbiota-mediated proinflammatory program which promotes the development of tumors with an immune signature that renders them responsive to the ICB therapy. CONCLUSIONS: Our findings reveal that C3aR may act as a previously unrecognized checkpoint to enhance antitumor immunity in CRC. C3aR can thus be exploited to overcome ICB resistance in a larger group of patients with CRC.

An Optimized Tissue Dissociation Protocol for Single-Cell RNA Sequencing Analysis of Fresh and Cultured Human Skin Biopsies.

We present an optimized dissociation protocol for preparing high-quality skin cell suspensions for in-depth single-cell RNA-sequencing (scRNA-seq) analysis of fresh and cultured human skin. Our protocol enabled the isolation of a consistently high number of highly viable skin cells from small freshly dissociated punch skin biopsies, which we use for scRNA-seq studies. We recapitulated not only the main cell populations of existing single-cell skin atlases, but also identified rare cell populations, such as mast cells. Furthermore, we effectively isolated highly viable single cells from ex vivo cultured skin biopsy fragments and generated a global single-cell map of the explanted human skin. The quality metrics of the generated scRNA-seq datasets were comparable between freshly dissociated and cultured skin. Overall, by enabling efficient cell isolation and comprehensive cell mapping, our skin dissociation-scRNA-seq workflow can greatly facilitate scRNA-seq discoveries across diverse human skin pathologies and ex vivo skin explant experimentations.

Phase I Trial Characterizing the Pharmacokinetic Profile of N-803, a Chimeric IL-15 Superagonist, in Healthy Volunteers.

The oncotherapeutic promise of IL-15, a potent immunostimulant, is limited by a short serum t 1/2 The fusion protein N-803 is a chimeric IL-15 superagonist that has a >20-fold longer in vivo t 1/2 versus IL-15. This phase 1 study characterized the pharmacokinetic (PK) profile and safety of N-803 after s.c. administration to healthy human volunteers. Volunteers received two doses of N-803, and after each dose, PK and safety were assessed for 9 d. The primary endpoint was the N-803 PK profile, the secondary endpoint was safety, and immune cell levels and immunogenicity were measures of interest. Serum N-803 concentrations peaked 4 h after administration and declined with a t 1/2 of ∼20 h. N-803 did not cause treatment-emergent serious adverse events (AEs) or grade ≥3 AEs. Injection site reactions, chills, and pyrexia were the most common AEs. Administration of N-803 was well tolerated and accompanied by proliferation of NK cells and CD8+ T cells and sustained increases in the number of NK cells. Our results suggest that N-803 administration can potentiate antitumor immunity.

Hand2 delineates mesothelium progenitors and is reactivated in mesothelioma.

The mesothelium lines body cavities and surrounds internal organs, widely contributing to homeostasis and regeneration. Mesothelium disruptions cause visceral anomalies and mesothelioma tumors. Nonetheless, the embryonic emergence of mesothelia remains incompletely understood. Here, we track mesothelial origins in the lateral plate mesoderm (LPM) using zebrafish. Single-cell transcriptomics uncovers a post-gastrulation gene expression signature centered on hand2 in distinct LPM progenitor cells. We map mesothelial progenitors to lateral-most, hand2-expressing LPM and confirm conservation in mouse. Time-lapse imaging of zebrafish hand2 reporter embryos captures mesothelium formation including pericardium, visceral, and parietal peritoneum. We find primordial germ cells migrate with the forming mesothelium as ventral migration boundary. Functionally, hand2 loss disrupts mesothelium formation with reduced progenitor cells and perturbed migration. In mouse and human mesothelioma, we document expression of LPM-associated transcription factors including Hand2, suggesting re-initiation of a developmental program. Our data connects mesothelium development to Hand2, expanding our understanding of mesothelial pathologies.

Mass cytometric and transcriptomic profiling of epithelial-mesenchymal transitions in human mammary cell lines.

Epithelial-mesenchymal transition (EMT) equips breast cancer cells for metastasis and treatment resistance. However, detection, inhibition, and elimination of EMT-undergoing cells is challenging due to the intrinsic heterogeneity of cancer cells and the phenotypic diversity of EMT programs. We comprehensively profiled EMT transition phenotypes in four non-cancerous human mammary epithelial cell lines using a flow cytometry surface marker screen, RNA sequencing, and mass cytometry. EMT was induced in the HMLE and MCF10A cell lines and in the HMLE-Twist-ER and HMLE-Snail-ER cell lines by prolonged exposure to TGFß1 or 4-hydroxytamoxifen, respectively. Each cell line exhibited a spectrum of EMT transition phenotypes, which we compared to the steady-state phenotypes of fifteen luminal, HER2-positive, and basal breast cancer cell lines. Our data provide multiparametric insights at single-cell level into the phenotypic diversity of EMT at different time points and in four human cellular models. These insights are valuable to better understand the complexity of EMT, to compare EMT transitions between the cellular models used here, and for the design of EMT time course experiments.

Disentangling tumorigenesis-associated DNA methylation changes in colorectal tissues from those associated with ageing.

Physiological ageing and tumorigenesis are both associated with epigenomic alterations in human tissue cells, the most extensively investigated of which entails de novo cytosine methylation (i.e., hypermethylation) within the CpG dinucleotides of CpG islands. Genomic regions that become hypermethylated during tumorigenesis are generally believed to overlap regions that acquire methylation in normal tissues as an effect of ageing. To define the extension of this overlap, we analysed the DNA methylomes of 48 large-bowel tissue samples taken from women of different ages during screening colonoscopy: 18 paired samples of normal and lesional tissues from donors harbouring a precancerous lesion and 12 samples of normal mucosa from tumour-free donors. Each sample was subjected to targeted, genome-wide bisulphite sequencing of ~2.5% of the genome, including all CpG islands. In terms of both its magnitude and extension along the chromatin, tumour-associated DNA hypermethylation in these regions was much more conspicuous than that observed in the normal mucosal samples from older (vs. younger) tumour-free donors. 83% of the ageing-associated hypermethylated regions (n = 2501) coincided with hypermethylated regions observed in tumour samples. However, 86% of the regions displaying hypermethylation in precancerous lesions (n = 16,772) showed no methylation changes in the ageing normal mucosa. The tumour-specificity of this latter hypermethylation was validated using published sets of data on DNA methylation in normal and neoplastic colon tissues. This extensive set of genomic regions displaying tumour-specific hypermethylation represents a rich vein of putative biomarkers for the early, non-invasive detection of colorectal tumours in women of all ages.