Psychrophilic proteases dramatically reduce single-cell RNA-seq artifacts: a molecular atlas of kidney development

A method for selective and efficient isolation of gray matter astrocytes from the spinal cord of adult mice

A growing body of evidence indicates intra- and inter-regional heterogeneity of astrocytes in the brain. However, because of a lack of an efficient method for isolating astrocytes from the spinal cord, little is known about how much spinal cord astrocytes are heterogeneous in adult mice. In this study, we developed a new method for isolating spinal astrocytes from adult mice using a cold-active protease from Bacillus licheniformis with an astrocyte cell surface antigen-2 (ACSA-2) antibody. Using fluorescence-activated cell sorting, isolated spinal ACSA-2+ cells were divided into two distinct populations, ACSA-2high and ACSA-2low. By analyzing the expression of cell-type marker genes, the ACSA-2high and ACSA-2low populations were identified as astrocytes and ependymal cells, respectively. Furthermore, ACSA-2high cells had mRNAs encoding genes that were abundantly expressed in the gray matter (GM) but not white matter astrocytes. By optimizing enzymatic isolation procedures, the yield of GM astrocytes also increased. Therefore, our newly established method enabled the selective and efficient isolation of GM astrocytes from the spinal cord of adult mice and may be useful for bulk- or single-cell RNA-sequencing under physiological and pathological conditions.

Single-cell transcriptomic analyses reveal distinct immune cell contributions to epithelial barrier dysfunction in checkpoint inhibitor colitis

Immune checkpoint inhibitor (ICI) therapy has revolutionized oncology, but treatments are limited by immune-related adverse events, including checkpoint inhibitor colitis (irColitis). Little is understood about the pathogenic mechanisms driving irColitis, which does not readily occur in model organisms, such as mice. To define molecular drivers of irColitis, we used single-cell multi-omics to profile approximately 300,000 cells from the colon mucosa and blood of 13 patients with cancer who developed irColitis (nine on anti-PD-1 or anti-CTLA-4 monotherapy and four on dual ICI therapy; most patients had skin or lung cancer), eight controls on ICI therapy and eight healthy controls. Patients with irColitis showed expanded mucosal Tregs, ITGAEHi CD8 tissue-resident memory T cells expressing CXCL13 and Th17 gene programs and recirculating ITGB2Hi CD8 T cells. Cytotoxic GNLYHi CD4 T cells, recirculating ITGB2Hi CD8 T cells and endothelial cells expressing hypoxia gene programs were further expanded in colitis associated with anti-PD-1/CTLA-4 therapy compared to anti-PD-1 therapy. Luminal epithelial cells in patients with irColitis expressed PCSK9, PD-L1 and interferon-induced signatures associated with apoptosis, increased cell turnover and malabsorption. Together, these data suggest roles for circulating T cells and epithelial-immune crosstalk critical to PD-1/CTLA-4-dependent tolerance and barrier function and identify potential therapeutic targets for irColitis.

Short-Term Disruption of TGFβ Signaling in Adult Mice Renders the Aorta Vulnerable to Hypertension-Induced Dissection

Hypertension and transient increases in blood pressure from extreme exertion are risk factors for aortic dissection in patients with age-related vascular degeneration or inherited connective tissue disorders. Yet, the common experimental model of angiotensin II-induced aortopathy in mice appears independent of high blood pressure as lesions do not occur in response to an alternative vasoconstrictor, norepinephrine, and are not prevented by co-treatment with a vasodilator, hydralazine. We investigated vasoconstrictor administration to adult mice 1 week after disruption of TGFβ signaling in smooth muscle cells. Norepinephrine increased blood pressure and induced aortic dissection by 7 days and even within 30 minutes that was rescued by hydralazine; results were similar with angiotensin II. Changes in regulatory contractile molecule expression were not of pathological significance. Rather, reduced synthesis of extracellular matrix yielded a vulnerable aortic phenotype by decreasing medial collagen, most dynamically type XVIII, and impairing cell-matrix adhesion. We conclude that transient and sustained increases in blood pressure cause dissection in aortas rendered vulnerable by inhibition of TGFβ-driven extracellular matrix production by smooth muscle cells. A corollary is that medial fibrosis, a frequent feature of medial degeneration, may afford some protection against aortic dissection.

Advancements in Single-Cell RNA Sequencing Research for Neurological Diseases

Neurological diseases are a major cause of the global burden of disease. Although the mechanisms of the occurrence and development of neurological diseases are not fully clear, most of them are associated with cells mediating neuroinflammation. Yet medications and other therapeutic options to improve treatment are still very limited. Single-cell RNA sequencing (scRNA-seq), as a delightfully potent breakthrough technology, not only identifies various cell types and response states but also uncovers cell-specific gene expression changes, gene regulatory networks, intercellular communication, and cellular movement trajectories, among others, in different cell types. In this review, we describe the technology of scRNA-seq in detail and discuss and summarize the application of scRNA-seq in exploring neurological diseases, elaborating the corresponding specific mechanisms of the diseases as well as providing a reliable basis for new therapeutic approaches. Finally, we affirm that scRNA-seq promotes the development of the neuroscience field and enables us to have a deeper cellular understanding of neurological diseases in the future, which provides strong support for the treatment of neurological diseases and the improvement of patients' prognosis.

scFSNN: a feature selection method based on neural network for single-cell RNA-seq data

While single-cell RNA sequencing (scRNA-seq) allows researchers to analyze gene expression in individual cells, its unique characteristics like over-dispersion, zero-inflation, high gene-gene correlation, and large data volume with many features pose challenges for most existing feature selection methods. In this paper, we present a feature selection method based on neural network (scFSNN) to solve classification problem for the scRNA-seq data. scFSNN is an embedded method that can automatically select features (genes) during model training, control the false discovery rate of selected features and adaptively determine the number of features to be eliminated. Extensive simulation and real data studies demonstrate its excellent feature selection ability and predictive performance.

Defining the T cell transcriptional landscape in pediatric liver transplant rejection at single cell resolution

Acute cellular rejection (ACR) affects >80% of pediatric liver transplant recipients within 5 years, and late ACR is associated with graft failure. Traditional anti-rejection therapy for late ACR is ineffective and has remained unchanged for six decades. Although CD8+ T cells promote late ACR, little has been done to define their specificity and gene expression. Here, we used single-cell sequencing and immune repertoire profiling (10X Genomics) on 30 cryopreserved 16G liver biopsies from 14 patients (5 pre-transplant or with no ACR, 9 with ACR). We identified expanded intragraft CD8+ T cell clonotypes (CD8EXP) and their gene expression profiles in response to anti-rejection treatment. Notably, we found that expanded CD8+ clonotypes (CD8EXP) bore markers of effector and CD56hiCD161- 'NK-like' T cells, retaining their clonotype identity and phenotype in subsequent biopsies from the same patients despite histologic ACR resolution. CD8EXP clonotypes localized to portal infiltrates during active ACR, and persisted in the lobule after histologic ACR resolution. CellPhoneDB analysis revealed differential crosstalk between KC and CD8EXP during late ACR, with activation of the LTB-LTBR pathway and downregulation of TGFß signaling. Therefore, persistently-detected intragraft CD8EXP clones remain active despite ACR treatment and may contribute to long-term allograft fibrosis and failure of operational tolerance.

Single-Cell Transcriptomic Profiling Identifies Molecular Phenotypes of Newborn Human Lung Cells

