Transcriptomic Patterns Research Articles

Deciphering transcriptome patterns in porcine mesenchymal stem cells promoting phenotypic maintenance and differentiation by key driver genes

Mesenchymal stem cells (MSC) are fibroblast-like non-hematopoietic cells with self-renewal and differentiation capacity, and thereby great potential in regeneration and wound healing. MSC populations are heterogeneous not only inherently, but also among different model species. In particular, porcine MSC serve as a frequently used resource for translational research, due to pigs’ distinctive closeness to human anatomy and physiology. However, information on gene expression profiles from porcine MSC and its dynamics during differentiation is sparse, especially with regard to cell surface and inner cell markers. In this study, we investigated the transcriptome of bone marrow-derived MSC and its differentiated cell types in a minipig breed for experimental research, known as Mini-LEWE, using bulk mRNA sequencing. Our data highlighted Rap1 signaling and downstream pathways PI3K-Akt and MAPK signaling as potential players for the maintenance of stemness of BM-MSC. In addition, we were able to link the process of differentiation to changes in the regulation of actin cytoskeleton. A total of 18 “BM-MSC differentiation driver markers” were identified, potentially promoting the process of differentiation into adipocytes, chondrocytes as well as osteocytes. Our results offer a new perspective on the molecular phenotype of porcine BM-MSC and the transcriptional responses in new differentiated progeny.

Transcriptomic profiling of 'Candidatus Liberibacter asiaticus' in different citrus tissues reveals novel insights into Huanglongbing pathogenesis.

'Candidatus Liberibacter asiaticus' (Las) is a gram-negative bacterial pathogen associated with citrus huanglongbing (HLB) or greening disease. Las is transmitted by the Asian citrus psyllid (ACP) where it colonizes the phloem tissue, resulting in substantial economic losses to citrus industry worldwide. Despite extensive efforts, effective management strategies against HLB remain elusive, necessitating a deeper understanding of the pathogen's biology. Las undergoes cell-to-cell movement through phloem flow and colonizes different tissues in which Las may have varying interactions with the host. Here, we investigate the transcriptomic landscape of Las in citrus seed coat vasculatures, enabling a complete gene expression profiling of Las genome and revealing unique transcriptomic patterns compared to previous studies using midrib tissues. Comparative transcriptomics between seed coat, midrib and ACP identified specific responses and metabolic states of Las in different host tissue. Two Las virulence factors that exhibit higher expression in seed coat can suppress callose deposition. Therefore, they may contribute to unclogging sieve plate pores during Las colonization in seed coat vasculature. Furthermore, analysis of regulatory elements uncovers a potential role of LuxR-type transcription factors in regulating Liberibacter effector gene expression during plant colonization. Together, this work provides novel insights into the pathogenesis of the devastating citrus HLB.

Shared and unique patterns of autonomous human endogenous retrovirus loci transcriptomes in CD14 + monocytes from individuals with physical trauma or infection with COVID-19

Since previous studies have suggested that the RNAs of human endogenous retrovirus (HERV) might be involved in regulating innate immunity, it is important to investigate the HERV transcriptome patterns in innate immune cell types such as CD14 + monocytes. Using single cell RNA-seq datasets from resting or stimulated PBMCs mapped to 3,220 known discrete autonomous proviral HERV loci, we found individual-specific variation in HERV transcriptomes between HERV loci in CD14 + monocytes. Analysis of paired datasets from the same individual that were cultured in vitro with LPS or without (i.e. control) revealed 36 HERV loci in CD14 + monocytes that were detected only after activation. To extend our analysis to in vivo activated CD14 + monocytes, we used two scRNA-seq datasets from studies that had demonstrated activation of circulating CD14 + monocytes in patients with physical trauma or patients hospitalized with COVID-19 infections. For direct comparison between the trauma and COVID-19 datasets, we first analyzed 1.625 billion sequence reads from a composite pangenome control of 21 normal individuals. Comparison of the sequence read depth of HERV loci in the trauma or COVID-19 samples to the pangenome control revealed that 39 loci in the COVID-19 and 11 HERV loci in the trauma samples were significantly different (Mann-Whitney U test), with 9 HERV loci shared between the COVID-19 and trauma datasets. The capacity to compare HERV loci transcriptome patterns in innate immune cells, like CD14 + monocytes, across different pathological conditions will lead to greater understanding of the physiological role of HERV expression in health and disease.

