Transcriptomic Fingerprint Research Articles

Spatial transcriptomic clocks reveal cell proximity effects in brain ageing.

Old age is associated with a decline in cognitive function and an increase in neurodegenerative disease risk1. Brain ageing is complex and is accompanied by many cellular changes2. Furthermore, the influence that aged cells have on neighbouring cells and how this contributes to tissue decline is unknown. More generally, the tools to systematically address this question in ageing tissues have not yet been developed. Here we generate a spatially resolved single-cell transcriptomics brain atlas of 4.2 million cells from 20 distinct ages across the adult lifespan and across two rejuvenating interventions-exercise and partial reprogramming. We build spatial ageing clocks, machine learning models trained on this spatial transcriptomics atlas, to identify spatial and cell-type-specific transcriptomic fingerprints of ageing, rejuvenation and disease, including for rare cell types. Using spatial ageing clocks and deep learning, we find that T cells, which increasingly infiltrate the brain with age, have a marked pro-ageing proximity effect on neighbouring cells. Surprisingly, neural stem cells have a strong pro-rejuvenating proximity effect on neighbouring cells. We also identify potential mediators of the pro-ageing effect of T cells and the pro-rejuvenating effect of neural stem cells on their neighbours. These results suggest that rare cell types can have a potent influence on their neighbours and could be targeted to counter tissue ageing. Spatial ageing clocks represent a useful tool for studying cell-cell interactions in spatial contexts and should allow scalable assessment of the efficacy of interventions for ageing and disease.

A Single Cell Transcriptomic Fingerprint of Stressed Premature, Imbalanced Differentiation of Embryonic Stem Cells.

Miscarriages cause a greater loss-of-life than cardiovascular diseases, but knowledge about environmentally induced miscarriages is limited. Cultured naïve pluripotent embryonic stem cells (ESC) differentiate into extra-embryonic endoderm/extraembryonic endoderm (XEN) or formative pluripotent ESC, during the period emulating maximal miscarriage of peri-implantation development. In previous reports using small marker sets, hyperosmotic sorbitol, or retinoic acid (RA) decreased naïve pluripotency and increased XEN by FACS quantitation. Bulk and single cell (sc)RNAseq analyses of two cultured ESC lines was done, corroborated by qPCR. Transcriptomic responses were analyzed of cultured ESC stressed by Sorbitol, with Leukemia inhibitory factor (LIF + ; stemness growth factor), RA without LIF to control for XEN induction, and compared with normal differentiation (LIF - , ND). Sorbitol and RA increase subpopulations of 2-cell embryo-like (2CEL) and XEN sub-lineages; primitive, parietal, and visceral endoderm (VE) cells and suppress formative pluripotency, imbalancing alternate lineage choices of initial naïve pluripotent cultured ESC compared with ND. Although bulk RNAseq and gene ontology (GO) group analyses suggest that stress induces anterior VE-head organizer and placental markers, scRNAseq reveals relatively few cells. But VE and placental markers/cells were in adjacent stressed cell clusters in the UMAP, like recent, normal UMAP of conceptuses. UMAPs show that dose-dependent stress overrides stemness to force premature lineage imbalance. Hyperosmotic stress, and other toxicological stresses, like drugs with active ingredient RA, may cause premature, lineage imbalance, resulting in miscarriages or birth defects.

Spatiotemporal transcriptomic profiling and modeling of mouse brain at single-cell resolution reveals cell proximity effects of aging and rejuvenation.

