UCSC Cell Browser: visualize your single-cell data

Maternal autoantibodies to the sodium potassium pump α1 subunit AT1A1 and fetal autoimmune congenital heart block: a case-control study

BACKGROUND

Fetal autoimmune congenital heart block is a rare but life-threatening condition that is difficult to predict. This study sought to identify a serological biomarker predictive of autoimmune congenital heart block in pregnancies at risk due to maternal systemic lupus erythematosus, Sjögren's disease, undifferentiated connective tissue disease, or a history of congenital heart block offspring.

METHODS

This case-control study analysed maternal blood samples from pregnancies affected and unaffected by autoimmune congenital heart block from two centres (The Hospital for Sick Children, Canada, and The University of Padua, Italy). Serum samples collected across varying gestational ages were used for biomarker discovery, verification, and validation. The key inclusion criterion was a positive clinical test for maternal anti-Sjögren's syndrome-related antigen A/Ro autoantibodies, and the key exclusion criterion was structural congenital heart disease associated with heart block. Cases were pregnancies that resulted in fetal or neonatal congenital heart block and controls were pregnancies that resulted in offspring with a normal heart rhythm. Two-dimensional western blotting was used to identify maternal autoantibodies targeting fetal cardiac proteins, using fetal heart tissue and stem-cell-derived cardiomyocytes. Findings were validated using commercial proteins. Sensitivity and specificity for predicting fetal heart block outcomes were assessed using receiver-operating characteristic curves. The primary study outcomes were the presence and specificity of anti-cardiac autoantibodies in autoimmune congenital heart block cases versus controls, the association between specific autoantibodies and congenital heart block development, and the predictive accuracy of identified autoantibodies. This study did not involve individuals with lived experience in the study design or implementation.

FINDINGS

Serum samples were collected between Jan 1, 2010, and Dec 31, 2020, in the discovery cohort (The Hospital for Sick Children), between Aug 1, 1999, and Dec 31, 2019 in the verification cohort (University Hospital of Padua), and between June 1, 2018, and Aug 31, 2023, in the validation cohort (The Hospital for Sick Children). Mean maternal age was 34·2 years (SD 5·2) in the discovery cohort, 33·5 years (4·6) in the verification cohort, and 32·8 years (4·3) in the validation cohort. Maternal serum samples from pregnancies affected by fetal autoimmune congenital heart block (cases; n=46) showed an expanded repertoire of anti-cardiac autoantibodies compared with unaffected pregnancies (controls; n=65). Mass spectrometry identified 11 potential cardiac protein targets for maternal autoantibodies and the presence of seven autoantibodies were confirmed in serum samples from affected pregnancies. Four targets were identified before the onset of congenital heart block (ANXA1, BIP, MYPC3, and AT1A1). Maternal autoantibodies against the Na+ and K+ ATPase α1 isoform (AT1A1) were detected as early as 7 weeks gestation and were present in all affected pregnancies of mothers with and without previous fetal heart block history (six [100%] of six in the discovery cohort, 22 [100%] of 22 in the verification cohort, and 16 [100%] of 16 in the validation cohort). Anti-AT1A1 autoantibodies were absent in all unaffected pregnancies in pregnant women with no fetal heart block history (54/65 [83%]).

INTERPRETATION

This study identifies a novel maternal autoantibody targeting the fetal AT1A1 cardiac protein as a potential biomarker for the early and accurate detection of fetal autoimmune congenital heart block in pregnant women who had not previously had an affected pregnancy. This advancement might improve the maternal-fetal management of at-risk pregnancies through enhanced surveillance and timely and informed therapeutic interventions.

FUNDING

Heart and Stroke Foundation of Canada and Canadian Institutes of Health Research.

De novo variants disrupt an LDB1-regulated transcriptional network in congenital ventriculomegaly

Congenital hydrocephalus, characterized by cerebral ventriculomegaly, is among the most common and least understood paediatric neurosurgical disorders. We have identified, in the largest assembled cerebral ventriculomegaly cohort (2697 parent-proband trios), an exome-wide significant enrichment of protein-altering de novo variants in LDB1 (P = 1.11 × 10-15). Eight unrelated patients with ventriculomegaly, developmental delay and dysmorphic features harboured loss-of-function de novo variants that truncate carboxy-terminal LIM interaction domain of LDB1, which regulates assembly of LIM homeodomain-containing transcriptional modulators. Integrative multiomic analyses suggest that LDB1 is a key transcriptional regulator in ventricular neuroprogenitors through its binding to LIM-homeodomain proteins, including SMARCC1 and ARID1B. Indeed, LIM-homeodomain-containing genes carry a disproportionate burden of protein-damaging de novo variants in our cohort, with SMARCC1 (P = 5.83 × 10-9) and ARID1B (P = 1.80 × 10-17) surpassing exome-wide significance thresholds. These data identify LBD1 as a novel neurodevelopmental disorder gene and suggest that an LDB1-regulated transcriptional programme is essential for human brain morphogenesis.

Integrated analysis of molecular atlases unveils modules driving developmental cell subtype specification in the human cortex

Human brain development requires generating diverse cell types, a process explored by single-cell transcriptomics. Through parallel meta-analyses of the human cortex in development (seven datasets) and adulthood (16 datasets), we generated over 500 gene co-expression networks that can describe mechanisms of cortical development, centering on peak stages of neurogenesis. These meta-modules show dynamic cell subtype specificities throughout cortical development, with several developmental meta-modules displaying spatiotemporal expression patterns that allude to potential roles in cell fate specification. We validated the expression of these modules in primary human cortical tissues. These include meta-module 20, a module elevated in FEZF2+ deep layer neurons that includes TSHZ3, a transcription factor associated with neurodevelopmental disorders. Human cortical chimeroid experiments validated that both FEZF2 and TSHZ3 are required to drive module 20 activity and deep layer neuron specification but through distinct modalities. These studies demonstrate how meta-atlases can engender further mechanistic analyses of cortical fate specification.

Monoallelic loss-of-function variants in GSK3B lead to autism and developmental delay

De novo variants adjacent to the canonical splicing sites or in the well-defined splicing-related regions are more likely to impair splicing but remain under-investigated in autism spectrum disorder (ASD). By analyzing large, recent ASD genome sequencing cohorts, we find a significant burden of de novo potential splicing-disrupting variants (PSDVs) in 5048 probands compared to 4090 unaffected siblings. We identified 55 genes with recurrent de novo PSDVs that were highly intolerant to variation. Forty-six of these genes have not been strongly implicated in ASD or other neurodevelopmental disorders previously, including GSK3B. Through international, multicenter collaborations, we assembled genotype and phenotype data for 15 individuals with GSK3B variants and identified common phenotypes including developmental delay, ASD, sleeping disturbance, and aggressive behavior. Using available single-cell transcriptomic data, we show that GSK3B is enriched in dorsal progenitors and intermediate forms of excitatory neurons in the developing brain. We showed that Gsk3b knockdown in mouse excitatory neurons interferes with dendrite arborization and spine maturation which could not be rescued by de novo missense variants identified from affected individuals. In summary, our findings suggest that PSDVs may play an important role in the genetic etiology of ASD and allow for the prioritization of new ASD candidate genes. Importantly, we show that genetic variation resulting in GSK3B loss-of-function can lead to a neurodevelopmental disorder with core features of ASD and developmental delay.

Multi-omics pan-cancer analysis of monocyte to macrophage differentiation-associated (MMD) and its significance in hepatocellular carcinoma

BackgroundMalignant tumors are serious diseases that endanger human health. Therefore, it is crucial to identify markers that facilitate tumor diagnosis and prognostic assessment.ObjectiveThis study analyzed the significance of Monocyte to macrophage differentiation-associated (MMD) in various tumors from multiple perspectives, to explore the possibility of using MMD as a novel tumor marker.MethodsUsing the R software, an examination of MMD levels was conducted across diverse human cancers and their influence on cancer outcomes. MMD methylation, mutations, and immune infiltration analyses of various tumors were performed. A Cox regression model was used to predict the survival rates of patients with hepatocellular carcinoma (HCC). Finally, MMD expression and function were validated in Hep-3B cells.ResultsMMD was aberrantly expressed in diverse tumors and can predict patient outcomes. Methylation and functional enrichment studies indicated possible function of MMD in tumor progression, whereas immune infiltration data suggested its involvement in tumor immune evasion. Cox regression analysis revealed that elevated MMD levels were independent predictors of HCC patient outcomes. The quantitative real-time polymerase chain reaction (qPCR) data demonstrated high MMD levels in Hep-3B cells, and its suppression impeded Hep-3B cell growth.ConclusionsMMD was abnormally expressed in various tumors and was closely associated with tumor prognosis. Thus, it had the potential to be used as a novel tumor marker.