Brain connectivity and transcriptomic similarity inform abnormal morphometric similarity patterns in first-episode, treatment-naïve major depressive disorder

Brain connectivity and transcriptomic similarity inform abnormal morphometric similarity patterns in first-episode, treatment-naïve major depressive disorder

Genome-wide identification and analysis of the HD-Zip transcription factor family in oats

IntroductionHD-Zip transcription factors are an important class of plant transcription factors involved in regulating plant growth and development as well as various stress responses. In order to explore the characteristics of oat HD-Zip transcription factors and their transcriptomic expression patterns under abiotic stress, this study identified members of the HD-Zip gene family in the oat genome through bioinformatics methods and analyzed their basic physicochemical properties, evolutionary relationships, conserved structural domains, gene duplication relationships, and expression profiles.Results74 HD-Zip gene sequences were identified in the oat genome, unevenly distributed in all chromosomes except the 4D chromosome. The 74 HD-Zip genes can be divided into four subfamilies (HD-Zip Ⅰ-Ⅳ), containing 30 (HD-Zip Ⅰ), 38 (HD-Zip Ⅱ), 4 (HD-Zip III), and 2 (HD-Zip IV) genes, respectively. The grouping of this study is completely consistent with the clustering results of family members in the oat HD-Zip gene phylogenetic tree, further supporting the reliability of the sequence grouping. In addition, there are significant differences in conserved motifs and gene lengths between subfamilies, but they are conserved within the same subfamily. Under drought and salt stress, the expression levels of a large number of oat HD-Zip genes were significantly induced or suppressed, indicating that these HD-Zip genes are likely to be involved in regulating the oat’s response to adversity.ConclusionsThe results of this study provide excellent candidate HD-Zip genes for the study of oat resistance, and provide a reference for oat gene improvement and genetic breeding.

Effects of simulated microgravity on colorectal cancer organoids growth and drug response

Cellular and molecular dynamics of human cells are constantly affected by gravity. Alteration of the gravitational force disturbs the cellular equilibrium, which might modify physiological and molecular characteristics. Nevertheless, biological responses of cancer cells to reduced gravitational force remains obscure. Here, we aimed to comprehend not only transcriptomic patterns but drug responses of colorectal cancer (CRC) under simulated microgravity. We established four organoids directly from CRC patients, and organoids cultured in 3D clinostat were subjected to genome wide expression profiling and drug library screening. Our observations revealed changes in cell morphology and an increase in cell viability under simulated microgravity compared to their static controls. Transcriptomic analysis highlighted a significant dysregulation in the TBC1D3 family of genes. The upregulation of cell proliferation observed under simulated microgravity conditions was further supported by enriched cell cycle processes, as evidenced by the functional clustering of mRNA expressions using cancer hallmark and gene ontology terms. Our drug screening results indicated an enhanced response rate to 5-FU under conditions of simulated microgravity, suggesting potential implications for cancer treatment strategies in simulated microgravity.

Distinct Cutibacterium acnes subspecies defendens strains classified by multi-omics dissection alleviate inflammatory skin lesions of a rosacea-like mouse model

IntroductionCutibacterium acnes (C. acnes) resides in various organs such as the skin, prostate, eye, nose, stomach, and intestine, indicating the possibility of extensive crosstalk between this bacterium and the human body. C. acnes strains are classified into three subspecies based on phylogenetics and distinguishable phenotypes. Among them, C. acnes subsp. defendens strains are characterized by anti-inflammatory features, raising expectations for their potential as future microbiome therapeutics. However, the heterogeneity of C. acnes subsp. defendens and its corresponding immunological functions have not been clearly addressed.MethodsThe genetic diversity of the strains was assessed using single- and multi-locus sequence typing. Their immune-modulatory functions were evaluated in vitro using 2D and 3D assays with immune and epithelial cells. The anti-inflammatory effects were further confirmed in vivo using a rosacea-like mouse model. Comparative genomic and transcriptomic analyses were conducted to uncover mechanisms underlying the immunosuppressive activity of the strains.ResultsWe demonstrated that the newly isolated C. acnes subsp. defendens strains, exhibiting phenotypic heterogeneity, are distinctly clustered using single- and multi-locus sequence typing methods. These strains showed strong immune-regulatory functions in immune and epithelial cell-based 2D and 3D in vitro assays. Furthermore, their anti-inflammatory role was functionally confirmed in vivo using a rosacea-like mouse model, where they alleviated skin lesions characterized by hyperplasia and dermal inflammation. Comparative transcriptomics revealed that these strains may exert their immunosuppressive effects through the enhanced expression of acnecins and transcriptional variation in envelope stress regulators (specifically the two-component systems, CesSR homologs). Additionally, we propose that these C. acnes type II strains produce anti-inflammatory metabolites or peptides smaller than 3 kDa, which are associated with elevated pyrimidine and reduced L-arginine biosynthesis.DiscussionThe newly isolated C. acnes subsp. defendens strains demonstrate significant anti-inflammatory properties both in vitro and in vivo, suggesting their potential as microbiome-based therapeutics. Their unique genomic and transcriptomic profiles, including the production of small bioactive compounds and specific transcriptomic patterns, underpin their immunosuppressive capabilities. These findings provide a foundation for developing novel treatments for inflammatory skin conditions, such as rosacea.