Old age is associated with a decline in cognitive function and an increase in neurodegenerative disease risk1. Brain aging is complex and accompanied by many cellular changes2-20. However, the influence that aged cells have on neighboring cells and how this contributes to tissue decline is unknown. More generally, the tools to systematically address this question in aging tissues have not yet been developed. Here, we generate spatiotemporal data at single-cell resolution for the mouse brain across lifespan, and we develop the first machine learning models based on spatial transcriptomics ('spatial aging clocks') to reveal cell proximity effects during brain aging and rejuvenation. We collect a single-cell spatial transcriptomics brain atlas of 4.2 million cells from 20 distinct ages and across two rejuvenating interventions-exercise and partial reprogramming. We identify spatial and cell type-specific transcriptomic fingerprints of aging, rejuvenation, and disease, including for rare cell types. Using spatial aging clocks and deep learning models, we find that T cells, which infiltrate the brain with age, have a striking pro-aging proximity effect on neighboring cells. Surprisingly, neural stem cells have a strong pro-rejuvenating effect on neighboring cells. By developing computational tools to identify mediators of these proximity effects, we find that pro-aging T cells trigger a local inflammatory response likely via interferon-γ whereas pro-rejuvenating neural stem cells impact the metabolism of neighboring cells possibly via growth factors (e.g. vascular endothelial growth factor) and extracellular vesicles, and we experimentally validate some of these predictions. These results suggest that rare cells can have a drastic influence on their neighbors and could be targeted to counter tissue aging. We anticipate that these spatial aging clocks will not only allow scalable assessment of the efficacy of interventions for aging and disease but also represent a new tool for studying cell-cell interactions in many spatial contexts.

Computational approaches identify a transcriptomic fingerprint of drug-induced structural cardiotoxicity

Structural cardiotoxicity (SCT) presents a high-impact risk that is poorly tolerated in drug discovery unless significant benefit is anticipated. Therefore, we aimed to improve the mechanistic understanding of SCT. First, we combined machine learning methods with a modified calcium transient assay in human-induced pluripotent stem cell-derived cardiomyocytes to identify nine parameters that could predict SCT. Next, we applied transcriptomic profiling to human cardiac microtissues exposed to structural and non-structural cardiotoxins. Fifty-two genes expressed across the three main cell types in the heart (cardiomyocytes, endothelial cells, and fibroblasts) were prioritised in differential expression and network clustering analyses and could be linked to known mechanisms of SCT. This transcriptomic fingerprint may prove useful for generating strategies to mitigate SCT risk in early drug discovery.Graphical

Differences in the neural correlates of schizophrenia with positive and negative formal thought disorder in patients with schizophrenia in the ENIGMA dataset

Formal thought disorder (FTD) is a clinical key factor in schizophrenia, but the neurobiological underpinnings remain unclear. In particular, the relationship between FTD symptom dimensions and patterns of regional brain volume loss in schizophrenia remains to be established in large cohorts. Even less is known about the cellular basis of FTD. Our study addresses these major obstacles by enrolling a large multi-site cohort acquired by the ENIGMA Schizophrenia Working Group (752 schizophrenia patients and 1256 controls), to unravel the neuroanatomy of FTD in schizophrenia and using virtual histology tools on implicated brain regions to investigate the cellular basis. Based on the findings of previous clinical and neuroimaging studies, we decided to separately explore positive, negative and total formal thought disorder. We used virtual histology tools to relate brain structural changes associated with FTD to cellular distributions in cortical regions. We identified distinct neural networks positive and negative FTD. Both networks encompassed fronto-occipito-amygdalar brain regions, but positive and negative FTD demonstrated a dissociation: negative FTD showed a relative sparing of orbitofrontal cortical thickness, while positive FTD also affected lateral temporal cortices. Virtual histology identified distinct transcriptomic fingerprints associated for both symptom dimensions. Negative FTD was linked to neuronal and astrocyte fingerprints, while positive FTD also showed associations with microglial cell types. These results provide an important step towards linking FTD to brain structural changes and their cellular underpinnings, providing an avenue for a better mechanistic understanding of this syndrome.

A Transcriptomic Fingerprint for Parkinson's Disease Found in Patients with Chronic Inflammatory Response Syndrome: Implications for Diagnosis, Treatment and Prevention

