Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “transcriptome sequencing”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 235 records · Page 13Linked to original sources
PubMed2026

Exploring the transcriptional crosstalk between adipose tissue and locoregional recurrence in breast cancer using independent component analysis.

Locoregional recurrence (LRR) poses a persistent clinical challenge in breast cancer, with emerging evidence implicating the tumor-associated adipose tissue in modulating recurrence risk. This study investigates shared transcriptional programs between adipose tissue and breast tumors and examines their association with disease-free survival (DFS), particularly in the context of reconstructive surgery where adipose tissue from different body compartments are commonly used. We analyzed bulk gene expression data from 5,691 breast tumors and 978 human adipose tissue samples from different body compartments using consensus-independent component analysis (c-ICA) to identify transcriptional components (TCs). Gene set enrichment analysis (GSEA) and copy number alteration profiling were used for biological annotation. Associations between TCs and DFS were evaluated through univariate Cox regression. Key findings were validated using spatial transcriptomic and single-cell RNA sequencing datasets. Among the 411 TCs identified, 332 showed biological enrichment, and 35 were significantly associated with DFS. Four DFS-associated TCs (TC257, TC350, TC371, TC400) were enriched for adipogenesis-related genes and exhibited heightened activity in high-grade, triple-negative tumors and in patients with elevated BMI. Notably, TC350 was highly active in adipose tissue from common reconstructive donor sites (abdomen, omentum, subcutis) but not in native breast adipose tissue. Spatial transcriptomic and single-cell analyses confirmed the increased activity of these adipogenesis-related TCs in tumor regions and adipose cells. TC350 included FABP4, a gene previously linked to poor prognosis in breast cancer and considered as a potential new therapeutic target. Adipose tissue-derived transcriptional programs influence breast cancer prognosis and this seems to differ by tissue origin. These findings generate a hypothesis that donor site selection for adipose tissue in reconstructive surgery may impact LRR risk through adipogenesis-associated mechanisms. Further research is warranted to elucidate the biological and clinical implications of adipose-tumor transcriptional interactions.

Humans↗
PubMed2025

Analysis of Leishbuviridae from Trypanosomatids.

Over the last decade, considerable progress has been made in unraveling RNA virus diversity. This has contributed to our understanding of the evolution of these viruses, which include emerging zoonotic human pathogens. Current success has been greatly facilitated by the development of next-generation sequencing platforms instrumental for meta-transcriptomic studies. However, due to the rapid evolution of RNA viruses, there are numerous "blind spots" waiting to be explored; one of those is the RNA virome of unicellular eukaryotes. Here, we present the pipeline, which has been successfully used to characterize various types of RNA viruses, including Leishbuviridae (Bunyaviricetes, Hareavirales) in the parasitic flagellates of the family Trypanosomatidae. The pipeline relies on axenic in vitro cell culture and double-stranded RNA enrichment, followed by direct RNA-sequencing. A detailed procedure description starting from the initial total RNA preparation to the final assembly of the viral segments is provided.

High-Throughput Nucleotide Sequencing↗
PubMed2026

Multimodal Analysis Reveals Aberrant Expression of SUMO2 and Its Significant Association With Key Mechanisms of Metabolic Pathways in Hepatocellular Carcinoma.