While animal model studies have extensively defined the mechanisms controlling cell diversity in the developing mammalian lung, there exists a significant knowledge gap with regards to late-stage human lung development. The NHLBI Molecular Atlas of Lung Development Program (LungMAP) seeks to fill this gap by creating a structural, cellular and molecular atlas of the human and mouse lung. Transcriptomic profiling at the single-cell level created a cellular atlas of newborn human lungs. Frozen single-cell isolates obtained from two newborn human lungs from the LungMAP Human Tissue Core Biorepository, were captured, and library preparation was completed on the Chromium 10X system. Data was analyzed in Seurat, and cellular annotation was performed using the ToppGene functional analysis tool. Transcriptional interrogation of 5500 newborn human lung cells identified distinct clusters representing multiple populations of epithelial, endothelial, fibroblasts, pericytes, smooth muscle, immune cells and their gene signatures. Computational integration of data from newborn human cells and with 32,000 cells from postnatal days 1 through 10 mouse lungs generated by the LungMAP Cincinnati Research Center facilitated the identification of distinct cellular lineages among all the major cell types. Integration of the newborn human and mouse cellular transcriptomes also demonstrated cell type-specific differences in maturation states of newborn human lung cells. Specifically, newborn human lung matrix fibroblasts could be separated into those representative of younger cells (n = 393), or older cells (n = 158). Cells with each molecular profile were spatially resolved within newborn human lung tissue. This is the first comprehensive molecular map of the cellular landscape of neonatal human lung, including biomarkers for cells at distinct states of maturity.

SOX9 switch links regeneration to fibrosis at the single-cell level in mammalian kidneys

The steps governing healing with or without fibrosis within the same microenvironment are unclear. After acute kidney injury (AKI), injured proximal tubular epithelial cells activate SOX9 for self-restoration. Using a multimodal approach for a head-to-head comparison of injury-induced SOX9 lineages, we identified a dynamic SOX9 switch in repairing epithelia. Lineages that regenerated epithelia silenced SOX9 and healed without fibrosis (SOX9on-off). By contrast, lineages with unrestored apicobasal polarity maintained SOX9 activity in sustained efforts to regenerate, which were identified as a SOX9on-on Cadherin6pos cell state. These reprogrammed cells generated substantial single-cell WNT activity to provoke a fibroproliferative response in adjacent fibroblasts, driving AKI to chronic kidney disease. Transplanted human kidneys displayed similar SOX9/CDH6/WNT2B responses. Thus, we have uncovered a sensor of epithelial repair status, the activity of which determines regeneration with or without fibrosis.

TMEM16A in smooth muscle cells acts as a pacemaker channel in the internal anal sphincter

Maintenance of fecal continence requires a continuous or basal tone of the internal anal sphincter (IAS). Paradoxically, the basal tone results largely from high-frequency rhythmic contractions of the IAS smooth muscle. However, the cellular and molecular mechanisms that initiate these contractions remain elusive. Here we show that the IAS contains multiple pacemakers. These pacemakers spontaneously generate propagating calcium waves that drive rhythmic contractions and establish the basal tone. These waves are myogenic and act independently of nerve, paracrine or autocrine signals. Using cell-specific gene knockout mice, we further found that TMEM16A Cl- channels in smooth muscle cells (but not in the interstitial cells of Cajal) are indispensable for pacemaking, rhythmic contractions, and basal tone. Our results identify TMEM16A in smooth muscle cells as a critical pacemaker channel that enables the IAS to contract rhythmically and continuously. This study provides cellular and molecular insights into fecal continence.

Embryonic Spinal Cord Innervation in Human Trunk Organogenesis Gastruloids: Cardiac Versus Enteric Customization and Beyond

Trunk-biased human gastruloids provide the ability to couple developmentally relevant spinal neurogenesis and organ morphogenesis via spatiotemporal self-organization events from derivatives of the three germ layers. The multi-lineage nature of gastruloids provides the full complexity of regulatory signaling cues that surpasses directed organoids and lays the foundation for an ex vivo self-evolving system. Here we detail two distinct protocols for trunk-biased gastruloids from an elongated, polarized structure with coordinated organ-specific neural patterning. Following an induction phase to caudalize iPSCs to trunk phenotype, divergent features of organogenesis and end-organ innervation yield separate models of enteric and cardiac nervous system formation. Both protocols are permissive to multi-lineage development and allow the study of neural integration events within a native, embryo-like context. We discuss the customizability of human gastruloids and the optimization of initial and extended conditions that maintain a permissive environment for multi-lineage differentiation and integration.

TFAP2 paralogs regulate midfacial development in part through a conserved ALX genetic pathway

Cranial neural crest development is governed by positional gene regulatory networks (GRNs). Fine-tuning of the GRN components underlies facial shape variation, yet how those networks in the midface are connected and activated remain poorly understood. Here, we show that concerted inactivation of Tfap2a and Tfap2b in the murine neural crest, even during the late migratory phase, results in a midfacial cleft and skeletal abnormalities. Bulk and single-cell RNA-seq profiling reveal that loss of both TFAP2 family members dysregulates numerous midface GRN components involved in midface morphogenesis, patterning and differentiation. Notably, Alx1, Alx3 and Alx4 (ALX) transcript levels are reduced, whereas ChIP-seq analyses suggest TFAP2 family members directly and positively regulate ALX gene expression. Tfap2a, Tfap2b and ALX co-expression in midfacial neural crest cells of both mouse and zebrafish implies conservation of this regulatory axis across vertebrates. Consistent with this notion, tfap2a zebrafish mutants present with abnormal alx3 expression patterns, Tfap2a binds ALX loci and tfap2a-alx3 genetic interactions are observed. Together, these data demonstrate TFAP2 paralogs regulate vertebrate midfacial development in part by activating expression of ALX transcription factor genes.

Comparative analysis of methods to reduce activation signature gene expression in PBMCs

Preserving the in vivo cell transcriptome is essential for accurate profiling, yet factors during cell isolation including time ex vivo and temperature induce artifactual gene expression, particularly in stress-responsive immune cells. In this study, we investigated two methods to mitigate ex vivo activation signature gene (ASG) expression in peripheral blood mononuclear cells (PBMCs): transcription and translation inhibitors (TTis) and cold temperatures during isolation. Comparative analysis of PBMCs isolated with TTis revealed reduced ASG expression. However, TTi treatment impaired responsiveness to LPS stimulation in subsequent in vitro experiments. In contrast, cold isolation methods also prevented ASG expression; up to a point where the addition of TTis during cold isolation offered minimal additional advantage. These findings highlight the importance of considering the advantages and drawbacks of different isolation methods to ensure accurate interpretation of PBMC transcriptomic profiles.

Protocol to dissociate epithelia from non-pregnant and pregnant mouse cervical tissue for single-cell RNA-sequencing

A challenge in studying cervical epithelial cell biology at the single-cell level is that differentiated subtypes, in particular mucus-secreting goblet cells, are sensitive to disassociating enzymes making isolation of all epithelial subpopulations difficult. Here we present a protocol to dissociate epithelia from non-pregnant and pregnant mouse cervical tissue for single-cell RNA-sequencing. We describe steps for harvesting cervices, preparing cervical tissue, dissociation of cervical cells, and viability checks. We then detail library preparation, sequencing, and procedure for data analysis. For complete details on the use and execution of this protocol, please refer to Cooley et al. (2023).1.

Prior viral infection primes cross-reactive CD8+ T cells that respond to mouse heart allografts

Introduction

Significant evidence suggests a connection between transplant rejection and the presence of high levels of pre-existing memory T cells. Viral infection can elicit viral-specific memory T cells that cross-react with allo-MHC capable of driving allograft rejection in mice. Despite these advances, and despite their critical role in transplant rejection, a systematic study of allo-reactive memory T cells, their specificities, and the role of cross-reactivity with viral antigens has not been performed.

Methods