A Transcriptomic Fingerprint for Parkinson's Disease Found in Patients with Chronic Inflammatory Response Syndrome: Implications for Diagnosis, Treatment and Prevention

Reduced Expression of Oligodendrocyte Linage-Enriched Transcripts During the Endoplasmic Reticulum Stress/Integrated Stress Response.

Endoplasmic reticulum (ER) stress in oligodendrocyte (OL) linage cells contributes to several CNS pathologies including traumatic spinal cord injury (SCI) and multiple sclerosis. Therefore, primary rat OL precursor cell (OPC) transcriptomes were analyzed using RNASeq after treatments with two ER stress-inducing drugs, thapsigargin (TG) or tunicamycin (TM). Gene ontology term (GO) enrichment showed that both drugs upregulated mRNAs associated with the general stress response. The GOs related to ER stress were only enriched for TM-upregulated mRNAs, suggesting greater ER stress selectivity of TM. Both TG and TM downregulated cell cycle/cell proliferation-associated transcripts, indicating the anti-proliferative effects of ER stress. Interestingly, many OL lineage-enriched mRNAs were downregulated, including those for transcription factors that drive OL identity such as Olig2. Moreover, ER stress-associated decreases of OL-specific gene expression were found in mature OLs from mouse models of white matter pathologies including contusive SCI, toxin-induced demyelination, and Alzheimer's disease-like neurodegeneration. Taken together, the disrupted transcriptomic fingerprint of OL lineage cells may facilitate myelin degeneration and/or dysfunction when pathological ER stress persists in OL lineage cells.

Transcriptomic and proteomic fingerprints induced by the fungicides difenoconazole and metalaxyl in zebrafish embryos

In this study, we applied OMICs analysis to identify substance-specific biomarker candidates, which may act as early indicators for specific ecotoxic modes of actions (MoA). Zebrafish embryos were exposed to two sublethal concentrations of difenoconazole and metalaxyl according to a modified protocol of the OECD test guideline No. 236. At the end of exposure, total RNA and protein were extracted, followed by transcriptomics and proteomics analysis. The analysis of significantly differentially expressed genes (DEGs) and differentially expressed proteins (DEPs) revealed a positive exposure-response correlation in all test concentrations for both fungicides. Similarly, also a positive correlation between the obtained transcriptome and proteome data was observed, highlighting the robustness of our approach. From the detected DEGs, candidate biomarkers specific for difenoconazole (apoa1b, gatm, mylpfb and acta1b) and metalaxyl (lgals2b, abat, fabp1b.1 and myh9a) were selected, and their biological functions were discussed to assess the predictive potential.

Isothiazolinone dysregulates the pattern of miRNA secretion: Endocrine implications for neurogenesis

Isothiazolinones are extensively used as preservatives and disinfectants in personal care products and household items. The unintended exposure of humans and animals to isothiazolinones has led to increasing concerns about their health hazards. The compound 4,5-Dichloro-2-n-octyl-4-isothiazolin-3-one (DCOIT), a representative isothiazolinone, can simultaneously induce endocrine disruption and neurotoxicity. However, the underlying mechanisms and linkages remain unclear. Our purpose was to elucidate the role of miRNAs as the signaling communicator during the crosstalk between endocrine and nervous systems in response to DCOIT stress. H295R cells were exposed to DCOIT, after which the alterations in intracellular miRNA composition, exosome secretory machinery, and extracellular miRNA composition were examined. Then, a PC12 cell line of neuronal differentiation potential was cultured with the extract of extracellular miRNAs from DCOIT-exposed H295R cell media to explore the functional implications in neurogenesis. The results showed that DCOIT exposure resulted in 349 differentially expressed miRNAs (DEMs) in H295R cells, which were closely related to the regulation of multiple endocrine pathways. In the media of H295R cells exposed to DCOIT, 66 DEMs were identified, showing distinct compositions compared to intracellular DEMs with only 2 common DEMs (e.g., novel-m0541-5p of inverse changes in the cell and medium). Functional annotation showed that extracellular DEMs were not only associated with sex endocrine synchronization, but were also implicated in nervous system development, morphogenesis, and tumor. Incubating PC12 cells with the extracellular exosomes (containing miRNAs) from DCOIT-exposed H295R cells significantly increased the neurite growth, promoted neuronal differentiation, and shaped the transcriptomic fingerprint, implying that miRNAs may communicate transduction of toxic information of DCOIT in endocrine system to neurons. Overall, the present findings provide novel insight into the endocrine disrupting and neural toxicity of DCOIT. The miRNAs have the potential to serve as the epigenetic mechanism of systems toxicology.

