Direct prediction of antimicrobial resistance in Pseudomonas aeruginosa by metagenomic next-generation sequencing.

Objective: Pseudomonas aeruginosa has strong drug resistance and can tolerate a variety of antibiotics, which is a major problem in the management of antibiotic-resistant infections. Direct prediction of multi-drug resistance (MDR) resistance phenotypes of P. aeruginosa isolates and clinical samples by genotype is helpful for timely antibiotic treatment. Methods: In the study, whole genome sequencing (WGS) data of 494 P. aeruginosa isolates were used to screen key anti-microbial resistance (AMR)-associated genes related to imipenem (IPM), meropenem (MEM), piperacillin/tazobactam (TZP), and levofloxacin (LVFX) resistance in P. aeruginosa by comparing genes with copy number differences between resistance and sensitive strains. Subsequently, for the direct prediction of the resistance of P. aeruginosa to four antibiotics by the AMR-associated features screened, we collected 74 P. aeruginosa positive sputum samples to sequence by metagenomics next-generation sequencing (mNGS), of which 1 sample with low quality was eliminated. Then, we constructed the resistance prediction model. Results: We identified 93, 88, 80, 140 AMR-associated features for IPM, MEM, TZP, and LVFX resistance in P. aeruginosa. The relative abundance of AMR-associated genes was obtained by matching mNGS and WGS data. The top 20 features with importance degree for IPM, MEM, TZP, and LVFX resistance were used to model, respectively. Then, we used the random forest algorithm to construct resistance prediction models of P. aeruginosa, in which the areas under the curves of the IPM, MEM, TZP, and LVFX resistance prediction models were all greater than 0.8, suggesting these resistance prediction models had good performance. Conclusion: In summary, mNGS can predict the resistance of P. aeruginosa by directly detecting AMR-associated genes, which provides a reference for rapid clinical detection of drug resistance of pathogenic bacteria.

Identifying important microbial biomarkers for the diagnosis of colon cancer using a random forest approach.

Colon cancer is one of the most common cancers, with 30-50 % of patients returning or metastasizing within 5 years of treatment. Increasingly, researchers have highlighted the influence of microbes on cancer malignant activity, while no studies have explored the relationship between colon cancer and the microbes in tumors. Here, we used tissue and blood samples from 67 colon cancer patients to identify pathogenic microorganisms associated with the diagnosis and prediction of colon cancer and evaluate the predictive performance of each pathogenic marker and its combination based on the next-generation sequencing data by using random forest algorithms. The results showed that we constructed a database of 13,187 pathogenic microorganisms associated with human disease and identified 2 pathogenic microorganisms (Synthetic.construct_32630 and Dicrocoelium.dendriticum_57078) associated with colon cancer diagnosis, and the constructed diagnostic prediction model performed well for tumor tissue samples and blood samples. In summary, for the first time, we provide new molecular markers for the diagnosis of colon cancer based on the expression of pathogenic microorganisms in order to provide a reference for improving the effective screening rate of colon cancer in clinical practice and ameliorating the personalized treatment of colon cancer patients.

Comparative analyses of the prognosis, tumor immune microenvironment, and drug treatment response between left-sided and right-sided colon cancer by integrating scRNA-seq and bulk RNA-seq data.