Here, we established a model to identify, isolate, and characterize cross-reactive T cells using Nur77 reporter mice (C57BL/6 background), which transiently express GFP exclusively upon TCR engagement. We infected Nur77 mice with lymphocytic choriomeningitis virus (LCMV-Armstrong) to generate a robust memory compartment, where quiescent LCMV-specific memory CD8+ T cells could be readily tracked with MHC tetramer staining. Then, we transplanted LCMV immune mice with allogeneic hearts and monitored expression of GFP within MHC-tetramer defined viral-specific T cells as an indicator of their ability to cross-react with alloantigens.

Results

Strikingly, prior LCMV infection significantly increased the kinetics and magnitude of rejection as well as CD8+ T cell recruitment into allogeneic, but not syngeneic, transplanted hearts, relative to non-infected controls. Interestingly, as early as day 1 after allogeneic heart transplant an average of ~8% of MHC-tetramer+ CD8+ T cells expressed GFP, in contrast to syngeneic heart transplants, where the frequency of viral-specific CD8+ T cells that were GFP+ was <1%. These data show that a significant percentage of viral-specific memory CD8+ T cells expressed T cell receptors that also recognized alloantigens in vivo. Notably, the frequency of cross-reactive CD8+ T cells differed depending upon the viral epitope. Further, TCR sequences derived from cross-reactive T cells harbored distinctive motifs that may provide insight into cross-reactivity and allo-specificity.

Discussion

In sum, we have established a mouse model to track viral-specific, allo-specific, and cross-reactive T cells; revealing that prior infection elicits substantial numbers of viral-specific T cells that cross-react to alloantigen, respond very early after transplant, and may promote rapid rejection.

Hedgehog signalling in Foxd1+ embryonic kidney stromal progenitors controls nephron formation via Cxcl12 and Wnt5a

Congenital anomalies of the kidney and urinary tract (CAKUT) are characterised by a spectrum of structural and histologic abnormalities and are the major cause of childhood kidney failure. During kidney morphogenesis, the formation of a critical number of nephrons is an embryonic process supported, in part, by signalling between nephrogenic precursors and Foxd1-positive stromal progenitor cells. Low nephron number and abnormal patterning of the stroma are signature pathological features among CAKUT phenotypes with decreased kidney function. Despite their critical contribution to CAKUT pathogenesis, the mechanisms that underlie a low nephron number and the functional contribution of a disorganised renal stroma to nephron number are both poorly defined. Here, we identify a primary pathogenic role for increased Hedgehog signalling in embryonic renal stroma in the genesis of congenital low nephron number. Pharmacologic activation of Hedgehog (Hh) signalling in human kidney organoid tissue decreased the number of nephrons and generated excess stroma. The mechanisms underlying these pathogenic effects were delineated in genetic mouse models in which Hh signalling was constitutively activated in a cell lineage-specific manner. Cre-mediated excision of Ptch1 in Foxd1+ stromal progenitor cells, but not in Six2+ nephrogenic precursor cells, generated kidney malformation, identifying the stroma as a driver of low nephron number. Single-cell RNA sequencing analysis identified Cxcl12 and Wnt5a as downstream targets of increased stromal Hh signalling, findings supported by analysis in human kidney organoids. In vivo deficiency of Cxcl12 or Wnt5a in mice with increased stromal Hh signalling improved nephron endowment. These results demonstrate that dysregulated Hh signalling in embryonic renal stromal cells inhibits nephron formation in a manner dependent on Cxcl12 and Wnt5a. © 2023 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Ezh2-dependent methylation in oral epithelia promotes secondary palatogenesis

BACKGROUND

In addition to genomic risk variants and environmental influences, increasing evidence suggests epigenetic modifications are important for orofacial development and their alterations can contribute to orofacial clefts. Ezh2 encodes a core catalytic component of the Polycomb repressive complex responsible for addition of methyl marks to Histone H3 as a mechanism of repressing target genes. The role of Ezh2 in orofacial clefts remains unknown.

AIMS

To investigate the epithelial role of Ezh2-dependent methylation in secondary palatogenesis.

METHODS

We used conditional gene-targeting methods to ablate Ezh2 in the surface ectoderm-derived oral epithelium of mouse embryos. We then performed single-cell RNA sequencing combined with immunofluorescence and RT-qPCR to investigate gene expression in conditional mutant palate. We also employed double knockout analyses of Ezh1 and Ezh2 to address if they have synergistic roles in palatogenesis.

RESULTS

We found that conditional inactivation of Ezh2 in oral epithelia results in partially penetrant cleft palate. Double knockout analyses revealed that another family member Ezh1 is dispensable in orofacial development, and it does not have synergistic roles with Ezh2 in palatogenesis. Histochemistry and single-cell RNA-seq analyses revealed dysregulation of cell cycle regulators in the palatal epithelia of Ezh2 mutant mouse embryos disrupts palatogenesis.

CONCLUSION

Ezh2-dependent histone H3K27 methylation represses expression of cell cycle regulator Cdkn1a and promotes proliferation in the epithelium of the developing palatal shelves. Loss of this regulation may perturb movement of the palatal shelves, causing a delay in palate elevation which may result in failure of the secondary palate to close altogether.

Single-Nucleus RNA-Sequencing in Brain Tissue

Neuroscience research greatly benefits from single-cell sequencing technologies, which can reveal transcriptional alterations on a cellular level. However, preparing single-cell suspensions is technically challenging, requires experience, and has several limitations that can influence the transcriptional readout. Performing sequencing of single nuclei instead of single cells alleviates several of the challenges of sample preparation and highlights acute nuclear transcription. Here, we provide a protocol to prepare a nuclei suspension for single-nucleus RNA-sequencing for cell type-specific transcriptional profiling of brain tissue using the 10x Genomics single-cell gene expression assay. Furthermore, we highlight important aspects to consider during experimental design and data analysis. © 2023 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: Preparation of single-nucleus suspension Basic Protocol 2: Preparation and sequencing of single-nucleus libraries for RNA-seq.

Pax Protein Depletion in Proximal Tubules Triggers Conserved Mechanisms of Resistance to Acute Ischemic Kidney Injury and Prevents Transition to Chronic Kidney Disease

Acute kidney injury (AKI) is a common condition that lacks effective treatments. In part this shortcoming is due to an incomplete understanding of the genetic mechanisms that control pathogenesis and recovery. Pax2 and Pax8 are homologous transcription factors with overlapping functions that are critical for kidney development and are re-activated in AKI. In this report, we examined the role of Pax2 and Pax8 in recovery from ischemic AKI. We found that Pax2 and Pax8 are upregulated after severe AKI and correlate with chronic injury. Surprisingly, we then discovered that proximal-tubule-selective deletion of Pax2 and Pax8 resulted in a less severe chronic injury phenotype. This effect was mediated by protection against the acute insult, similar to preconditioning. Prior to injury, Pax2 and Pax8 mutant mice develop a unique subpopulation of S3 proximal tubule cells that display features usually seen only in acute or chronic injury. The expression signature of these cells was strongly enriched with genes associated with other mechanisms of protection against ischemic AKI including caloric restriction, hypoxic preconditioning, and female sex. Taken together, our results identify a novel role for Pax2 and Pax8 in mature proximal tubules that regulates critical genes and pathways involved in both injury response and protection from ischemic AKI.

TRANSLATIONAL STATEMENT

Identifying the molecular and genetic regulators unique to the nephron that dictate vulnerability to injury and regenerative potential could lead to new therapeutic targets to treat ischemic kidney injury. Pax2 and Pax8 are two homologous nephron-specific transcription factors that are critical for kidney development and physiology. Here we report that proximal-tubule-selective depletion of Pax2 and Pax8 protects against both acute and chronic injury and induces an expression profile in the S3 proximal tubule with common features shared among diverse conditions that protect against ischemia. These findings highlight a new role for Pax proteins as potential therapeutic targets to treat AKI.

Smarca4 deficiency induces Pttg1 oncogene upregulation and hyperproliferation of tubular and interstitial cells during kidney development