Pathological autoantibody internalisation in myositis

ObjectivesAutoantibodies targeting intracellular proteins are common in various autoimmune diseases. In the context of myositis, the pathologic significance of these autoantibodies has been questioned due to the assumption that autoantibodies...

Transcriptome dynamics in developing leaves from C3 and C4 Flaveria species.

C4 species have evolved more than 60 times independently from C3 ancestors. This multiple and parallel evolution of the complex C4 trait suggests common underlying evolutionary mechanisms, which could be identified by comparative analysis of closely related C3 and C4 species. Efficient C4 function depends on a distinctive leaf anatomy that is characterised by enlarged, chloroplast-rich bundle sheath cells and narrow vein spacing. To elucidate the molecular mechanisms that generate the Kranz anatomy, we analysed a developmental series of leaves from the C4 plant Flaveria bidentis and the closely related C3 species Flaveria robusta by comparing anatomies and transcriptomes. Vascular density measurements of all nine leaf developmental stages identified three leaf anatomical zones whose proportions vary with respect to the developmental stage. We then deconvoluted the transcriptome datasets using non-negative matrix factorisation, which identified four distinct transcriptome patterns in the growing leaves of both species. By integrating the leaf anatomy and transcriptome data, we were able to correlate the different transcriptional profiles with different developmental zones in the leaves. These comparisons revealed an important role for auxin metabolism, in particular auxin homeostasis (conjugation and deconjugation), in establishing the high vein density typical of C4 species.

Inhibition of TGF-β signaling in bone marrow endothelial cells promotes hematopoietic recovery in acute myeloid leukemia patients

Although it is an effective treatment for acute myeloid leukemia (AML), chemotherapy leads to myelosuppression and poor hematopoietic reconstruction. Hematopoiesis is regulated by bone marrow (BM) endothelial cells (ECs), and BM ECs are dysfunctional in acute leukemia patients with poor hematopoietic reconstitution after allogenic hematopoietic stem cell transplantation. Thus, it is crucial to explore the underlying mechanism of EC impairment and establish strategies for targeted therapy. TGF-β signaling was found to be upregulated in ECs from AML patients in complete remission (CR ECs) and led to CR EC damage. Administration of a TGF-β inhibitor rescued the dysfunction of ECs caused by TGF-β1 expression in vitro, especially their hematopoiesis-supporting ability. Moreover, inhibition of TGF-β expression repaired the BM EC damage triggered by chemotherapy in both AML patients in vitro and in an AML-CR murine model, and restored normal hematopoiesis without promoting AML progression. Mechanistically, our data reveal alterations in the transcriptomic pattern of damaged BM ECs, accompanied by the overexpression of downstream molecules TGF-βRI, pSmad2/3, and functional genes related to adhesion, angiogenesis suppression and pro-apoptosis. Collectively, our findings reveal for the first time that the activation of TGF-β signaling leads to BM EC dysfunction and poor hematopoietic reconstitution. Targeting TGF-β represents a potential therapeutic strategy to promote multilineage hematopoiesis, thereby benefiting more cancer patients who suffer from myelosuppression after chemotherapy.

IL23R-specific CAR Tregs for the treatment of Crohn's disease.

