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 1,045 records · Page 58Linked to original sources
PubMed2025

Identification of two biological subgroups of complex regional pain syndrome type 1 by transcriptomic profiling of skin and blood in women.

BACKGROUND: Patients with Complex Regional Pain Syndrome (CRPS) present prolonged, debilitating pain and functional impairment. Treatments are not disease-modifying due to the poorly understood underlying pathomechanisms. This study aimed to identify the molecular signatures of potential CRPS type 1 subgroups. METHODS: Twelve women with CRPS type 1 were included. Demographics and pain questionnaires were recorded. Skin biopsies of the affected and non-affected limbs (n&#x2009;=&#x2009;6&#x2009;+&#x2009;6) and peripheral blood (n&#x2009;=&#x2009;11) were collected. RNA sequencing was performed on skin and peripheral blood mononuclear cells (PBMCs). Twenty cytokines were quantified in blood plasma (n&#x2009;=&#x2009;12). RESULTS: Cluster analysis of the affected skin identified two CRPS subgroups (SG). SG1 exhibited increased gene expression related to epidermal development, metabolic processes, and a greater abundance of keratinocytes. SG2 showed enhanced transcriptomic changes in inflammatory, immune, and fibrotic processes, along with higher abundance of fibroblasts, macrophages, and endothelial cells. PBMCs transcriptomics revealed the same SG1/SG2 clusters and highlighted a stronger inflammatory response in the blood of SG1, suggesting distinct tissue-specific immune responses for the subgroups. Interleukin-1 receptor antagonist (IL-1RA) levels were higher in the blood plasma of SG1 (FDR&#x2009;=&#x2009;0.01), consistent with its encoding gene IL1RN expression in PBMCs (log2 FC&#x2009;=&#x2009;1.10, P&#x2009;<&#x2009;0.001) and affected skin (log2 FC&#x2009;=&#x2009;0.88, P&#x2009;=&#x2009;0.006). Subgroups did not differ in demographic or clinical parameters but correlations among clinical factors varied between them. CONCLUSIONS: This study identified two potential biological subgroups of CRPS type 1 in women through skin and blood transcriptomic profiling, advancing the understanding of this condition. This could facilitate the development of targeted treatments for CRPS type 1.

Humans↗
PubMed2026

Dietary Titanium Dioxide (E171) Alters the Colon Transcriptome-Evidence From a Human Dietary Intervention Study.

Food-grade titanium dioxide (E171) has been removed from the European food market due to concerns about its potential genotoxicity, yet human evidence remains limited. Given its continuous use outside the European Union, clarifying its potential hazard is essential for human risk assessment. In a randomized crossover study, 31 adults consumed 2&#xa0;mg/kg body weight/day of E171 for 2 weeks, with biospecimens collected after both control and exposure periods. Oral intake resulted in elevated titanium concentrations in feces (CTRL: 6.3&#xa0;mg/kg; E171: 377&#xa0;mg/kg). Particle characterization showed a median size of 228&#xa0;nm, with 9% <100&#xa0;nm. Systemic oxidative stress increased, evidenced by increased superoxide radical levels in whole blood (p = 0.057). Transcriptomic profiling of colon biopsies revealed activation of 73 pathways linked to oxidative stress, cellular metabolism, and colorectal cancer-associated processes. The observed transcriptional changes were consistent with findings of its proposed mode of action. This human intervention study characterizes the effects of dietary E171 following human gastrointestinal digestion and provides evidence that exposure to dietary E171 induces systemic oxidative stress and transcriptional changes in the colon. These findings strengthen regulatory concerns about E171 as a food additive and its potentially harmful effects on human gastrointestinal health.

Humans↗
PubMed2026

An immune exhaustion signature predicts prognosis and identifies patients with diffuse large B-cell lymphoma (DLBCL) who derive preferential benefit from chimeric antigen receptor (CAR)-T cell therapy.

BACKGROUND: The tumor microenvironment (TME) is a key determinant of prognosis in diffuse large B-cell lymphoma (DLBCL). While T-cell exhaustion is implicated in therapeutic failure, its precise molecular hallmarks and utility for predicting response to modern immunotherapies, such as chimeric antigen receptor (CAR)-T cell therapy, remain unclear. METHODS: We performed an integrative analysis of transcriptomic and clinical data from multiple DLBCL cohorts (The Cancer Genome Atlas [TCGA], GSE181063, GSE10846, GSE248835, GSE182434). We used unsupervised clustering, exploratory analysis of single-cell RNA sequencing data, and the least absolute shrinkage and selection operator for variable selection (LASSO-Cox) regression to characterize the exhausted TME, construct a prognostic model, and evaluate its predictive value for CAR-T cell therapy. The model's dynamic behavior was assessed in a proof-of-concept longitudinal cohort of patients treated with the T-cell-engaging bispecific antibody glofitamab. RESULTS: We identified a "high-exhaustion" subtype associated with significantly poorer overall survival (OS; log-rank P = 0.016). Based on this, we developed a five-gene immune exhaustion-Related Prognostic Score (IERPS) that served as a robust independent predictor of poor OS across multiple cohorts. Critically, in a cohort of 256 relapsed/refractory patients, the IERPS was strongly prognostic for event-free survival (EFS) in the standard-of-care (SOC) arm (HR = 2.02, 95% confidence interval [95% CI]: 1.07-3.81, P = 0.029) but lost prognostic significance in the CAR-T arm (HR = 0.70, 95 % CI: 0.35-1.40, P = 0.314). This significant interaction suggests that CAR-T cell therapy may abrogate the poor prognosis associated with a high IERPS. Biologically, exploratory single-cell analysis (n = 4 samples) defined the high-IERPS state by hallmarks of classical T-cell exhaustion, and a descriptive case study showed the score dynamically tracked clinical response to glofitamab. CONCLUSIONS: A state of active T-cell exhaustion and a suppressive TME drive the adverse immune phenotype in DLBCL. Our IERPS model captures this dysfunctional state, acting as a powerful prognostic tool and, more importantly, as a potential predictive biomarker to identify high-risk patients who appear to overcome their inherently poor prognosis through CAR-T cell therapy.

Biomarkers↗
PubMed2026

Paired analysis of primary adenoid cystic carcinoma and derived cell lines reveals a mesenchymal and stem-like shift associated with therapy resistance.

Adenoid cystic carcinoma (ACC) is a salivary gland malignancy characterized by slow but persistent growth, frequent local recurrence, and late metastatic progression. Patients with unresectable, recurrent, or metastatic disease have limited therapeutic options. Efforts to identify effective therapeutic targets have been hindered by the limited availability of well-characterized ACC models. In this study, we established 11 ACC cell lines and performed RNA sequencing of nine cell lines and their matched primary tumors to evaluate the preservation and evolution of molecular and lineage-associated characteristics during cell line establishment. Comparative transcriptomic analysis revealed reduced epithelial and luminal differentiation programs in the cell lines, accompanied by enrichment of myoepithelial, EMT-, and cancer stem cell-associated transcriptional programs. Digital deconvolution and single-sample gene set enrichment analysis supported enrichment of hybrid EMT/stem-like states during in vitro propagation, while comparison with publicly available primary-recurrent ACC data demonstrated partial preservation of recurrence-associated plasticity and invasion programs. Protein-level validation of representative epithelial, myoepithelial, EMT, and stemness markers supported the major transcriptomic changes. In addition, a cell line with a higher stemness signature showed reduced sensitivity to cisplatin. Together, these findings indicate that ACC cell line establishment is associated with transcriptional reprogramming and enrichment of plastic, EMT/stem-like states while retaining selected ACC lineage characteristics. These models provide experimentally tractable platforms for investigating ACC progression, therapeutic response, and mechanisms of treatment resistance.

Adenoid cystic carcinoma↗
PubMed2026

A ribozyme ligase that requires a 3' terminal phosphate on its RNA substrate.

Ribozymes likely played essential roles in catalyzing metabolic processes and facilitating genome replication in primordial RNA-based life. In vitro evolution has allowed us to expand the biochemical capabilities of RNA, especially new ribozyme chemistries. Here, we report the serendipitous discovery of ribozyme ligases that catalyze the attack of the 2'-hydroxyl group of an RNA substrate on its own 5'-triphosphate group, but only when the substrate possesses a 3'-phosphate vicinal to its nucleophilic 2'-hydroxyl group. The ligases' requirement for a 3'-phosphate group on its substrate resembles enzymatic mechanisms found in protein-based RNA repair pathways. We propose that ribozyme-catalyzed ligation of 3'-phosphorylated RNA could have provided pathways for RNA repair in primordial cells. We demonstrate that these ribozymes ligate specifically to 3'-phosphorylated RNA present in a heterogeneous mixture of cellular RNAs. We further show that these ribozymes can capture cleaved RNAs with 3'-phosphate and 2'-3'-cyclic phosphate termini, enabling us to selectively amplify the captured RNAs. These results demonstrate their potential utility as enrichment reagents for profiling RNA cleavage products in transcriptomics studies. Our findings not only report a new catalytic reactivity in RNA but also provide insights into ribozyme evolution, primordial RNA repair, and potential applications in RNA sequencing.

RNA, Catalytic↗
PubMed2026

A-to-I RNA editing remodels 5'-UTR initiation codons to tune translational output.

A-to-I RNA editing is a prevalent post-transcriptional modification in higher eukaryotes that converts adenosine to inosine within RNA molecules. Because inosine is interpreted as guanosine during translation, editing can alter codon identity and potentially influence translation initiation signals. Here, we examined whether A-to-I editing within the 5' untranslated region (5'-UTR) can remodel upstream initiation codons and thereby tune downstream translation. Using luciferase-based reporter systems, we show that AUA-to-AUI editing generates an initiation-competent inosine-containing codon, whereas AUG-to-IUG editing markedly attenuates initiation and can relieve uORF-mediated repression. Quantitative in vitro and cellular assays establish the initiation hierarchy AUA&#x2009;<&#x2009;AUI&#x2009;<&#x2009;AUG, with IUG exhibiting strongly reduced initiation efficiency. Importantly, AUI-mediated upstream initiation did not behave like a canonical AUG-initiated uORF in the tested contexts; its effect on downstream ORF translation was modest and context-dependent. Transcriptome-wide bioinformatic analysis identified endogenous human transcripts whose 5'-UTRs harbor editing sites compatible with initiation-codon gain or attenuation. Reporter validation using native 5'-UTR sequences supports the possibility that editing-dependent initiation-codon remodeling can tune translational output in living cells, particularly through AUG-to-IUG-mediated derepression. Together, these findings establish a reporter-based framework in which A-to-I editing can remodel 5'-UTR initiation codons, while highlighting the need for endogenous protein-level and native-locus validation to determine physiological relevance.