Kidney formation and nephrogenesis are controlled by precise spatiotemporal gene expression programs, which are coordinately regulated by cell-cycle, cell type-specific transcription factors and epigenetic/chromatin regulators. However, the roles of epigenetic/chromatin regulators in kidney development and disease remain poorly understood. In this study, we investigated the impact of deleting the chromatin remodeling factor Smarca4 (Brg1), a human Wilms tumor-associated gene, in Wnt4-expressing cells. Smarca4 deficiency led to severe tubular defects and a shortened medulla. Through unbiased single-cell RNA sequencing analyses, we identified multiple types of Wnt4 Cre-labeled interstitial cells, along with nephron-related cells. Smarca4 deficiency increased interstitial cells but markedly reduced tubular cells, resulting in cells with mixed identity and elevated expression of cell-cycle regulators and genes associated with extracellular matrix and epithelial-to-mesenchymal transition/fibrosis. We found that Smarca4 loss induced a significant upregulation of the oncogene Pttg1 and hyperproliferation of Wnt4 Cre-labeled cells. These changes in the cellular state could hinder the cellular transition into characteristic tubular structures, eventually leading to fibrosis. In conclusion, our findings shed light on novel cell types and genes associated with Wnt4 Cre-labeled cells and highlight the critical role of Smarca4 in regulating tubular cell differentiation and the expression of the cancer-causing gene Pttg1 in the kidney. These findings may provide valuable insights into potential therapeutic strategies for renal cell carcinoma resulting from SMARCA4 deficiency.

scAce: an adaptive embedding and clustering method for single-cell gene expression data

MOTIVATION

Since the development of single-cell RNA sequencing (scRNA-seq) technologies, clustering analysis of single-cell gene expression data has been an essential tool for distinguishing cell types and identifying novel cell types. Even though many methods have been available for scRNA-seq clustering analysis, the majority of them are constrained by the requirement on predetermined cluster numbers or the dependence on selected initial cluster assignment.

RESULTS

In this article, we propose an adaptive embedding and clustering method named scAce, which constructs a variational autoencoder to simultaneously learn cell embeddings and cluster assignments. In the scAce method, we develop an adaptive cluster merging approach which achieves improved clustering results without the need to estimate the number of clusters in advance. In addition, scAce provides an option to perform clustering enhancement, which can update and enhance cluster assignments based on previous clustering results from other methods. Based on computational analysis of both simulated and real datasets, we demonstrate that scAce outperforms state-of-the-art clustering methods for scRNA-seq data, and achieves better clustering accuracy and robustness.

AVAILABILITY AND IMPLEMENTATION

The scAce package is implemented in python 3.8 and is freely available from https://github.com/sldyns/scAce.

Single cell genomics as a transformative approach for aquaculture research and innovation

Single cell genomics encompasses a suite of rapidly maturing technologies that measure the molecular profiles of individual cells within target samples. These approaches provide a large up-step in biological information compared to long-established 'bulk' methods that profile the average molecular profiles of all cells in a sample, and have led to transformative advances in understanding of cellular biology, particularly in humans and model organisms. The application of single cell genomics is fast expanding to non-model taxa, including aquaculture species, where numerous research applications are underway with many more envisaged. In this review, we highlight the potential transformative applications of single cell genomics in aquaculture research, considering barriers and potential solutions to the broad uptake of these technologies. Focusing on single cell transcriptomics, we outline considerations for experimental design, including the essential requirement to obtain high quality cells/nuclei for sequencing in ectothermic aquatic species. We further outline data analysis and bioinformatics considerations, tailored to studies with the under-characterized genomes of aquaculture species, where our knowledge of cellular heterogeneity and cell marker genes is immature. Overall, this review offers a useful source of knowledge for researchers aiming to apply single cell genomics to address biological challenges faced by the global aquaculture sector though an improved understanding of cell biology.

Spinal cord injury regulates circular RNA expression in axons

Introduction

Neurons transport mRNA and translational machinery to axons for local translation. After spinal cord injury (SCI), de novo translation is assumed to enable neurorepair. Knowledge of the identity of axonal mRNAs that participate in neurorepair after SCI is limited. We sought to identify and understand how axonal RNAs play a role in axonal regeneration.

Methods

We obtained preparations enriched in axonal mRNAs from control and SCI rats by digesting spinal cord tissue with cold-active protease (CAP). The digested samples were then centrifuged to obtain a supernatant that was used to identify mRNA expression. We identified differentially expressed genes (DEGS) after SCI and mapped them to various biological processes. We validated the DEGs by RT-qPCR and RNA-scope.

Results

The supernatant fraction was highly enriched for mRNA from axons. Using Gene Ontology, the second most significant pathway for all DEGs was axonogenesis. Among the DEGs was Rims2, which is predominately a circular RNA (circRNA) in the CNS. We show that Rims2 RNA within spinal cord axons is circular. We found an additional 200 putative circRNAs in the axonal-enriched fraction. Knockdown in primary rat cortical neurons of the RNA editing enzyme ADAR1, which inhibits formation of circRNAs, significantly increased axonal outgrowth and increased the expression of circRims2. Using Rims2 as a prototype we used Circular RNA Interactome to predict miRNAs that bind to circRims2 also bind to the 3'UTR of GAP-43, PTEN or CREB1, all known regulators of axonal outgrowth. Axonally-translated GAP-43 supports axonal elongation and we detect GAP-43 mRNA in the rat axons by RNAscope.

Discussion

By enriching for axonal RNA, we detect SCI induced DEGs, including circRNA such as Rims2. Ablation of ADAR1, the enzyme that regulates circRNA formation, promotes axonal outgrowth of cortical neurons. We developed a pathway model using Circular RNA Interactome that indicates that Rims2 through miRNAs can regulate the axonal translation GAP-43 to regulate axonal regeneration. We conclude that axonal regulatory pathways will play a role in neurorepair.

Opportunities and challenges in the application of single-cell and spatial transcriptomics in plants

Single-cell and spatial transcriptomics have diverted researchers' attention from the multicellular level to the single-cell level and spatial information. Single-cell transcriptomes provide insights into the transcriptome at the single-cell level, whereas spatial transcriptomes help preserve spatial information. Although these two omics technologies are helpful and mature, further research is needed to ensure their widespread applicability in plant studies. Reviewing recent research on plant single-cell or spatial transcriptomics, we compared the different experimental methods used in various plants. The limitations and challenges are clear for both single-cell and spatial transcriptomic analyses, such as the lack of applicability, spatial information, or high resolution. Subsequently, we put forth further applications, such as cross-species analysis of roots at the single-cell level and the idea that single-cell transcriptome analysis needs to be combined with other omics analyses to achieve superiority over individual omics analyses. Overall, the results of this review suggest that combining single-cell transcriptomics, spatial transcriptomics, and spatial element distribution can provide a promising research direction, particularly for plant research.

Dynamic thresholding and tissue dissociation optimization for CITE-seq identifies differential surface protein abundance in metastatic melanoma

Multi-omics profiling by CITE-seq bridges the RNA-protein gap in single-cell analysis but has been largely applied to liquid biopsies. Applying CITE-seq to clinically relevant solid biopsies to characterize healthy tissue and the tumor microenvironment is an essential next step in single-cell translational studies. In this study, gating of cell populations based on their transcriptome signatures for use in cell type-specific ridge plots allowed identification of positive antibody signals and setting of manual thresholds. Next, we compare five skin dissociation protocols by taking into account dissociation efficiency, captured cell type heterogeneity and recovered surface proteome. To assess the effect of enzymatic digestion on transcriptome and epitope expression in immune cell populations, we analyze peripheral blood mononuclear cells (PBMCs) with and without dissociation. To further assess the RNA-protein gap, RNA-protein we perform codetection and correlation analyses on thresholded protein values. Finally, in a proof-of-concept study, using protein abundance analysis on selected surface markers in a cohort of healthy skin, primary, and metastatic melanoma we identify CD56 surface marker expression on metastatic melanoma cells, which was further confirmed by multiplex immunohistochemistry. This work provides practical guidelines for processing and analysis of clinically relevant solid tissue biopsies for biomarker discovery.