Regulatory T cells (Tregs) are key regulators in maintaining tissue homeostasis. Disrupted immune homeostasis is associated with Crohn's disease (CD) pathogenesis. Thus, Treg therapy represents a promising long-acting treatment to restore immune balance in the diseased intestine. CAR (Chimeric Antigen Receptor) T-cell therapy has revolutionized cancer treatment. This innovative approach also provides the opportunity to improve therapy for CD. By targeting a disease-relevant protein, Interleukin-23 receptor (IL23R), we engineered Tregs expressing IL23R-CAR for treating active CD. Intestinal IL23R expression from active CD was verified by immunohistochemical analysis. Phenotypic and functional characteristics of IL23R-CAR Tregs were assessed using in vitro assays and their migration capacity was monitored in a xenograft tumor model. Transcriptomic and proteomic analyses were performed to associate molecular profiles with IL23R-CAR Treg activation against colon biopsy-derived cells from active CD patients. Our study showed that IL23R-CAR displayed negligible tonic signalling and strong signal-to-noise ratio. IL23R-CAR Tregs maintained regulatory phenotype during in vitro expansion, even when chronically exposed to proinflammatory cytokines and target antigen. IL23R engagement on IL23R-CAR Tregs triggered CAR-specific activation and significantly enhanced their suppressive activity. Also, IL23R-CAR Tregs migrated to IL23R-expressing tissue in humanized mice. Finally, IL23R-CAR Tregs elicited a specific activation against colon biopsy-derived cells from active CD, suggesting an efficient CAR engagement in active CD. Molecular profiling of CD patient biopsies also revealed transcriptomic and proteomic patterns associated with IL23R-CAR activation. Overall, our results demonstrate that IL23R-CAR Tregs represent a promising therapy for active CD.

On discovery of novel hub genes for ER+ and TN breast cancer types through RNA seq data analyses and classification models

Breast cancer (BC) is a malignant neoplasm which is classified into various types defined by underlying molecular factors such as estrogen receptor positive (ER+), progesterone receptor positive (PR+), human epidermal growth factor positive (HER2+) and triple negative (TNBC). Early detection of ER+ and TNBC is crucial in the choice of diagnosis and appropriate treatment strategy. Here we report the key genes associated to ER+ and TNBC using RNA-Seq analysis and machine learning models. Three ER+ and TNBC RNA seq datasets comprising 164 patients in-toto were selected for standard NGS hierarchical data processing and data analyses protocols. Enrichment pathway analysis and network analysis was done and finally top hub genes were identified. To come with a reliable classifier which could distinguish the distinct transcriptome patterns associated to ER+ and TNBC, ML models were built employing Naïve Bayes, SVM and kNN. 1730 common DEG’s exhibiting significant logFC values with 0.05 p-value threshold were identified. A list of top ten hub genes were screened on the basis of maximal clique centrality (MCC) which included CDC20, CDK1, BUB1, AURKA, CDCA8, RRM2, TTK, CENPF, CEP55 and NDC80.These genes were found to be involved in crucial cell cycle pathways. k-Nearest Neighbor (kNN) model was observed to be best classifier with accuracy 84%, specificity 66% and sensitivity 95% to differentiate between ER+ and TNBC RNA-Seq transcriptomes. Our screened list of 10 hub genes can thus help unearth novel molecular signatures implicated in ER+ and TNBC onset, prognosis and design of novel protocols for breast cancer diagnostics and therapeutics.

Individual variation in the emergence of anterior-to-posterior neural fates from human pluripotent stem cells

Variability between human pluripotent stem cell (hPSC) lines remains a challenge and opportunity in biomedicine. In this study, hPSC lines from multiple donors were differentiated toward neuroectoderm and mesendoderm lineages. We revealed dynamic transcriptomic patterns that delineate the emergence of these lineages, which were conserved across lines, along with individual line-specific transcriptional signatures that were invariant throughout differentiation. These transcriptomic signatures predicted an antagonism between SOX21-driven forebrain fates and retinoic acid-induced hindbrain fates. Replicate lines and paired adult tissue demonstrated the stability of these line-specific transcriptomic traits. We show that this transcriptomic variation in lineage bias had both genetic and epigenetic origins, aligned with the anterior-to-posterior structure of early mammalian development, and was present across a large collection of hPSC lines. These findings contribute to developing systematic analyses of PSCs to define the origin and consequences of variation in the early events orchestrating individual human development.