Transcriptomic insights into the roles of the transcription factors Clr1, Clr2 and Clr4 in lignocellulose degradation of the thermophilic fungal platform Thermothelomyces thermophilus.

Introduction: Thermothelomyces thermophilus, formerly known as Myceliophthora thermophila, is used in industry to produce lignocellulolytic enzymes and heterologous proteins. However, the transcriptional network driving the expression of these proteins remains elusive. As a first step to systematically uncover this network, we investigated growth, protein secretion, and transcriptomic fingerprints of strains deficient in the cellulolytic transcriptional regulators Clr1, Clr2, and Clr4, respectively. Methods: The genes encoding Clr1, Clr2, and Clr4 were individually deleted using split marker or the CRISPR/Cas12a technology and the resulting strains as well as the parental strain were cultivated in bioreactors under chemostat conditions using glucose as the carbon source. During steady state conditions, cellulose was added instead of glucose to study the genetic and cellular responses in all four strains to the shift in carbon source availability. Results: Notably, the clr1 and clr2 deletion strains were unable to continue to grow on cellulose, demonstrating a key role of both regulators in cellulose catabolism. Their transcriptomic fingerprints uncovered not only a lack of cellulase gene expression but also reduced expression of genes predicted to encode hemicellulases, pectinases, and esterases. In contrast, the growth of the clr4 deletion strain was very similar compared to the parental strain. However, a much stronger expression of cellulases, hemicellulases, pectinases, and esterases was observed. Discussion: The data gained in this study suggest that both transcriptional regulators Clr1 and Clr2 activate the expression of genes predicted to encode cellulases as well as hemicellulases, pectinases, and esterases. They further suggest that Clr1 controls the basal expression of cellulases and initiates the main lignocellulolytic response to cellulose via induction of clr2 expression. In contrast, Clr4 seems to act as a repressor of the lignocellulolytic response presumably via controlling clr2 expression. Comparative transcriptomics in all four strains revealed potentially new regulators in carbohydrate catabolism and lignocellulolytic enzyme expression that define a candidate gene list for future analyses.

Antibody-mediated neutralization of galectin-3 as a strategy for the treatment of systemic sclerosis

Systemic sclerosis (SSc) is an autoimmune, inflammatory and fibrotic disease with limited treatment options. Developing new therapies is therefore crucial to address patient needs. To this end, we focused on galectin-3 (Gal-3), a lectin known to be associated with several pathological processes seen in SSc. Using RNA sequencing of whole-blood samples in a cross-sectional cohort of 249 patients with SSc, Gal-3 and its interactants defined a strong transcriptomic fingerprint associated with disease severity, pulmonary and cardiac malfunctions, neutrophilia and lymphopenia. We developed new Gal-3 neutralizing monoclonal antibodies (mAb), which were then evaluated in a mouse model of hypochlorous acid (HOCl)-induced SSc. We show that two of these antibodies, D11 and E07, reduced pathological skin thickening, lung and skin collagen deposition, pulmonary macrophage content, and plasma interleukin-5 and -6 levels. Moreover, E07 changed the transcriptional profiles of HOCl-treated mice, resulting in a gene expression pattern that resembled that of control mice. Similarly, pathological pathways engaged in patients with SSc were counteracted by E07 in mice. Collectively, these findings demonstrate the translational potential of Gal-3 blockade as a therapeutic option for SSc.

