Doublet identification in single-cell sequencing data using scDblFinder

Single-cell RNA sequencing of mouse lower respiratory tract epithelial cells: A meta-analysis

The respiratory system is a vital component of our body, essential for both oxygen uptake and immune defense. Knowledge of cellular composition and function in different parts of the respiratory tract provides the basis for a better understanding of the pathological processes involved in various diseases such as chronic respiratory diseases and cancer. Single-cell RNA sequencing (scRNA-seq) is a proficient approach for the identification and transcriptional characterization of cellular phenotypes. Although the mouse is an essential tool for the study of lung development, regeneration, and disease, a scRNA-seq mouse atlas of the lung in which all epithelial cell types are included and annotated systematically is lacking. Here, we established a single-cell transcriptome landscape of the mouse lower respiratory tract by performing a meta-analysis of seven different studies in which mouse lungs and trachea were analyzed by droplet and/or plate-based scRNA-seq technologies. We provide information on the best markers for each epithelial cell type, propose surface markers for the isolation of viable cells, harmonized the annotation of cell types, and compare the mouse single-cell transcriptomes with human scRNA-seq data of the lung.

Activation of Angiopoietin-Tie2 Signaling Protects the Kidney from Ischemic Injury by Modulation of Endothelial-Specific Pathways

SIGNIFICANCE STATEMENT

Ischemia-reperfusion AKI (IR-AKI) is common and causes significant morbidity. Effective treatments are lacking. However, preclinical studies suggest that inhibition of angiopoietin-Tie2 vascular signaling promotes injury, whereas activation of Tie2 is protective. We show that kidney ischemia leads to increased levels of the endothelial-specific phosphatase vascular endothelial protein tyrosine phosphatase (VE-PTP; PTPRB), which inactivates Tie2. Activation of Tie2 through VE-PTP deletion, or delivery of a novel angiopoietin mimetic (Hepta-ANG1), abrogated IR-AKI in mice. Single-cell RNAseq analysis showed Tie2 activation promotes increased Entpd1 expression, downregulation of FOXO1 target genes in the kidney vasculature, and emergence of a new subpopulation of glomerular endothelial cells. Our data provide a molecular basis and identify a candidate therapeutic to improve endothelial integrity and kidney function after IR-AKI.

BACKGROUND

Ischemia-reperfusion AKI (IR-AKI) is estimated to affect 2%-7% of all hospitalized patients. The significant morbidity and mortality associated with AKI indicates urgent need for effective treatments. Previous studies have shown activation of the vascular angiopoietin-Tie2 tyrosine kinase signaling pathway abrogates ischemia-reperfusion injury (IRI). We extended previous studies to (1) determine the molecular mechanism(s) underlying kidney injury and protection related to decreased or increased activation of Tie2, respectively, and (2) to test the hypothesis that deletion of the Tie2 inhibitory phosphatase vascular endothelial protein tyrosine phosphatase (VE-PTP) or injection of a new angiopoietin mimetic protects the kidney from IRI by common molecular mechanism(s).

METHODS

Bilateral IR-AKI was performed in VE-PTP wild-type or knockout mice and in C57BL/6J mice treated with Hepta-ANG1 or vehicle. Histologic, immunostaining, and single-cell RNA sequencing analyses were performed.

RESULTS

The phosphatase VE-PTP, which negatively regulates the angiopoietin-Tie2 pathway, was upregulated in kidney endothelial cells after IRI, and genetic deletion of VE-PTP in mice protected the kidney from IR-AKI. Injection of Hepta-ANG1 potently activated Tie2 and protected the mouse kidney from IRI. Single-cell RNAseq analysis of kidneys from Hepta-ANG1-treated and vehicle-treated mice identified endothelial-specific gene signatures and emergence of a new glomerular endothelial subpopulation associated with improved kidney function. Overlap was found between endothelial-specific genes upregulated by Hepta-ANG1 treatment and those downregulated in HUVECs with constitutive FOXO1 activation, including Entpd1 / ENTPD1 that modulates purinergic receptor signaling.

CONCLUSIONS

Our data support a key role of the endothelium in the development of IR-AKI, introduce Hepta-ANG1 as a putative new therapeutic biologic, and report a model to explain how IRI reduces Tie2 signaling and how Tie2 activation protects the kidney.

PODCAST