RNA Editing↗
PubMed2005

Transcriptome analysis of ozone-responsive genes in leaves of European beech (Fagus sylvatica L.).

Suppression subtractive hybridization (SSH) was performed to isolate cDNAs representing genes that are differentially expressed in leaves of Fagus sylvatica upon ozone exposure. 1248 expressed sequence tags (ESTs) were obtained from 2 subtractive libraries containing early and late ozone-responsive genes. Sequences of 1139 clones (91 %) matched the EBI/NCBI database entries. For 578 clones, no putative function could be assigned. Most abundant transcripts were O-methyltransferases, representing 7 % of all sequenced clones. ESTs were organized into 12 functional categories according to the MIPS database. Among them, 12 % (early)/15 % (late) were associated with disease and defence, 19/11 % with cell structure, 4/10 % with signal transduction, and 9/6 % with transcription. The expression pattern of selected ESTs (ribulose-1,5-bisphosphate carboxylase/oxygenase small subunit [rbcS], WRKY-type transcription factor, ultraviolet-B-repressible protein, aquaporine, glutathione S-transferase, catalase, caffeic acid O-methyltransferase, and pathogenesis-related protein 1 [PR1]) was analysed by quantitative real-time RT-PCR (qRT-PCR) which confirmed changed transcript levels upon ozone treatment of European beech saplings. The ESTs characterized will contribute to a better understanding of forest tree genomics and also to a comparison of ozone-responsive genes in woody and herbaceous plants.