Unsupervised clustering of brain tumor leukocytes to reveal atypical lymphoid cell subsets.

25 Background: The brain tumor immune microenvironment (TiME) is characterised by suppression of tumor-infiltrating leukocytes. The lymphocyte niche activity is particularly downregulated within this TiME milieu. Methods: Here, using data-driven approaches in flow cytometry and RNA sequencing analysis, we aimed at dissecting leukocyte cell types infiltrating human brain tumor TiME. Results: Unsupervised clustering provided T lymphocyte subsets including double-negative (CD3+CD4- CD8-) T cells (DNT). Our DNT phenotype was validated with an independent dataset, suggesting DNT to be compatible with γδ T cell phenotype. In our cohort, myeloid immune cell frequencies were associated with malignancy type as previously reported, while a gliosarcoma notably presented unexpectedly high lymphocyte frequencies. Consistently, double-positive (CD4+ CD8+) T cells (DPT) had higher frequency in the gliosarcoma than other tumors. The flow cytometry results were supported by transcriptome patterns and the activity of immune-related genes in these tumors. Conclusions: Our findings feature the benefit of data-driven approaches for brain TiME lymphocyte analysis, and show that primary brain tumors can entail atypical lymphocyte subsets like DNTs and DPT cells. Unbiased lymphocyte characterisation can potentially outline novel targets for immune checkpoint blockade therapy.

Adaptive Evolution and Functional Differentiation of Testis-Expressed Genes in Theria.

Gene expression patterns differ in different tissues, and the expression pattern of genes in the mammalian testis is known to be extremely variable in different species. To clarify how the testis transcriptomic pattern has evolved in particular species, we examined the evolution of the adult testis transcriptome in Theria using 10 species: two marsupials (opossum and Tasmanian devil), six eutherian (placental) mammals (human, chimpanzee, bonobo, gorilla, rhesus macaque, and mouse), and two outgroup species (platypus and chicken). We show that 22 testis-expressed genes are marsupial-specific, suggesting their acquisition in the stem lineage of marsupials after the divergence from eutherians. Despite the time length of the eutherian stem lineage being similar to that of the marsupial lineage, acquisition of testis-expressed genes was not found in the stem lineage of eutherians; rather, their expression patterns differed by species, suggesting rapid gene evolution in the eutherian ancestors. Fifteen testis-expressed genes are therian-specific, and for three of these genes, the evolutionary tempo is markedly faster in eutherians than in marsupials. Our phylogenetic analysis of Rho GTPase-activating protein 28 (ARHGAP28) suggests the adaptive evolution of this gene in the eutherians, probably together with the expression pattern differentiation.

Light has a principal role in the Arabidopsis transcriptomic response to the spaceflight environment

The Characterizing Arabidopsis Root Attractions (CARA) spaceflight experiment provides comparative transcriptome analyses of plants grown in both light and dark conditions within the same spaceflight. CARA compared three genotypes of Arabidopsis grown in ambient light and in the dark on board the International Space Station (ISS); Col-0, Ws, and phyD, a phytochrome D mutant in the Col-0 background. In all genotypes, leaves responded to spaceflight with a higher number of differentially expressed genes (DEGs) than root tips, and each genotype displayed distinct light / dark transcriptomic patterns that were unique to the spaceflight environment. The Col-0 leaves exhibited a substantial dichotomy, with ten-times as many spaceflight DEGs exhibited in light-grown plants versus dark-grown plants. Although the total number of DEGs in phyD leaves is not very different from Col-0, phyD altered the manner in which light-grown leaves respond to spaceflight, and many genes associated with the physiological adaptation of Col-0 to spaceflight were not represented. This result is in contrast to root tips, where a previous CARA study showed that phyD substantially reduced the number of DEGs. There were few DEGs, but a series of space-altered gene categories, common to genotypes and lighting conditions. This commonality indicates that key spaceflight genes are associated with signal transduction for light, defense, and oxidative stress responses. However, these key signaling pathways enriched from DEGs showed opposite regulatory direction in response to spaceflight under light and dark conditions, suggesting a complex interaction between light as a signal, and light-signaling genes in acclimation to spaceflight.