This article contains a podcast at https://dts.podtrac.com/redirect.mp3/www.asn-online.org/media/podcast/JASN/2023_05_23_JSN_Ang_EP23_052323.mp3.

Disrupted memory T cell expansion in HIV-exposed uninfected infants is preceded by premature skewing of T cell receptor clonality

While preventing vertical HIV transmission has been very successful, the increasing number of HIV-exposed uninfected infants (iHEU) experience an elevated risk to infections compared to HIV-unexposed and uninfected infants (iHUU). Immune developmental differences between iHEU and iHUU remains poorly understood and here we present a longitudinal multimodal analysis of infant immune ontogeny that highlights the impact of HIV/ARV exposure. Using mass cytometry, we show alterations and differences in the emergence of NK cell populations and T cell memory differentiation between iHEU and iHUU. Specific NK cells observed at birth were also predictive of acellular pertussis and rotavirus vaccine-induced IgG and IgA responses, respectively, at 3 and 9 months of life. T cell receptor Vβ clonotypic diversity was significantly and persistently lower in iHEU preceding the expansion of T cell memory. Our findings show that HIV/ARV exposure disrupts innate and adaptive immunity from birth which may underlie relative vulnerability to infections.

Molecular and cognitive signatures of ageing partially restored through synthetic delivery of IL2 to the brain

Cognitive decline is a common pathological outcome during aging, with an ill-defined molecular and cellular basis. In recent years, the concept of inflammaging, defined as a low-grade inflammation increasing with age, has emerged. Infiltrating T cells accumulate in the brain with age and may contribute to the amplification of inflammatory cascades and disruptions to the neurogenic niche observed with age. Recently, a small resident population of regulatory T cells has been identified in the brain, and the capacity of IL2-mediated expansion of this population to counter neuroinflammatory disease has been demonstrated. Here, we test a brain-specific IL2 delivery system for the prevention of neurological decline in aging mice. We identify the molecular hallmarks of aging in the brain glial compartments and identify partial restoration of this signature through IL2 treatment. At a behavioral level, brain IL2 delivery prevented the age-induced defect in spatial learning, without improving the general decline in motor skill or arousal. These results identify immune modulation as a potential path to preserving cognitive function for healthy aging.

Heterogeneity of Repolarization and Cell-Cell Variability of Cardiomyocyte Remodeling Within the Myocardial Infarction Border Zone Contribute to Arrhythmia Susceptibility

BACKGROUND

After myocardial infarction, the infarct border zone (BZ) is the dominant source of life-threatening arrhythmias, where fibrosis and abnormal repolarization create a substrate for reentry. We examined whether repolarization abnormalities are heterogeneous within the BZ in vivo and could be related to heterogeneous cardiomyocyte remodeling.

METHODS

Myocardial infarction was induced in domestic pigs by 120-minute ischemia followed by reperfusion. After 1 month, remodeling was assessed by magnetic resonance imaging, and electroanatomical mapping was performed to determine the spatial distribution of activation-recovery intervals. Cardiomyocytes were isolated and tissue samples collected from the BZ and remote regions. Optical recording allowed assessment of action potential duration (di-8-ANEPPS, stimulation at 1 Hz, 37 °C) of large cardiomyocyte populations while gene expression in cardiomyocytes was determined by single nuclear RNA sequencing.

RESULTS

In vivo, activation-recovery intervals in the BZ tended to be longer than in remote with increased spatial heterogeneity evidenced by a greater local SD (3.5±1.3 ms versus remote: 2.0±0.5 ms, P=0.036, npigs=5). Increased activation-recovery interval heterogeneity correlated with enhanced arrhythmia susceptibility. Cellular population studies (ncells=635-862 cells per region) demonstrated greater heterogeneity of action potential duration in the BZ (SD, 105.9±17.0 ms versus remote: 73.9±8.6 ms; P=0.001; npigs=6), which correlated with heterogeneity of activation-recovery interval in vivo. Cell-cell gene expression heterogeneity in the BZ was evidenced by increased Euclidean distances between nuclei of the BZ (12.1 [9.2-15.0] versus 10.6 [7.5-11.6] in remote; P<0.0001). Differentially expressed genes characterizing BZ cardiomyocyte remodeling included hypertrophy-related and ion channel-related genes with high cell-cell variability of expression. These gene expression changes were driven by stress-responsive TFs (transcription factors). In addition, heterogeneity of left ventricular wall thickness was greater in the BZ than in remote.