Europe↗
PubMed2026

Plasticity in a bacterial global regulatory switch that drives a shift in antibiotic resistance and virulence.

Antibiotic resistance and expression of virulence factors impact the outcome of infection by Pseudomonas aeruginosa. Pathogenesis is often modelled using the PAO1 reference strain but laboratory lineages vary in the sequence and activity of MexT, a global regulator impacting virulence, biofilm formation, and ciprofloxacin resistance. We defined the impact of active versus inactive MexT in PAO1 and observed transcriptomic changes affecting the expression of ~900 genes. Phenotyping revealed altered metabolism, antibiotic resistance, and virulence, resulting in striking variation across a 'single' model organism. We propose that antibiotic resistance promotes plasticity in mexT accounting for variation across lineages. We introduced antibiotic resistance into clinical P. aeruginosa isolates and observed mutations in mexT when selective pressure was removed, supporting the proposed evolutionary pathway. Overall, we have demonstrated the transcriptomic basis of MexT as a phenotypic switch in PAO1 and implicated antibiotic resistance as a cause of changes in mexT. Furthermore, MexS/MexT-regulated efflux is implicated in the antibiotic stress response and virulence, helping identify the mechanisms for rapid phenotypic switching through mexT and confirming that PAO1 is unlike most isolates. Improved understanding of the regulatory changes linked to antibiotic resistance is particularly relevant to P. aeruginosa where cycles of antibiotic treatment are common.

