Epididymis cell atlas in a patient with a sex development disorder and a novel NR5A1 gene mutation / 亚洲男科学杂志(英文版)

This study aims to characterize the cell atlas of the epididymis derived from a 46,XY disorders of sex development (DSD) patient with a novel heterozygous mutation of the nuclear receptor subfamily 5 group A member 1 (NR5A1) gene. Next-generation sequencing found a heterozygous c.124C>G mutation in NR5A1 that resulted in a p.Q42E missense mutation in the conserved DNA-binding domain of NR5A1. The patient demonstrated feminization of external genitalia and Tanner stage 1 breast development. The surgical procedure revealed a morphologically normal epididymis and vas deferens but a dysplastic testis. Microfluidic-based single-cell RNA sequencing (scRNA-seq) analysis found that the fibroblast cells were significantly increased (approximately 46.5%), whereas the number of main epididymal epithelial cells (approximately 9.2%), such as principal cells and basal cells, was dramatically decreased. Bioinformatics analysis of cell-cell communications and gene regulatory networks at the single-cell level inferred that epididymal epithelial cell loss and fibroblast occupation are associated with the epithelial-to-mesenchymal transition (EMT) process. The present study provides a cell atlas of the epididymis of a patient with 46,XY DSD and serves as an important resource for understanding the pathophysiology of DSD.

Single-cell analysis reveals an Angpt4-initiated EPDC-EC-CM cellular coordination cascade during heart regeneration

Mammals exhibit limited heart regeneration ability, which can lead to heart failure after myocardial infarction. In contrast, zebrafish exhibit remarkable cardiac regeneration capacity. Several cell types and signaling pathways have been reported to participate in this process. However, a comprehensive analysis of how different cells and signals interact and coordinate to regulate cardiac regeneration is unavailable. We collected major cardiac cell types from zebrafish and performed high-precision single-cell transcriptome analyses during both development and post-injury regeneration. We revealed the cellular heterogeneity as well as the molecular progress of cardiomyocytes during these processes, and identified a subtype of atrial cardiomyocyte exhibiting a stem-like state which may transdifferentiate into ventricular cardiomyocytes during regeneration. Furthermore, we identified a regeneration-induced cell (RIC) population in the epicardium-derived cells (EPDC), and demonstrated Angiopoietin 4 (Angpt4) as a specific regulator of heart regeneration. angpt4 expression is specifically and transiently activated in RIC, which initiates a signaling cascade from EPDC to endocardium through the Tie2-MAPK pathway, and further induces activation of cathepsin K in cardiomyocytes through RA signaling. Loss of angpt4 leads to defects in scar tissue resolution and cardiomyocyte proliferation, while overexpression of angpt4 accelerates regeneration. Furthermore, we found that ANGPT4 could enhance proliferation of neonatal rat cardiomyocytes, and promote cardiac repair in mice after myocardial infarction, indicating that the function of Angpt4 is conserved in mammals. Our study provides a mechanistic understanding of heart regeneration at single-cell precision, identifies Angpt4 as a key regulator of cardiomyocyte proliferation and regeneration, and offers a novel therapeutic target for improved recovery after human heart injuries.

PathogenTrack and Yeskit: tools for identifying intracellular pathogens from single-cell RNA-sequencing datasets as illustrated by application to COVID-19 / 医学前沿

Pathogenic microbes can induce cellular dysfunction, immune response, and cause infectious disease and other diseases including cancers. However, the cellular distributions of pathogens and their impact on host cells remain rarely explored due to the limited methods. Taking advantage of single-cell RNA-sequencing (scRNA-seq) analysis, we can assess the transcriptomic features at the single-cell level. Still, the tools used to interpret pathogens (such as viruses, bacteria, and fungi) at the single-cell level remain to be explored. Here, we introduced PathogenTrack, a python-based computational pipeline that uses unmapped scRNA-seq data to identify intracellular pathogens at the single-cell level. In addition, we established an R package named Yeskit to import, integrate, analyze, and interpret pathogen abundance and transcriptomic features in host cells. Robustness of these tools has been tested on various real and simulated scRNA-seq datasets. PathogenTrack is competitive to the state-of-the-art tools such as Viral-Track, and the first tools for identifying bacteria at the single-cell level. Using the raw data of bronchoalveolar lavage fluid samples (BALF) from COVID-19 patients in the SRA database, we found the SARS-CoV-2 virus exists in multiple cell types including epithelial cells and macrophages. SARS-CoV-2-positive neutrophils showed increased expression of genes related to type I interferon pathway and antigen presenting module. Additionally, we observed the Haemophilus parahaemolyticus in some macrophage and epithelial cells, indicating a co-infection of the bacterium in some severe cases of COVID-19. The PathogenTrack pipeline and the Yeskit package are publicly available at GitHub.

Single-cell RNA Sequencing Reveals Thoracolumbar Vertebra Heterogeneity and Rib-genesis in Pigs / 基因组蛋白质组与生物信息学报·英文版

Development of thoracolumbar vertebra (TLV) and rib primordium (RP) is a common evolutionary feature across vertebrates, although whole-organism analysis of the expression dynamics of TLV- and RP-related genes has been lacking. Here, we investigated the single-cell transcriptome landscape of thoracic vertebra (TV), lumbar vertebra (LV), and RP cells from a pig embryo at 27 days post-fertilization (dpf) and identified six cell types with distinct gene expression signatures. In-depth dissection of the gene expression dynamics and RNA velocity revealed a coupled process of osteogenesis and angiogenesis during TLV and RP development. Further analysis of cell type-specific and strand-specific expression uncovered the extremely high level of HOXA10 3'-UTR sequence specific to osteoblasts of LV cells, which may function as anti-HOXA10-antisense by counteracting the HOXA10-antisense effect to determine TLV transition. Thus, this work provides a valuable resource for understanding embryonic osteogenesis and angiogenesis underlying vertebrate TLV and RP development at the cell type-specific resolution, which serves as a comprehensive view on the transcriptional profile of animal embryo development.