A contamination focused approach for optimizing the single-cell RNA-seq experiment

Droplet-based single-cell RNA-seq (scRNA-seq) data are plagued by ambient contaminations caused by nucleic acid material released by dead and dying cells. This material is mixed into the buffer and is co-encapsulated with cells, leading to a lower signal-to-noise ratio. Although there exist computational methods to remove ambient contaminations post-hoc, the reliability of algorithms in generating high-quality data from low-quality sources remains uncertain. Here, we assess data quality before data filtering by a set of quantitative, contamination-based metrics that assess data quality more effectively than standard metrics. Through a series of controlled experiments, we report improvements that can minimize ambient contamination outside of tissue dissociation, via cell fixation, improved cell loading, microfluidic dilution, and nuclei versus cell preparation; many of these parameters are inaccessible on commercial platforms. We provide end-users with insights on factors that can guide their decision-making regarding optimizations that minimize ambient contamination, and metrics to assess data quality.

Liver in infections: a single-cell and spatial transcriptomics perspective

The liver is an immune organ that plays a vital role in the detection, capture, and clearance of pathogens and foreign antigens that invade the human body. During acute and chronic infections, the liver transforms from a tolerant to an active immune state. The defence mechanism of the liver mainly depends on a complicated network of intrahepatic and translocated immune cells and non-immune cells. Therefore, a comprehensive liver cell atlas in both healthy and diseased states is needed for new therapeutic target development and disease intervention improvement. With the development of high-throughput single-cell technology, we can now decipher heterogeneity, differentiation, and intercellular communication at the single-cell level in sophisticated organs and complicated diseases. In this concise review, we aimed to summarise the advancement of emerging high-throughput single-cell technologies and re-define our understanding of liver function towards infections, including hepatitis B virus, hepatitis C virus, Plasmodium, schistosomiasis, endotoxemia, and corona virus disease 2019 (COVID-19). We also unravel previously unknown pathogenic pathways and disease mechanisms for the development of new therapeutic targets. As high-throughput single-cell technologies mature, their integration into spatial transcriptomics, multiomics, and clinical data analysis will aid in patient stratification and in developing effective treatment plans for patients with or without liver injury due to infectious diseases.

TFAP2 paralogs regulate midfacial development in part through a conserved ALX genetic pathway

Cranial neural crest development is governed by positional gene regulatory networks (GRNs). Fine-tuning of the GRN components underly facial shape variation, yet how those in the midface are connected and activated remain poorly understood. Here, we show that concerted inactivation of Tfap2a and Tfap2b in the murine neural crest even during the late migratory phase results in a midfacial cleft and skeletal abnormalities. Bulk and single-cell RNA-seq profiling reveal that loss of both Tfap2 members dysregulated numerous midface GRN components involved in midface fusion, patterning, and differentiation. Notably, Alx1/3/4 (Alx) transcript levels are reduced, while ChIP-seq analyses suggest TFAP2 directly and positively regulates Alx gene expression. TFAP2 and ALX co-expression in midfacial neural crest cells of both mouse and zebrafish further implies conservation of this regulatory axis across vertebrates. Consistent with this notion, tfap2a mutant zebrafish present abnormal alx3 expression patterns, and the two genes display a genetic interaction in this species. Together, these data demonstrate a critical role for TFAP2 in regulating vertebrate midfacial development in part through ALX transcription factor gene expression.

Dissociation protocols used for sarcoma tissues bias the transcriptome observed in single-cell and single-nucleus RNA sequencing

BACKGROUND

Single-cell RNA-seq has emerged as an innovative technology used to study complex tissues and characterize cell types, states, and lineages at a single-cell level. Classification of bulk tumors by their individual cellular constituents has also created new opportunities to generate single-cell atlases for many organs, cancers, and developmental models. Despite the tremendous promise of this technology, recent evidence studying epithelial tissues and diverse carcinomas suggests the methods used for tissue processing, cell disaggregation, and preservation can significantly bias gene expression and alter the observed cell types. To determine whether sarcomas - tumors of mesenchymal origin - are subject to the same technical artifacts, we profiled patient-derived tumor explants (PDXs) propagated from three aggressive subtypes: osteosarcoma (OS), Ewing sarcoma (ES), desmoplastic small round cell tumor (DSRCT). Given the rarity of these sarcoma subtypes, we explored whether single-nuclei RNA-seq from more widely available archival frozen specimens could accurately be identified by gene expression signatures linked to tissue phenotype or pathognomonic fusion proteins.

RESULTS

We systematically assessed dissociation methods across different sarcoma subtypes. We compared gene expression from single-cell and single-nucleus RNA-sequencing of 125,831 whole-cells and nuclei from ES, DSRCT, and OS PDXs. We detected warm dissociation artifacts in single-cell samples and gene length bias in single-nucleus samples. Classic sarcoma gene signatures were observed regardless of the dissociation method. In addition, we showed that dissociation method biases could be computationally corrected.

CONCLUSIONS

We highlighted transcriptional biases, including warm dissociation and gene-length biases, introduced by the dissociation method for various sarcoma subtypes. This work is the first to characterize how the dissociation methods used for sc/snRNA-seq may affect the interpretation of the molecular features in sarcoma PDXs.

Single cell transcriptomic analysis of HPV16-infected epithelium identifies a keratinocyte subpopulation implicated in cancer

Persistent HPV16 infection is a major cause of the global cancer burden. The viral life cycle is dependent on the differentiation program of stratified squamous epithelium, but the landscape of keratinocyte subpopulations which support distinct phases of the viral life cycle has yet to be elucidated. Here, single cell RNA sequencing of HPV16 infected compared to uninfected organoids identifies twelve distinct keratinocyte populations, with a subset mapped to reconstruct their respective 3D geography in stratified squamous epithelium. Instead of conventional terminally differentiated cells, an HPV-reprogrammed keratinocyte subpopulation (HIDDEN cells) forms the surface compartment and requires overexpression of the ELF3/ESE-1 transcription factor. HIDDEN cells are detected throughout stages of human carcinogenesis including primary human cervical intraepithelial neoplasias and HPV positive head and neck cancers, and a possible role in promoting viral carcinogenesis is supported by TCGA analyses. Single cell transcriptome information on HPV-infected versus uninfected epithelium will enable broader studies of the role of individual keratinocyte subpopulations in tumor virus infection and cancer evolution.

Multi-omic single cell sequencing: Overview and opportunities for kidney disease therapeutic development

Single cell sequencing technologies have rapidly advanced in the last decade and are increasingly applied to gain unprecedented insights by deconstructing complex biology to its fundamental unit, the individual cell. First developed for measurement of gene expression, single cell sequencing approaches have evolved to allow simultaneous profiling of multiple additional features, including chromatin accessibility within the nucleus and protein expression at the cell surface. These multi-omic approaches can now further be applied to cells in situ, capturing the spatial context within which their biology occurs. To extract insights from these complex datasets, new computational tools have facilitated the integration of information across different data types and the use of machine learning approaches. Here, we summarize current experimental and computational methods for generation and integration of single cell multi-omic datasets. We focus on opportunities for multi-omic single cell sequencing to augment therapeutic development for kidney disease, including applications for biomarkers, disease stratification and target identification.

B cell-dependent subtypes and treatment-based immune correlates to survival in stage 3 and 4 lung adenocarcinomas