70 Circulating leptin and immune-related transcriptomic patterns in clear cell renal cell carcinoma

Abstract Background Obesity has been associated with better outcomes in localized and metastatic clear cell renal cell carcinoma (ccRCC), as well as improved survival on immunotherapy treatment. However, the mechanisms underlying this association remain unclear. Leptin is an adipokine secreted by adipocytes, and circulating levels of leptin are higher in individuals with obesity. Additionally, leptin promotes inflammation and angiogenesis and has been proposed as pro-tumorigenic in colorectal and breast cancer. In this study, we examined how circulating leptin levels relate to both tumor and perinephric fat transcriptomic patterns in a cohort of ccRCC patients by evaluating the associations with body mass index (BMI), sex, and immune-related gene expression patterns. Methods We conducted a retrospective cohort study of 92 treatment-naïve ccRCC patients undergoing nephrectomy at Memorial Sloan Kettering Cancer Center. Available data included circulating leptin from fasting blood samples, clinical characteristics from medical records, and, in a subset of patients, RNA sequencing from tumor and perinephric fat specimens. We analyzed differentially expressed genes (DEG) according to circulating leptin levels using Gene Set Enrichment Analysis (GSEA) to describe pathways that were enriched in patients with high levels of leptin. We performed analysis for the whole cohort and stratified by sex based on differences in leptin levels. Results From the 92 patients with available peripheral leptin measurements, 51 and 44 had available tumor and perinephric fat RNA sequencing data, respectively. Of these, 64 (70%) were male, the median age was 60 years old, and most tumors were low grade and stage. Leptin distribution was significantly different in males than females, with circulating leptin levels of 7 ng/ml (IQR 4-14) and 22 ng/ml (IQR 9-52), respectively. Higher BMI was associated with higher leptin levels, with a correlation coefficient of 0.63 (p<0.001) in males and 0.77 (p<0.001) in females. GSEA of tumor DEGs by circulating leptin, showed upregulation of pathways related to adaptive immune activity in patients with higher leptin levels across sexes, although this association was attenuated in females. Strikingly, in the perinephric fat there were stark differences in opposite directions in female and male specimens, with females showing significantly enriched transcription of genes associated with B cell signaling, humoral and adaptive immune responses. Conclusions Higher leptin levels were associated with a modest increase of immune-related gene expression in the tumors in both males and females, but significantly different directional changes were observed in the perinephric fat. Circulating leptin may be involved in peritumoral immune responses which may link host factors to sex related tumor outcomes. Further studies should aim to address the relationship between leptin activity and the tumor and fat microenvironment. DOD CDMRP Funding: yes

Barramundi (Lates calcarifer) rare coloration patterns: a multiomics approach to understand the "panda" phenotype.

The barramundi (Lates calcarifer), a significant aquaculture species, typically displays silver to bronze coloration. However, attention is now drawn to rare variants like the "panda" phenotype, characterized by blotch-like patterns of black (PB) and golden (PG) patches. This phenotype presents an opportunity to explore the molecular mechanisms underlying color variations in teleosts. Unlike stable color patterns in many fish, the "panda" variant demonstrates phenotypic plasticity, responding dynamically to unknown cues. We propose a complex interplay of genetic factors and epigenetic modifications, focusing on DNA methylation. Through a multiomics approach, we analyze transcriptomic and methylation patterns between PB and PG patches. Our study reveals differential gene expression related to melanosome trafficking and chromatophore differentiation. Although the specific gene responsible for the PB-PG difference remains elusive, candidate genes like asip1, asip2, mlph, and mreg have been identified. Methylation emerges as a potential contributor to the "panda" phenotype, with changes in gene promoters like hand2 and dynamin possibly influencing coloration. This research lays the groundwork for further exploration into rare barramundi color patterns, enhancing our understanding of color diversity in teleosts. Additionally, it underscores the "panda" phenotype's potential as a model for studying adult skin coloration.

Dynamics of peripheral blood inflammatory index predict tumor pathological response and survival among patients with locally advanced non-small cell lung cancer who underwent neoadjuvant immunochemotherapy: a multi-cohort retrospective study.