BACKGROUND: In this study, we compared the prognosis, tumor immune microenvironment (TIM), and drug treatment response between left-sided (LCC) and right-sided (RCC) colon cancer to predict outcomes in patients with LCC and RCC. METHODS: Based on identified differentially expressed genes and using single-cell RNA sequencing data, we constructed and validated a prognostic model for LCC and RCC patients in the TCGA-COAD cohort and GSE103479 cohort. Moreover, we compared the differences of TIM characteristics and drug treatment response between LCC and RCC patients. RESULTS: We constructed and validated a five-gene prognostic model for LCC patients and a four-gene prognostic model for RCC patients, and both showed excellent performance. The RCC patients with higher risk scores were significantly associated with greater metastasis (P = 2.6×10-5), N stage (P = 0.012), advanced pathological stage (P = 1.4×10-4), and more stable microsatellite status (P = 0.007) but not T stage (P = 0.200). For LCC patients, the risk scores were not significantly associated with tumor stage and microsatellite status (P > 0.05). Additionally, immune infiltration by CD8 and regulatory T cells and M0, M1, and M2 macrophages differed significantly between LCC and RCC patients (P < 0.05). APC and TP53 mutations were significantly more common in LCC patients (P < 0.05). In contrast, KRAS, SYNE1, and MUC16 mutations were significantly more common in RCC patients (P < 0.05). In addition, tumor mutation burden values were significantly higher in RCC patients than in LCC patients (P = 5.9×10-8). Moreover, the expression of immune checkpoint targets was significantly higher in RCC patients than in LCC patients (P < 0.05), indicating that RCC patients maybe more sensitive to immunotherapy. However, LCC and RCC patients did not differ significantly in their sensitivity to eight selected chemicals or target drugs (P > 0.05). The average half-maximal inhibitory concentrations for camptothecin, teniposide, vinorelbine, and mitoxantrone were significantly lower in low-risk than in high-risk RCC patients (P < 0.05), indicating that the lower risk score of RCC patients, the more sensitive they were to these four drugs. CONCLUSIONS: We investigated the differences in prognosis, TIM, and drug treatment response between LCC and RCC patients, which may contribute to accurate colon cancer prognosis and treatment of colon cancer.

Identification of m6A-related lncRNAs as prognostic signature within colon tumor immune microenvironment.

BACKGROUND: The current study probed prognosis-related potential for m6A-related lncRNAs signatures within colon tumor immune microenvironment (TIM). METHODS: After downloading transcriptomic datasets for colon cancer (CC) patients from The Cancer Genome Atlas (TCGA), they were divided, in a 1:1 ratio, within training or test datasets. m6A-related lncRNAs were then scrutinized across such dataset using Pearson correlation assessment before generating a m6A-related lncRNAs prognosis-related model using the training dataset. The latter was then validated with the test and the whole dataset. In addition, we compared the differences of TIM and the estimated IC50 of drug response between the high- and low-risk groups. RESULTS: Overall survival (OS) resulted as linked with 11 m6A-related lncRNAs, while within the developed prognosis-related model, areas-under-curves were as follows: within training dataset, values at 3-, 4-, and 5-years were 0.777, 0.819, and 0.805, accordingly, and for test one, they were 0.697, 0.682, and 0.706, respectively. Finally, the values for the whole dataset were 0.675 (3-year), 0.682 (4-years), and 0.679 (5-years), accordingly. Moreover, CC cases categorized within low-risk cohort demonstrated enhanced OS (p < .0001), lower metastasis (p = 2e-06) and lower T stage (p = .0067), more instability for microsatellite status (p = .012), and downregulation for PD-L1, PD-1, CTLA-4, LAG3, and HAVCR2 (p < .05). In addition, risk scorings were significantly linked to the degree of infiltrative intensity for CD8 and CD4 (memory resting) T-cells, T-regulatory (Tregs), and Mast cells triggering (p < .05). Patients with low infiltrative propensity for CD4 T-cells also had better OS (p = .016). Moreover, six representative drugs were found to be sensitive for treating CC patients. CONCLUSION: A robust m6A-related prognostic model with great performances was developed before exploring the TIM characteristics and its potential therapeutic drugs, which might improve the prognosis and therapeutic efficacy.

Unusual Talaromyces marneffei and Pneumocystis jirovecii coinfection in a child with a STAT1 mutation: A case report and literature review.

Talaromyces marneffei and Pneumocystis jirovecii are the common opportunistic pathogens in immunodeficient patients. There have been no reports of T. marneffei and P. jirovecii coinfection in immunodeficient children. Signal transducer and activator of transcription 1 (STAT1) is a key transcription factor in immune responses. STAT1 mutations are predominately associated with chronic mucocutaneous candidiasis and invasive mycosis. We report a 1-year-2-month-old boy diagnosed with severe laryngitis and pneumonia caused by T. marneffei and P. jirovecii coinfection, which was confirmed by smear, culture, polymerase chain reaction and metagenome next-generation sequencing of bronchoalveolar lavage fluid. He has a known STAT1 mutation at amino acid 274 in the coiled-coil domain of STAT1 according to whole exome sequencing. Based on the pathogen results, itraconazole and trimethoprim-sulfamethoxazole were administered. This patient's condition improved, and he was discharged after two weeks of targeted therapy. In the one-year follow-up, the boy remained symptom-free without recurrence.