CONCLUSIONS

Heterogeneous cardiomyocyte remodeling in the BZ is driven by uniquely altered gene expression, related to heterogeneity in the local microenvironment, and translates to heterogeneous repolarization and arrhythmia vulnerability in vivo.

Inulin diet uncovers complex diet-microbiota-immune cell interactions remodeling the gut epithelium

BACKGROUND

The continuous proliferation of intestinal stem cells followed by their tightly regulated differentiation to epithelial cells is essential for the maintenance of the gut epithelial barrier and its functions. How these processes are tuned by diet and gut microbiome is an important, but poorly understood question. Dietary soluble fibers, such as inulin, are known for their ability to impact the gut bacterial community and gut epithelium, and their consumption has been usually associated with health improvement in mice and humans. In this study, we tested the hypothesis that inulin consumption modifies the composition of colonic bacteria and this impacts intestinal stem cells functions, thus affecting the epithelial structure.

METHODS

Mice were fed with a diet containing 5% of the insoluble fiber cellulose or the same diet enriched with an additional 10% of inulin. Using a combination of histochemistry, host cell transcriptomics, 16S microbiome analysis, germ-free, gnotobiotic, and genetically modified mouse models, we analyzed the impact of inulin intake on the colonic epithelium, intestinal bacteria, and the local immune compartment.

RESULTS

We show that the consumption of inulin diet alters the colon epithelium by increasing the proliferation of intestinal stem cells, leading to deeper crypts and longer colons. This effect was dependent on the inulin-altered gut microbiota, as no modulations were observed in animals deprived of microbiota, nor in mice fed cellulose-enriched diets. We also describe the pivotal role of γδ T lymphocytes and IL-22 in this microenvironment, as the inulin diet failed to induce epithelium remodeling in mice lacking this T cell population or cytokine, highlighting their importance in the diet-microbiota-epithelium-immune system crosstalk.

CONCLUSION

This study indicates that the intake of inulin affects the activity of intestinal stem cells and drives a homeostatic remodeling of the colon epithelium, an effect that requires the gut microbiota, γδ T cells, and the presence of IL-22. Our study indicates complex cross kingdom and cross cell type interactions involved in the adaptation of the colon epithelium to the luminal environment in steady state. Video Abstract.

Integrated single-nucleus and spatial transcriptomics captures transitional states in soybean nodule maturation

Legumes form symbiosis with rhizobium leading to the development of nitrogen-fixing nodules. By integrating single-nucleus and spatial transcriptomics, we established a cell atlas of soybean nodules and roots. In central infected zones of nodules, we found that uninfected cells specialize into functionally distinct subgroups during nodule development, and revealed a transitional subtype of infected cells with enriched nodulation-related genes. Overall, our results provide a single-cell perspective for understanding rhizobium-legume symbiosis.

Gata3 is detrimental to regulatory T cell function in autoimmune diabetes

Regulatory T cells (Tregs) protect against autoimmunity. In type 1 diabetes (T1D), Tregs slow the progression of beta cell autoimmunity within pancreatic islets. Increasing the potency or frequency of Tregs can prevent diabetes, as evidenced by studies in the nonobese diabetic (NOD) mouse model for T1D. We report herein that a significant proportion of islets Tregs in NOD mice express Gata3. The expression of Gata3 was correlated with the presence of IL-33, a cytokine known to induce and expand Gata3+ Tregs. Despite significantly increasing the frequency of Tregs in the pancreas, exogenous IL-33 was not protective. Based on these data, we hypothesized that Gata3 is deleterious to Treg function in autoimmune diabetes. To test this notion, we generated NOD mice with a Treg-specific deletion of Gata3. We found that deleting Gata3 in Tregs strongly protected against diabetes. Disease protection was associated with a shift of islet Tregs toward a suppressive CXCR3+Foxp3+ population. Our results suggest that islet Gata3+ Tregs are maladaptive and that this Treg subpopulation compromises the regulation of islet autoimmunity, contributing to diabetes onset.

Clonal associations of lymphocyte subsets and functional states revealed by single cell antigen receptor profiling of T and B cells in rheumatoid arthritis synovium