Modeling forebrain regional development and connectivity by human brain organoids

The forebrain is one of the most important brain structures for modern human existence, which houses the uniquely sophisticated social and cognitive functions that distinguish our species. Therefore, modeling the forebrain development by using human cells is especially critical for our understanding of the intricacies of human development and devising treatments for related diseases. Recent advancements in brain organoid fields have offered unprecedented tools to investigate forebrain development from studies on specific regions to exploring tract formation and connectivity between different regions of the forebrain. In this review, we discuss the developmental biology of the forebrain and diverse methods for modeling its development by using organoids.

Methylation patterns of the nasal epigenome of hospitalized SARS-CoV-2 positive patients reveal insights into molecular mechanisms of COVID-19

BACKGROUND

Coronavirus disease 2019 (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has varied presentations from asymptomatic to death. Efforts to identify factors responsible for differential COVID-19 severity include but are not limited to genome wide association studies (GWAS) and transcriptomic analysis. More recently, variability in host epigenomic profiles have garnered attention, providing links to disease severity. However, whole epigenome analysis of the respiratory tract, the target tissue of SARS-CoV-2, remains ill-defined.

RESULTS

We interrogated the nasal methylome to identify pathophysiologic drivers in COVID-19 severity through whole genome bisulfite sequencing (WGBS) of nasal samples from COVID-19 positive individuals with severe and mild presentation of disease. We noted differential DNA methylation in intergenic regions and low methylated regions (LMRs), demonstrating the importance of distal regulatory elements in gene regulation in COVID-19 illness. Additionally, we demonstrated differential methylation of pathways implicated in immune cell recruitment and function, and the inflammatory response. We found significant hypermethylation of the FUT4 promoter implicating impaired neutrophil adhesion in severe disease. We also identified hypermethylation of ELF5 binding sites suggesting downregulation of ELF5 targets in the nasal cavity as a factor in COVID-19 phenotypic variability.

CONCLUSIONS

This study demonstrated DNA methylation as a marker of the immune response to SARS-CoV-2 infection, with enhancer-like elements playing significant roles. It is difficult to discern whether this differential methylation is a predisposing factor to severe COVID-19, or if methylation differences occur in response to disease severity. These differences in the nasal methylome may contribute to disease severity, or conversely, the nasal immune system may respond to severe infection through differential immune cell recruitment and immune function, and through differential regulation of the inflammatory response.

Microglia regulate myelin clearance and cholesterol metabolism after demyelination via interferon regulatory factor 5

Interferon regulatory factor 5 (IRF5) is a transcription factor that plays a role in orchestrating innate immune responses, particularly in response to viral infections. Notably, IRF5 has been identified as a microglia risk gene linked to multiple sclerosis (MS), but its specific role in MS pathogenesis remains unclear. Through the use of Irf5-/- mice, our study uncovers a non-canonical function of IRF5 in MS recovery. Irf5-/- mice exhibited increased damage in an experimental autoimmune encephalomyelitis (EAE) model and demonstrated impaired oligodendrocyte recruitment into the lesion core following lysolecithin-induced demyelination. Transcriptomic and lipidomic analyses revealed that IRF5 has a role in microglia-mediated myelin phagocytosis, lipid metabolism, and cholesterol homeostasis. Indeed, Irf5-/- microglia phagocytose myelin, but myelin debris is not adequately degraded, leading to an accumulation of lipid droplets, cholesterol esters, and cholesterol crystals within demyelinating lesions. This abnormal buildup can hinder remyelination processes. Importantly, treatments that promote cholesterol transport were found to reduce lipid droplet accumulation and mitigate the exacerbated damage in Irf5-/- mice with EAE. Altogether, our study identified the antiviral transcription factor IRF5 as a key transcriptional regulator of lipid degradation and cholesterol homeostasis and suggest that loss of IRF5 function leads to pathogenic lipid accumulation in microglia, thereby obstructing remyelination. These data and the fact that Irf5 polymorphisms are significantly associated with MS, highlight IRF5 as a potential therapeutic target to promote regenerative responses.

scMEDAL for the interpretable analysis of single-cell transcriptomics data with batch effect visualization using a deep mixed effects autoencoder

scRNA-seq data has the potential to provide new insights into cellular heterogeneity and data acquisition; however, a major challenge is unraveling confounding from technical and biological batch effects. Existing batch correction algorithms suppress and discard these effects, rather than quantifying and modeling them. Here, we present scMEDAL, a framework for s ingle- c ell M ixed E ffects D eep A utoencoder L earning, which separately models batch-invariant and batch-specific effects using two complementary autoencoder networks. One network is trained through adversarial learning to capture a batch-invariant representation, while a Bayesian autoencoder learns a batch-specific representation. Comprehensive evaluations spanning conditions (e.g., autism, leukemia, and cardiovascular), cell types, and technical and biological effects demonstrate that scMEDAL suppresses batch effects while modeling batch-specific variation, enhancing accuracy and interpretability. Unlike prior approaches, the framework's fixed- and random-effects autoencoders enable retrospective analyses, including predicting a cell's expression as if it had been acquired in a different batch via genomap projections at the cellular level, revealing the impact of biological (e.g., diagnosis) and technical (e.g., acquisition) effects. By combining scMEDAL's batch-agnostic and batch-specific latent spaces, it enables more accurate predictions of disease status, donor group, and cell type, making scMEDAL a valuable framework for gaining deeper insight into data acquisition and cellular heterogeneity.

scMEDAL for the interpretable analysis of single-cell transcriptomics data with batch effect visualization using a deep mixed effects autoencoder

scRNA-seq data has the potential to provide new insights into cellular heterogeneity and data acquisition; however, a major challenge is unraveling confounding from technical and biological batch effects. Existing batch correction algorithms suppress and discard these effects, rather than quantifying and modeling them. Here, we present scMEDAL, a framework for single-cell Mixed Effects Deep Autoencoder Learning, which separately models batch-invariant and batch-specific effects using two complementary autoencoder networks. One network is trained through adversarial learning to capture a batch-invariant representation, while a Bayesian autoencoder learns a batch-specific representation. Comprehensive evaluations spanning conditions (e.g., autism, leukemia, and cardiovascular), cell types, and technical and biological effects demonstrate that scMEDAL suppresses batch effects while modeling batch-specific variation, enhancing accuracy and interpretability. Unlike prior approaches, the framework's fixed- and random-effects autoencoders enable retrospective analyses, including predicting a cell's expression as if it had been acquired in a different batch via genomap projections at the cellular level, revealing the impact of biological (e.g., diagnosis) and technical (e.g., acquisition) effects. By combining scMEDAL's batch-agnostic and batch-specific latent spaces, it enables more accurate predictions of disease status, donor group, and cell type, making scMEDAL a valuable framework for gaining deeper insight into data acquisition and cellular heterogeneity.

Single-Cell Analysis of Sex and Gender Differences in the Human Brain During Development and Disease

Sex and gender (SG) differences in the human brain are of interest to society and science as numerous processes are impacted by them, including brain development, behavior, and diseases. By collecting publicly available single-cell data from the in-utero to elderly age in healthy, Alzheimer's disease and multiple sclerosis samples, we identified and characterized SG-biased genes in ten brain cell types across 9 age and disease groups. Sex and gender differences in the transcriptome were present throughout the lifespan and across all cell types. Although there was limited overlap among SG-biased genes across different age and disease groups, we observed significant functional overlap. Female-biased genes are consistently enriched for brain-related processes, while male-biased genes are enriched for metabolic pathways. Additionally, mitochondrial genes showed a consistent female bias across cell types. We also found that androgen response elements (not estrogen) were significantly enriched in both male- and female-biased genes, and thymosin hormone targets being consistently enriched only in male-biased genes. We systematically characterised SG differences in brain development and brain-related disorders at a single-cell level, by analysing a total of publicly available 419,885 single nuclei from 161 human brain samples (72 females, 89 males). The significant enrichment of androgen (not estrogen) response elements in both male- and female-biased genes suggests that androgens are important regulators likely establishing these SG differences. Finally, we provide full characterization of SG-biased genes at different thresholds for the scientific community as a web resource.

Mosquito Cell Atlas: A single-nucleus transcriptomic atlas of the adult Aedes aegypti mosquito

The female mosquito's remarkable ability to hunt humans and transmit pathogens relies on her unique biology. Here, we present the Mosquito Cell Atlas (MCA), a comprehensive single-nucleus RNA sequencing dataset of more than 367,000 nuclei from 19 dissected tissues of adult female and male Aedes aegypti, providing cellular-level resolution of mosquito biology. We identify novel cell types and expand our understanding of sensory neuron organization of chemoreceptors to all sensory tissues. Our analysis uncovers male-specific cells and sexually dimorphic gene expression in the antenna and brain. In female mosquitoes, we find that glial cells in the brain, rather than neurons, undergo the most extensive transcriptional changes following blood feeding. Our findings provide insights into the cellular basis of mosquito behavior and sexual dimorphism. The MCA aims to serve as a resource for the vector biology community, enabling systematic investigation of cell-type specific expression across all mosquito tissues.

Characterization of canine tumor-infiltrating leukocyte transcriptomic signatures reveals conserved expression patterns with human osteosarcoma

Immune cells play key roles in host responses to malignant tumors. The selective pressure that immune cells elicit on tumors promotes immune escape, while tumor-associated modulation of immune cells creates an environment favorable to tumor growth and progression. In this study we used publicly available single-cell RNA sequencing (scRNA-seq) data from the translationally relevant canine osteosarcoma (OS) model to compare tumor-infiltrating immune cells to circulating leukocytes. Through computational analysis we investigated the differences in cell type proportions and how the OS TME impacted infiltrating immune cell transcriptomic profiles relative to circulating leukocytes. Differential abundance analysis revealed increased proportions of follicular helper T cells, regulatory T cells, and mature regulatory dendritic cells (mregDCs) in the OS TME. Differential gene expression analysis identified exhaustion markers (LAG3, HAVCR2, PDCD1) to be upregulated in CD4 and CD8 T cells within the OS TME. Comparisons of B cell gene expression profiles revealed an enrichment of protein processing and endoplasmic reticulum pathways, suggesting infiltrating B cells were activated following tumor infiltration. Gene expression changes within myeloid cells identified increased expression of immune suppressive molecules (CD274, OSM, MSR1) in the OS TME, indicating the TME skews myeloid cells toward an immunosuppressive phenotype. Comparisons to human literature and analysis of human scRNA-seq data revealed conserved transcriptomic responses to tumor infiltration, while also identifying species differences. Overall, the analysis presented here provides new insights into how the OS TME impacts the transcriptional programs of major immune cell populations in dogs and acts as a resource for comparative immuno-oncology research.

paraCell: a novel software tool for the interactive analysis and visualization of standard and dual host-parasite single-cell RNA-seq data

Advances in sequencing technology have led to a dramatic increase in the number of single-cell transcriptomic datasets. In the field of parasitology, these datasets typically describe the gene expression patterns of a given parasite species at the single-cell level under experimental conditions, in specific hosts or tissues, or at different life cycle stages. However, while this wealth of available data represents a significant resource, analysing these datasets often requires expert computational skills, preventing a considerable proportion of the parasitology community from meaningfully integrating existing single-cell data into their work. Here, we present paraCell, a novel software tool that allows the user to visualize and analyse pre-loaded single-cell data without requiring any programming ability. The source code is free to allow remote installation. On our web server, we demonstrated how to visualize and re-analyse published Plasmodium and Trypanosoma datasets. We have also generated Toxoplasma-mouse and Theileria-cow scRNA-seq datasets to highlight the functionality of paraCell for pathogen-host interaction. The analysis of the data highlights the impact of the host interferon-γ response and gene expression profiles associated with disease susceptibility by these intracellular parasites, respectively.

Meteorin-like protein/METRNL/Interleukin-41 ameliorates atopic dermatitis-like inflammation

BACKGROUND

Meteorin-like protein (METRNL)/Interleukin-41 (IL-41) is a novel immune-secreted cytokine/myokine involved in several inflammatory diseases. However, how METRNL exerts its regulatory properties on skin inflammation remains elusive. This study aims to elucidate the functionality and regulatory mechanism of METRNL in atopic dermatitis (AD).

METHODS

METRNL levels were determined in skin and serum samples from patients with AD and subsequently verified in the vitamin D3 analogue MC903-induced AD-like mice model. The cellular target of METRNL activity was identified by multiplex immunostaining, single-cell RNA-seq and RNA-seq.

RESULTS

METRNL was significantly upregulated in lesions and serum of patients with dermatitis compared to healthy controls (p <.05). Following repeated MC903 exposure, AD model mice displayed elevated levels of METRNL in both ears and serum. Administration of recombinant murine METRNL protein (rmMETRNL) ameliorated allergic skin inflammation and hallmarks of AD in mice, whereas blocking of METRNL signaling led to the opposite. METRNL enhanced β-Catenin activation, limited the expression of Th2-related molecules that attract the accumulation of Arginase-1 (Arg1)hi macrophages, dendritic cells, and activated mast cells.

CONCLUSIONS

METRNL can bind to KIT receptor and subsequently alleviate the allergic inflammation of AD by inhibiting the expansion of immune cells, and downregulating inflammatory gene expression by regulating the level of active WNT pathway molecule β-Catenin.

ciBAR1 loss in mice causes laterality defects, pancreatic degeneration, and altered glucose tolerance

Bin/Amphiphysin/Rvs (BAR) domains are highly conserved domains found in all eukaryotes. BAR domain proteins form crescent-shaped dimers that sense and sculpt curved lipid membranes and play key roles in various cellular processes. However, their functions in mammalian development are poorly understood. We previously demonstrated that Chibby1-interacting BAR domain-containing 1 (ciBAR1, formerly known as FAM92A) localizes to the ciliary base and plays a critical role in ciliogenesis. Here, we report ciliopathy phenotypes of ciBAR1-KO mice. We found that ∼28% of ciBAR1-KO mice show embryonic lethality because of randomized left-right asymmetry; the rest survive into adulthood with no gross morphological abnormalities. Histological assessments of ciliated tissues revealed exocrine pancreatic lesions. Although overall endocrine islet morphology appeared to be normal, ciBAR1-KO mice showed impaired glucose tolerance. Examination of ductal and islet cilia revealed that cilia number and length were significantly reduced in ciBAR1-KO pancreata. ciBAR1-KO MEFs also exhibited ciliary defects. Our findings indicate that ciBAR1 plays a critical role in ciliogenesis depending on the tissue and cell type in mice.

Manipulating the EphB4-ephrinB2 axis to reduce metastasis in HNSCC

The EphB4-ephrinB2 signaling axis has been heavily implicated in metastasis across numerous cancer types. Our emerging understanding of the dichotomous roles that EphB4 and ephrinB2 play in head and neck squamous cell carcinoma (HNSCC) poses a significant challenge to rational drug design. We find that EphB4 knockdown in cancer cells enhances metastasis in preclinical HNSCC models by augmenting immunosuppressive cells like T regulatory cells (Tregs) within the tumor microenvironment. EphB4 inhibition in cancer cells also amplifies their ability to metastasize through increased expression of genes associated with hallmark pathways of metastasis along with classical and non-classical epithelial-mesenchymal transition. In contrast, vascular ephrinB2 knockout coupled with radiation therapy (RT) enhances anti-tumor immunity, reduces Treg accumulation into the tumor, and decreases metastasis. Notably, targeting the EphB4-ephrinB2 signaling axis with the engineered ligands ephrinB2-Fc-His and Fc-TNYL-RAW-GS reduces local tumor growth and distant metastasis in a preclinical model of HNSCC. Our data suggests that targeted inhibition of vascular ephrinB2 while avoiding inhibition of EphB4 in cancer cells could be a promising strategy to mitigate HNSCC metastasis.

Human neural rosettes secrete bioactive extracellular vesicles enriched in neuronal and glial cellular components

Extracellular vesicles (EVs) play a critical role in the development of neural cells in the central nervous system (CNS). Human neural rosettes (hNRs) are radial cell structures that assemble from induced pluripotent stem cells (hiPSCs) and recapitulate some stages of neural tube morphogenesis. Here we show that hiPSCs and hNRs secrete EVs (hiPSC-EVs and hNR-EVs) with distinctive protein cargoes. Remarkably, hNR-EVs carry neuronal and glial cellular components involved in human CNS development. Importantly, hNR-EVs stimulate stem cells to change their cellular morphology and promote neurite growth in human and murine neurons with a significant dysregulation of SOX2 levels. This transcription factor modulates both neural differentiation and pluripotency. Interestingly, these effects were inhibited by antibodies against an unexpected neuroglial cargo of hNR-EVs: the major proteolipid protein (PLP). These findings show that hNRs secrete bioactive EVs containing neural components and might contribute as trophic factors during human neurodevelopment.

MAPbrain: a multi-omics atlas of the primate brain

The brain is the central hub of the entire nervous system. Its development is a lifelong process guided by a genetic blueprint. Understanding how genes influence brain development is critical for deciphering the formation of human cognitive functions and the underlying mechanisms of neurological disorders. Recent advances in multi-omics techniques have now made it possible to explore these aspects comprehensively. However, integrating and analyzing extensive multi-omics data presents significant challenges. Here, we introduced MAPbrain (http://bigdata.ibp.ac.cn/mapBRAIN/), a multi-omics atlas of the primate brain. This repository integrates and normalizes both our own lab's published data and publicly available multi-omics data, encompassing 21 million brain cells from 38 key brain regions and 436 sub-regions across embryonic and adult stages, with 164 time points in humans and non-human primates. MAPbrain offers a unique, robust, and interactive platform that includes transcriptomics, epigenomics, and spatial transcriptomics data, facilitating a comprehensive exploration of brain development. The platform enables the exploration of cell type- and time point-specific markers, gene expression comparison between brain regions and species, joint analyses across transcriptome and epigenome, and navigation of cell types across species, brain regions, and development stages. Additionally, MAPbrain provides an online integration module for users to navigate and analyze their own data within the platform.

Accelerating Single-Cell Sequencing Data Analysis with SciDAP: A User-Friendly Approach

Single-cell (sc) RNA, ATAC, and Multiome sequencing became powerful tools for uncovering biological and disease mechanisms. Unfortunately, manual analysis of sc data presents multiple challenges due to large data volumes and complexity of configuration parameters. This complexity, as well as not being able to reproduce a computational environment, affects the reproducibility of analysis results. The Scientific Data Analysis Platform ( https://SciDAP.com ) allows biologists without computational expertise to analyze sequencing-based data using portable and reproducible pipelines written in Common Workflow Language (CWL). Our suite of computational pipelines addresses the most common needs in scRNA-Seq, scATAC-Seq and scMultiome data analysis. When executed on SciDAP, it offers a user-friendly alternative to manual data processing, eliminating the need for coding expertise. In this protocol, we describe the use of SciDAP to analyze scMultiome data. Similar approaches can be used for analysis of scRNA-Seq, scATAC-Seq and scVDJ-Seq datasets.

Considerations for building and using integrated single-cell atlases

The rapid adoption of single-cell technologies has created an opportunity to build single-cell 'atlases' integrating diverse datasets across many laboratories. Such atlases can serve as a reference for analyzing and interpreting current and future data. However, it has become apparent that atlasing approaches differ, and the impact of these differences are often unclear. Here we review the current atlasing literature and present considerations for building and using atlases. Importantly, we find that no one-size-fits-all protocol for atlas building exists, but rather we discuss context-specific considerations and workflows, including atlas conceptualization, data collection, curation and integration, atlas evaluation and atlas sharing. We further highlight the benefits of integrated atlases for analyses of new datasets and deriving biological insights beyond what is possible from individual datasets. Our overview of current practices and associated recommendations will improve the quality of atlases to come, facilitating the shift to a unified, reference-based understanding of single-cell biology.

Further evidence of biallelic NAV3 variants associated with recessive neurodevelopmental disorder with dysmorphism, developmental delay, intellectual disability, and behavioral abnormalities

Neuron navigators (NAVs) are cytoskeleton-associated proteins well known for their role in axonal guidance, neuronal migration, and neurite growth necessary for neurodevelopment. Neuron navigator 3 (NAV3) is one of the three NAV proteins highly expressed in the embryonic and adult brain. However, the role of the NAV3 gene in human disease is not well-studied. Recently, five bi-allelic and three mono-allelic variants in NAV3 were reported in 12 individuals from eight unrelated families with neurodevelopmental disorder (NDD). Here, we report five patients from three unrelated consanguineous families segregating autosomal recessive NDD. Patients have symptoms of dysmorphism, intellectual disability, developmental delay, and behavioral abnormalities. Exome sequencing (ES) was performed on two affected individuals from one large family, and one affected individual from each of the other two families. ES revealed two homozygous nonsense c.6325C > T; p.(Gln2109Ter) and c.6577C > T; p.(Arg2193Ter) and a homozygous splice site (c.243 + 1G > T) variants in the NAV3 (NM_001024383.2). Analysis of single-cell sequencing datasets from embryonic and young adult human brains revealed that NAV3 is highly expressed in the excitatory neurons, inhibitory neurons, and microglia, consistent with its role in neurodevelopment. In conclusion, in this study, we further validate biallelic protein truncating variants in NAV3 as a cause of NDD, expanding the spectrum of pathogenic variants in this newly discovered NDD gene.

Alleviating batch effects in cell type deconvolution with SCCAF-D

Cell type deconvolution methods can impute cell proportions from bulk transcriptomics data, revealing changes in disease progression or organ development. But benchmarking studies often use simulated bulk data from the same source as the reference, which limits its application scenarios. This study examines batch effects in deconvolution and introduces SCCAF-D, a computational workflow that ensures a Pearson Correlation Coefficient above 0.75 across simulated and real bulk data for various tissue types. Applied to non-alcoholic fatty liver disease, SCCAF-D unveils meaningful insights into changes in cell proportions during disease progression.

Perivascular Macrophages Convert Physical Wound Signals Into Rapid Vascular Responses

Leukocytes detect distant wounds within seconds to minutes, which is essential for effective pathogen defense, tissue healing, and regeneration. Blood vessels must detect distant wounds just as rapidly to initiate local leukocyte extravasation, but the mechanism behind this immediate vascular response remains unclear. Using high-speed imaging of live zebrafish larvae, we investigated how blood vessels achieve rapid wound detection. We monitored two hallmark vascular responses: vessel dilation and serum exudation. Our experiments-including genetic, pharmacologic, and osmotic perturbations, along with chemogenetic leukocyte depletion-revealed that the cPla2 nuclear shape sensing pathway in perivascular macrophages converts a fast (~50 μm/s) osmotic wound signal into a vessel-permeabilizing, 5-lipoxygenase (Alox5a) derived lipid within seconds of injury. These findings demonstrate that perivascular macrophages act as physicochemical relays, bridging osmotic wound signals and vascular responses. By uncovering this novel type of communication, we provide new insights into the coordination of immune and vascular responses to injury.

Hypoxia-inducible factor 2α promotes pathogenic polarization of stem-like Th2 cells via modulation of phospholipid metabolism

T helper 2 (Th2) cells orchestrate immunity against parasite infection and promote tissue repair but promote pathology in asthma and tissue fibrosis. Here, we examined the mechanisms driving pathogenic differentiation of Th2 cells. Single-cell analyses of CD4+ T cells from asthma and chronic rhinosinusitis patients revealed high expression of the hypoxia-inducible factor (HIF)2α in Th2 cells. In mice, HIF2α deficiency impaired Th2 differentiation and alleviated asthmatic inflammation. Single-cell and lineage tracing approaches delineated a differentiation trajectory from TCF1+Ly108+ stem-like Th2 cells to the ST2+CD25+ pathogenic progeny, depending on a HIF2α-GATA3 circuit that modulated phospholipid metabolism and T cell receptor (TCR)-phosphatidylinositol 3-kinase (PI3K)-protein kinase B (AKT) activation via transcriptional regulation of the inositol polyphosphate multikinase (IPMK). Overexpression of IPMK in HIF2α-deficient cells promoted Phosphatidylinositol (3,4,5)-trisphosphate (PIP3) synthesis and pathogenic Th2 cell differentiation, whereas pharmacological inhibition of HIF2α impaired pathogenic differentiation of Th2 cells and mitigated airway inflammation. Our findings provide insight into the contextual cues that promote Th2-mediated pathology and suggest HIF2α as a therapeutic target in asthma.

3D genome topology distinguishes molecular subgroups of medulloblastoma

Four main medulloblastoma (MB) molecular subtypes have been identified based on transcriptional, DNA methylation, and genetic profiles. However, it is currently not known whether 3D genome architecture differs between MB subtypes. To address this question, we performed in situ Hi-C to reconstruct the 3D genome architecture of MB subtypes. In total, we generated Hi-C and matching transcriptome data for 28 surgical specimens and Hi-C data for one patient-derived xenograft. The average resolution of the Hi-C maps was 6,833 bp. Using these data, we found that insulation scores of topologically associating domains (TADs) were effective at distinguishing MB molecular subgroups. TAD insulation score differences between subtypes were globally not associated with differential gene expression, although we identified few exceptions near genes expressed in the lineages of origin of specific MB subtypes. Our study therefore supports the notion that TAD insulation scores can distinguish MB subtypes independently of their transcriptional differences.

Comprehensive Analysis of the Tumour Immune Microenvironment in Canine Urothelial Carcinoma Reveals Immunosuppressive Mechanisms Induced by the COX-Prostanoid Cascade

A comprehensive understanding of the tumour immune microenvironment (TIME) is essential for advancing precision medicine and identifying potential therapeutic targets. This study focused on canine urothelial carcinoma (cUC) recognised for its high sensitivity to cyclooxygenase (COX) inhibitors. Using immunohistochemical techniques, we quantified the infiltration of seven immune cell populations within cUC tumour tissue to identify clinicopathological features that characterise the TIME in cUC. Our results revealed several notable factors, including the significantly higher levels of CD3+ T cells and CD8+ T cells within tumour cell nests in cases treated with preoperative COX inhibitors compared to untreated cases. Based on the immunohistochemistry data, we further performed a comparative analysis using publicly available RNA-seq data from untreated cUC tissues (n = 29) and normal bladder tissues (n = 4) to explore the link between COX-prostanoid pathways and the immune response to tumours. We observed increased expression of COX-2, microsomal prostaglandin E2 synthase-1 (mPGES-1) and mPGES-2 in cUC tissues. However, only mPGES-2 showed a negative correlation with the cytotoxic T-cell (CTL)-related genes CD8A and granzyme B (GZMB). In addition, a broader analysis of solid tumours using The Cancer Genome Atlas (TCGA) database revealed similar patterns in several human tumours, suggesting a common mechanism in dogs and humans. Our results suggest that the COX-2/mPGES-2 pathway may act as a cross-species tumour-intrinsic factor that weakens anti-tumour immunity, and that COX inhibitors may convert TIME from a 'cold tumour' to a 'hot tumour' state by counteracting COX/mPGES-2-mediated immunosuppression.

Building a FAIR data ecosystem for incorporating single-cell transcriptomics data into agricultural genome to phenome research

Introduction

The agriculture genomics community has numerous data submission standards available, but the standards for describing and storing single-cell (SC, e.g., scRNA- seq) data are comparatively underdeveloped.

Methods

To bridge this gap, we leveraged recent advancements in human genomics infrastructure, such as the integration of the Human Cell Atlas Data Portal with Terra, a secure, scalable, open-source platform for biomedical researchers to access data, run analysis tools, and collaborate. In parallel, the Single Cell Expression Atlas at EMBL-EBI offers a comprehensive data ingestion portal for high-throughput sequencing datasets, including plants, protists, and animals (including humans). Developing data tools connecting these resources would offer significant advantages to the agricultural genomics community. The FAANG data portal at EMBL-EBI emphasizes delivering rich metadata and highly accurate and reliable annotation of farmed animals but is not computationally linked to either of these resources.

Results

Herein, we describe a pilot-scale project that determines whether the current FAANG metadata standards for livestock can be used to ingest scRNA-seq datasets into Terra in a manner consistent with HCA Data Portal standards. Importantly, rich scRNA-seq metadata can now be brokered through the FAANG data portal using a semi-automated process, thereby avoiding the need for substantial expert curation. We have further extended the functionality of this tool so that validated and ingested SC files within the HCA Data Portal are transferred to Terra for further analysis. In addition, we verified data ingestion into Terra, hosted on Azure, and demonstrated the use of a workflow to analyze the first ingested porcine scRNA-seq dataset. Additionally, we have also developed prototype tools to visualize the output of scRNA-seq analyses on genome browsers to compare gene expression patterns across tissues and cell populations. This JBrowse tool now features distinct tracks, showcasing PBMC scRNA-seq alongside two bulk RNA-seq experiments.

Discussion

We intend to further build upon these existing tools to construct a scientist-friendly data resource and analytical ecosystem based on Findable, Accessible, Interoperable, and Reusable (FAIR) SC principles to facilitate SC-level genomic analysis through data ingestion, storage, retrieval, re-use, visualization, and comparative annotation across agricultural species.

Patient iPSC models reveal glia-intrinsic phenotypes in multiple sclerosis

Multiple sclerosis (MS) is an inflammatory and neurodegenerative disease of the central nervous system (CNS), resulting in neurological disability that worsens over time. While progress has been made in defining the immune system's role in MS pathophysiology, the contribution of intrinsic CNS cell dysfunction remains unclear. Here, we generated a collection of induced pluripotent stem cell (iPSC) lines from people with MS spanning diverse clinical subtypes and differentiated them into glia-enriched cultures. Using single-cell transcriptomic profiling and orthogonal analyses, we observed several distinguishing characteristics of MS cultures pointing to glia-intrinsic disease mechanisms. We found that primary progressive MS-derived cultures contained fewer oligodendrocytes. Moreover, MS-derived oligodendrocyte lineage cells and astrocytes showed increased expression of immune and inflammatory genes, matching those of glia from MS postmortem brains. Thus, iPSC-derived MS models provide a unique platform for dissecting glial contributions to disease phenotypes independent of the peripheral immune system and identify potential glia-specific targets for therapeutic intervention.

Exploring CD26-/lo subpopulations of lymphocytes in asthma phenotype and severity: A novel CD4+ T cell subset expressing archetypical granulocyte proteins

BACKGROUND

Asthma pathology may induce changes in naïve/memory lymphocyte proportions assessable through the evaluation of surface CD26 (dipeptidyl peptidase 4/DPP4) levels. Our aim was to investigate the association of asthma phenotype/severity with the relative frequency of CD26-/lo, CD26int and CD26hi subsets within different lymphocyte populations.

METHODS

The proportion of CD26-/lo, CD26int and CD26hi subsets within CD4+ effector T cells (Teff), total CD4- lymphocytes, γδ-T cells, NK cells and NKT cells was measured in peripheral blood samples from healthy (N = 30) and asthma (N = 119) donors with different phenotypes/severities by flow cytometry. We performed K-means clustering analysis and further characterised the CD4+CD26-/lo Teff cell subset by LC-MS/MS and immunofluorescence.

RESULTS

Cluster analysis including clinical and flow cytometry data resulted in four groups, two of them with opposite inflammatory profiles (neutrophilic vs. eosinophilic). Neutrophilic asthma presented reduced CD4-CD26hi cells, which negatively correlated with systemic inflammation. Eosinophilic asthma displayed a general expansion of CD26-/lo subsets. Specifically, CD4+CD26-/lo Teff expansion was confirmed in asthma, especially in atopic patients. Proteomic characterisation of this subset with a TEM/TEMRA phenotype revealed upregulated levels of innate (e.g. MPO and RNASE2) and cytoskeleton/extracellular matrix (e.g. MMP9 and ACTN1) proteins. Immunofluorescence assays confirmed the presence of atypical proteins for CD4+ T cells, and an enrichment in 'flower-like' nuclei and MMP9/RNASE2 levels in CD4+CD26-/lo Teff compared to CD4+ T lymphocytes.

CONCLUSION

There is an association between CD26 levels in different lymphocyte subsets and asthma phenotype/severity. CD4+CD26-/loTEMRA cells expressing innate proteins specific to eosinophils/neutrophils could be determinant in sustaining long-term inflammation in adult allergic asthma.

An open-source framework for end-to-end analysis of electronic health record data

With progressive digitalization of healthcare systems worldwide, large-scale collection of electronic health records (EHRs) has become commonplace. However, an extensible framework for comprehensive exploratory analysis that accounts for data heterogeneity is missing. Here we introduce ehrapy, a modular open-source Python framework designed for exploratory analysis of heterogeneous epidemiology and EHR data. ehrapy incorporates a series of analytical steps, from data extraction and quality control to the generation of low-dimensional representations. Complemented by rich statistical modules, ehrapy facilitates associating patients with disease states, differential comparison between patient clusters, survival analysis, trajectory inference, causal inference and more. Leveraging ontologies, ehrapy further enables data sharing and training EHR deep learning models, paving the way for foundational models in biomedical research. We demonstrate ehrapy's features in six distinct examples. We applied ehrapy to stratify patients affected by unspecified pneumonia into finer-grained phenotypes. Furthermore, we reveal biomarkers for significant differences in survival among these groups. Additionally, we quantify medication-class effects of pneumonia medications on length of stay. We further leveraged ehrapy to analyze cardiovascular risks across different data modalities. We reconstructed disease state trajectories in patients with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) based on imaging data. Finally, we conducted a case study to demonstrate how ehrapy can detect and mitigate biases in EHR data. ehrapy, thus, provides a framework that we envision will standardize analysis pipelines on EHR data and serve as a cornerstone for the community.

FGF8-mediated gene regulation affects regional identity in human cerebral organoids

The morphogen FGF8 establishes graded positional cues imparting regional cellular responses via modulation of early target genes. The roles of FGF signaling and its effector genes remain poorly characterized in human experimental models mimicking early fetal telencephalic development. We used hiPSC-derived cerebral organoids as an in vitro platform to investigate the effect of FGF8 signaling on neural identity and differentiation. We found that FGF8 treatment increases cellular heterogeneity, leading to distinct telencephalic and mesencephalic-like domains that co-develop in multi-regional organoids. Within telencephalic regions, FGF8 affects the anteroposterior and dorsoventral identity of neural progenitors and the balance between GABAergic and glutamatergic neurons, thus impacting spontaneous neuronal network activity. Moreover, FGF8 efficiently modulates key regulators responsible for several human neurodevelopmental disorders. Overall, our results show that FGF8 signaling is directly involved in both regional patterning and cellular diversity in human cerebral organoids and in modulating genes associated with normal and pathological neural development.

Neuronal SLC39A8 deficiency impairs cerebellar development by altering manganese homeostasis

Solute carrier family 39, member 8 (SLC39A8), is a transmembrane transporter that mediates the cellular uptake of zinc, iron, and manganese (Mn). Human genetic studies document the involvement of SLC39A8 in Mn homeostasis, brain development, and function. However, the role and pathophysiological mechanisms of SLC39A8 in the central nervous system remain elusive. We generated Slc39a8 neuron-specific knockout (Slc39a8-NSKO) mice to study SLC39A8 function in neurons. The Slc39a8-NSKO mice displayed markedly decreased Mn levels in the whole brain and brain regions, especially the cerebellum. Radiotracer studies using 54Mn revealed that Slc39a8-NSKO mice had impaired brain uptake of Mn. Slc39a8-NSKO cerebellums exhibited morphological defects and abnormal dendritic arborization of Purkinje cells. Reduced neurogenesis and increased apoptotic cell death occurred in the cerebellar external granular layer of Slc39a8-NSKO mice. Brain Mn deficiency in Slc39a8-NSKO mice was associated with motor dysfunction. Unbiased RNA-Seq analysis revealed downregulation of key pathways relevant to neurodevelopment and synaptic plasticity, including cAMP signaling pathway genes. We further demonstrated that Slc39a8 was required for the optimal transcriptional response to the cAMP-mediated signaling pathway. In summary, our study highlighted the essential roles of SLC39A8 in brain Mn uptake and cerebellum development and functions.

The single-cell transcriptomic atlas iPain identifies senescence of nociceptors as a therapeutical target for chronic pain treatment

Chronic pain remains a significant medical challenge with complex underlying mechanisms, and an urgent need for new treatments. Our research built and utilized the iPain single-cell atlas to study chronic pain progression in dorsal root and trigeminal ganglia. We discovered that senescence of a small subset of pain-sensing neurons may be a driver of chronic pain. This mechanism was observed in animal models after nerve injury and in human patients diagnosed with chronic pain or diabetic painful neuropathy. Notably, treatment with senolytics, drugs that remove senescent cells, reversed pain symptoms in mice post-injury. These findings highlight the role of cellular senescence in chronic pain development, demonstrate the therapeutic potential of senolytic treatments, and underscore the value of the iPain atlas for future pain research.

Protocol to search for genetic factors related to severe COVID-19 by analyzing publicly available genome-wide association studies

COVID-19 casualties vary among different ancestral groups due to a variety of factors. Here, we present a protocol for analyzing publicly available genome-wide association studies (GWASs) to search for ancestry-specific genetic factors related to severe COVID-19. We describe steps for downloading and comparing two COVID-19 GWASs, calculating expression quantitative trait loci, and single-cell gene expression analysis. We demonstrate this approach using GWASs from Host Genetics Initiative; however, it is applicable to other databases such as the UK Biobank. For complete details on the use and execution of this protocol, please refer to Cheng et al.1.

CelltypeR: A flow cytometry pipeline to characterize single cells from brain organoids

Motivated by the cellular heterogeneity in complex tissues, particularly in brain and induced pluripotent stem cell (iPSC)-derived brain models, we developed a complete workflow to reproducibly characterize cell types in complex tissues. Our approach combines a flow cytometry (FC) antibody panel with our computational pipeline CelltypeR, enabling dataset aligning, unsupervised clustering optimization, cell type annotating, and statistical comparisons. Applied to human iPSC derived midbrain organoids, it successfully identified the major brain cell types. We performed fluorescence-activated cell sorting of CelltypeR-defined astrocytes, radial glia, and neurons, exploring transcriptional states by single-cell RNA sequencing. Among the sorted neurons, we identified subgroups of dopamine neurons: one reminiscent of substantia nigra cells most vulnerable in Parkinson's disease. Finally, we used our workflow to track cell types across a time course of organoid differentiation. Overall, our adaptable analysis framework provides a generalizable method for reproducibly identifying cell types across FC datasets in complex tissues.

Case report: Artificial thymic organoids facilitate clinical decisions for a patient with a TP63 variant and severe persistent T cell lymphopenia

Pathogenic variants in the transcription factor TP63 are associated with clinically overlapping syndromes including ectrodactyly-ectodermal dysplasia clefting (EEC) and ankyloblepharon-ectodermal defects-cleft lip/palate (AEC). T cell lymphopenia has rarely been described in individuals with TP63 variants and the cause of the T cell defect is unclear. Here, we present a case of a female infant born with TP63-related syndrome and profound T cell lymphopenia, first uncovered through newborn screening. Flow cytometry analysis revealed low CD4+ naïve T cells and nearly absent CD8+ T cells with intact B and NK cell compartments. A de novo heterozygous pathogenic variant c.1040 G>A (C347Y) in exon 8 of TP63 was identified. An artificial thymic organoid system, to assess the intrinsic ability of the patient's hematopoietic cells to develop into T cells, was performed twice using separate peripheral blood samples. Ex vivo T cell differentiation was evident with the artificial organoid system, suggesting that a thymic stromal cell defect may be the cause of the T cell lymphopenia. Consistent with this, interrogation of publicly available data indicated that TP63 expression in the human thymus is restricted to thymic epithelial cells. Based on these data, congenital athymia was suspected and the patient received an allogenic cultured thymus tissue implant (CTTI). This is the first report of suspected congenital athymia and attempted treatment with CTTI associated with TP63 variant. At 9 months post-implant, peripheral lymphocyte analysis revealed measurable T cell receptor excision circles and presence of CD4+ recent thymic emigrants suggestive of early thymopoiesis. She will continue regular monitoring to ensure restoration of T cell immunity.

A brain cell atlas integrating single-cell transcriptomes across human brain regions

While single-cell technologies have greatly advanced our comprehension of human brain cell types and functions, studies including large numbers of donors and multiple brain regions are needed to extend our understanding of brain cell heterogeneity. Integrating atlas-level single-cell data presents a chance to reveal rare cell types and cellular heterogeneity across brain regions. Here we present the Brain Cell Atlas, a comprehensive reference atlas of brain cells, by assembling single-cell data from 70 human and 103 mouse studies of the brain throughout major developmental stages across brain regions, covering over 26.3 million cells or nuclei from both healthy and diseased tissues. Using machine-learning based algorithms, the Brain Cell Atlas provides a consensus cell type annotation, and it showcases the identification of putative neural progenitor cells and a cell subpopulation of PCDH9high microglia in the human brain. We demonstrate the gene regulatory difference of PCDH9high microglia between hippocampus and prefrontal cortex and elucidate the cell-cell communication network. The Brain Cell Atlas presents an atlas-level integrative resource for comparing brain cells in different environments and conditions within the Human Cell Atlas.

Single-nucleus atlas of the Artemia female reproductive system suggests germline repression of the Z chromosome

Our understanding of the molecular pathways that regulate oogenesis and define cellular identity in the Arthropod female reproductive system and the extent of their conservation is currently very limited. This is due to the focus on model systems, including Drosophila and Daphnia, which do not reflect the observed diversity of morphologies, reproductive modes, and sex chromosome systems. We use single-nucleus RNA and ATAC sequencing to produce a comprehensive single nucleus atlas of the adult Artemia franciscana female reproductive system. We map our data to the Fly Cell Atlas single-nucleus dataset of the Drosophila melanogaster ovary, shedding light on the conserved regulatory programs between the two distantly related Arthropod species. We identify the major cell types known to be present in the Artemia ovary, including germ cells, follicle cells, and ovarian muscle cells. Additionally, we use the germ cells to explore gene regulation and expression of the Z chromosome during meiosis, highlighting its unique regulatory dynamics and allowing us to explore the presence of meiotic sex chromosome silencing in this group.

Single-cell multiregion dissection of Alzheimer's disease

Alzheimer's disease is the leading cause of dementia worldwide, but the cellular pathways that underlie its pathological progression across brain regions remain poorly understood1-3. Here we report a single-cell transcriptomic atlas of six different brain regions in the aged human brain, covering 1.3 million cells from 283 post-mortem human brain samples across 48 individuals with and without Alzheimer's disease. We identify 76 cell types, including region-specific subtypes of astrocytes and excitatory neurons and an inhibitory interneuron population unique to the thalamus and distinct from canonical inhibitory subclasses. We identify vulnerable populations of excitatory and inhibitory neurons that are depleted in specific brain regions in Alzheimer's disease, and provide evidence that the Reelin signalling pathway is involved in modulating the vulnerability of these neurons. We develop a scalable method for discovering gene modules, which we use to identify cell-type-specific and region-specific modules that are altered in Alzheimer's disease and to annotate transcriptomic differences associated with diverse pathological variables. We identify an astrocyte program that is associated with cognitive resilience to Alzheimer's disease pathology, tying choline metabolism and polyamine biosynthesis in astrocytes to preserved cognitive function late in life. Together, our study develops a regional atlas of the ageing human brain and provides insights into cellular vulnerability, response and resilience to Alzheimer's disease pathology.

Manipulating the EphB4-ephrinB2 axis to reduce metastasis in HNSCC

The EphB4-ephrinB2 signaling axis has been heavily implicated in metastasis across numerous cancer types. Our emerging understanding of the dichotomous roles that EphB4 and ephrinB2 play in head and neck squamous cell carcinoma (HNSCC) poses a significant challenge to rational drug design. We find that EphB4 knockdown in cancer cells enhances metastasis in preclinical HNSCC models by augmenting immunosuppressive cells like T regulatory cells (Tregs) within the tumor microenvironment. EphB4 inhibition in cancer cells also amplifies their ability to metastasize through increased expression of genes associated with epithelial mesenchymal transition and hallmark pathways of metastasis. In contrast, vascular ephrinB2 knockout coupled with radiation therapy (RT) enhances anti-tumor immunity, reduces Treg accumulation into the tumor, and decreases metastasis. Notably, targeting the EphB4-ephrinB2 signaling axis with the engineered EphB4 ligands EFNB2-Fc-His and Fc-TNYL-RAW-GS reduces local tumor growth and distant metastasis in a preclinical model of HNSCC. Our data suggest that targeted inhibition of vascular ephrinB2 while avoiding inhibition of EphB4 in cancer cells could be a promising strategy to mitigate HNSCC metastasis.

SANTO: a coarse-to-fine alignment and stitching method for spatial omics

With the flourishing of spatial omics technologies, alignment and stitching of slices becomes indispensable to decipher a holistic view of 3D molecular profile. However, existing alignment and stitching methods are unpractical to process large-scale and image-based spatial omics dataset due to extreme time consumption and unsatisfactory accuracy. Here we propose SANTO, a coarse-to-fine method targeting alignment and stitching tasks for spatial omics. SANTO firstly rapidly supplies reasonable spatial positions of two slices and identifies the overlap region. Then, SANTO refines the positions of two slices by considering spatial and omics patterns. Comprehensive experiments demonstrate the superior performance of SANTO over existing methods. Specifically, SANTO stitches cross-platform slices for breast cancer samples, enabling integration of complementary features to synergistically explore tumor microenvironment. SANTO is then applied to 3D-to-3D spatiotemporal alignment to study development of mouse embryo. Furthermore, SANTO enables cross-modality alignment of spatial transcriptomic and epigenomic data to understand complementary interactions.

An Adipo-Pulmonary Axis Mediated by FABP4 Hormone Defines a Therapeutic Target Against Obesity-Induced Airway Disease

Obesity-related airway disease is a clinical condition without a clear description and effective treatment. Here, we define this pathology and its unique properties, which differ from classic asthma phenotypes, and identify a novel adipo-pulmonary axis mediated by FABP4 hormone as a critical mediator of obesity-induced airway disease. Through detailed analysis of murine models and human samples, we elucidate the dysregulated lipid metabolism and immunometabolic responses within obese lungs, particularly highlighting the stress response activation and downregulation of surfactant-related genes, notably SftpC. We demonstrate that FABP4 deficiency mitigates these alterations, demonstrating a key role in obesity-induced airway disease pathogenesis. Importantly, we identify adipose tissue as the source of FABP4 hormone in the bronchoalveolar space and describe strong regulation in the context of human obesity, particularly among women. Finally, our exploration of antibody-mediated targeting of circulating FABP4 unveils a novel therapeutic avenue, addressing a pressing unmet need in managing obesity-related airway disease. These findings not only define the presence of a critical adipo-pulmonary endocrine link but also present FABP4 as a therapeutic target for managing this unique airway disease that we refer to as fatty lung disease associated with obesity.

One Sentence Summary

Investigating FABP4's pivotal role in obesity-driven airway disease, this study unveils an adipo-pulmonary axis with potential therapeutic implications.

Perimenopausal and Menopausal Mammary Glands In A 4-Vinylcyclohexene Diepoxide Mouse Model

As both perimenopausal and menopausal periods are recognized critical windows of susceptibility for breast carcinogenesis, development of a physiologically relevant model has been warranted. The traditional ovariectomy model causes instant removal of the entire hormonal repertoire produced by the ovary, which does not accurately approximate human natural menopause with gradual transition. Here, we characterized the mammary glands of 4-vinylcyclohexene diepoxide (VCD)-treated animals at different time points, revealing that the model can provide the mammary glands with both perimenopausal and menopausal states. The perimenopausal gland showed moderate regression in ductal structure with no responsiveness to external hormones, while the menopausal gland showed severe regression with hypersensitivity to hormones. Leveraging the findings on the VCD model, effects of a major endocrine disruptor (polybrominated diphenyl ethers, PBDEs) on the mammary gland were examined during and after menopausal transition, with the two exposure modes; low-dose, chronic (environmental) and high-dose, subacute (experimental). All conditions of PBDE exposure did not augment or compromise the macroscopic ductal reorganization resulting from menopausal transition and/or hormonal treatments. Single-cell RNA sequencing revealed that the experimental PBDE exposure during the post-menopausal period caused specific transcriptomic changes in the non-epithelial compartment such as Errfi1 upregulation in fibroblasts. The environmental PBDE exposure resulted in similar transcriptomic changes to a lesser extent. In summary, the VCD mouse model provides both perimenopausal and menopausal windows of susceptibility for the breast cancer research community. PBDEs, including all tested models, may affect the post-menopausal gland including impacts on the non-epithelial compartments.

Lysosomal endonuclease RNase T2 and PLD exonucleases cooperatively generate RNA ligands for TLR7 activation

Toll-like receptor 7 (TLR7) is essential for recognition of RNA viruses and initiation of antiviral immunity. TLR7 contains two ligand-binding pockets that recognize different RNA degradation products: pocket 1 recognizes guanosine, while pocket 2 coordinates pyrimidine-rich RNA fragments. We found that the endonuclease RNase T2, along with 5' exonucleases PLD3 and PLD4, collaboratively generate the ligands for TLR7. Specifically, RNase T2 generated guanosine 2',3'-cyclic monophosphate-terminated RNA fragments. PLD exonuclease activity further released the terminal 2',3'-cyclic guanosine monophosphate (2',3'-cGMP) to engage pocket 1 and was also needed to generate RNA fragments for pocket 2. Loss-of-function studies in cell lines and primary cells confirmed the critical requirement for PLD activity. Biochemical and structural studies showed that PLD enzymes form homodimers with two ligand-binding sites important for activity. Previously identified disease-associated PLD mutants failed to form stable dimers. Together, our data provide a mechanistic basis for the detection of RNA fragments by TLR7.

Variants of NAV3, a neuronal morphogenesis protein, cause intellectual disability, developmental delay, and microcephaly

Microtubule associated proteins (MAPs) are widely expressed in the central nervous system, and have established roles in cell proliferation, myelination, neurite formation, axon specification, outgrowth, dendrite, and synapse formation. We report eleven individuals from seven families harboring predicted pathogenic biallelic, de novo, and heterozygous variants in the NAV3 gene, which encodes the microtubule positive tip protein neuron navigator 3 (NAV3). All affected individuals have intellectual disability (ID), microcephaly, skeletal deformities, ocular anomalies, and behavioral issues. In mouse brain, Nav3 is expressed throughout the nervous system, with more prominent signatures in postmitotic, excitatory, inhibiting, and sensory neurons. When overexpressed in HEK293T and COS7 cells, pathogenic variants impaired NAV3 ability to stabilize microtubules. Further, knocking-down nav3 in zebrafish led to severe morphological defects, microcephaly, impaired neuronal growth, and behavioral impairment, which were rescued with co-injection of WT NAV3 mRNA and not by transcripts encoding the pathogenic variants. Our findings establish the role of NAV3 in neurodevelopmental disorders, and reveal its involvement in neuronal morphogenesis, and neuromuscular responses.

Reduced Nephrin Tyrosine Phosphorylation Enhances Insulin Secretion and Increases Glucose Tolerance With Age

BACKGROUND

Nephrin is a transmembrane protein with well-established signaling roles in kidney podocytes, and a smaller set of secretory functions in pancreatic β cells are implicated in diabetes. Nephrin signaling is mediated in part through its 3 cytoplasmic YDxV motifs, which can be tyrosine phosphorylated by high glucose and β cell injuries. Although in vitro studies demonstrate these phosphorylated motifs can regulate β cell vesicle trafficking and insulin release, in vivo evidence of their role in this cell type remains to be determined.

METHODS

To further explore the role of nephrin YDxV phosphorylation in β cells, we used a mouse line with tyrosine to phenylalanine substitutions at each YDxV motif (nephrin-Y3F) to inhibit phosphorylation. We assessed islet function via primary islet glucose-stimulated insulin secretion assays and oral glucose tolerance tests.

RESULTS

Nephrin-Y3F mice successfully developed pancreatic endocrine and exocrine tissues with minimal structural differences. Unexpectedly, male and female nephrin-Y3F mice showed elevated insulin secretion, with a stronger increase observed in male mice. At 8 months of age, no differences in glucose tolerance were observed between wild-type (WT) and nephrin-Y3F mice. However, aged nephrin-Y3F mice (16 months of age) demonstrated more rapid glucose clearance compared to WT controls.

CONCLUSION

Taken together, loss of nephrin YDxV phosphorylation does not alter baseline islet function. Instead, our data suggest a mechanism linking impaired nephrin YDxV phosphorylation to improved islet secretory ability with age. Targeting nephrin phosphorylation could provide novel therapeutic opportunities to improve β cell function.

Zebrafish do not have calprotectin

The protein heterodimer calprotectin and its component proteins, S100A8 and S100A9, play important antibacterial and proinflammatory roles in the mammalian innate immune response. Gaining mechanistic insights into the regulation and biological function of calprotectin will help facilitate patient diagnostics and therapy and further our understanding of the host-microbe interface. Recent literature has identified zebrafish s100a10b as zebrafish calprotectin based on sequence similarity, genomic context, and transcriptional upregulation during the immune response to bacterial infections. The field would benefit from expanding the breadth of calprotectin studies into a zebrafish innate immunity model. Here, we carefully evaluated the possibility that zebrafish possess a calprotectin. We found that zebrafish do not possess an ortholog of mammalian S100A8 or S100A9. We then identified four zebrafish s100 proteins- including s100a10b-that are expressed in immune cells and upregulated during the immune response. We recombinantly expressed and purified these proteins and measured the antimicrobial and proinflammatory characteristics. We found that none of the zebrafish proteins exhibited activity comparable to mammalian calprotectin. Our work demonstrates conclusively that zebrafish have no ortholog of calprotectin, and the most plausible candidate proteins have not convergently evolved similar functions.

Systematic single-cell analysis reveals dynamic control of transposable element activity orchestrating the endothelial-to-hematopoietic transition

BACKGROUND

The endothelial-to-hematopoietic transition (EHT) process during definitive hematopoiesis is highly conserved in vertebrates. Stage-specific expression of transposable elements (TEs) has been detected during zebrafish EHT and may promote hematopoietic stem cell (HSC) formation by activating inflammatory signaling. However, little is known about how TEs contribute to the EHT process in human and mouse.

RESULTS

We reconstructed the single-cell EHT trajectories of human and mouse and resolved the dynamic expression patterns of TEs during EHT. Most TEs presented a transient co-upregulation pattern along the conserved EHT trajectories, coinciding with the temporal relaxation of epigenetic silencing systems. TE products can be sensed by multiple pattern recognition receptors, triggering inflammatory signaling to facilitate HSC emergence. Interestingly, we observed that hypoxia-related signals were enriched in cells with higher TE expression. Furthermore, we constructed the hematopoietic cis-regulatory network of accessible TEs and identified potential TE-derived enhancers that may boost the expression of specific EHT marker genes.

CONCLUSIONS

Our study provides a systematic vision of how TEs are dynamically controlled to promote the hematopoietic fate decisions through transcriptional and cis-regulatory networks, and pre-train the immunity of nascent HSCs.