Static tumor features before initiating anti-tumor treatment were insufficient to distinguish responding from non-responding tumors under the selective pressure of immuno-therapy. Herein we investigated the longitudinal dynamics of peripheral blood inflammatory indexes (dPBI) and its value in predicting major pathological response (MPR) in non-small cell lung cancer (NSCLC). A total of 147 patients with NSCLC who underwent neoadjuvant immunochemotherapy were retrospectively reviewed as training cohort, and 26 NSCLC patients from a phase II trial were included as validation cohort. Peripheral blood inflammatory indexes were collected at baseline and as posttreatment status; their dynamics were calculated as their posttreatment values minus their baseline level. Least absolute shrinkage and selection operator algorithm was utilized to screen out predictors for MPR, and a MPR score was integrated. We constructed a model incorporating this MPR score and clinical predictors for predicting MPR and evaluated its predictive capacity via the area under the curve (AUC) of the receiver operating characteristic and calibration curves. Furthermore, we sought to interpret this MPR score in the context of micro-RNA transcriptomic analysis in plasma exosomes for 12 paired samples (baseline and posttreatment) obtained from the training cohort. Longitudinal dynamics of monocyte-lymphocyte ratio, platelet-to-lymphocyte ratio, platelet-to-albumin ratio, and prognostic nutritional index were screened out as significant indicators for MPR and a MPR score was integrated, which was further identified as an independent predictor of MPR. Then, we constructed a predictive model incorporating MPR score, histology, and differentiated degree, which discriminated MPR and non-MPR patients well in both the training and validation cohorts with an AUC value of 0.803 and 0.817, respectively. Furthermore, micro-RNA transcriptomic analysis revealed the association between our MPR score and immune regulation pathways. A significantly better event-free survival was seen in subpopulations with a high MPR score. Our findings suggested that dPBI reflected responses to neoadjuvant immuno-chemotherapy for NSCLC. The MPR score, a non-invasive biomarker integrating their dynamics, captured the miRNA transcriptomic pattern in the tumor microenvironment and distinguished MPR from non-MPR for neoadjuvant immunochemotherapy, which could support the clinical decisions on the utilization of immune checkpoint inhibitor-based treatments in NSCLC patients.

Identification of novel immune-related signatures for keloid diagnosis and treatment: insights from integrated bulk RNA-seq and scRNA-seq analysis

BackgroundKeloid is a disease characterized by proliferation of fibrous tissue after the healing of skin tissue, which seriously affects the daily life of patients. However, the clinical treatment of keloids still has limitations, that is, it is not effective in controlling keloids, resulting in a high recurrence rate. Thus, it is urgent to identify new signatures to improve the diagnosis and treatment of keloids.MethodBulk RNA seq and scRNA seq data were downloaded from the GEO database. First, we used WGCNA and MEGENA to co-identify keloid/immune-related DEGs. Subsequently, we used three machine learning algorithms (Randomforest, SVM-RFE, and LASSO) to identify hub immune-related genes of keloid (KHIGs) and investigated the heterogeneous expression of KHIGs during fibroblast subpopulation differentiation using scRNA-seq. Finally, we used HE and Masson staining, quantitative reverse transcription-PCR, western blotting, immunohistochemical, and Immunofluorescent assay to investigate the dysregulated expression and the mechanism of retinoic acid in keloids.ResultsIn the present study, we identified PTGFR, RBP5, and LIF as KHIGs and validated their diagnostic performance. Subsequently, we constructed a novel artificial neural network molecular diagnostic model based on the transcriptome pattern of KHIGs, which is expected to break through the current dilemma faced by molecular diagnosis of keloids in the clinic. Meanwhile, the constructed IG score can also effectively predict keloid risk, which provides a new strategy for keloid prevention. Additionally, we observed that KHIGs were also heterogeneously expressed in the constructed differentiation trajectories of fibroblast subtypes, which may affect the differentiation of fibroblast subtypes and thus lead to dysregulation of the immune microenvironment in keloids. Finally, we found that retinoic acid may treat or alleviate keloids by inhibiting RBP5 to differentiate pro-inflammatory fibroblasts (PIF) to mesenchymal fibroblasts (MF), which further reduces collagen secretion.ConclusionIn summary, the present study provides novel immune signatures (PTGFR, RBP5, and LIF) for keloid diagnosis and treatment, and identifies retinoic acid as potential anti-keloid drugs. More importantly, we provide a new perspective for understanding the interactions between different fibroblast subtypes in keloids and the remodeling of their immune microenvironment.