BACKGROUND: Hepatocellular carcinoma (HCC) is the third leading cause of cancer-related deaths worldwide. However, the role of small ubiquitin-like modifier 2 (SUMO2), a core member of the small ubiquitin-like modifier (SUMO) family, regarding its expression patterns and metabolism-related functions in HCC remains inadequately understood. METHODS: A multidimensional analytical framework was applied, integrating immunohistochemistry (153 HCC vs. 21 non-HCC samples), proteomics (159 paired samples), bulk transcriptomics (3240 HCC vs. 2267 non-HCC samples), single-cell RNA sequencing (RNA-seq) (10 HCC vs. 8 non-HCC samples), spatial transcriptomics, and external CRISPR/Cas9 functional genomics data. Systematic analyses included standardized mean difference (SMD), pathway enrichment, pseudotime trajectory inference, in silico knockout, cell-cell communication, metabolic flux scoring, immune infiltration, clinical correlation, drug sensitivity prediction, and molecular docking. RESULTS: At the protein level, immunohistochemistry (nuclear positivity) and external proteomic data collectively demonstrated consistent SUMO2 overexpression in HCC. Consistent upregulation was also observed at the mRNA level across large-scale cohorts. Single-cell RNA-seq and spatial transcriptomics localized SUMO2 enrichment to malignant hepatocytes and tumor-dominant regions. CRISPR-mediated SUMO2 knockout suppressed proliferation in multiple HCC cell lines. Mechanistically, high SUMO2 expression was significantly associated with metabolic reprogramming involving glycolysis/gluconeogenesis, pyruvate metabolism, and the tricarboxylic acid cycle. SUMO2-high malignant hepatocyte subpopulations exhibited enhanced activity of the macrophage migration inhibitory factor signaling axis and enhanced iron-sensor interactions. Further, the immune infiltration analysis revealed a negative correlation between SUMO2 expression and M1 macrophages and a positive correlation between follicular helper T cells and regulatory T cells. Clinically, elevated SUMO2 levels were found to be associated with adverse prognostic features. Furthermore, high SUMO2 expression was associated with increased sensitivity to dasatinib, and molecular docking simulations predicted potential binding between SUMO2 and dasatinib, with a Vina score of -8.5 kcal/mol. CONCLUSIONS: SUMO2 is aberrantly expressed at the protein, mRNA, single-cell, and spatial transcriptomic levels in HCC and is significantly associated with metabolic reprogramming and altered migration inhibitory factor (MIF)-mediated intercellular communication, suggesting its potential as a novel biomarker for diagnosis and treatment.

Humans↗
PubMed2026

Comprehensive Microbiome Analyses for Regenerative Endodontic Therapy.

INTRODUCTION: Comprehensive microbiome analyses include the study of all microbial taxa, including bacteria, archaea, viruses, and fungi, as well as their functional activities and antibiotic resistance gene expression. Regenerative endodontic therapy presents a clinical situation where the most effective antimicrobial approaches are needed in order to ensure clinical success. METHODS AND RESULTS: In this paper, different contemporary technologies for the identification of endodontic microorganisms, such as with next generation sequencing, and their functional characterization, such as with whole genome sequencing, are described. The role of transcriptomics, as well as resistome analysis, are also discussed. Furthermore, the manner in which all this work and knowledge could be incorporated into clinical endodontics in general, and regenerative endodontic therapy as a special treatment, is outlined. CONCLUSIONS: Comprehensive microbiome analysis can lead to the development of more effective and personalized antimicrobial treatment.

Endodontics↗
PubMed2026

Malignant epithelial states drive immune dysfunction in ampulla of Vater carcinoma.

BACKGROUND: Ampulla of Vater (AoV) carcinoma is a rare malignancy arising at the junction of intestinal and pancreatobiliary epithelium. Its heterogeneous clinical behavior and histological diversity have hindered therapeutic advances, and the cellular basis of this heterogeneity remains unclear. We aimed to construct a single-cell transcriptomic atlas of AoV carcinoma, with a focus on identifying epithelial subtypes and their interactions with the tumor microenvironment (TME). METHODS: We performed single-cell RNA sequencing on eight primary AoV tumors and four matched normal tissues. Comprehensive clustering and transcriptomic analyses identified cell-type composition, epithelial heterogeneity, and tumor-immune interactions. Findings were validated using deconvolution of bulk RNA-seq data from 62 AoV carcinoma patients. Results Malignant epithelial cells were categorized into four distinct subtypes: Int-Wnt, PB-KRAS, Int-Hypoxia, and Cycling stage. PB-KRAS cells exhibited stem-like transcriptional programs and high genomic instability. Deconvolution analysis of bulk RNA-seq data from the independent AoV cohort revealed that enrichment of the PB-KRAS subtype correlated with tumor recurrence and poor survival. Our immune profiling analysis discovered a significant association between PB-KRAS subtype and GZMK+ CD8+ T cells, which are in a pre-dysfunctional state, alongside SPP1+ macrophages exhibiting immunosuppressive traits. Spatial transcriptome data further supports the immunosuppressive natures of TME around PB-KRAS subtype malignant epithelial cells in AoV carcinoma. CONCLUSIONS: Our study presents a single-cell atlas of AoV carcinoma, highlighting the molecular diversity of malignant epithelium and its association with the immune microenvironment. The PB-KRAS subtype emerges as a stem-like, immunosuppressive tumor state associated with poor prognosis, providing insights for future therapeutic targeting.