Lung cancer is the leading cause of cancer-related deaths worldwide. Surgery and chemoradiation are the standard of care in early stages of non-small cell lung cancer (NSCLC), while immunotherapy is the standard of care in late-stage NSCLC. The immune composition of the tumor microenvironment (TME) is recognized as an indicator for responsiveness to immunotherapy, although much remains unknown about its role in responsiveness to surgery or chemoradiation. In this pilot study, we characterized the NSCLC TME using mass cytometry (CyTOF) and bulk RNA sequencing (RNA-Seq) with deconvolution of RNA-Seq being performed by Kassandra, a recently published deconvolution tool. Stratification of patients based on the intratumoral abundance of B cells identified that the B-cell rich patient group had increased expression of CXCL13 and greater abundance of PD1+ CD8 T cells. The presence of B cells and PD1+ CD8 T cells correlated positively with the presence of intratumoral tertiary lymphoid structures (TLS). We then assessed the predictive and prognostic utility of these cell types and TLS within publicly available stage 3 and 4 lung adenocarcinoma (LUAD) RNA-Seq datasets. As previously described by others, pre-treatment expression of intratumoral 12-chemokine TLS gene signature is associated with progression free survival (PFS) in patients who receive treatment with immune checkpoint inhibitors (ICI). Notably and unexpectedly pre-treatment percentages of intratumoral B cells are associated with PFS in patients who receive surgery, chemotherapy, or radiation. Further studies to confirm these findings would allow for more effective patient selection for both ICI and non-ICI treatments.

Type I interferon responsive microglia shape cortical development and behavior

Microglia are brain resident phagocytes that can engulf synaptic components and extracellular matrix as well as whole neurons. However, whether there are unique molecular mechanisms that regulate these distinct phagocytic states is unknown. Here we define a molecularly distinct microglial subset whose function is to engulf neurons in the developing brain. We transcriptomically identified a cluster of Type I interferon (IFN-I) responsive microglia that expanded 20-fold in the postnatal day 5 somatosensory cortex after partial whisker deprivation, a stressor that accelerates neural circuit remodeling. In situ, IFN-I responsive microglia were highly phagocytic and actively engulfed whole neurons. Conditional deletion of IFN-I signaling (Ifnar1fl/fl) in microglia but not neurons resulted in dysmorphic microglia with stalled phagocytosis and an accumulation of neurons with double strand DNA breaks, a marker of cell stress. Conversely, exogenous IFN-I was sufficient to drive neuronal engulfment by microglia and restrict the accumulation of damaged neurons. IFN-I deficient mice had excess excitatory neurons in the developing somatosensory cortex as well as tactile hypersensitivity to whisker stimulation. These data define a molecular mechanism through which microglia engulf neurons during a critical window of brain development. More broadly, they reveal key homeostatic roles of a canonical antiviral signaling pathway in brain development.

Protocol to dissociate healthy and infected murine- and hamster-derived lung tissue for single-cell transcriptome analysis

In infectious disease research, single-cell RNA sequencing allows dissection of host-pathogen interactions. As a prerequisite, we provide a protocol to transform solid and complex organs such as lungs into representative diverse, viable single-cell suspensions. Our protocol describes performance of vascular perfusion, pneumonectomy, enzymatic digestion, and mechanical dissociation of lung tissue, as well as red blood cell lysis and counting of isolated cells. A challenge remains, however, to further increase the proportion of pulmonary endothelial cells without compromising on viability. For complete details on the use and execution of this protocol, please refer to Nouailles et al. (2021),1 Wyler et al. (2022),2 and Ebenig et al. (2022).3.

scGCL: an imputation method for scRNA-seq data based on graph contrastive learning

MOTIVATION

Single-cell RNA-sequencing (scRNA-seq) is widely used to reveal cellular heterogeneity, complex disease mechanisms and cell differentiation processes. Due to high sparsity and complex gene expression patterns, scRNA-seq data present a large number of dropout events, affecting downstream tasks such as cell clustering and pseudo-time analysis. Restoring the expression levels of genes is essential for reducing technical noise and facilitating downstream analysis. However, existing scRNA-seq data imputation methods ignore the topological structure information of scRNA-seq data and cannot comprehensively utilize the relationships between cells.

RESULTS

Here, we propose a single-cell Graph Contrastive Learning method for scRNA-seq data imputation, named scGCL, which integrates graph contrastive learning and Zero-inflated Negative Binomial (ZINB) distribution to estimate dropout values. scGCL summarizes global and local semantic information through contrastive learning and selects positive samples to enhance the representation of target nodes. To capture the global probability distribution, scGCL introduces an autoencoder based on the ZINB distribution, which reconstructs the scRNA-seq data based on the prior distribution. Through extensive experiments, we verify that scGCL outperforms existing state-of-the-art imputation methods in clustering performance and gene imputation on 14 scRNA-seq datasets. Further, we find that scGCL can enhance the expression patterns of specific genes in Alzheimer's disease datasets.

AVAILABILITY AND IMPLEMENTATION

The code and data of scGCL are available on Github: https://github.com/zehaoxiong123/scGCL.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Growth and differentiation of human induced pluripotent stem cell (hiPSC)-derived kidney organoids using fully synthetic peptide hydrogels

Human induced pluripotent stem cell (hiPSC)-derived kidney organoids have prospective applications ranging from basic disease modelling to personalised medicine. However, there remains a necessity to refine the biophysical and biochemical parameters that govern kidney organoid formation. Differentiation within fully-controllable and physiologically relevant 3D growth environments will be critical to improving organoid reproducibility and maturation. Here, we matured hiPSC-derived kidney organoids within fully synthetic self-assembling peptide hydrogels (SAPHs) of variable stiffness (storage modulus, G'). The resulting organoids contained complex structures comparable to those differentiated within the animal-derived matrix, Matrigel. Single-cell RNA sequencing (scRNA-seq) was then used to compare organoids matured within SAPHs to those grown within Matrigel or at the air-liquid interface. A total of 13,179 cells were analysed, revealing 14 distinct clusters. Organoid compositional analysis revealed a larger proportion of nephron cell types within Transwell-derived organoids, while SAPH-derived organoids were enriched for stromal-associated cell populations. Notably, differentiation within a higher G' SAPH generated podocytes with more mature gene expression profiles. Additionally, maturation within a 3D microenvironment significantly reduced the derivation of off-target cell types, which are a known limitation of current kidney organoid protocols. This work demonstrates the utility of synthetic peptide-based hydrogels with a defined stiffness, as a minimally complex microenvironment for the selected differentiation of kidney organoids.

Single cell transcriptomic analysis of renal allograft rejection reveals novel insights into intragraft TCR clonality

Bulk analysis of renal allograft biopsies (rBx) identified RNA transcripts associated with acute cellular rejection (ACR); however, these lacked cellular context critical to mechanistic understanding. We performed combined single cell RNA transcriptomic and TCRα/β sequencing on rBx from patients with ACR under differing immunosuppression (IS): tacrolimus, iscalimab, and belatacept. TCR analysis revealed a highly restricted CD8 ⁺ T cell clonal expansion (CD8 _EXP ), independent of HLA mismatch or IS type. Subcloning of TCRα/β cDNAs from CD8 _EXP into Jurkat76 cells (TCR ^-/- ) conferred alloreactivity by mixed lymphocyte reaction. scRNAseq analysis of CD8 _EXP revealed effector, memory, and exhausted phenotypes that were influenced by IS type. Successful anti-rejection treatment decreased, but did not eliminate, CD8 _EXP , while CD8 _EXP were maintained during treatment-refractory rejection. Finally, most rBx-derived CD8 _EXP were also observed in matching urine samples. Overall, our data define the clonal CD8 ⁺ T cell response to ACR, providing novel insights to improve detection, assessment, and treatment of rejection.

Gene length is a pivotal feature to explain disparities in transcript capture between single transcriptome techniques

The scale and capability of single-cell and single-nucleus RNA-sequencing technologies are rapidly growing, enabling key discoveries and large-scale cell mapping operations. However, studies directly comparing technical differences between single-cell and single-nucleus RNA sequencing are still lacking. Here, we compared three paired single-cell and single-nucleus transcriptomes from three different organs (Heart, Lung and Kidney). Differently from previous studies that focused on cell classification, we explored disparities in the transcriptome output of whole cells relative to the nucleus. We found that the major cell clusters could be recovered by either technique from matched samples, but at different proportions. In 2/3 datasets (kidney and lung) we detected clusters exclusively present with single-nucleus RNA sequencing. In all three organ groups, we found that genomic and gene structural characteristics such as gene length and exon content significantly differed between the two techniques. Genes recovered with the single-nucleus RNA sequencing technique had longer sequence lengths and larger exon counts, whereas single-cell RNA sequencing captured short genes at higher rates. Furthermore, we found that when compared to the whole host genome (mouse for kidney and lung datasets and human for the heart dataset), single transcriptomes obtained with either technique skewed from the expected proportions in several points: a) coding sequence length, b) transcript length and c) genomic span; and d) distribution of genes based on exons counts. Interestingly, the top-100 DEG between the two techniques returned distinctive GO terms. Hence, the type of single transcriptome technique used affected the outcome of downstream analysis. In summary, our data revealed both techniques present disparities in RNA capture. Moreover, the biased RNA capture affected the calculations of basic cellular parameters, raising pivotal points about the limitations and advantages of either single transcriptome techniques.