A single-cell RNA-seq analysis unravels the heterogeneity of primary cultured human corneal endothelial cells

The cornea is a transparent and avascular tissue located in front of the eye. Its inner surface is lined by a monolayer of corneal endothelial cells (CECs), which maintain the cornea transparency. CECs remain arrested in a non-proliferative state and damage to these cells can compromise their function leading to corneal opacity. The primary culture of donor-derived CECs is a promising cell therapy. It confers the potential to treat multiple patients from a single donor, alleviating the global donor shortage. Nevertheless, this approach has limitations preventing its adoption, particularly culture protocols allow limited expansion of CECs and there is a lack of clear parameters to identify therapy-grade CECs. To address this limitation, a better understanding of the molecular changes arising from the primary culture of CECs is required. Using single-cell RNA sequencing on primary cultured CECs, we identify their variable transcriptomic fingerprint at the single cell level, provide a pseudo-temporal reconstruction of the changes arising from primary culture, and suggest markers to assess the quality of primary CEC cultures. This research depicts a deep transcriptomic understanding of the cellular heterogeneity arising from the primary expansion of CECs and sets the basis for further improvement of culture protocols and therapies.

Microisolation of Spatially Characterized Single Populations of Neurons for RNA Sequencing from Mouse and Postmortem Human Brain Tissues.

Single-cell and single-population RNA sequencing (RNA-seq) is a rapidly evolving new field of intense investigation. Recent studies indicate unique transcriptomic profiles are derived based on the spatial localization of neurons within circuits and regions. Individual neuronal subtypes can have vastly different transcriptomic fingerprints, well beyond the basic excitatory neuron and inhibitory neuron designations. To study single-population gene expression profiles of spatially characterized neurons, we have developed a methodology combining laser capture microdissection (LCM), RNA purification of single populations of neurons, and subsequent library preparation for downstream applications, including RNA-seq. LCM provides the benefit of isolating single neurons characterized by morphology or via transmitter-identified and/or receptor immunoreactivity and enables spatial localization within the sample. We utilize unfixed human postmortem and mouse brain tissue that is frozen to preserve RNA quality in order to isolate the desired neurons of interest. Microisolated neurons are then pooled for RNA purification utilizing as few as 250 individual neurons from a tissue section, precluding extraneous nonspecific tissue contaminants. Library preparation is performed from picogram RNA quantities extracted from LCM-captured neurons. Single-population RNA-seq analysis demonstrates that microisolated neurons from both postmortem human and mouse brain tissues are viable for transcriptomic profiling, including differential gene expression assessment and bioinformatic pathway inquiry.

Transcriptomic fingerprint of neurons and glial cells isolated and differentiated from olfactory neuroepithelium of subjects with schizophrenia

Transcriptomic fingerprint of neurons and glial cells isolated and differentiated from olfactory neuroepithelium of subjects with schizophrenia

Transcriptomic Interaction between Young Fecal Transplantation and Perfluorobutanesulfonate in Aged Zebrafish Gonads.