Rheumatoid arthritis (RA) is an autoimmune disease initiated by antigen-specific T cells and B cells, which promote synovial inflammation through a complex set of interactions with innate immune and stromal cells. To better understand the phenotypes and clonal relationships of synovial T and B cells, we performed single-cell RNA and repertoire sequencing on paired synovial tissue and peripheral blood samples from 12 donors with seropositive RA ranging from early to chronic disease. Paired transcriptomic-repertoire analyses highlighted 3 clonally distinct CD4 T cells populations that were enriched in RA synovium: T peripheral helper (Tph) and T follicular helper (Tfh) cells, CCL5+ T cells, and T regulatory cells (Tregs). Among these cells, Tph cells showed a unique transcriptomic signature of recent T cell receptor (TCR) activation, and clonally expanded Tph cells expressed an elevated transcriptomic effector signature compared to non-expanded Tph cells. CD8 T cells showed higher oligoclonality than CD4 T cells, and the largest CD8 T cell clones in synovium were highly enriched in GZMK+ cells. TCR analyses revealed CD8 T cells with likely viral-reactive TCRs distributed across transcriptomic clusters and definitively identified MAIT cells in synovium, which showed transcriptomic features of TCR activation. Among B cells, non-naive B cells including age-associated B cells (ABC), NR4A1+ activated B cells, and plasma cells, were enriched in synovium and had higher somatic hypermutation rates compared to blood B cells. Synovial B cells demonstrated substantial clonal expansion, with ABC, memory, and activated B cells clonally linked to synovial plasma cells. Together, these results reveal clonal relationships between functionally distinct lymphocyte populations that infiltrate RA synovium.

An integrated single cell and spatial transcriptomic map of human white adipose tissue

To date, single-cell studies of human white adipose tissue (WAT) have been based on small cohort sizes and no cellular consensus nomenclature exists. Herein, we performed a comprehensive meta-analysis of publicly available and newly generated single-cell, single-nucleus, and spatial transcriptomic results from human subcutaneous, omental, and perivascular WAT. Our high-resolution map is built on data from ten studies and allowed us to robustly identify >60 subpopulations of adipocytes, fibroblast and adipogenic progenitors, vascular, and immune cells. Using these results, we deconvolved spatial and bulk transcriptomic data from nine additional cohorts to provide spatial and clinical dimensions to the map. This identified cell-cell interactions as well as relationships between specific cell subtypes and insulin resistance, dyslipidemia, adipocyte volume, and lipolysis upon long-term weight changes. Altogether, our meta-map provides a rich resource defining the cellular and microarchitectural landscape of human WAT and describes the associations between specific cell types and metabolic states.

Placental cell type deconvolution reveals that cell proportions drive preeclampsia gene expression differences

The placenta mediates adverse pregnancy outcomes, including preeclampsia, which is characterized by gestational hypertension and proteinuria. Placental cell type heterogeneity in preeclampsia is not well-understood and limits mechanistic interpretation of bulk gene expression measures. We generated single-cell RNA-sequencing samples for integration with existing data to create the largest deconvolution reference of 19 fetal and 8 maternal cell types from placental villous tissue (n = 9 biological replicates) at term (n = 40,494 cells). We deconvoluted eight published microarray case-control studies of preeclampsia (n = 173 controls, 157 cases). Preeclampsia was associated with excess extravillous trophoblasts and fewer mesenchymal and Hofbauer cells. Adjustment for cellular composition reduced preeclampsia-associated differentially expressed genes (log2 fold-change cutoff = 0.1, FDR < 0.05) from 1154 to 0, whereas downregulation of mitochondrial biogenesis, aerobic respiration, and ribosome biogenesis were robust to cell type adjustment, suggesting direct changes to these pathways. Cellular composition mediated a substantial proportion of the association between preeclampsia and FLT1 (37.8%, 95% CI [27.5%, 48.8%]), LEP (34.5%, 95% CI [26.0%, 44.9%]), and ENG (34.5%, 95% CI [25.0%, 45.3%]) overexpression. Our findings indicate substantial placental cellular heterogeneity in preeclampsia contributes to previously observed bulk gene expression differences. This deconvolution reference lays the groundwork for cellular heterogeneity-aware investigation into placental dysfunction and adverse birth outcomes.

Single nucleus multiomics identifies ZEB1 and MAFB as candidate regulators of Alzheimer's disease-specific cis-regulatory elements