Preterm birth leads to a decreased number of differentiated podocytes and accelerated podocyte differentiation

Preterm birth was previously identified as a high-risk factor for the long-term development of chronic kidney disease. However, the detailed pattern of podocyte (PD) changes caused by preterm birth and the potential mechanism underlying this process have not been well clarified. In present study, a rat model of preterm birth was established by delivery of pups 2 days early and podometric methods were applied to identify the changes in PDs number caused by preterm birth. In addition, single-cell RNA sequencing (scRNA-seq) and subsequent bioinformatic analysis were performed in the preterm rat kidney to explore the possible mechanism caused by preterm birth. As results, when the kidney completely finished nephrogenesis at the age of 3 weeks, a reduction in the total number of differentiated PDs in kidney sections was detected. In addition, 20 distinct clusters and 12 different cell types were identified after scRNA-seq in preterm rats (postnatal day 2) and full-term rats (postnatal day 0). The numbers of PDs and most types of inherent kidney cells were decreased in the preterm birth model. In addition, 177 genes were upregulated while 82 genes were downregulated in the PDs of full-term rats compared with those of preterm rats. Further functional GO analysis revealed that ribosome-related genes were enriched in PDs from full-term rats, and kidney development-related genes were enriched in PDs from preterm rats. Moreover, known PD-specific and PD precursor genes were highly expressed in PDs from preterm rats, and pseudotemporal analysis showed that PDs were present earlier in preterm rats than in full-term rats. In conclusion, the present study showed that preterm birth could cause a reduction in the number of differentiated PDs and accelerate the differentiation of PDs.

Fishing for Developmental Regulatory Regions: Zebrafish Tissue-Specific ATAC-seq

Interactions between transcription factors and regulatory DNA can be described by gene regulatory networks. These networks provide a systems-level view of embryonic tissue development. Here, we describe a protocol for the isolation, identification, and experimental manipulation of tissue-specific cis-regulatory elements during zebrafish embryonic development using low-input ATAC-seq. With the methods described, genome-wide assessments of regulatory DNA in small populations of developing tissues can be identified, allowing for the construction of gene regulatory networks.

Single-cell multiomics reveals the complexity of TGFβ signalling to chromatin in iPSC-derived kidney organoids

TGFβ1 plays a regulatory role in the determination of renal cell fate and the progression of renal fibrosis. Here we show an association between SMAD3 and the histone methyltransferase, EZH2, during cell differentiation; ChIP-seq revealed that SMAD3 and EZH2 co-occupy the genome in iPSCs and in iPSC-derived nephron progenitors. Through integration of single cell gene expression and epigenome profiling, we identified de novo ACTA2^+ve/POSTN^+ve myofibroblasts in kidney organoids treated with TGFβ1, characterised by increased SMAD3-dependent cis chromatin accessibility and gene expression associated with fibroblast activation. We have identified fibrosis-associated regulons characterised by enrichment of SMAD3, AP1, the ETS family of transcription factors, and NUAK1, CREB3L1, and RARG, corresponding to enriched motifs at accessible loci identified by scATACseq. Treatment with the EZH2 specific inhibitor GSK343, blocked SMAD3-dependent cis co-accessibility and inhibited myofibroblast activation. This mechanism, through which TGFβ signals directly to chromatin, represents a critical determinant of fibrotic, differentiated states.

Characterization of alternative mRNA splicing in cultured cell populations representing progressive stages of human fetal kidney development

Nephrons are the functional units of the kidney. During kidney development, cells from the cap mesenchyme-a transient kidney-specific progenitor state-undergo a mesenchymal to epithelial transition (MET) and subsequently differentiate into the various epithelial cell types that create the tubular structures of the nephron. Faults in this transition can lead to a pediatric malignancy of the kidney called Wilms' tumor that mimics normal kidney development. While human kidney development has been characterized at the gene expression level, a comprehensive characterization of alternative splicing is lacking. Therefore, in this study, we performed RNA sequencing on cell populations representing early, intermediate, and late developmental stages of the human fetal kidney, as well as three blastemal-predominant Wilms' tumor patient-derived xenografts. Using this newly generated RNAseq data, we identified a set of transcripts that are alternatively spliced between the different developmental stages. Moreover, we found that cells from the earliest developmental stage have a mesenchymal splice-isoform profile that is similar to that of blastemal-predominant Wilms' tumor xenografts. RNA binding motif enrichment analysis suggests that the mRNA binding proteins ESRP1, ESRP2, RBFOX2, and QKI regulate alternative mRNA splicing during human kidney development. These findings illuminate new molecular mechanisms involved in human kidney development and pediatric kidney cancer.

Early Molecular Events Mediating Loss of Aquaporin-2 during Ureteral Obstruction in Rats

BACKGROUND

Ureteral obstruction is marked by disappearance of the vasopressin-dependent water channel aquaporin-2 (AQP2) in the renal collecting duct and polyuria upon reversal. Most studies of unilateral ureteral obstruction (UUO) models have examined late time points, obscuring the early signals that trigger loss of AQP2.

METHODS

We performed RNA-Seq on microdissected rat cortical collecting ducts (CCDs) to identify early signaling pathways after establishment of UUO.

RESULTS

Vasopressin V2 receptor (AVPR2) mRNA was decreased 3 hours after UUO, identifying one cause of AQP2 loss. Collecting duct principal cell differentiation markers were lost, including many not regulated by vasopressin. Immediate early genes in CCDs were widely induced 3 hours after UUO, including Myc, Atf3, and Fos (confirmed at the protein level). Simultaneously, expression of NF-κB signaling response genes known to repress Aqp2 increased. RNA-Seq for CCDs at an even earlier time point (30 minutes) showed widespread mRNA loss, indicating a "stunned" profile. Immunocytochemical labeling of markers of mRNA-degrading P-bodies DDX6 and 4E-T indicated an increase in P-body formation within 30 minutes.

CONCLUSIONS

Immediately after establishment of UUO, collecting ducts manifest a stunned state with broad disappearance of mRNAs. Within 3 hours, there is upregulation of immediate early and inflammatory genes and disappearance of the V2 vasopressin receptor, resulting in loss of AQP2 (confirmed by lipopolysaccharide administration). The inflammatory response seen rapidly after UUO establishment may be relevant to both UUO-induced polyuria and long-term development of fibrosis in UUO kidneys.