antibiotic resistance↗
PubMed2005

Identification of the testis-specific protein 10 (TSGA10) as serologically defined tumour-associated antigen in primary cutaneous T-cell lymphoma.

BACKGROUND: The number of identified tumour-associated antigens for cutaneous lymphoma is still very restricted, which limits the elucidation of the tumour immunology of these malignancies and the development of specific immunotherapies and immunodiagnostics. OBJECTIVES: To identify new serologically defined antigens associated with cutaneous lymphoma. METHODS: A phage expression library of the human testis transcriptom was established and immunoscreened with sera from 100 patients with cutaneous lymphoma and nine with parapsoriasis, and 81 age-matched control donors. Positive expression clones were sequenced to identify the respective antigen. RESULTS: The testis-specific protein 10 (TSGA10) was identified as an antigen recognized by sera of two patients with Mycosis fungoides but not by sera from healthy donors. By reverse transcription-polymerase chain reaction analysis, TSGA10 was found expressed in all cutaneous lymphoma samples tested, various tumour cell lines, testis, peripheral blood mononuclear cells, skin, isolated lymphocytes, keratinocytes and fibroblasts. TSGA10 overexpression had previously been reported for other cancers. CONCLUSIONS: TSGA10 is a new tumour-associated antigen of cutaneous lymphoma.

5' Untranslated Regions↗
PubMed2026

An integrated single-cell and spatial transcriptomic atlas of thyroid cancer progression identifies prognostic fibroblast subpopulations.

Although well-differentiated thyroid carcinoma (WDTC) is characterized by a robust treatment response, aggressive subtypes, such as anaplastic thyroid carcinoma (ATC), remain highly lethal. To understand thyroid cancer evolution in both children and adults, we analyzed single-cell transcriptomes of 423,733 cells from 81 samples and spatially resolved key tumor and microenvironment populations across 28 tumors with spatial transcriptomics, including rare and unique composite WDTC/ATC tumors and pediatric diffuse sclerosing thyroid carcinomas. Additionally, we identified gene signatures of stromal cell populations in 5 large thyroid cancer bulk RNA-sequencing cohorts. Through this multi-institutional effort, we defined a population of POSTN+ myofibroblast cancer-associated fibroblasts (myCAFs) that are intimately associated with invasive tumor cells and correlate with poor prognosis, lymph node metastasis, and disease progression in thyroid carcinoma. We also revealed a population of inflammatory CAFs that are distant to tumor cells and are found in the inflammatory stromal microenvironment of autoimmune thyroiditis. Together, our study provides spatial profiling of thyroid cancer evolution in samples with mixed WDTC/ATC histopathology and identifies a prognostic myCAF subtype with potential clinical utility in predicting aggressive disease in both children and adults.

Humans↗
PubMed2026

CoxFormer enables spatial omics inference with multimodal generative modeling.

Gene co-expression maps transcriptome-wide gene-gene relationships, yet high-quality estimates cover less than half the genome. Meanwhile, spatial omics either profiles restricted in situ panels or lacks cellular resolution. Extending co-expression transcriptome-wide could overcome these limitations by inferring unassayed gene expression at subcellular resolution. Here we show that CoxFormer integrates literature-derived gene knowledge with co-expression networks from bulk tissues and large-scale single-cell atlases to learn 512-dimensional representations for 32,016 human genes. These embeddings capture functional gene relationships and serve as a generative prior for spatial inference across platforms and modalities. Without requiring a matched single-cell RNA-sequencing reference, CoxFormer supports four applications beyond measured genes: histology-based expression imputation, gene activity prediction from chromatin accessibility, subcellular super-resolution inference, and pathological region detection. Together, CoxFormer extends gene embedding from gene- and cell-level tasks to whole-transcriptome spatial inference, providing a unified framework for biological analysis beyond the limited gene coverage of current spatial omics technologies.

Humans↗
PubMed2004

A method for gene expression analysis by oligonucleotide arrays from minute biological materials.