Cell type-specific transcriptional differences between brain tissues from donors with Alzheimer's disease (AD) and unaffected controls have been well documented, but few studies have rigorously interrogated the regulatory mechanisms responsible for these alterations. We performed single nucleus multiomics (snRNA-seq plus snATAC-seq) on 105,332 nuclei isolated from cortical tissues from 7 AD and 8 unaffected donors to identify candidate cis-regulatory elements (CREs) involved in AD-associated transcriptional changes. We detected 319,861 significant correlations, or links, between gene expression and cell type-specific transposase accessible regions enriched for active CREs. Among these, 40,831 were unique to AD tissues. Validation experiments confirmed the activity of many regions, including several candidate regulators of APP expression. We identified ZEB1 and MAFB as candidate transcription factors playing important roles in AD-specific gene regulation in neurons and microglia, respectively. Microglia links were globally enriched for heritability of AD risk and previously identified active regulatory regions.

Preparation of noninfectious scRNAseq samples from SARS-CoV-2-infected epithelial cells

Coronavirus disease (COVID-19) is an infectious disease caused by the SARS coronavirus 2 (SARS-CoV-2) virus. Direct assessment, detection, and quantitative analysis using high throughput methods like single-cell RNA sequencing (scRNAseq) is imperative to understanding the host response to SARS-CoV-2. One barrier to studying SARS-CoV-2 in the laboratory setting is the requirement to process virus-infected cell cultures, and potentially infectious materials derived therefrom, under Biosafety Level 3 (BSL-3) containment. However, there are only 190 BSL3 laboratory facilities registered with the U.S. Federal Select Agent Program, as of 2020, and only a subset of these are outfitted with the equipment needed to perform high-throughput molecular assays. Here, we describe a method for preparing non-hazardous RNA samples from SARS-CoV-2 infected cells, that enables scRNAseq analyses to be conducted safely in a BSL2 facility-thereby making molecular assays of SARS-CoV-2 cells accessible to a much larger community of researchers. Briefly, we infected African green monkey kidney epithelial cells (Vero-E6) with SARS-CoV-2 for 96 hours, trypsin-dissociated the cells, and inactivated them with methanol-acetone in a single-cell suspension. Fixed cells were tested for the presence of infectious SARS-CoV-2 virions using the Tissue Culture Infectious Dose Assay (TCID50), and also tested for viability using flow cytometry. We then tested the dissociation and methanol-acetone inactivation method on primary human lung epithelial cells that had been differentiated on an air-liquid interface. Finally, we performed scRNAseq quality control analysis on the resulting cell populations to evaluate the effects of our virus inactivation and sample preparation protocol on the quality of the cDNA produced. We found that methanol-acetone inactivated SARS-CoV-2, fixed the lung epithelial cells, and could be used to obtain noninfectious, high-quality cDNA libraries. This methodology makes investigating SARS-CoV-2, and related high-containment RNA viruses at a single-cell level more accessible to an expanded community of researchers.

scONE-seq: A single-cell multi-omics method enables simultaneous dissection of phenotype and genotype heterogeneity from frozen tumors

Single-cell multi-omics can provide a unique perspective on tumor cellular heterogeneity. Most previous single-cell whole-genome RNA sequencing (scWGS-RNA-seq) methods demonstrate utility with intact cells from fresh samples. Among them, many are not applicable to frozen samples that cannot produce intact single-cell suspensions. We have developed scONE-seq, a versatile scWGS-RNA-seq method that amplifies single-cell DNA and RNA without separating them from each other and hence is compatible with frozen biobanked samples. We benchmarked scONE-seq against existing methods using fresh and frozen samples to demonstrate its performance in various aspects. We identified a unique transcriptionally normal-like tumor clone by analyzing a 2-year frozen astrocytoma sample, demonstrating that performing single-cell multi-omics interrogation on biobanked tissue by scONE-seq could enable previously unidentified discoveries in tumor biology.

Single-cell transcriptomic analysis in two patients with rare systemic autoinflammatory diseases treated with anti-TNF therapy

Systemic autoinflammatory diseases (SAIDs) are a group of rare diseases characterized by recurrent or continuous inflammation, typically accompanied by genetic variants. Good responses to anti-TNF therapy were observed in SAIDs patients. However, the mechanisms underlying the disease flare and the response to TNF blocking therapy have not been fully elucidated. Here, single-cell RNA sequencing technology was used to describe the transcriptomic profile of PBMCs and PMNs in two SAID patients both before and after anti-TNF treatment. Interferon responses were involved in the disease flare. After anti-TNF therapy, clinical symptoms were alleviated while TNF and IL-1 were unexpectedly increased, indicating that these inflammatory cytokines are not positively correlated with disease activity. Trajectory analysis showed that inhibition of macrophage differentiation, rather than reduction of the inflammatory cytokines, as the potential mechanism of anti-TNF treatment response in SAIDs.