Single-cell RNA sequencing in atherosclerosis: Mechanism and precision medicine

Atherosclerosis is the pathological basis of various vascular diseases, including those with high mortality, such as myocardial infarction and stroke. However, its pathogenesis is complex and has not been fully elucidated yet. Over the past few years, single-cell RNA sequencing (scRNA-seq) has been developed and widely used in many biological fields to reveal biological mechanisms at the cellular level and solve the problems of cellular heterogeneity that cannot be solved using bulk RNA sequencing. In this review, we briefly summarize the existing scRNA-seq technologies and focus on their application in atherosclerosis research to provide insights into the occurrence, development and treatment of atherosclerosis.

Practical Considerations for Single-Cell Genomics

The single-cell revolution in the field of genomics is in full bloom, with clever new molecular biology tricks appearing regularly that allow researchers to explore new modalities or scale up their projects to millions of cells and beyond. Techniques abound to measure RNA expression, DNA alterations, protein abundance, chromatin accessibility, and more, all with single-cell resolution and often in combination. Despite such a rapidly changing technology landscape, there are several fundamental principles that are applicable to the majority of experimental workflows to help users avoid pitfalls and exploit the advantages of the chosen platform. In this overview article, we describe a variety of popular single-cell genomics technologies and address some common questions pertaining to study design, sample preparation, quality control, and sequencing strategy. As the majority of relevant publications currently revolve around single-cell RNA-seq, we will prioritize this genomics modality in our discussion. © 2022 Wiley Periodicals LLC.

Single-cell analysis of endometriosis reveals a coordinated transcriptional programme driving immunotolerance and angiogenesis across eutopic and ectopic tissues

Endometriosis is characterized by the growth of endometrial-like tissue outside the uterus. It affects many women during their reproductive age, causing years of pelvic pain and potential infertility. Its pathophysiology remains largely unknown, which limits early diagnosis and treatment. We characterized peritoneal and ovarian lesions at single-cell transcriptome resolution and compared them to matched eutopic endometrium, unaffected endometrium and organoids derived from these tissues, generating data on over 122,000 cells across 14 individuals. We spatially localized many of the cell types using imaging mass cytometry. We identify a perivascular mural cell specific to the peritoneal lesions, with dual roles in angiogenesis promotion and immune cell trafficking. We define an immunotolerant peritoneal niche, fundamental differences in eutopic endometrium and between lesion microenvironments and an unreported progenitor-like epithelial cell subpopulation. Altogether, this study provides a holistic view of the endometriosis microenvironment that represents a comprehensive cell atlas of the disease in individuals undergoing hormonal treatment, providing essential information for future therapeutics and diagnostics.

Single-cell transcriptomic signatures and gene regulatory networks modulated by Wls in mammalian midline facial formation and clefts

Formation of highly unique and complex facial structures is controlled by genetic programs that are responsible for the precise coordination of three-dimensional tissue morphogenesis. However, the underlying mechanisms governing these processes remain poorly understood. We combined mouse genetic and genomic approaches to define the mechanisms underlying normal and defective midfacial morphogenesis. Conditional inactivation of the Wnt secretion protein Wls in Pax3-expressing lineage cells disrupted frontonasal primordial patterning, cell survival and directional outgrowth, resulting in altered facial structures, including midfacial hypoplasia and midline facial clefts. Single-cell RNA sequencing revealed unique transcriptomic atlases of mesenchymal subpopulations in the midfacial primordia, which are disrupted in the conditional Wls mutants. Differentially expressed genes and cis-regulatory sequence analyses uncovered that Wls modulates and integrates a core gene regulatory network, consisting of key midfacial regulatory transcription factors (including Msx1, Pax3 and Pax7) and their downstream targets (including Wnt, Shh, Tgfβ and retinoic acid signaling components), in a mesenchymal subpopulation of the medial nasal prominences that is responsible for midline facial formation and fusion. These results reveal fundamental mechanisms underlying mammalian midfacial morphogenesis and related defects at single-cell resolution.

Tissue-resident memory CD8+ T cells possess unique transcriptional, epigenetic and functional adaptations to different tissue environments

Tissue-resident memory T cells (T_RM cells) provide protective immunity, but the contributions of specific tissue environments to T_RM cell differentiation and homeostasis are not well understood. In the present study, the diversity of gene expression and genome accessibility by mouse CD8⁺ T_RM cells from distinct organs that responded to viral infection revealed both shared and tissue-specific transcriptional and epigenetic signatures. T_RM cells in the intestine and salivary glands expressed transforming growth factor (TGF)-β-induced genes and were maintained by ongoing TGF-β signaling, whereas those in the fat, kidney and liver were not. Constructing transcriptional-regulatory networks identified the transcriptional repressor Hic1 as a critical regulator of T_RM cell differentiation in the small intestine and showed that Hic1 overexpression enhanced T_RM cell differentiation and protection from infection. Provision of a framework for understanding how CD8⁺ T_RM cells adapt to distinct tissue environments, and identification of tissue-specific transcriptional regulators mediating these adaptations, inform strategies to boost protective memory responses at sites most vulnerable to infection.

Single-Cell versus Single-Nucleus: Transcriptome differences in murine kidney after ischemia-reperfusion injury

The kidney is a complex organ, which consists of multiple components with highly diverse cell types. A detailed understanding of these cell types in health and disease is crucial for future development of preventive and curative treatment strategies. In recent years, single cell RNA sequencing (scRNAseq) and single nucleus RNA sequencing (snRNAseq) technology has opened up completely new possibilities in investigating the variety of renal cell populations in physiological and pathological states. Here, we systematically assess differences between scRNAseq and snRNAseq approaches in transcriptome analysis of murine kidneys after ischemia reperfusion injury. We included tissues from control kidneys and from kidneys harvested one week after mild (17 minutes clamping time) and severe (27 minutes clamping time) transient unilateral ischemia. Our findings reveal important methodological differences in the discovery of inflammatory cells, tubular cells, and other specialized cell types. While the scRNAseq approach is advantageous for investigating immune cells, the snRNAseq approach allows superior insight into healthy and damaged tubular cells. Apart from differences in the quantitative discovery rate, we found important qualitative discrepancies in the captured transcriptomes with crucial consequences for the interpretation of cell states and molecular functions. Together, we provide an overview of method-dependent differences between scRNAseq and snRNAseq results from identical post-ischemic kidney tissues. Our results highlight the importance of choosing the right approach for specific research questions.

Visual barcodes for clonal-multiplexing of live microscopy-based assays

While multiplexing samples using DNA barcoding revolutionized the pace of biomedical discovery, multiplexing of live imaging-based applications has been limited by the number of fluorescent proteins that can be deconvoluted using common microscopy equipment. To address this limitation, we develop visual barcodes that discriminate the clonal identity of single cells by different fluorescent proteins that are targeted to specific subcellular locations. We demonstrate that deconvolution of these barcodes is highly accurate and robust to many cellular perturbations. We then use visual barcodes to generate ‘Signalome’ cell-lines by mixing 12 clones of different live reporters into a single population, allowing simultaneous monitoring of the activity in 12 branches of signaling, at clonal resolution, over time. Using the ‘Signalome’ we identify two distinct clusters of signaling pathways that balance growth and proliferation, emphasizing the importance of growth homeostasis as a central organizing principle in cancer signaling. The ability to multiplex samples in live imaging applications, both in vitro and in vivo may allow better high-content characterization of complex biological systems.

Multiplex analyses of samples allow understanding complex processes in cancer initiation, progression and therapy response. Here, the authors present a fluorescence imaging-based visual barcode for livecell clonal-multiplexing which allows identifying signalling pathways clusters in response to different chemotherapy compounds.