The transfer of young fecal microbiota has been found to significantly refresh the reproductive endocrine system and effectively ameliorate the toxicity of perfluorobutanesulfonate (PFBS) in aged zebrafish recipients. However, the mechanisms underlying the antagonistic action of young fecal microbiota against the reproductive endocrine toxicity of PFBS remain largely unknown. In this study, the aged zebrafish were transplanted with feces from young donors and then exposed to PFBS for 14 days. After exposure, the shift in the transcriptomic fingerprint of the gonads was profiled by using high-throughput sequencing, aiming to provide mechanistic clues into the interactive mode of action between young fecal transplantation and PFBS's innate toxicity. The results showed that the gene transcription pattern associated with protein and lipid synthesis in the gonads of the aged individuals was quite different from the young counterparts. It was intriguing that the transplantation of young feces established a youth-like transcriptomic phenotype in the elderly recipients, thus attenuating the functional decline and maintaining a healthy aging state of the gonads. A sex specificity response was clearly observed. Compared to the aged females, more metabolic pathways (e.g., glycine, serine, and threonine metabolism; glyoxylate and dicarboxylate metabolism; pyrimidine metabolism) were significantly enriched in aged males receiving young feces transplants. PFBS dramatically altered the transcriptome of aged testes, while a much milder effect was observable in aged ovaries. Accordingly, a suite of biological processes related to germ cell proliferation were disrupted by PFBS in aged males, including the ECM-receptor interaction, retinol metabolism, and folate biosynthesis. In aged ovaries exposed to PFBS, mainly the fatty acid and arginine biosynthesis pathway was significantly affected. However, these transcriptomic disorders caused by PFBS were largely mitigated in aged gonads by transferring young feces. Overall, the present findings highlighted the potential of young fecal transplantation to prevent the functional compromise of gonads resulting from aging and PFBS.

Distinct immune signature predicts progression of vestibular schwannoma and unveils a possible viral etiology

BackgroundThe management of sub-totally resected sporadic vestibular schwannoma (VS) may include observation, re-resection or irradiation. Identifying the optimal choice can be difficult due to the disease’s variable progression rate.We aimed to define an immune signature and associated transcriptomic fingerprint characteristic of rapidly-progressing VS to elucidate the underpinnings of rapidly progressing VS and identify a prognostic model for determining rate of progression.MethodsWe used multiplex immunofluorescence to characterize the immune microenvironment in 17 patients with sporadic VS treated with subtotal surgical resection alone. Transcriptomic analysis revealed differentially-expressed genes and dysregulated pathways when comparing rapidly-progressing VS to slowly or non-progressing VS.ResultsRapidly progressing VS was distinctly enriched in CD4+, CD8+, CD20+, and CD68+ immune cells. RNA data indicated the upregulation of anti-viral innate immune response and T-cell senescence. K − Top Scoring Pair analysis identified 6 pairs of immunosenescence-related genes (CD38-KDR, CD22-STAT5A, APCS-CXCR6, MADCAM1-MPL, IL6-NFATC3, and CXCL2-TLR6) that had high sensitivity (100%) and specificity (78%) for identifying rapid VS progression.ConclusionRapid progression of residual vestibular schwannoma following subtotal surgical resection has an underlying immune etiology that may be virally originating; and despite an abundant adaptive immune response, T-cell immunosenescence may be associated with rapid progression of VS. These findings provide a rationale for clinical trials evaluating immunotherapy in patients with rapidly progressing VS.

A conserved transcriptional fingerprint of multi-neurotransmitter neurons necessary for social behavior

BackgroundAn essential determinant of a neuron’s functionality is its neurotransmitter phenotype. We previously identified a defined subpopulation of cholinergic neurons required for social orienting behavior in zebrafish.ResultsWe transcriptionally profiled these neurons and discovered that they are capable of synthesizing both acetylcholine and GABA. We also established a constellation of transcription factors and neurotransmitter markers that can be used as a “transcriptomic fingerprint” to recognize a homologous neuronal population in another vertebrate.ConclusionOur results suggest that this transcriptomic fingerprint and the cholinergic-GABAergic neuronal subtype that it defines are evolutionarily conserved.

Gene Expression Signatures Reveal Common Virus Infection Pathways in Target Tissues of Type 1 Diabetes, Hashimoto's Thyroiditis, and Celiac Disease.