Targeting VEGF-A/VEGFR2 Y949 Signaling-Mediated Vascular Permeability Alleviates Hypoxic Pulmonary Hypertension

BACKGROUND

Pulmonary hypertension (PH) is associated with increased expression of VEGF-A (vascular endothelial growth factor A) and its receptor, VEGFR2 (vascular endothelial growth factor 2), but whether and how activation of VEGF-A signal participates in the pathogenesis of PH is unclear.

METHODS

VEGF-A/VEGFR2 signal activation and VEGFR2 Y949-dependent vascular leak were investigated in lung samples from patients with PH and mice exposed to hypoxia. To study their mechanistic roles in hypoxic PH, we examined right ventricle systolic pressure, right ventricular hypertrophy, and pulmonary vasculopathy in mutant mice carrying knock-in of phenylalanine that replaced the tyrosine at residual 949 of VEGFR2 (Vefgr2^Y949F) and mice with conditional endothelial deletion of Vegfr2 after chronic hypoxia exposure.

RESULTS

We show that PH leads to excessive pulmonary vascular leak in both patients and hypoxic mice, and this is because of an overactivated VEGF-A/VEGFR2 Y949 signaling axis. In the context of hypoxic PH, activation of Yes1 and c-Src and subsequent VE-cadherin phosphorylation in endothelial cells are involved in VEGFR2 Y949-induced vascular permeability. Abolishing VEGFR2 Y949 signaling by Vefgr2^Y949F point mutation was sufficient to prevent pulmonary vascular permeability and inhibit macrophage infiltration and Rac1 activation in smooth muscle cells under hypoxia exposure, thereby leading to alleviated PH manifestations, including muscularization of distal pulmonary arterioles, elevated right ventricle systolic pressure, and right ventricular hypertrophy. It is important that we found that VEGFR2 Y949 signaling in myeloid cells including macrophages was trivial and dispensable for hypoxia-induced vascular abnormalities and PH. In contrast with selective blockage of VEGFR2 Y949 signaling, disruption of the entire VEGFR2 signaling by conditional endothelial deletion of Vegfr2 promotes the development of PH.

CONCLUSIONS

Our results support the notion that VEGF-A/VEGFR2 Y949-dependent vascular permeability is an important determinant in the pathogenesis of PH and might serve as an attractive therapeutic target pathway for this disease.

A single-cell transcriptomic atlas of the human ciliary body

The ciliary body critically contributes to the ocular physiology with multiple responsibilities in the production of aqueous humor, vision accommodation and intraocular immunity. Comparatively little work, however, has revealed the single-cell molecular taxonomy of the human ciliary body required for studying these functionalities. In this study, we report a comprehensive atlas of the cellular and molecular components of human ciliary body as well as their interactions using single-cell RNA sequencing (scRNAseq). Cluster analysis of the transcriptome of 14,563 individual ciliary cells from the eyes of 3 human donors identified 14 distinct cell types, including the ciliary epithelium, smooth muscle, vascular endothelial cell, immune cell and other stromal cell populations. Cell-type discriminative gene markers were also revealed. Unique gene expression patterns essential for ciliary epithelium-mediated aqueous humor inflow and ciliary smooth muscle contractility were identified. Importantly, we discovered the transitional states that probably contribute to the transition of ciliary macrophage into retina microglia and verified no lymphatics in the ciliary body. Moreover, the utilization of CellPhoneDB allowed us to systemically infer cell–cell interactions among diverse ciliary cells including those that potentially participate in the pathogenesis of glaucoma and uveitis. Altogether, these new findings provide insights into the regulation of intraocular pressure, accommodation reflex and immune homeostasis under physiological and pathological conditions.

Transcription and splicing regulation by NLRC5 shape the interferon response in human pancreatic β cells

