PubMed HealthSearch

SEARCH · PubMed Health

Results for “transcriptomics”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

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

At least 19 recordsLinked to original sources

Dual-transcriptomic analysis of human nasal transcriptome and microbiome reveals host-bacteria associations in symptomatic respiratory infection.

BACKGROUND: The human nasopharynx is colonized by a diverse community of commensal microbiota linked to many respiratory diseases, yet their associations with the host remain unclear. RESULTS: In this study, we introduced a dual-transcriptomics analysis strategy, which can characterize the host transcriptome and microbiome from nasal samples simultaneously. We applied this workflow to a local SARS-CoV-2 cohort with 76 asymptomatic infected patients, among whom 52 (68.42%) developed symptomatic infection during a 1-week follow-up period. Nasal swabs were collected from all 76 patients at enrollment and from 73 patients at one-week later follow-up. We detected a median of 8.94% reads that did not map to the human genome across all 149 samples, among which around half (median 49.68%) were successfully mapped to microbiome genome. Meta-transcriptomic analysis detected significantly higher SARS-related coronavirus loads in samples from the symptomatic group at enrollment (P&#x2009;=&#x2009;0.004), and both groups showed decreased loads one week later (symptomatic, P&#x2009;=&#x2009;0.001; asymptomatic, P&#x2009;=&#x2009;0.035). Compared with benchmarking 16&#xa0;S rRNA sequencing on 53 samples, our computational strategy showed high correlation of relative abundance in all top 20 genera (median Rho&#x2009;=&#x2009;0.90, Pmax < 0.001). A total of 670 bacteria species were identified to show a relative abundance&#x2009;&#x2265;&#x2009;0.01% in at least 10% samples. Differential abundance analysis identified 76 species (DASs) from six phyla with significantly decreased abundance in samples from the symptomatic group (log2(fold change or FC) < -1 and adjusted P&#x2009;<&#x2009;0.05) compared to the asymptomatic group at enrollment. Integrating these symptom-associated DASs with host's gene expression using an expression quantitative trait bacteria (eQTB) model, we found 45 symptom-associated DASs identified at enrollment were significantly associated with one to 14 genes (adjusted P&#x2009;<&#x2009;0.05). GSEA showed a series of symptom-associated DASs were significantly correlated with pathways related to olfactory function, keratinocyte differentiation, and DNA methylation. CONCLUSIONS: In summary, our dual-transcriptomic analysis strategy effectively characterized host-microbiome associations, offering insights into microbial contributions to respiratory diseases.

Humans

Spatial transcriptomics-aided localization for single-cell transcriptomics with STALocator.

Single-cell RNA-sequencing (scRNA-seq) techniques can measure gene expression at single-cell resolution but lack spatial information. Spatial transcriptomics (ST) techniques simultaneously provide gene expression data and spatial information. However, the data quality of the spatial resolution or gene coverage is still much lower than the quality of the single-cell transcriptomics data. To this end, we develop a ST-Aided Locator for single-cell transcriptomics (STALocator) to localize single cells to corresponding ST data. Applications on simulated data showed that STALocator performed better than other localization methods. When applied to the human brain and squamous cell carcinoma data, STALocator could robustly reconstruct the relative spatial organization of critical cell populations. Moreover, STALocator could enhance gene expression patterns for Slide-seqV2 data and predict genome-wide gene expression data for fluorescence in situ hybridization (FISH) and Xenium data, leading to the identification of more spatially variable genes and more biologically relevant Gene Ontology (GO) terms compared with the raw data. A record of this paper's transparent peer review process is included in the supplemental information.

Single-Cell Analysis

Transcriptome-wide association analysis of Alzheimer's disease: construction and clinical validation of transcriptomic risk scores.

Early identification of individuals at high risk for Alzheimer's disease (AD) is crucial for disease prevention and intervention. This study aims to develop AD-specific transcriptomic risk scores (TRSs) through multi-tissue transcriptome-wide association study (TWAS) and to evaluate its clinical utility in AD diagnosis and risk prediction. Using GWAS summary statistics combined with expression quantitative trait loci (eQTL) data from 14 tissues, a multi-tissue TWAS approach was applied to identify AD-associated genes. Peripheral blood RNA expression data from the ADNI and GEO databases were used to construct the AD-specific TRSs. The associations of TRSs with AD pathological features and cognitive function were assessed in two independent cohorts. Furthermore, the diagnostic performance, differential diagnostic capability, and risk prediction efficiency of TRSs were evaluated. The TWAS identified 131 genes significantly associated with AD. The TRSs were significantly elevated in patients with AD and mild cognitive impairment (MCI) compared to cognitively normal (CN) individuals, and showed significant correlations with AD pathological markers and cognitive performance. When combined with APOE4 status, the TRSs demonstrated robust diagnostic ability for AD and MCI. When combined with age, the TRSs showed good diagnostic performance in distinguishing AD from frontotemporal dementia (FTD) (AUC&#x2009;=&#x2009;0.86). Additionally, the TRSs effectively predicted the risk of progression to AD in non-AD individuals (HR&#x2009;=&#x2009;1.74). The AD-specific TRSs developed in this study shows promising clinical utility in AD diagnosis, differential diagnosis, and risk prediction, providing valuable translational medical evidence for early screening and precision prevention of Alzheimer's disease.

Humans

Coordinated inflammatory macrophage and vascular smooth muscle cell remodeling signatures in human atherosclerosis: An integrative single-cell and bulk transcriptomic analysis.

Atherosclerotic plaque progression is shaped by coordinated inflammatory and remodeling programs involving immune cells and vascular wall cells. Inflammatory macrophage activation and vascular smooth muscle cell (VSMC) phenotypic remodeling are central features of human atherosclerosis, but their transcriptomic relationships during plaque progression remain incompletely characterized. This study integrated single-cell and bulk transcriptomic datasets to examine highly inflammatory macrophage states, VSMC remodeling-related transcriptional programs, and candidate ligand-receptor expression patterns in human atherosclerotic plaques. Human atherosclerotic plaque single-cell RNA sequencing data from GSE260657 and bulk transcriptomic data from GSE28829 were analyzed. After quality control, 7628 cells were retained for single-cell analysis. Major cell types were annotated using canonical markers, followed by reclustering of macrophages and VSMC-related cells. Functional module scoring, differential expression analysis, Gene Ontology biological process enrichment, and Kyoto Encyclopedia of Genes and Genomes pathway analyses were performed to characterize macrophage transcriptional states. Slingshot was applied to infer VSMC pseudotime ordering. CellChat and NicheNet were used to prioritize candidate ligand-receptor expression patterns and ligand-associated VSMC target gene programs. External bulk transcriptomic analysis was performed to examine whether single-cell-derived inflammatory and remodeling signatures were represented at the tissue-transcriptome level during plaque progression. Macrophage reclustering identified a highly inflammatory macrophage state characterized by prominent inflammatory activation, cytokine-response, and stress-response features. Genes upregulated in this population were enriched in pathways related to tumor necrosis factor (TNF) response, nuclear factor kappa B signaling, leukocyte activation, cytokine signaling, lipid and atherosclerosis, toll-like receptor signaling, and inflammasome-associated inflammation. VSMC reclustering revealed contractile VSMCs, PTHLH+&#xa0;synthetic VSMCs, KRT7+ VSMC-like cells, interferon-responsive VSMCs, pericyte-like mural cells, and osteogenic/modulated VSMCs. Pseudotime analysis showed a broad contractile-to-osteogenic/modulated transcriptional continuum accompanied by increased expression of remodeling-associated genes and selected inflammatory or remodeling-associated receptor genes. CellChat and NicheNet analyses prioritized candidate ligand-receptor and ligand-associated target gene expression patterns involving SPP1-CD44, TNF-TNFRSF1A, IL1B-IL1R1/IL1RAP, MIF-ACKR3, PDGFB-PDGFRB, and FN1-SDC1/ITGB1. In GSE28829, inflammatory macrophage-, osteogenic/modulated VSMC-, candidate ligand-receptor expression-, SPP1-CD44 candidate axis-, and NicheNet-prioritized target program-related signatures were more prominent in advanced plaques and were positively correlated with each other. This integrative transcriptomic analysis identified a highly inflammatory macrophage state and a VSMC remodeling continuum in human atherosclerotic plaques. Candidate ligand-receptor and ligand-associated target gene expression patterns linked inflammatory macrophage activation with osteogenic/modulated VSMC remodeling at the computational level. External bulk data further showed coordinated enrichment of inflammatory and remodeling signatures in advanced plaques. These findings provide a descriptive and hypothesis-generating transcriptomic framework for understanding inflammatory macrophage activation and VSMC remodeling in human atherosclerosis.

atherosclerosis

Evolutionary Reorganization of Transcriptomic Architecture Across a UVB Tolerance Gradient in Fish.

Environmental stressors such as ultraviolet radiation impose strong selective pressures on organisms, yet how adaptation to such stressors shapes transcriptomic responses at the network level remains poorly understood. Although stratospheric ozone is recovering globally, substantial regional variation in UV exposure persists, particularly in high-altitude environments where extreme UV levels can occur. Here, we compared three fish models representing distinct biological responses to UVB exposure: wild-type zebrafish (Danio rerio), a melanin-deficient zebrafish mutant (nacre) lacking a major protective mechanism against UVB damage, and the high-altitude Andean killifish Orestias ascotanensis, a species naturally exposed to extreme UVB radiation. Together, these models define a gradient spanning physiological protection, impaired protection, and evolutionary adaptation to UVB stress. Using RNA-seq and protein-protein interaction networks, we show that transcriptomic responses differ markedly across this gradient. Wild-type and nacre zebrafish exhibited relatively limited transcriptomic changes (&#x223c;2%-2.4% of genes changing), whereas O. ascotanensis displayed a large-scale and highly coordinated response (&#x223c;21.6% of genes changing) characterized by functionally specialized networks enriched in DNA repair pathways. These differences involved not only transcriptomic magnitude but also marked reorganization of transcriptomic architecture. Integration with positive selection analyses revealed that positively selected genes were concentrated within highly interconnected regions of transcriptomic networks, consistent with adaptation involving network reorganization. Furthermore, ortholog-based analyses suggest that adaptive responses involve differential reorganization of a conserved functional background. Together, our results support a model in which adaptation to environmental stress is associated with the reorganization of conserved transcriptomic networks across physiological and evolutionary contexts, providing a systems-level perspective on the molecular basis of adaptation.

UVB radiation

Transcriptomic responses of gill and intestinal tissues in Nile tilapia (Oreochromis niloticus) to bacterial infection following sequential nanoimmersion and hydrogel-based multivalent vaccination.

Bacterial pathogens, including Flavobacterium oreochromis, Aeromonas veronii, Streptococcus agalactiae, and Edwardsiella tarda, represent major infectious threats to Nile tilapia (Oreochromis niloticus). A multivalent vaccination strategy integrating cationic nanoemulsion immersion with oral hydrogel boosters was developed to investigate tissue-specific immune responses at the transcriptomic level. Gill tissues were collected following immersion challenge and intestinal tissues following intraperitoneal injection challenge, reflecting the physiologically relevant infection biology of each pathogen and the mechanistic rationale of each delivery platform. RNA sequencing (RNA-seq) generated high-quality datasets (mapping rate&#xa0;>&#xa0;81.64%) with strong concordance to quantitative real-time PCR (qRT-PCR) validation (r&#xa0;=&#xa0;0.83). Comparative transcriptomic analysis revealed distinct yet complementary immune signatures between tissues. Gill transcriptomes were enriched in phagosome, focal adhesion, extracellular matrix-receptor interaction (ECM-receptor interaction), and cytokine-cytokine receptor interaction pathways, accompanied by increased expression of major histocompatibility complex class I/II (MHC class I/II), mannose receptor, &#x3b1;V&#x3b2;3 integrin, and calnexin, indicating innate activation, enhanced phagocytic capacity, epithelial barrier reinforcement, and adaptive immune coordination. Intestinal transcriptomes showed predominant enrichment of adaptive immune pathways, including the intestinal immune network for immunoglobulin (Ig) production, Forkhead box O (FoxO) signaling, and mitogen-activated protein kinase (MAPK) signaling, with increased expression of T-cell receptor (TCR), inducible T-cell co-stimulator ligand (ICOS-L), C-X-C chemokine receptor type 4 (CXCR4), and polymeric immunoglobulin receptor (pIgR), reflecting T and B cell coordination, lymphocyte trafficking, and mucosal immunoglobulin transport, alongside innate engagement through phagosome pathway enrichment. Shared upregulation of MHC class II, B-cell receptor (BCR) signaling, integrin alpha M (ITGAM), and immunoglobulin-associated components across both tissues suggests coordinated mucosal immune activation through a conserved immune module, warranting direct experimental validation. Collectively, these findings provide transcriptomic evidence that this vaccination strategy elicits an integrated, tissue-specialized immune response, advancing mechanistic understanding of gill and intestinal immunity in vaccine-induced protection of teleost fish.

Animals

Dynamic effects of short-term storage temperature and duration on the transcriptome and functional pathways of umbilical cord blood-derived NK cells.

OBJECTIVE: This study aims to elucidate the dynamic impact of different storage conditions (temperature and time) on the function of natural killer (NK) cells derived from umbilical cord blood (UCB) at the transcriptome level, providing a theoretical basis for optimizing the standardized post-collection processing protocol of UCB in clinical settings. METHODS: Four healthy full-term UCB samples were collected and assigned to a fresh control group (0H) and experimental groups stored at 4&#xa0;&#xb0;C (4C) or 25&#xa0;&#xb0;C (RT) for 24&#xa0;h (24H) and 72&#xa0;h (72H). Umbilical cord blood mononuclear cells (CBMCs) were isolated and expanded in vitro to derive NK cells. Using RNA sequencing (RNA-seq) technology, combined with principal component analysis (PCA), screening of differentially expressed genes (DEGs), Gene Ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses, the effects of storage conditions on the transcriptome of NK cells were comprehensively analyzed. RESULTS: The results indicated that storage time was the primary factor causing transcriptomic differences. Compared to the 0H group, storage for 24&#xa0;h induced a limited number of Differentially Expressed Genes (DEGs) (only 11 in the 24H-RT group), indicating overall transcriptome stability. In contrast, storage for 72&#xa0;h triggered profound transcriptomic reprogramming, with DEGs significantly enriched in immune activation-related pathways. Specifically, at the 72-h time point, storage at 4&#xa0;&#xb0;C compared to 25&#xa0;&#xb0;C prevented the downregulation of key metabolic pathways such as 'response to decreased oxygen levels' and 'regulation of leukocyte cell-cell adhesion'. CONCLUSION: CBMC-derived NK cells can tolerate up to 24&#xa0;h of ex vivo storage while maintaining transcriptome and core functional stability. However, extending the ex vivo storage period to 72&#xa0;h leads to significant reprogramming in the expression of immune and metabolism-related genes. Notably, the degree of metabolic suppression under refrigerated conditions at 4&#xa0;&#xb0;C was less pronounced than that at 25&#xa0;&#xb0;C room temperature storage. This study recommends that the ex vivo storage time of UCB should ideally be limited to within 24&#xa0;h, with room temperature storage being a viable option. If prolonged storage is necessary, refrigeration at 4&#xa0;&#xb0;C is advised.

Humans

Exogenous ABA enhances cold tolerance of Rhododendron yedoense var. poukhanense under subzero temperature: integrating physiology, transcriptome, and proteome.

Low temperature limits the growth and ornamental value of evergreen shrubs. Rhododendron yedoense var. poukhanense, an important ornamental shrub from Northeast China, frequently suffers freezing damage during winter. While exogenous abscisic acid (ABA) enhances cold tolerance in many plants, its molecular mechanisms at subzero temperatures remain poorly understood in non-model species lacking chromosome-level reference genomes. This study investigated the effects of exogenous ABA on freezing tolerance in R. yedoense var. poukhanense at -4&#xa0;&#xb0;C using an integrated physiological, transcriptomic, and proteomic approach. Cutting seedlings were subjected to four treatments: CK (22&#xb0;C control), A (22&#xb0;C + ABA), LT (-4&#xb0;C), and ALT (-4&#xb0;C + ABA). Photosynthetic pigments, osmotic regulation substances, antioxidant enzyme activities, and malondialdehyde (MDA) content were measured. Transcriptome sequencing and quantitative proteomics were performed, and transcriptome data were validated by quantitative real-time PCR (qRT-PCR) of 15 selected genes. ABA pretreatment reduced visible cold injury severity, partially preserved photosynthetic pigments, decreased MDA content by 28.7%, and promoted recovery of catalase (+43.6%), superoxide dismutase (+31.1%), and peroxidase (+20.0%) activities under freezing stress. Transcriptome analysis revealed 8, 444 differentially expressed genes (DEGs) in LT versus CK and 6, 481 DEGs in ALT versus CK, representing a 23% reduction in transcriptional reprogramming scope attributable to ABA priming. The ALT versus LT comparison identified only 1, 690 additional DEGs, indicating that most cold-responsive genes were pre-activated during the ABA priming phase. Proteome analysis identified 1, 461 differentially expressed proteins (DEPs) in ALT versus CK. Integrated analysis revealed extensive post-transcriptional regulation, with transcript-protein concordance of only 1.0-4.1%, and co-enriched Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways in both omics layers. qRT-PCR validation confirmed high reliability of the transcriptome data (R2&#xa0;=&#xa0;0.8500). These findings demonstrate that exogenous ABA enhances freezing tolerance through multi-layered molecular regulation encompassing transcriptional buffering, translational reprogramming, and functional reallocation from photosynthesis to stress protection. This study provides the first integrated physiology-transcriptome-proteome framework for ABA-mediated freezing tolerance in an evergreen ornamental shrub and offers theoretical support for ABA-based winter protection strategies.

Rhododendron yedoense var. Poukhanense

Comprehensive transcriptomics and proteomics analysis of neointima formation in human saphenous vein: implications for bypass graft disease.

Human saphenous veins (SVs) are widely used as grafts in coronary artery bypass (CABG) surgery but often fail due to neointima formation. Little is known, however, regarding the cellular, transcriptomic, and proteomic dynamics of neointima formation in human veins. Here, we performed transcriptomics and proteomics analysis in an ex vivo tissue culture model of neointima formation in human SVs procured for CABG surgery. Histological examination demonstrated significant elastin degradation and neointima formation (indicated by increased neointima area and neointima-to-media ratio) in SVs subjected to tissue culture. Analysis of data from 72 patients suggests that the progression of SV remodeling and neointima formation differs according to sex and body mass index, which is negatively associated with neointima formation in males only. RNA sequencing demonstrated upregulation of proinflammatory and proliferation-related genes during neointima formation and identified novel processes, including increased cellular stress and DNA damage responses, reflecting tissue trauma associated with vein harvesting. Proteomic analysis identified upregulated extracellular matrix-related and coagulation/thrombosis proteins and downregulated metabolic proteins. Spatial transcriptomics, used to infer regionally enriched gene expression, suggested dynamic alterations in fibroblast and vascular smooth muscle cell (VSMC) states during neointima formation. Specifically, we identified the emergence of HES1+ and matrix metalloproteinase 2- and 14-positive (MMP2+/MMP14+) expression in VSMCs and fibroblasts, respectively, during neointima formation. Furthermore, our data suggest that MIR647, identified through screening, maintains VSMC contractile gene expression. Our findings suggest dynamic transcriptomic and proteomic changes during neointima formation in human veins and provide useful mechanistic information for the pathogenesis of SV graft disease.NEW & NOTEWORTHY Using multiomics and spatial transcriptomics, we uncover dynamic molecular and cellular changes driving neointima proliferation in human saphenous veins, the most common conduit for bypass surgery. Our study highlights sex- and body mass index-associated differences, novel fibroblast and smooth muscle cell states, and a role for microRNA-647 in preserving vascular contractile phenotype. These findings provide new insight into the mechanisms of vein graft failure and may guide future strategies to improve coronary bypass outcomes.

Humans

A contextual activity score (CAS) for inferring ADAR-associated transcriptional activity across RNA-seq, single-cell, and spatial transcriptomics.

BACKGROUND AND OBJECTIVE: Adenosine-to-inosine RNA editing, catalyzed by Adenosine Deaminases Acting on RNA (ADARs), is a widespread modification involved in neural function, immune regulation, and cancer. The Alu Editing Index (AEI) is the standard metric to estimate ADAR activity but requires raw sequencing reads and is poorly suited for single-cell and spatial transcriptomic data. This study aimed to develop an alternative framework for inferring ADAR-associated transcriptional activity from gene expression data across diverse transcriptomic technologies. METHODS: We developed the Contextual Activity Score (CAS), a framework based on transcriptional signatures from ADAR perturbation experiments. Context-specific signatures were generated for human neurons, mouse neurons, and cancer models to infer ADAR1 and ADAR2 activity. CAS was computed from normalized gene expression matrices using regulon-based enrichment analysis. Performance was evaluated by comparing with the Alu Editing Index across bulk RNA sequencing datasets, simulated sequencing depths, and library preparation protocols. RESULTS: CAS showed strong concordance with the Alu Editing Index across multiple datasets, while remaining robust to reduced sequencing depth and different library protocols. Unlike the Alu Editing Index, CAS can be applied to single-cell and spatial transcriptomic data and enables the independent assessment of ADAR2 activity. In cancer and neuronal contexts, CAS captured biologically meaningful variations in ADAR-associated transcriptional activity at sample, cell-type, and spatial levels. CONCLUSION: CAS provides a scalable approach applicable across multiple RNA-seq protocols for estimating ADAR-associated transcriptional activity using gene expression data. This method, implemented in an open-source R package for broad adoption, expands the ability to study ADAR-associated transcriptional activity across transcriptomic modalities where direct editing quantification is challenging, such as single-cell and spatial transcriptomics.

Adenosine Deaminase

Preliminary Exploration on Melatonin-Mediated Protective Effects in Intracranial Aneurysms: Transcriptomic, Proteomic, and Metabolomic Profiling of Cerebral Vascular Tissues Combined with in vivo Animal Experiments.

BACKGROUND: Intracranial aneurysm (IA) is a life-threatening cerebrovascular disease with unclear molecular mechanisms and limited drug treatment. Our previous research has shown that melatonin (MLT) has potential protective effects in IA, but its mechanism remains unclear. The purpose of this study is to explore the pathological mechanism of IA and the therapeutic mechanism of MLT by integrating transcriptomic, proteomic and metabolomic analyses. METHODS: In this study, mouse models of IA were successfully established by combining elastase injection with angiotensin II infusion. C57BL/6 mice were divided into control, IA model, IA model+MLT, and IA model+nimodipine groups. The pathological conditions were evaluated by hematoxylin-eosin (HE) staining, TUNEL staining, and scanning electron microscopy. Transcriptomic (n=3 for each group), proteomic (n=3 for each group), and metabolomic (n=6 for each group) analyses were performed based on cerebral vascular tissue samples. The screening thresholds for differentially expressed genes and differentially expressed proteins were P <0.05 and fold change >1.5 and fold change <0.667. The screening criteria for differential metabolites were variable importance for the projection (VIP)> 1.0, fold change >1.2 and fold change <0.833, and P <0.05. RESULTS: MLT alleviated brain tissue damage, vascular endothelial damage, structural disruption, and apoptosis in IA mice. Transcriptomic, proteomic and metabolomic analyses identified numerous differential molecules. Functional annotation revealed that these molecules may be involved in biological pathways and processes such as immune inflammation, vascular remodeling, extracellular matrix remodeling, neuropeptide activity, oxidative stress and metabolic pathways, thereby regulating the occurrence and development of IA or mediating the therapeutic effects of MLT. Furthermore, transcriptomic and proteomic analyses also suggest that there may be extensive post-transcriptional, translational and post-translational regulatory events in the progression of IA and the therapeutic effects of MLT. Integrated transcriptomic and proteomic analyses suggest that Npy may be a key molecule in regulating IA progression and mediating MLT therapeutic effects, and its potential value is further supported by our immunohistochemical validation results. CONCLUSION: Multi-omics integrative analysis preliminarily revealed that the potential mechanisms of MLT may involve the regulation of inflammatory response, vascular remodeling, extracellular matrix remodeling, neuropeptide activity, oxidative stress, metabolic pathways, and post-transcriptional/translational regulation.

Animals

Identification of transcriptome SNPs between Xiphophorus lines and species for assessing allele specific gene expression within F&#x2081; interspecies hybrids.

Variations in gene expression are essential for the evolution of novel phenotypes and for speciation. Studying allelic specific gene expression (ASGE) within interspecies hybrids provides a unique opportunity to reveal underlying mechanisms of genetic variation. Using Xiphophorus interspecies hybrid fishes and high-throughput next generation sequencing technology, we were able to assess variations between two closely related vertebrate species, Xiphophorus maculatus and Xiphophorus couchianus, and their F(1) interspecies hybrids. We constructed transcriptome-wide SNP polymorphism sets between two highly inbred X. maculatus lines (JP 163 A and B), and between X. maculatus and a second species, X. couchianus. The X. maculatus JP 163 A and B parental lines have been separated in the laboratory for &#x2248;70 years and we were able to identify SNPs at a resolution of 1 SNP per 49 kb of transcriptome. In contrast, SNP polymorphisms between X. couchianus and X. maculatus species, which diverged &#x2248;5-10 million years ago, were identified about every 700 bp. Using 6524 transcripts with identified SNPs between the two parental species (X. maculatus and X. couchianus), we mapped RNA-seq reads to determine ASGE within F(1) interspecies hybrids. We developed an in silico X. couchianus transcriptome by replacing 90,788 SNP bases for X. maculatus transcriptome with the consensus X. couchianus SNP bases and provide evidence that this procedure overcomes read mapping biases. Employment of the in silico reference transcriptome and tolerating 5 mismatches during read mapping allow direct assessment of ASGE in the F(1) interspecies hybrids. Overall, these results show that Xiphophorus is a tractable vertebrate experimental model to investigate how genetic variations that occur during speciation may affect gene interactions and the regulation of gene expression.

Alleles

Genomic and Transcriptomic Profiling of Radiation-Resistant, Locally Recurrent Prostate Cancer.

PURPOSE: The biology of locally radiorecurrent prostate cancer (LRR-PCa) is poorly understood. METHODS AND MATERIALS: We sought to explore the genomic and transcriptomic landscape of LRR-PCa with targeted DNA sequencing and RNA expression analysis from 41 biopsy-proven LRR-PCa tumors from 36 unique patients who had a recurrence at a median interval of 84 months (IQR, 70-124 months). Genomic alteration frequencies and transcriptomic data were compared between the LRR-PCa cohort and treatment-na&#xef;ve patients from the Cancer Genome Atlas (genomic; n = 496) and Gleason grade-at-recurrence-matched patients from the Decipher Genomics Resource for Intelligent Discovery (transcriptomic; n = 22,320). RESULTS: Twenty-five patients (69%) had pathologic upgrading at recurrence (17% vs 64% with Gleason grade 4-5 disease; P < .001). The LRR-PCa cohort demonstrated significantly greater single-nucleotide variations in 29 genes known to be associated with prostate cancer, including several associated with increased aggressiveness and DNA repair: FAT1 (58.5% vs 1.0%), RAD51B (36.6% vs 0.4%), POLQ (34.1% vs 1.4%), KMT2C (34.1% vs 4.9%), BRCA2 (29.3% vs 1.8%), ATRX (26.8% vs 0.8%), and BRCA1 (24.4% vs 0.4%) (Pvalues < .001 for all). The LRR-PCa cohort had a significantly higher Decipher score (median, 0.80 vs 0.66; P = .05) and demonstrated significantly greater basal subtype based on PAM50 (56% vs 20%; P < .001) and lower androgen receptor activity (61% for LRR vs 9%; P < .001). CONCLUSIONS: Overall, these results suggest that LRR-PCa has a distinct genomic and transcriptomic landscape from de novo prostate cancer. Specifically, LRR-PCa has an enrichment in SNVs in genes associated with tumor aggressiveness and/or DNA repair, has higher Decipher scores, a more basal subtype, and has transcriptomic evidence of lower androgen receptor activity and loss of tumor suppressor genes.

Humans

A novel reusable transcriptome-wide association study workflow used to map key genes linked to important cattle traits.

Transcriptome-wide association studies (TWAS) are a powerful approach for studying the genes underlying complex traits by directly integrating GWAS and gene expression datasets. In cattle, they have been previously applied to identify genes driving fertility, milk production, and health. However, these studies have also highlighted several challenges, from difficulties in reproducing these complex analyses to limitations from poor genotype calls, especially when called directly from RNA sequencing data. To address these and other challenges, for the H2020 BovReg Project, we have developed a streamlined, species-agnostic, and reusable Nextflow TWAS workflow to integrate transcriptomic and GWAS summary statistic datasets. Our workflow first generates accurate genotype calls and gene expression prediction models from transcriptomic datasets and then applies these tools to impute gene expression levels into GWAS cohorts, enabling the association of genes with traits of interest. We explore optimal strategies for calling genetic variants directly from transcriptomic data and illustrate that using imputation approaches specifically designed for low-pass sequencing data can improve variant calling over previously adopted methods. We demonstrate the utility of our TWAS workflow by applying it to both novel and publicly available GWAS cohorts for cattle, detecting novel gene-trait associations for complex traits. Using a new transcriptome annotation of the cattle genome generated for the BovReg project we also illustrate how previously un-assayable associations can be detected. The results and the workflow we present, provide a new resource for the community and contribute to a better understanding of the molecular drivers of complex traits in cattle with the goal of eventually leveraging this information in future breeding decisions.

Animals

Species-specific transcriptomic changes upon respiratory syncytial virus infection in cotton rats.

The cotton rat (Sigmodon) is the gold standard pre-clinical small animal model for respiratory viral pathogens, especially for respiratory syncytial virus (RSV). However, without a reference genome or a published transcriptome, studies requiring gene expression analysis in cotton rats are severely limited. The aims of this study were to generate a comprehensive transcriptome from multiple tissues of two species of cotton rats that are commonly used as animal models (Sigmodon fulviventer and Sigmodon hispidus), and to compare and contrast gene expression changes and immune responses to RSV infection between the two species. Transcriptomes were assembled from lung, spleen, kidney, heart, and intestines for each species with a contig N50&#x2009;>&#x2009;1600. Annotation of contigs generated nearly 120,000 gene annotations for each species. The transcriptomes of S. fulviventer and S. hispidus were then used to assess immune response to RSV infection. We identified 238 unique genes that are significantly differentially expressed, including several genes implicated in RSV infection (e.g., Mx2, I27L2, LY6E, Viperin, Keratin 6A, ISG15, CXCL10, CXCL11, IRF9) as well as novel genes that have not previously described in RSV research (LG3BP, SYWC, ABEC1, IIGP1, CREB1). This study presents two comprehensive transcriptome references as resources for future gene expression analysis studies in the cotton rat model, as well as provides gene sequences for mechanistic characterization of molecular pathways. Overall, our results provide generalizable insights into the effect of host genetics on host-virus interactions, as well as identify new host therapeutic targets for RSV treatment and prevention.

Animals

Enhancing and accelerating cell type deconvolution of large-scale spatial transcriptomics slices with dual network model.

MOTIVATION: Cell type deconvolution deciphers spatial distribution of mRNA transcripts at single cell level by integrating single-cell RNA sequencing (scRNA-seq) and spatial transcriptomics data to infer mixture of cell types of spots in slices. Current algorithms are criticized for neglecting connection between scRNA-seq and spatial transcriptomics data, as well as time-consuming, hampering their application to large-scale datasets. RESULTS: In this study, we propose a joint learning nonnegative matrix factorization algorithm for fast cell type deconvolution (aka jMF2D), which integrates scRNA-seq and spatial transcriptomics data with network models. To bridge scRNA-seq and spatial transcriptomics data, jMF2D jointly learns cell type similarity network to enhance quality of signatures of cell types, thereby promoting accuracy and efficiency of deconvolution. Experiments demonstrate that jMF2D outperforms state-of-the-art baselines in terms of accuracy by saving about 90% running time on various datasets generated by different platforms. Furthermore, it can also facilitates the identification of spatial domains and bio-marker genes, providing an efficient and effective model for analyzing spatial transcriptomics data. AVAILABILITY AND IMPLEMENTATION: The software is coded using python, and is free available for academic https://github.com/xkmaxidian/jMF2D.

Algorithms

Endurance training attenuates acute exercise-induced monocyte transcriptomic responses in adolescents: sex-specific molecular adaptations.

Circulating monocytes contribute to atherogenesis and vascular dysfunction. Although clinical cardiovascular disease presents in adulthood, its biological origins often begin in youth and differ substantially between females and males. Twelve males and nine females (13-17 yr) completed an acute exercise protocol consisting of 10, 2-min cycling bouts at 70% of maximal work rate interspersed with 1-min rest intervals, performed before and after an 8-wk supervised endurance training intervention with brief supplementary strength work (60 min/session, 3 sessions/wk). Blood was collected before and immediately after each exercise challenge. Peripheral blood monocytes were isolated, and whole transcriptome RNA sequencing (RNA-seq) was performed. Before training, acute exercise induced a markedly greater monocyte transcriptomic response in females compared with males [5,135 vs. 567 differentially expressed transcripts, false discovery rate (FDR) < 0.1]. Pathway analyses identified vascular function-related pathways in males, whereas females showed enrichment of pathways related to adipose tissue cross talk and oxidative metabolism. Following training, the acute transcriptomic response was markedly attenuated in both sexes (165 transcripts in males and 94 in females, FDR < 0.1), representing an &#x223c;98% reduction in females and a &#x223c;70% reduction in males relative to pre-training responses. These findings reveal sex-specific monocyte responses to acute exercise in youth and suggest that endurance exercise alters immune transcriptional responsiveness in pathways relevant to vascular and cardiovascular health.NEW & NOTEWORTHY Acute exercise induced a markedly greater monocyte transcriptomic response in female adolescents than in males. Females showed activation of pathways related to adipose tissue signaling and oxidative metabolism, whereas males exhibited vascular-function pathways. Following exercise training, the monocyte transcriptomic response to acute exercise was substantially attenuated in both sexes. These findings identify sex-specific immune transcriptional responses to exercise during adolescence with potential implications for cardiovascular health.

Humans

Spatiotemporal transcriptomic analysis during cold ischemic injury to the murine kidney reveals compartment-specific changes.

BACKGROUND: Kidney transplantation is the preferred treatment strategy for end-stage kidney disease. Deceased donor kidneys usually undergo cold storage until kidney transplantation, leading to cold ischemia injury that may contribute to poor graft outcomes. However, the molecular characterization of potential mechanisms of cold ischemia injury remains incomplete. RESULTS: To bridge this knowledge gap, we leverage 10x&#x2009;Visium spatial transcriptomic technology to perform full transcriptome profiling of murine kidneys subject to varying durations of cold ischemia typical in a deceased donor kidney transplant setting. We develop a computational workflow to identify and compare spatiotemporal transcriptomic changes that accompany the injury pathophysiology in a tissue compartment-specific manner. We identify proportional enrichment of oxidative phosphorylation (OXPHOS) genes with increasing duration of cold ischemia injury within the oxygen-lean inner medulla region, suggestive of atypical metabolic presentation. This is distinct in cold ischemia injury tissue compared to warm ischemia-reperfusion kidney injury tissue. Spatiotemporal trends are validated by qPCR and immunofluorescence in a larger cohort of mice. CONCLUSIONS: Altogether, our spatiotemporal transcriptomic analysis identifies coordinated molecular changes within metabolic pathways such as OXPHOS deep within the cold ischemic kidney, highlighting the need for increased attention to the inner medulla and potential opportunities for new insights beyond those available from superficial biopsy-focused tissue examination.

Animals