Ferroptosis-related genes are considered as potential targets for CPAP treatment of obstructive sleep apnea.

Obstructive sleep apnea (OSA) is a common syndrome characterized by upper airway dysfunction during sleep. Continuous positive airway pressure (CPAP) is the most frequently utilized non-surgical treatment for OSA. Ferroptosis play a crucial role in the physiological diseases caused by chronic intermittent hypoxia, but its involvement in the development of OSA and the exact mechanisms have incompletely elucidated. GSE75097 microarray dataset was used to identify differentially expressed genes between OSA patients and CPAP-treated OSA patients. Subsequently, Gene Ontology (GO) annotation, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway, STRING database, and FerrDb database were conducted to analyze the biological functions of differentially expressed genes and screen ferroptosis-related genes. Finally, GSE135917 dataset employed for validation. There were 1,540 differentially expressed genes between OSA patients and CPAP-treated OSA patients. These differentially expressed genes were significantly enriched in the regulation of interleukin-1-mediated signaling pathway and ferroptosis-related signaling pathway. Subsequently, 13 ferroptosis-related genes (DRD5, TSC22D3, TFAP2A, STMN1, DDIT3, MYCN, ELAVL1, JUN, DUSP1, MIB1, PSAT1, LCE2C, and MIR27A) were identified from the interaction between differentially expressed genes and FerrDb database, which are regarded as the potential targets of CPAP-treated OSA. These ferroptosis-related genes were mainly involved in cell proliferation and apoptosis and MAPK signaling pathway. Furthermore, DRD5 and TFAP2A were downregulated in OSA patients, which showed good diagnostic properties for OSA, but these abnormal signatures are not reversed with short-term effective CPAP therapy. In summary, the identification of 13 ferroptosis-related genes as potential targets for the CPAP treatment of OSA provides valuable insights into the development of novel, reliable, and accurate therapeutic options.

Exploring immune-related signatures for predicting immunotherapeutic responsiveness, prognosis, and diagnosis of patients with colon cancer.

The present study focused on identifying the immune-related signatures and exploring their performance in predicting the prognosis, immunotherapeutic responsiveness, and diagnosis of patients with colon cancer. Firstly, the immunotherapeutic response-related differential expressed genes (DEGs) were identified by comparing responders and non-responders from an anti-PD-L1 cohort using the edgeR R package. Then, the immunotherapeutic response related DEGs was intersected with immune-related genes (IRGs) to obtain the immunotherapeutic response and immune-related genes (IRIGs). Then, an immunotherapeutic response and immune-related risk score (IRIRScore) model consisting of 6 IRIGs was constructed using the univariable Cox regression analysis and multivariate Cox regression analysis based on the COAD cohort from the cancer genome atlas (TCGA) database, which was further validated in two independent gene expression omnibus database (GEO) datasets (GSE39582 and GSE17536) and anti-PD-L1 cohort. A nomogram with good accuracy was established based on the immune-related signatures and clinical factors (C-index = 0.75). In the training dataset and GSE39582, higher IRIRScore was significantly associated with higher TMN and advanced pathological stages. Based on the anti-PD-L1 cohort, patients who were sensitive to immunotherapy had significantly lower risk score than non-responders. Furthermore, we explored the immunotherapy-related signatures based on the training dataset. Kaplan-Meier curve revealed a high level of T cells regulatory (Tregs) was significantly related to poor overall survival (OS), while a high level of T cells CD4 memory resting was significantly related to better OS. Besides, the TMB value of patients in the high-risk group was significantly higher than those in a low-risk group. Moreover, patients in the high-risk group had significantly higher expression levels of immune checkpoint inhibitors. In addition, the immune-related signatures were applied to establish prediction models using the random forest algorithm. Among them, TDGF1 and NRG1 revealed excellent diagnostic predictive performance (AUC >0.8). In conclusion, the current findings provide new insights into immune-related immunotherapeutic responsiveness, prognosis, and diagnosis of colon cancer.