Ampulla of Vater carcinoma↗
PubMed2025

Multi-omics analyses reveal DjTcf4 critical for proper timing of differentiation in planarian regeneration.

The blastema is key to forming complete tissues in regenerating Dugesia japonica (D. japonica). However, the dynamic changes in cellular compositions and transcription landscapes in blastema during regeneration are understudied. Here, through genome reannotation, 3D spatial transcriptome construction, single-cell RNA sequencing (scRNA-seq), and single-cell assay for transposase-accessible chromatin sequencing (scATAC-seq) analyses of changes in gene expression and chromatin structures, we delineate key transcription factors regulating the developmental trajectories of major cell clusters in the regenerating head. Importantly, we find that the T cell factor 4 (DjTcf4)-positive cells highly accumulate at wound areas, and its gene network is critical for the proper timing of development during regeneration in multiple progenitor cells. Depletion of DjTcf4 and its target genes leads to singular eye and/or dull tail phenotypes and delays regeneration. Taken together, we build multi-omics atlases in D. japonica and reveal the noncanonical function of the DjTcf4 network in developmental pattern formation, laying a foundation for studies of regeneration in D. japonica.

Animals↗
PubMed2026

Nested co-expression network analysis identifies compact gene clusters in a black box.

MOTIVATION: Digital analysis of biological systems requires methods capable of identifying both broad and nested gene modules reflecting complex biological processes. Existing transcriptomic methods often miss compact gene sets corresponding to subprocesses in specialized cell types, limiting insights into functional heterogeneity. RESULTS: We present Nested-WGCNA, a two-stage unsupervised network analysis algorithm designed to identify coarse-grained and fine-grained gene modules. Applied to bulk RNA-Seq data, Nested-WGCNA reveals stable modules reproducible across datasets. When validated against scRNA-Seq data, these modules correspond to both major and minor immune cell subtypes. Application to immunotherapy response datasets uncovers predictive and prognostic biomarkers, highlighting its utility in treatment stratification and biomarker discovery. AVAILABILITY: The NestedWGCNA source code and analysis pipeline are available on GitHub (https://github.com/ilyada/NestedWGCNA) and archived on Zenodo (https://doi.org/10.5281/zenodo.18959244).

Algorithms↗
PubMed2026

Single-cell and spatial transcriptomics define a progenitor subpopulation and fibroinflammatory niche at the leading edge of parathyroid carcinoma.

Parathyroid carcinoma (PC) is a rare but clinically aggressive endocrine malignancy with limited treatment options and a poorly defined tumor microenvironment (TME). To elucidate its cellular heterogeneity and spatial architecture, we integrated single-cell and spatial transcriptomic profiling with whole-exome sequencing and multiplex immunohistochemistry on eight parathyroid neoplasm specimens, including PC, parathyroid adenoma, and atypical parathyroid tumor. We identified a distinct progenitor-like endocrine subpopulation (Ca-1) enriched in CDC73-mutant PC, exhibiting stem-like properties, elevated cell cycle activity, and pronounced genomic instability. Spatial mapping revealed that Ca-1 cells preferentially localize at the leading edge, forming a fibroinflammatory niche characterized by the enrichment of inflammatory cancer-associated fibroblasts (iCAFs) and SPP1+ macrophages. Within this niche, the dipeptidyl peptidase 4 (DPP4) is selectively expressed in Ca-1 cells and iCAFs, implicating a potential paracrine axis driving stromal remodeling and immunosuppression. These findings suggest that a spatially organized ecosystem may promote PC progression through TME remodeling and highlight the DPP4-CXCL2 axis as a candidate pathway for future investigation in aggressive parathyroid neoplasms.

Humans↗
PubMed2026

Multimodal computational framework resolves B cell maturation in autoimmunity and ageing.

Identification of the origin of pathogenic immune cells is crucial for therapeutic interventions and diagnosis but pseudotime methods struggle to trace immune cells accurately. Current trajectory inference methods for B cell development and response in health and disease either ignore or underutilize antigen receptor sequence information, limiting their ability to resolve developmental pathways, particularly for pathogenic populations. Widely used methods such as Monocle 3 reconstruct developmental paths from transcriptomic similarity alone, discarding the features from immune receptors. Dandelion has combined the immune receptor features with transcriptomics but it struggles to simulate the trajectory path of B cells. Here we present ClonoTrace, a computational framework that integrates BCR sequence features with transcriptomic trajectory inference through gated fusion of multimodal embeddings. In fetal B cell development and germinal centre development, ClonoTrace demonstrates closer concordance with the canonical reference ordering than Monocle 3 and Dandelion. Applied to systemic lupus erythematosus, ClonoTrace indicates a memory B cell extrafollicular maturation route alongside the naïve B cell route, accompanied by induction of ZEB2 with a concomitant decline of BACH2 along the trajectory, as a candidate alternative route to pathogenic double negative 2 B cells (DN2) in systemic lupus erythematosus (SLE) patients. In healthy ageing, ClonoTrace resolved three candidate age-related B cell maturation routes, from naïve, IgM+ memory and switched-memory B cells, each passing through a DN2-associated transcriptional state that is ordered before age-associated B cells along the inferred trajectory. ClonoTrace's fate probability algorithm indicated that IgM+ memory B cell to ABC transition as the leading candidate age-associated transition, which may be distinct from SLE DN2 maturation. ClonoTrace provides a generalizable framework for receptor-informed trajectory inference, describing candidate developmental routes of pathogenic B cell populations in autoimmunity and ageing.

Humans↗
PubMed2026

Antennal transcriptome analysis of chemosensory proteins in the raspberry weevil, Aegorhinus superciliosus (Coleoptera: Curculionidae).

Aegorhinus superciliosus (Coleoptera: Curculionidae) is a polyphagous pest of economic importance in southern Chile, the chemical ecology of which remains poorly characterized. Across insect species, chemosensory proteins, including odorant receptors (ORs), gustatory receptors (GRs), ionotropic receptors (IRs), odorant-binding proteins (OBPs), chemosensory proteins (CSPs), and sensory neuron membrane proteins (SNMPs), mediate the detection of chemical cues involved in host selection, reproduction, and other ecologically relevant behaviors. In this study, the antennal transcriptome of adult A. superciliosus was sequenced and analyzed using a de novo RNA-seq approach. Three independent biological replicates per sex were used for RNA-seq, and the same number of independent biological replicates was used for RT-qPCR validation; sequencing yielded 147,409,936 high-quality reads after quality filtering. A total of 112 candidate chemosensory genes were identified, comprising 43 ORs, 34 OBPs, 10 CSPs, 18 IRs, 5 GRs, and 2 SNMPs. Phylogenetic analyses assigned these candidate proteins to established clades, providing a comparative framework for functional inference for ORs and OBPs. Sex- and tissue-biased expression analyses revealed that several ORs, including AsupOR4, AsupOR19, and AsupOBP13, exhibit antennal enrichment and sex-specific expression patterns. Notably, AsupOR19 and AsupOBP13 displayed strong female-biased expression. In addition, transcripts of selected ORs and OBPs were detected in non-antennal tissues, such as the rostrum and legs, suggesting potential functional versatility beyond canonical olfaction. Together, these findings represent the first molecular identification of the chemosensory repertoire of A. superciliosus. This study establishes a foundation for reverse chemical ecology approaches aimed at identifying behaviorally active volatile organic compounds (VOCs) toward environmentally sustainable strategies for integrated pest management.

Animals↗
PubMed2025

CountASAP: a lightweight, easy to use python package for processing ASAPseq data.

BACKGROUND: Declining sequencing costs coupled with the increasing availability of easy-to-use kits for the isolation of DNA and RNA transcripts from single cells have driven a rapid proliferation of studies centered around genomic and transcriptomic data. Simultaneously, a wealth of new techniques have been developed that utilize single cell technologies to interrogate a broad range of cell-biological processes. One recently developed technique, transposase-accessible chromatin with sequencing (ATAC) with select antigen profiling by sequencing (ASAPseq), provides a combination of chromatin accessibility assessments with measurements of cell-surface marker expression levels. While software exists for the characterization of these datasets, there currently exists no tool explicitly designed to reformat ASAP surface marker FASTQ data into a count matrix which can then be used for these downstream analyses. RESULTS: To address this lack of a dedicated tool for ASAPseq data processing, we created CountASAP, an easy-to-use Python package purposefully designed to transform FASTQ files from ASAP experiments into count matrices compatible with commonly-used downstream bioinformatic analysis packages. CountASAP takes advantage of the independence of the relevant data structures to perform fully parallelized matches of each sequenced read to user-supplied input ASAP oligos and unique cell-identifier sequences. We directly compare the performance and user-friendliness of CountASAP to existing tools using similarly-structured data from a more common sequencing experiment: cellular indexing of transcriptomes and epitopes by sequencing (CITEseq). Further benchmarking against existing tools helps to identify proper defaults for CountASAP and assess the agreement of outputs from all tested software. A final test using a novel ASAPseq dataset provides evidence that CountASAP can generate biologically meaningful results that correlate well with paired chromatin accessibility data. CONCLUSIONS: CountASAP shows good agreement with existing, well-tested data processing tools in the analysis of similarly-structured benchmarking data. CountASAP runs efficiently on a standard laptop, has user-friendly documentation, a one-step installation, and represents the first and only tool designed specifically for the processing of ASAPseq data.

Software↗
PubMed2025

An Immunosenescent CD8+ T Cell Subset in Patients with Axial Spondyloarthritis and Psoriatic Arthritis Links Spontaneous Motility to Telomere Shortening and Dysfunction.

OBJECTIVE: A pathogenetic role of CD8+ T lymphocytes in radiographic axial spondyloarthritis (r-axSpA) and other spondyloarthritis (SpA) is sustained by genome-wide association studies and by the expansion of public T cell clonotypes in the target tissues. This study investigates the migration of CD8+ T cells along with their phenotype and functions in patients with r-axSpA and psoriatic arthritis (PsA). METHODS: Peripheral blood CD8+ and CD4+ T cells were isolated from patients with r-axSpA (n = 128), PsA (n = 60), and rheumatoid arthritis (RA) (n = 74) and healthy donors (HDs) (n = 79). Transwell migration assay was performed in the presence of different chemokines. CD8+ T cell immunoprofiling and effector functions were assessed by multiparametric flow cytometry. Transcriptome signature was evaluated by RNA sequencing analysis, whereas telomere length and dysfunction were measured by reverse transcriptase-polymerase chain reaction and immunofluorescence-fluorescence in situ hybridization, respectively. RESULTS: A significantly higher number of CD8+ T cells migrating in the absence of chemokine stimuli was found in patients with SpA compared with HDs and patients with RA. This subset, producing cytotoxic (granzyme B, perforin, granulysin) and proinflammatory molecules (tumor necrosis factor), was significantly enriched in terminally differentiated (CCR7-CD45RA+) and senescent (CD28-CD57+) cells having a gene expression profile characterized by cytolytic signature and natural killer markers. Remarkably, these spontaneously migrating CD8+ T cells showed DNA damage response activation, telomere shortening, and dysfunction. CONCLUSION: These data describe a terminally differentiated CD8+ T cell subset with a senescent and cytotoxic/proinflammatory profile and an intrinsic invasive potential enriched in patients with SpA that represents a possible player in disease pathogenesis.

Humans↗
PubMed2026

Contact hypersensitivity promotes hair regeneration through SPP1-secreting macrophages.

Allergic contact dermatitis, or contact hypersensitivity (CHS), is a pathological adaptive immune response that paradoxically induces hair regeneration, yet its underlying mechanisms remain unclear. We integrated high-resolution spatial transcriptomics and single-cell RNA sequencing to map the intricate interactions between immune cells, stroma, and hair follicles during CHS-induced hair growth in mice. Among all immunocytes, macrophages underwent the most prominent compositional and functional remodeling. We resolved five transcriptionally distinct macrophage subsets, with contact hypersensitivity driving a shift from homeostatic antigen-presenting cells toward a pro-inflammatory CD14+SPP1+ population. Trajectory analysis revealed divergent differentiation paths under homeostatic versus allergic conditions, highlighting the plasticity of skin macrophages. Mechanistically, CD14+SPP1+ macrophages secreted SPP1 (osteopontin), which engaged CD44 on hair follicle stem cells to activate PI3K-AKT signaling and trigger their proliferation. Notably, canonical pro-inflammatory cytokine signaling through TNF-α and IL-1 was dispensable for this process, underscoring the specificity of the SPP1-CD44 axis in immune-mediated hair regeneration. These findings reveal a macrophage-dependent mechanism of immune-mediated hair regeneration, offering therapeutic insights into immune-stem cell crosstalk.

Journal Article↗
PubMed2026

Transmitochondrial pigs reveal causal effects of mitochondrial DNA on backfat thickness via nuclear epigenetic reprogramming.

Mitochondrial DNA (mtDNA) polymorphisms have been associated with production traits in farm animals, including backfat thickness in pigs, yet direct in vivo evidence establishing a causal link between specific mtDNA haplotypes and fat deposition remains limited. In this study, we generated transmitochondrial pigs (mitopigs) by combining the Dapulian nuclear genome with Wuzhishan mtDNA via somatic cell nuclear transfer, introducing 23 mtDNA mutations relative to controls. Mitopigs exhibited significantly increased backfat thickness at 5 months, a difference that persisted in their offspring, without significant differences in body weight, body size, or litter size. Fibroblasts derived from mitopigs exhibited reduced mtDNA copy numbers, decreased expression of mitochondrial biogenesis genes (PPARA, PPARGC1A, RRM2B, and LRPPRC), impaired mitochondrial respiration, elevated reactive oxygen species (ROS), and upregulated adipogenic transcription factors (CEBPA, CEBPB, and PPARG). Consistent with these fibroblast findings, backfat tissue of mitopigs showed corresponding upregulation of adipogenic transcription factors and downregulation of mitochondrial biogenesis genes. Integrated transcriptomic and whole-genome bisulfite sequencing (WGBS) analyses revealed nuclear transcriptional reprogramming that was closely associated with differential DNA methylation, predominantly affecting mitochondrial function and lipid metabolism pathways. Mitopig fibroblasts also showed a pro-inflammatory response to lipopolysaccharide stimulation, with elevated expression of IL-12, NOS2, RELA, and TNF-α. Our findings provide direct in vivo evidence that mtDNA variants regulate adiposity in pigs through mitochondrial dysfunction, oxidative stress, and nuclear epigenetic modulation, highlighting the potential for incorporating mtDNA haplotype information into pig breeding programs as a complementary strategy to nuclear genomic selection.

Adipogenesis↗
PubMed2024

Implications of noncoding regulatory functions in the development of insulinomas.

Insulinomas are rare neuroendocrine tumors arising from pancreatic β cells, characterized by aberrant proliferation and altered insulin secretion, leading to glucose homeostasis failure. With the aim of uncovering the role of noncoding regulatory regions and their aberrations in the development of these tumors, we coupled epigenetic and transcriptome profiling with whole-genome sequencing. As a result, we unraveled somatic mutations associated with changes in regulatory functions. Critically, these regions impact insulin secretion, tumor development, and epigenetic modifying genes, including polycomb complex components. Chromatin remodeling is apparent in insulinoma-selective domains shared across patients, containing a specific set of regulatory sequences dominated by the SOX17 binding motif. Moreover, many of these regions are H3K27me3 repressed in β cells, suggesting that tumoral transition involves derepression of polycomb-targeted domains. Our work provides a compendium of aberrant cis-regulatory elements affecting the function and fate of β cells in their progression to insulinomas and a framework to identify coding and noncoding driver mutations.

Humans↗
PubMed2025

UHRF1 restricts HCoV-229E infection through epigenetic silencing of the viral receptor APN.

The emergence of SARS-CoV-2 has posed significant threats to global health, particularly for the older population. Similarly, common human coronaviruses, such as HCoV-229E, which typically cause mild cold-like symptoms, can lead to severe diseases, underscoring the need to understand virus-host interactions and identify host factors contributing to viral pathogenesis and disease progression. In this study, we perform a genome-wide CRISPR knockout screen using HCoV-229E and identify UHRF1 as a potent restriction factor. Mechanistically, UHRF1 suppresses HCoV-229E infection by downregulating the expression of its cell entry receptor, APN, through promoter hypermethylation. Focused CRISPR activation screens of UHRF1-downregulated genes confirm the critical role of APN in HCoV-229E infection and identify additional genes (e.g., SIGLEC1, PLAC8, and heparan sulfate biosynthesis genes) contributing to the restrictive functions of UHRF1. Transcriptomic and single-cell RNA sequencing analysis reveal that UHRF1 expression decreases with age, negatively correlating with increased APN expression. This age-related decline in UHRF1 is validated in primary alveolar macrophages from elderly individuals, which exhibit heightened susceptibility to HCoV-229E compared to those from younger individuals. Our findings highlight UHRF1 as a key age-related host defense factor against coronavirus and provide insights into the epigenetic regulation of viral entry receptors.

Animals↗
PubMed2025

The offonome reveals on and off states of gene expression near the detection limit of RNA-seq.

RNA-seq, widely used for gene expression profiling, provides nucleotide level genome coverage and summary gene expression values. Generally, low-expressed genes are ignored due to their unfavorable signal-to-noise ratio, however, these genes may offer crucial information, such as detecting rare cells in bulk tissues. In this study, we applied an approach that transforms the expression levels of low-expressed genes into a robust dichotomized on/off state by leveraging similarities in transcript coverage shape. Applied to three human cancer cohorts from the Cancer Genome Atlas (TCGA), chosen based on tissue morphology and anatomic site, we identified genes, the "offonome" near the detection limit, consistently or occasionally off across samples. Genes in the offonome spectrum proved useful for supervised and unsupervised applications, including characterizing oncogenic pathways, and identifying rare populations of cells in bulk tissue. Interrogating the offonome is relevant to bulk tumor analyses like TCGA, potentially expediting gene investigation in low-input situations like single cell RNA-seq.

Humans↗
PubMed2026

Spatially guided in vivo single-cell functional genomics of postnatal heart.

Understanding how spatial organization and cell-cell interactions shape gene regulatory programs is central to decoding tissue development and function. The transition at birth, marked by increased circulatory demands and rapid tissue growth, requires precise spatiotemporal coordination of cardiac maturation. In this study, we generated a high-resolution spatial and temporal atlas of the postnatal mouse heart by integrating single-nucleus RNA sequencing with image-based spatial transcriptomics. This framework revealed dynamic cellular interactions, niche-specific signaling and transcriptional programs guiding cardiomyocyte maturation. To functionally test prioritized regulators in vivo and at scale, we developed PIP-seq (probe-based indel-detectable Perturb-seq), a high-throughput platform that detects single guide RNA identity, infers gene editing and profiles transcription from fixed nuclei. Applying PIP-seq to the developing postnatal heart, we identified 21 previously uncharacterized regulators of cardiomyocyte maturation, including genes essential for sarcomere assembly, metabolic reprogramming and electrophysiological transitions. Together, our findings define how microenvironmental signals and intrinsic gene programs cooperate to guide heart maturation and establish a broadly applicable framework for functional genomics in complex tissues.

Animals↗