Type 1 diabetes (T1D) patients are at heightened risk for other autoimmune disorders, particularly Hashimoto’s thyroiditis (HT) and celiac disease (CD). Recent evidence suggests that target tissues of autoimmune diseases engage in a harmful dialogue with the immune system. However, it is unclear whether shared mechanisms drive similar molecular signatures at the target tissues among T1D, HT, and CD. In our current study, microarray datasets were obtained and mined to identify gene signatures from disease-specific targeted tissues including the pancreas, thyroid, and intestine from individuals with T1D, HT, and CD, as well as their matched controls. Further, the threshold-free algorithm rank-rank hypergeometric overlap analysis (RRHO) was used to compare the genomic signatures of the target tissues of the three autoimmune diseases. Next, promising drugs that could potentially reverse the observed signatures in patients with two or more autoimmune disorders were identified using the cloud-based CLUE software platform. Finally, microarray data of auto-antibody positive individuals but not diagnosed with T1D and single cell sequencing data of patients with T1D and HT were used to validate the shared transcriptomic fingerprint. Our findings revealed significant common gene expression changes in target tissues of the three autoimmune diseases studied, many of which are associated with virus infections, including influenza A, human T-lymphotropic virus type 1, and herpes simplex infection. These findings support the importance of common environmental factors in the pathogenesis of T1D, HT, and CD.

Differences and Similarities between the Lung Transcriptomic Profiles of COVID-19, COPD, and IPF Patients: A Meta-Analysis Study of Pathophysiological Signaling Pathways.

Coronavirus disease 2019 (COVID-19) is a pandemic respiratory disease associated with high morbidity and mortality. Although many patients recover, long-term sequelae after infection have become increasingly recognized and concerning. Among other sequelae, the available data indicate that many patients who recover from COVID-19 could develop fibrotic abnormalities over time. To understand the basic pathophysiology underlying the development of long-term pulmonary fibrosis in COVID-19, as well as the higher mortality rates in patients with pre-existing lung diseases, we compared the transcriptomic fingerprints among patients with COVID-19, idiopathic pulmonary fibrosis (IPF), and chronic obstructive pulmonary disease (COPD) using interactomic analysis. Patients who died of COVID-19 shared some of the molecular biological processes triggered in patients with IPF, such as those related to immune response, airway remodeling, and wound healing, which could explain the radiological images seen in some patients after discharge. However, other aspects of this transcriptomic profile did not resemble the profile associated with irreversible fibrotic processes in IPF. Our mathematical approach instead showed that the molecular processes that were altered in COVID-19 patients more closely resembled those observed in COPD. These data indicate that patients with COPD, who have overcome COVID-19, might experience a faster decline in lung function that will undoubtedly affect global health.

Transcriptomic fingerprint of bacterial infection in lower extremity ulcers

Clinicians and researchers utilize subjective, clinical classification systems to stratify lower extremity ulcer infections for treatment and research. The purpose of this study was to examine whether these clinical classifications are reflected in the ulcer's transcriptome. RNA sequencing (RNA‐seq) was performed on biopsies from clinically infected lower extremity ulcers (n = 44). Resulting sequences were aligned to the host reference genome to create a transcriptome profile. Differential gene expression analysis and gene ontology (GO) enrichment analysis were performed between ulcer severities as well as between sample groups identified by k‐means clustering. Lastly, a support vector classifier was trained to estimate clinical infection score or k‐means cluster based on a subset of genes. Clinical infection severity did not explain the major sources of variability among the samples and samples with the same clinical classification demonstrated high inter‐sample variability. High proportions of bacterial RNA were identified in some samples, which resulted in a strong effect on transcription and increased expression of genes associated with immune response and inflammation. K‐means clustering identified two clusters of samples, one of which contained all of the samples with high levels of bacterial RNA. A support vector classifier identified a fingerprint of 20 genes, including immune‐associated genes such as CXCL8, GADD45B, and HILPDA, which accurately identified samples with signs of infection via cross‐validation. This study identified a unique, host‐transcriptome signature in the presence of infecting bacteria, often incongruent with clinical infection‐severity classifications. This suggests that stratification of infection status based on a transcriptomic fingerprint may be useful as an objective classification method to classify infection severity, as well as a tool for studying host–pathogen interactions.