Exploring Immune-Related Prognostic Signatures in the Tumor Microenvironment of Colon Cancer.

Background: Colon cancer is a common malignant tumor with poor prognosis. The aim of this study is to explore the immune-related prognostic signatures and the tumor immune microenvironment of colon cancer. Methods: The mRNA expression data of TCGA-COAD from the UCSC Xena platform and the list of immune-related genes (IRGs) from the ImmPort database were used to identify immune-related differentially expressed genes (DEGs). Then, we constructed an immune-related risk score prognostic model and validated its predictive performance in the test dataset, the whole dataset, and two independent GEO datasets. In addition, we explored the differences in tumor-infiltrating immune cell types, tumor mutation burden (TMB), microsatellite status, and expression levels of immune checkpoints and their ligands between the high-risk and low-risk score groups. Moreover, the potential value of the identified immune-related signature with respect to immunotherapy was investigated based on an immunotherapeutic cohort (Imvigor210) treated with an anti-PD-L1 agent. Results: Seven immune-related DEGs were identified as prognostic signatures. The areas under the curves (AUCs) of the constructed risk score model for overall survival (OS) were calculated (training dataset: 0.780 at 3 years, 0.801 at 4 years, and 0.766 at 5 years; test dataset: 0.642 at 3 years, 0.647 at 4 years, and 0.629 at 5 years; and the whole dataset: 0.642 at 3 years, 0.647 at 4 years, and 0.629 at 5 years). In the high-risk score group of the whole dataset, patients had worse OS, higher TMN stages, advanced pathological stages, and a higher TP53 mutation rate (p < 0.05). In addition, a high level of resting NK cells or M0 macrophages, and high TMB were significantly related to poor OS (p < 0.05). Also, we observed that high-risk score patients had a high expression level of PD-L1, PD-1, and CTLA-4 (p < 0.05). The patients with high-risk scores demonstrated worse prognosis than those with low-risk scores in multiple datasets (GSE39582: p = 0.0023; GSE17536: p = 0.0008; immunotherapeutic cohort without platinum treatment: p = 0.0014; immunotherapeutic cohort with platinum treatment: p = 0.0027). Conclusion: We developed a robust immune-related prognostic signature that performed great in multiple cohorts and explored the characteristics of the tumor immune microenvironment of colon cancer patients, which may give suggestions for the prognosis and immunotherapy in the future.

Construction of a novel methylation-related prognostic model for colorectal cancer based on microsatellite status.

The present study aimed to construct a novel methylation-related prognostic model based on microsatellite status that may enhance the prognosis of colorectal cancer (CRC) from methylation and microsatellite status perspective. DNA methylation and mRNA expression data with clinical information were downloaded from The Cancer Genome Atlas (TCGA) data set. The samples were divided into microsatellite stability and microsatellite instability group, and CIBERSORT was used to assess the immune cell infiltration characteristics. After identifying the differentially methylated genes and differentially expression genes using R packages, the methylation-driven genes were further identified. Prognostic genes that were used to establish the methylation-related risk score model were generated by the univariate and multivariate Cox regression model. Finally, we established and evaluated the methylation-related prognostic model for CRC patients. A total of 69 MDGs were obtained and three of these genes (MIOX, TH, DKFZP434K028) were selected to construct the prognostic model. Patients in the low-risk score group had a conspicuously better overall survival than those in the high-risk score group (p < .0001). The area under the receiver operating characteristic curve for this model was 0.689 at 3 years, 0.674 at 4 years, and 0.658 at 5 years. The Wilcoxon test showed that higher risk score was associated with higher T stage (p = .01), N stages (p = .0028), metastasis (p = .013), and advanced pathological stage (p = .0013). However, the more instability of microsatellite status, the lower risk score of CRC patients (p = .0048). Our constructed methylation-related prognostic model based on microsatellite status presents potential significance in assessing recurrence risk stratification, tumor staging, and immunotherapy for CRC patients.