Gene expression profiling has been widely used in identifying differentially expressed genes. One of the most popular formats is oligonucleotide array. A limitation of oligonucleotide arrays is the requirement of relatively large amounts of biological starting materials for gene expression analysis. We have developed a simple method for gene expression profiling from very small amounts of biological material by combining exponential (PCR) and linear (T7 RNA polymerase) amplification. By modifying the widely used SMART protocol, we combined T7 promoter ligation and PCR amplification in one step and generated around 0.5 microg of PCRcDNA from 30 ng of total RNA in a single PCR. The PCRcDNA was in vitro transcribed by T7 RNA polymerase to generate complementary RNA (cRNA), which then was used to hybridize Affymetrix GeneChips. Our results demonstrated a linear correlation between the PCR amplification and the conventional linear amplification in gene expression ratios of individual transcript species between two different RNA preparations. The method was further validated by TaqMan reactions. Moreover, both linear and PCR methods showed some inherent bias as to which transcripts were detected, suggesting that using both in parallel may provide a more comprehensive coverage of the transcriptome present in a given sample.

Base Sequence↗
PubMed2002

Control analysis of DNA microarray expression data.

DNA microarrays produce large amounts of data. Complex changes in gene expression are revealed; sometimes thousands of mRNAs change between experiments. Here we apply modular regulation analysis to microarray data to reveal and quantify the mRNA changes that are important for cellular responses. The mRNAs are sorted into clusters. How strongly a perturbation alters each cluster is multiplied by how strongly each cluster affects an output, to obtain coefficients that describe how much of the change in the output is transmitted through each mRNA cluster. An example published dataset is analysed to reveal that the response ('relative fitness') of yeast to 2-deoxy-D-glucose is not transmitted by a single mRNA cluster, but instead many clusters contribute to the overall response. The method is applicable to microarray, transcriptome, proteome and metabolome data.

Deoxyglucose↗
PubMed2026

Metabolome-driven rhizosphere microbiome assembly determining the health of medicinal herb (Angelica sinensis) against root rot.

BACKGROUND: The rhizosphere-associated microbiota plays a crucial role in plant responses to disease stress. Plant secondary metabolites are recognized as crucial mediators in the assembly of rhizosphere microbial communities, particularly by enhancing the colonization of beneficial microorganisms. Despite this recognized importance, a deeper understanding of how such metabolome-driven microbiome assembly specifically determines plant resistance against soil-borne diseases is still lacking. RESULTS: Here, we focused on the widely planted medicinal plant Angelica sinensis and demonstrated that root rot-diseased rhizosphere soils (DRS) exhibited a higher relative abundance of Fusarium and a lower relative abundance of Streptomyces compared to healthy rhizosphere soils (HRS). Shotgun metagenomic sequencing revealed that metabolism-associated genes, particularly those related to steroid degradation, are significantly enriched in HRS samples. Subsequent genome and functional gene analysis of Streptomyces revealed that the steroid degradation-related genes are associated with rhizosphere colonization in hosts. Rhizosphere Streptomyces S15 directly antagonized Fusarium and enhanced the root resistance of A. sinensis. Comparative metabolomics showed that A. sinensis plants from HRS secreted more lipid and lipid-like molecules than those from DRS, especially sterol lipids and long-chain fatty acids, which promoted the growth of Streptomyces S15 isolates. Transcriptome analysis validated that the lipid hormones are essential for sporulation, biofilm formation, and streptomycin biosynthesis of S15 strain. Finally, exogenous application of synbiotics (lipid prebiotics and S15) to A. sinensis resulted in the enrichment of S15-homologous Streptomyces amplicon sequence variant (ASV), further establishing beneficial bacterial communities in Fusarium-stressed rhizospheres. CONCLUSIONS: Our study proposes that A. sinensis recruits steroid-metabolizing Streptomyces species by exuding key lipid compounds (i.e., methyl jasmonate and brassinolide) to combat Fusarium root rot. This study provides novel insights into using functional synbiotics as a promising strategy for manipulating plant-microbiome interactions to promote sustainable agriculture. Video Abstract.

Rhizosphere↗
PubMed2006

Hypoxia modifies the transcriptome of primary human monocytes: modulation of novel immune-related genes and identification of CC-chemokine ligand 20 as a new hypoxia-inducible gene.