IFNα is a key regulator of the dialogue between pancreatic β cells and the immune system in early type 1 diabetes (T1D). IFNα up-regulates HLA class I expression in human β cells, fostering autoantigen presentation to the immune system. We observed by bulk and single-cell RNA sequencing that exposure of human induced pluripotent-derived islet-like cells to IFNα induces expression of HLA class I and of other genes involved in antigen presentation, including the transcriptional activator NLRC5. We next evaluated the global role of NLRC5 in human insulin-producing EndoC-βH1 and human islet cells by RNA sequencing and targeted gene/protein determination. NLRC5 regulates expression of HLA class I, antigen presentation-related genes, and chemokines. NLRC5 also mediates the effects of IFNα on alternative splicing, a generator of β cell neoantigens, suggesting that it is a central player of the effects of IFNα on β cells that contribute to trigger and amplify autoimmunity in T1D.

Single-cell RNA and protein profiling of immune cells from the mouse brain and its border tissues

Brain-immune cross-talk and neuroinflammation critically shape brain physiology in health and disease. A detailed understanding of the brain immune landscape is essential for developing new treatments for neurological disorders. Single-cell technologies offer an unbiased assessment of the heterogeneity, dynamics and functions of immune cells. Here we provide a protocol that outlines all the steps involved in performing single-cell multi-omic analysis of the brain immune compartment. This includes a step-by-step description on how to microdissect the border regions of the mouse brain, together with dissociation protocols tailored to each of these tissues. These combine a high yield with minimal dissociation-induced gene expression changes. Next, we outline the steps involved for high-dimensional flow cytometry and droplet-based single-cell RNA sequencing via the 10x Genomics platform, which can be combined with cellular indexing of transcriptomes and epitopes by sequencing (CITE-seq) and offers a higher throughput than plate-based methods. Importantly, we detail how to implement CITE-seq with large antibody panels to obtain unbiased protein-expression screening coupled to transcriptome analysis. Finally, we describe the main steps involved in the analysis and interpretation of the data. This optimized workflow allows for a detailed assessment of immune cell heterogeneity and activation in the whole brain or specific border regions, at RNA and protein level. The wet lab workflow can be completed by properly trained researchers (with basic proficiency in cell and molecular biology) and takes between 6 and 11 h, depending on the chosen procedures. The computational analysis requires a background in bioinformatics and programming in R.

Cell-type-specific cis-eQTLs in eight human brain cell types identify novel risk genes for psychiatric and neurological disorders

To date, most expression quantitative trait loci (eQTL) studies, which investigate how genetic variants contribute to gene expression, have been performed in heterogeneous brain tissues rather than specific cell types. In this study, we performed an eQTL analysis using single-nuclei RNA sequencing from 192 individuals in eight brain cell types derived from the prefrontal cortex, temporal cortex and white matter. We identified 7,607 eGenes, a substantial fraction (46%, 3,537/7,607) of which show cell-type-specific effects, with strongest effects in microglia. Cell-type-level eQTLs affected more constrained genes and had larger effect sizes than tissue-level eQTLs. Integration of brain cell type eQTLs with genome-wide association studies (GWAS) revealed novel relationships between expression and disease risk for neuropsychiatric and neurodegenerative diseases. For most GWAS loci, a single gene co-localized in a single cell type, providing new clues into disease etiology. Our findings demonstrate substantial contrast in genetic regulation of gene expression among brain cell types and reveal potential mechanisms by which disease risk genes influence brain disorders.

Spatial multi-omic map of human myocardial infarction

Myocardial infarction is a leading cause of death worldwide1. Although advances have been made in acute treatment, an incomplete understanding of remodelling processes has limited the effectiveness of therapies to reduce late-stage mortality2. Here we generate an integrative high-resolution map of human cardiac remodelling after myocardial infarction using single-cell gene expression, chromatin accessibility and spatial transcriptomic profiling of multiple physiological zones at distinct time points in myocardium from patients with myocardial infarction and controls. Multi-modal data integration enabled us to evaluate cardiac cell-type compositions at increased resolution, yielding insights into changes of the cardiac transcriptome and epigenome through the identification of distinct tissue structures of injury, repair and remodelling. We identified and validated disease-specific cardiac cell states of major cell types and analysed them in their spatial context, evaluating their dependency on other cell types. Our data elucidate the molecular principles of human myocardial tissue organization, recapitulating a gradual cardiomyocyte and myeloid continuum following ischaemic injury. In sum, our study provides an integrative molecular map of human myocardial infarction, represents an essential reference for the field and paves the way for advanced mechanistic and therapeutic studies of cardiac disease.