Peripheral blood monocytes migrate to and accumulate in hypoxic areas of inflammatory and tumor lesions. To characterize the molecular bases underlying monocyte functions within a hypoxic microenvironment, we investigated the transcriptional profile induced by hypoxia in primary human monocytes using high-density oligonucleotide microarrays. Profound changes in the gene expression pattern were detected following 16 h exposure to 1% O(2), with 536 and 677 sequences showing at least a 1.5-fold increase and decrease, respectively. Validation of this analysis was provided by quantitative RT-PCR confirmation of expression differences of selected genes. Among modulated genes, 74 were known hypoxia-responsive genes, whereas the majority were new genes whose responsiveness to hypoxia had not been previously described. The hypoxic transcriptome was characterized by the modulation of a significant cluster of genes with immunological relevance. These included scavenger receptors (CD163, STAB1, C1qR1, MSR1, MARCO, TLR7), immunoregulatory, costimulatory, and adhesion molecules (CD32, CD64, CD69, CD89, CMRF-35H, ITGB5, LAIR1, LIR9), chemokines/cytokines and receptors (CCL23, CCL15, CCL8, CCR1, CCR2, RDC1, IL-23A, IL-6ST). Furthermore, we provided conclusive evidence of hypoxic induction of CCL20, a chemoattractant for immature dendritic cells, activated/memory T lymphocytes, and naive B cells. CCL20 mRNA up-regulation was paralleled by increased protein expression and secretion. This study represents the first transcriptome analysis of hypoxic primary human monocytes, which provides novel insights into monocyte functional behavior within ischemic/hypoxic tissues. CCL20 up-regulation by hypoxia may constitute an important mechanism to promote recruitment of specific leukocyte subsets at pathological sites and may have implications for the pathogenesis of chronic inflammatory diseases.

Cell Hypoxia↗
PubMedDate not supplied

Genomic annotation and transcriptome analysis of the zebrafish (Danio rerio) hox complex with description of a novel member, hox b 13a.

The zebrafish (Danio rerio) is an important model in evolutionary developmental biology, and its study is being revolutionized by the zebrafish genome project. Sequencing is at an advanced stage, but annotation is largely the result of in silico analyses. We have performed genomic annotation, comparative genomics, and transcriptional analysis using microarrays of the hox homeobox-containing transcription factors. These genes have important roles in specifying the body plan. Candidate sequences were located in version Z v 4 of the Ensembl genome database by TBLASTN searching with Danio and other vertebrate published Hox protein sequences. Homologies were confirmed by alignment with reference sequences, and by the relative position of genes along each cluster. RT-PCR using adult Tübingen cDNA was used to confirm annotations, to check the genomic sequence and to confirm expression in vivo. Our RT-PCR and microarray data show that all 49 hox genes are expressed in adult zebrafish. Significant expression for all known hox genes could be detected in our microarray analysis. We also find significant expression of hox 8 paralogs and hox b 7 a in the anti-sense direction. A novel gene, D. rerio hox b 13 a, was identified, and a preliminary characterization by in situ hybridization showed expression at 24 hpf at the tip of the developing tail. We are currently characterizing this gene at the functional level. We argue that the oligo design for microarrays can be greatly enhanced by the availability of genomic sequences.

Animals↗
PubMed2026

Benchmarking urinary cell transcriptomes for noninvasive differentiation of BK polyomavirus-associated nephropathy from T cell-mediated rejection.

BK polyomavirus-associated nephropathy (BKVN) adversely impacts kidney allograft survival and often mimics acute T cell-mediated rejection (TCMR), confounding diagnosis and management. To address this conundrum, we performed unbiased RNA sequencing of urinary cells matched to biopsies classified as BKVN with intragraft inflammation (BKVN-P), BKVN without inflammation (BKVN-N), TCMR, or no rejection (NR). BKVN-N displayed dominant host DNA replication, cell cycle, and repair programs, while BKVN-P samples exhibited expansive innate immune activation, antigen presentation, chemokine upregulation, and epithelial injury. Both BKVN subtypes shared signatures of T cell exhaustion and mature and tolerogenic dendritic cell activation but differed in immune orientation - Th1 predominance in BKVN-N versus Treg and CD8 enrichment in BKVN-P. Compared with TCMR samples, BKVN-P lacked robust TCR/CD28 signaling and was enriched for viral and innate modules; BKVN-N lacked alloimmune activation. B cell exhaustion characterized BKVN-N, while BKVN-P displayed robust B cell activation with metabolic downregulation. A ratiometric urinary cell biomarker, CXCL10 mRNA/CD3E mRNA, distinguished both BKVN subtypes from TCMR with diagnostic accuracy, replicated by quantitative reverse transcription PCR for clinical translation, and confirmed in an independent cohort. These findings demonstrate the utility of urinary cell transcriptomics for resolving viral injury from alloimmunity, enabling precision diagnostics and targeted immunomodulation in kidney transplantation.

Humans↗