N6-methyladenosine-related non-coding RNAs are potential prognostic and immunotherapeutic responsiveness biomarkers for bladder cancer

RNA N6-methyladenosine modifications in urological cancers: from mechanism to application

The N6-methyladenosine (m6A) modification is the most common modification of messenger RNAs in eukaryotes and has crucial roles in multiple cancers, including in urological malignancies such as renal cell carcinoma, bladder cancer and prostate cancer. The m6A RNA modification is controlled by three types of regulators, including methyltransferases (writers), demethylases (erasers) and RNA-binding proteins (readers), which are responsible for gene regulation at the post-transcriptional level. This Review summarizes the current evidence indicating that aberrant or dysregulated m6A modification is associated with urological cancer development, progression and prognosis. The complex and context-dependent effects of dysregulated m6A modifications in urological cancers are described, along with the potential for aberrantly expressed m6A regulators to provide valuable diagnostic and prognostic biomarkers as well as new therapeutic targets.

Multi-omic profiling reveals potential biomarkers of hepatocellular carcinoma prognosis and therapy response among mitochondria-associated cell death genes in the context of 3P medicine

Background

Cancer cell growth, metastasis, and drug resistance are major challenges in treating liver hepatocellular carcinoma (LIHC). However, the lack of comprehensive and reliable models hamper the effectiveness of the predictive, preventive, and personalized medicine (PPPM/3PM) strategy in managing LIHC.

Methods

Leveraging seven distinct patterns of mitochondrial cell death (MCD), we conducted a multi-omic screening of MCD-related genes. A novel machine learning framework was developed, integrating 10 machine learning algorithms with 67 different combinations to establish a consensus mitochondrial cell death index (MCDI). This index underwent rigorous evaluation across training, validation, and in-house clinical cohorts. A comprehensive multi-omics analysis encompassing bulk, single-cell, and spatial transcriptomics was employed to achieve a deeper insight into the constructed signature. The response of risk subgroups to immunotherapy and targeted therapy was evaluated and validated. RT-qPCR, western blotting, and immunohistochemical staining were utilized for findings validation.

Results

Nine critical differentially expressed MCD-related genes were identified in LIHC. A consensus MCDI was constructed based on a 67-combination machine learning computational framework, demonstrating outstanding performance in predicting prognosis and clinical translation. MCDI correlated with immune infiltration, Tumor Immune Dysfunction and Exclusion (TIDE) score and sorafenib sensitivity. Findings were validated experimentally. Moreover, we identified PAK1IP1 as the most important gene for predicting LIHC prognosis and validated its potential as an indicator of prognosis and sorafenib response in our in-house clinical cohorts.

Conclusion

This study developed a novel predictive model for LIHC, namely MCDI. Incorporating MCDI into the PPPM framework will enhance clinical decision-making processes and optimize individualized treatment strategies for LIHC patients.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1007/s13167-024-00362-8.

Enhancing m6A modification of lncRNA through METTL3 and RBM15 to promote malignant progression in bladder cancer

Objective

Bladder cancer is one of the most prominent malignancies affecting the urinary tract, characterized by a poor prognosis. Our previous research has underscored the pivotal role of m6A methylation in the progression of bladder cancer. Nevertheless, the precise relationship between N6-methyladenosine (m6A) regulation of long non-coding RNA (lncRNA) and bladder cancer remains elusive.

Methods

This study harnessed sequencing data and clinical records from 408 bladder cancer patients in the TCGA database. Employing R software, we conducted bioinformatics analysis to establish an m6A-lncRNA co-expression network. Analyzing the differences between high and low-risk groups, particularly at the immunological level, and subsequently investigating the primary regulatory factors of these lncRNA, validating the findings through experiments, and exploring their specific cellular functions.

Results

We identified 50 m6A-related lncRNA with prognostic significance through univariate Cox regression analysis. In parallel, we employed a LASSO-Cox regression model to pinpoint 11 lncRNA and calculate risk scores for bladder cancer patients. Based on the median risk score, patients were categorized into low-risk and high-risk groups. The high-risk cohort exhibited notably lower survival rates than their low-risk counterparts. Further analysis pointed to RBM15 and METTL3 as potential master regulators of these m6A-lncRNA. Experimental findings also shed light on the upregulated expression of METTlL3 and RBM15 in bladder cancer, where they contributed to the malignant progression of tumors. The experimental findings demonstrated a significant upregulation of METTL3 and RBM15 in bladder cancer specimens, implicating their contributory role in the oncogenic progression. Knockdown of METTL3 and RBM15 resulted in a marked attenuation of tumor cell proliferation, invasion, and migration, which was concomitant with a downregulation in the cellular m6A methylation status. Moreover, these results revealed that RBM15 and METTL3 function in a synergistic capacity, positing their involvement in cancer promotion via the upregulation of m6A modifications in long non-coding RNAs. Additionally, this study successfully developed an N-methyl-N-nitrosourea (MNU)-induced rat model of in situ bladder carcinoma, confirming the elevated expression of RBM15 and METTL3, which paralleled the overexpression of m6A-related- lncRNAs observed in bladder cancer cell lines. This congruence underscores the potential utility of these molecular markers in in vivo models that mirror human malignancies.

Conclusion

This study not only offers novel molecular targets,but also enriches the research on m6A modification in bladder cancer, thereby facilitating its clinical translation.

Comprehensive analysis of the interaction of antigen presentation during anti-tumour immunity and establishment of AIDPS systems in ovarian cancer

There are hundreds of prognostic models for ovarian cancer. These genes are based on different gene classes, and there are many ways to construct the models. Therefore, this paper aims to build the most stable prognostic evaluation system known to date through 101 machine learning strategies. We combined 101 algorithm combinations with 10 machine learning algorithms to create antigen presentation-associated genetic markers (AIDPS) with outstanding precision and steady performance. The inclusive set of algorithms comprises the elastic network (Enet), Ridge, stepwise Cox, Lasso, generalized enhanced regression model (GBM), random survival forest (RSF), supervised principal component (SuperPC), Cox partial least squares regression (plsRcox), survival support vector machine (Survival-SVM). Then, in the train cohort, the prediction model was fitted using a leave-one cross-validation (LOOCV) technique, which involved 101 different possible combinations of prognostic genes. Seven validation data sets (GSE26193, GSE26712, GSE30161, GSE63885, GSE9891, GSE140082 and ICGC_OV_AU) were compared and analysed, and the C-index was calculated. Finally, we collected 32 published ovarian cancer prognostic models (including mRNA and lncRNA). All data sets and prognostic models were subjected to a univariate Cox regression analysis, and the C-index was calculated to demonstrate that the antigen presentation process should be the core criterion for evaluating ovarian cancer prognosis. In a univariate Cox regression analysis, 22 prognostic genes were identified based on the expression profiles of 283 genes involved in antigen presentation and the intersection of genes (p < 0.05). AIDPS were developed by our machine learning-based integration method, which was applied to these 22 genes. One hundred and one prediction models are fitted using the LOOCV framework, and the C-index is calculated for each model across all validation sets. Interestingly, RSF + Lasso was the best model overall since it had the greatest average C-index and the highest C-index of any combination of models tested on the validated data sets. In comparing external cohorts, we found that the C-index correlated AIDPS method using the RSF + Lasso method in 101 prediction models was in contrast to other features. Notably, AIDPS outperformed the vast majority of models across all data sets. Antigen-presenting anti-tumour immune pathways can be used as a representative gene set of ovarian cancer to track the prognosis of patients with cancer. The antigen-presenting model obtained by the RSF + Lasso method has the best C-INDEX, which plays a key role in developing antigen-presenting targeted drugs in ovarian cancer and improving the treatment outcome of patients.

PD-L1-related IncRNAs are associated with malignant characteristics and immune microenvironment in glioma

BACKGROUND

The expression of long non-coding RNA (lncRNA) can function as diagnostic and therapeutic biomarker for tumors. This research explores the role of PD-L1-related lncRNAs in affecting malignant characteristics and the immune microenvironment of glioma.

METHODS

Downloading gene expression profiles and clinicopathological information of glioma from TCGA and CGGA databases, 6 PD-L1-related lncRNAs were identified through correlation analysis, Cox and LASSO regression analysis, establishing the risk score model based on them. Bioinformatics analysis and cell experiments in vitro were adopted to verify the effects of LINC01271 on glioma.

RESULTS

Risk scores based on 6 PD-L1-related lncRNAs (AL355974.3, LINC01271, AC011899.3, MIR4500HG, LINC02594, AL357055.3) can reflect malignant characteristics and immunotherapy response of glioma. Patients with high LINC01271 expression had a worse prognosis, a higher abundance of M1 subtype macrophages in the immune microenvironment, and a higher degree of tumor malignancy. Experiments in vitro confirmed its positive regulatory effect on the proliferation and migration of glioma cells.

CONCLUSIONS

The risk score model based on 6 PD-L1-related lncRNAs can reflect the malignant characteristics and prognosis of glioma. LINC01271 can independently be used as a new target for prognosis evaluation and therapy.

Comprehensive analyses of fatty acid metabolism-related lncRNA for ovarian cancer patients

Ovarian cancer (OC) is a disease with difficult early diagnosis and treatment and poor prognosis. OC data profiles were downloaded from The Cancer Genome Atlas. Eight key fatty acid metabolism-related long non-coding RNAs (lncRNAs) were finally screened for building a risk scoring model by univariate/ multifactor and least absolute shrinkage and selection operator (LASSO) Cox regression. To make this risk scoring model more applicable to clinical work, we established a nomogram containing the clinical characteristics of OC patients after confirming that the model has good reliability and validity and the ability to distinguish patient prognosis. To further explore how these key lncRNAs are involved in OC progression, we explored their relationship with LUAD immune signatures and tumor drug resistance. The structure shows that the risk scoring model established based on these 8 fatty acid metabolism-related lncRNAs has good reliability and validity and can better predict the prognosis of patients with different risks of OC, and LINC00861in these key RNAs may be a hub gene that affects the progression of OC and closely related to the sensitivity of current OC chemotherapy drugs. In addition, combined with immune signature analysis, we found that patients in the high-risk group are in a state of immunosuppression, and Tfh cells may play an important role in it. We innovatively established a prognostic prediction model with excellent reliability and validity from the perspective of OC fatty acid metabolism reprogramming and lncRNA regulation and found new molecular/cellular targets for future OC treatment.

Underexplored reciprocity between genome-wide methylation status and long non-coding RNA expression reflected in breast cancer research: potential impacts for the disease management in the framework of 3P medicine

Breast cancer (BC) is the most common female malignancy reaching a pandemic scale worldwide. A comprehensive interplay between genetic alterations and shifted epigenetic regions synergistically leads to disease development and progression into metastatic BC. DNA and histones methylations, as the most studied epigenetic modifications, represent frequent and early events in the process of carcinogenesis. To this end, long non-coding RNAs (lncRNAs) are recognized as potent epigenetic modulators in pathomechanisms of BC by contributing to the regulation of DNA, RNA, and histones' methylation. In turn, the methylation status of DNA, RNA, and histones can affect the level of lncRNAs expression demonstrating the reciprocity of mechanisms involved. Furthermore, lncRNAs might undergo methylation in response to actual medical conditions such as tumor development and treated malignancies. The reciprocity between genome-wide methylation status and long non-coding RNA expression levels in BC remains largely unexplored. Since the bio/medical research in the area is, per evidence, strongly fragmented, the relevance of this reciprocity for BC development and progression has not yet been systematically analyzed. Contextually, the article aims at:consolidating the accumulated knowledge on both-the genome-wide methylation status and corresponding lncRNA expression patterns in BC andhighlighting the potential benefits of this consolidated multi-professional approach for advanced BC management. Based on a big data analysis and machine learning for individualized data interpretation, the proposed approach demonstrates a great potential to promote predictive diagnostics and targeted prevention in the cost-effective primary healthcare (sub-optimal health conditions and protection against the health-to-disease transition) as well as advanced treatment algorithms tailored to the individualized patient profiles in secondary BC care (effective protection against metastatic disease). Clinically relevant examples are provided, including mitochondrial health control and epigenetic regulatory mechanisms involved.

Multi-omics identification of an immunogenic cell death-related signature for clear cell renal cell carcinoma in the context of 3P medicine and based on a 101-combination machine learning computational framework

Background

Clear cell renal cell carcinoma (ccRCC) is a prevalent urological malignancy associated with a high mortality rate. The lack of a reliable prognostic biomarker undermines the efficacy of its predictive, preventive, and personalized medicine (PPPM/3PM) approach. Immunogenic cell death (ICD) is a specific type of programmed cell death that is tightly associated with anti-cancer immunity. However, the role of ICD in ccRCC remains unclear.

Methods

Based on AddModuleScore, single-sample gene set enrichment analysis (ssGSEA), and weighted gene co-expression network (WGCNA) analyses, ICD-related genes were screened at both the single-cell and bulk transcriptome levels. We developed a novel machine learning framework that incorporated 10 machine learning algorithms and their 101 combinations to construct a consensus immunogenic cell death-related signature (ICDRS). ICDRS was evaluated in the training, internal validation, and external validation sets. An ICDRS-integrated nomogram was constructed to provide a quantitative tool for predicting prognosis in clinical practice. Multi-omics analysis was performed, including genome, single-cell transcriptome, and bulk transcriptome, to gain a more comprehensive understanding of the prognosis signature. We evaluated the response of risk subgroups to immunotherapy and screened drugs that target specific risk subgroups for personalized medicine. Finally, the expression of ICD-related genes was validated by qRT-PCR.

Results

We identified 131 ICD-related genes at both the single-cell and bulk transcriptome levels, of which 39 were associated with overall survival (OS). A consensus ICDRS was constructed based on a 101-combination machine learning computational framework, demonstrating outstanding performance in predicting prognosis and clinical translation. ICDRS can also be used to predict the occurrence, development, and metastasis of ccRCC. Multivariate analysis verified it as an independent prognostic factor for OS, progression-free survival (PFS), and disease-specific survival (DSS) of ccRCC. The ICDRS-integrated nomogram provided a quantitative tool in clinical practice. Moreover, we observed distinct biological functions, mutation landscapes, and immune cell infiltration in the tumor microenvironment between the high- and low-risk groups. Notably, the immunophenoscore (IPS) score showed a significant difference between risk subgroups, suggesting a better response to immunotherapy in the high-risk group. Potential drugs targeting specific risk subgroups were also identified.

Conclusion

Our study constructed an immunogenic cell death-related signature that can serve as a promising tool for prognosis prediction, targeted prevention, and personalized medicine in ccRCC. Incorporating ICD into the PPPM framework will provide a unique opportunity for clinical intelligence and new management approaches.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13167-023-00327-3.

Establishment of a novel lysosomal signature for the diagnosis of gastric cancer with in-vitro and in-situ validation

Background

Gastric cancer (GC) represents a malignancy with a multi-factorial combination of genetic, environmental, and microbial factors. Targeting lysosomes presents significant potential in the treatment of numerous diseases, while lysosome-related genetic markers for early GC detection have not yet been established, despite implementing this process by assembling artificial intelligence algorithms would greatly break through its value in translational medicine, particularly for immunotherapy.

Methods

To this end, this study, by utilizing the transcriptomic as well as single cell data and integrating 20 mainstream machine-learning (ML) algorithms. We optimized an AI-based predictor for GC diagnosis. Then, the reliability of the model was initially confirmed by the results of enrichment analyses currently in use. And the immunological implications of the genes comprising the predictor was explored and response of GC patients were evaluated to immunotherapy and chemotherapy. Further, we performed systematic laboratory work to evaluate the build-up of the central genes, both at the expression stage and at the functional aspect, by which we could also demonstrate the reliability of the model to guide cancer immunotherapy.

Results

Eight lysosomal-related genes were selected for predictive model construction based on the inclusion of RMSE as a reference standard and RF algorithm for ranking, namely ADRB2, KCNE2, MYO7A, IFI30, LAMP3, TPP1, HPS4, and NEU4. Taking into account accuracy, precision, recall, and F1 measurements, a preliminary determination of our study was carried out by means of applying the extra tree and random forest algorithms, incorporating the ROC-AUC value as a consideration, the Extra Tree model seems to be the optimal option with the AUC value of 0.92. The superiority of diagnostic signature is also reflected in the analysis of immune features.

Conclusion

In summary, this study is the first to integrate around 20 mainstream ML algorithms to construct an AI-based diagnostic predictor for gastric cancer based on lysosomal-related genes. This model will facilitate the accurate prediction of early gastric cancer incidence and the subsequent risk assessment or precise individualized immunotherapy, thus improving the survival prognosis of GC patients.

Global burden of cancers attributable to tobacco smoking, 1990-2019: an ecological study

Aim and background

Identifying risk factors for cancer initiation and progression is the cornerstone of the preventive approach to cancer management and control (EPMA J. 4(1):6, 2013). Tobacco smoking is a well-recognized risk factor for initiation and spread of several cancers. The predictive, preventive, and personalized medicine (PPPM) approach to cancer management and control focuses on smoking cessation as an essential cancer prevention strategy. Towards this end, this study examines the temporal patterns of cancer burden due to tobacco smoking in the last three decades at global, regional, and national levels.

Data and methods

The data pertaining to the burden of 16 cancers attributable to tobacco smoking at global, regional, and national levels were procured from the Global Burden of Disease 2019 Study. Two main indicators, deaths and disability-adjusted life years (DALYs), were used to describe the burden of cancers attributable to tobacco smoking. The socio-economic development of countries was measured using the socio-demographic index (SDI).

Results

Globally, deaths due to neoplasms caused by tobacco smoking increased from 1.5 million in 1990 to 2.5 million in 2019, whereas the age-standardized mortality rate (ASMR) decreased from 39.8/100,000 to 30.6/100,000 and the age-standardized DALY rate (ASDALR) decreased from 948.9/100,000 to 677.3/100,000 between 1990 and 2019. Males accounted for approximately 80% of global deaths and DALYs in 2019. Populous regions of Asia and a few regions of Europe account for the largest absolute burden, whereas countries in Europe and America have the highest age-standardized rates of cancers due to tobacco smoking. In 8 out of 21 regions, there were more than 100,000 deaths due to cancers attributable to tobacco smoking led by East Asia, followed by Western Europe in 2019. The regions of Sub-Saharan Africa (except southern region) had one of the lowest absolute counts of deaths, DALYs, and age-standardized rates. In 2019, tracheal, bronchus, and lung (TBL), esophageal, stomach, colorectal, and pancreatic cancer were the top 5 neoplasms attributable to tobacco smoking, with different burdens in regions as per their development status. The ASMR and ASDALR of neoplasms due to tobacco smoking were positively correlated with SDI, with pairwise correlation coefficient of 0.55 and 0.52, respectively.

Conclusion

As a preventive tool, tobacco smoking cessation has the biggest potential among all risk factors for preventing millions of cancer deaths every year. Cancer burden due to tobacco smoking is found to be higher in males and is positively associated with socio-economic development of countries. As tobacco smoking begins mostly at younger ages and the epidemic is unfolding in several parts of the world, more accelerated efforts are required towards tobacco cessation and preventing youth from entering this addiction. The PPPM approach to medicine suggests that not only personalized and precision medicine must be provided to cancer patients afflicted by tobacco smoking but personalized and targeted preventive solutions must be provided to prevent initiation and progression of smoking.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13167-022-00308-y.

PXDNL activates the motility of urothelial bladder carcinoma cells through the Wnt/β-catenin pathway and has a prognostic value

AIMS

Although aberrant expression of peroxidasin-like (PXDNL) has been associated with carcinogenesis, its potential role in the Urothelial Carcinoma of the Bladder (UCB) remains unknown. The present study aimed to explore the role of PXDNL in UCB carcinogenesis and its potential clinical value.

MAIN METHODS

Based on The Cancer Genome Atlas (TCGA) data, bioinformatics was used to explore the potential clinical value of PXDNL. Wound healing and Transwell invasion assays were employed for the purpose of assessing the cell motility, while the Western Blotting experiments were utilized for investigating the protein expression pattern of PXDNL in UCB and investigating the Epithelial-to-Mesenchymal Transition (EMT) and Wnt/β-catenin pathways for understanding the probable mechanisms involved.

KEY FINDS

PXDNL mRNA was overexpressed in UCB tissues and indicated a poor prognosis. High PXDNL mRNA levels were also associated with advanced clinicopathological features and were regarded as independent prognostic factors for UCB. However, PXDNL showed a weak correlation with immune cell infiltration in UCB. In addition, the findings of the study verified that the existing form of the PXDNL protein was 57-kDa and it was upregulated in the UCB cell lines and tissue samples. Furthermore, silencing PXDNL inhibited, while overexpressing PXDNL promoted EMT and motility of UCB cells in vitro. Mechanistic studies showed that PXDNL activated UCB cell motility via the Wnt/β-catenin pathway.

SIGNIFICANCE

The results reveal a novel molecular target that could be further explored for developing preventive, predictive, and individualized treatment strategies for UCB.

Identification of cuproptosis-related long noncoding RNA signature for predicting prognosis and immunotherapy response in bladder cancer

Bladder cancer (BC) is the most common malignant tumour of the urinary system and one of the leading causes of cancer-related death. Cuproptosis is a novel form of programmed cell death, and its mechanism in tumours remains unclear. This study aimed to establish the prognostic signatures of cuproptosis-related lncRNAs and determine their clinical prognostic value. RNA sequencing data from The Cancer Genome Atlas were used to detect the expression levels of cuproptosis-related genes in BC. Cuproptosis-related lncRNAs linked to survival were identified using co-expression and univariate Cox regression. Furthermore, consensus cluster analysis divided the lncRNAs into two subtypes. Subsequently, we established a signature model consisting of seven cuproptosis-related lncRNAs (AC073534.2, AC021321.1, HYI-AS1, PPP1R26-AS1, AC010328.1, AC012568.1 and MIR4435-2Hg) using least absolute shrinkage and selection operator regression. Survival analysis based on risk score showed that the overall survival and progression-free survival of patients in the high-risk group were worse than those in the low-risk group. Multivariate Cox analysis demonstrated the independent prognostic potential of this signature model for patients with BC. Moreover, age and clinical stage were also significantly correlated with prognosis. The constructed nomogram plots revealed good predictive power for the prognosis of patients with BC and were validated using calibration plots. Additionally, enrichment analysis, Single sample gene set enrichment analysis and immune infiltration abundance analysis revealed significant differences in immune infiltration between the two risk groups, with high levels of immune cell subset infiltrations observed in the high-risk group accompanied by various immune pathway activation. Moreover, almost all the immune checkpoint genes showed high expression levels in the high-risk group. Moreover, TIDE analysis suggested that the high-risk group was more responsive to immunotherapy. Finally, eight drugs with low IC50 values were screened, which may prove to be beneficial for patients in the high-risk group.

Novel Pyroptosis-Related Gene Signatures Identified as the Prognostic Biomarkers for Bladder Carcinoma

Background

Bladder carcinoma (BLCA) is a common malignant tumor with high morbidity and mortality in the urinary system. Pyroptosis is a pattern of programmed cell death that is closely associated with progression of tumors. Therefore, it is significant to probe the expression of pyroptosis-related genes (PRGs) in BLCA.

Methods

The differentially expressed genes in normal and BLCA tissues were first obtained from the Cancer Genome Atlas (TCGA) database analysis, as well as PRGs from the National Center for Biotechnology Information (NCBI) database, intersecting to obtain differentially expressed pyroptosis-related genes (DEPRGs) in BLCA. With the construction of a prognostic model of pyroptosis by regression analysis, we derived and validated key genes, which were ascertained as a separate prognostic marker by individual prognostic and clinical relevance analysis. In addition, we gained six immune cells from the Tumor Immune Evaluation Resource (TIMER) website and analyzed the relationship between pyroptosis prognostic genes and immune infiltration.

Result

Our results revealed that 31 DEPRGs were available by comparing normal and BLCA tissues with |log2 (fold change, FC)| &gt; 0.5 and FDR &lt;0.05. Four key genes (CRTAC1, GSDMB, AIM2, and FOXO3) derived from the pyroptosis prognostic model were experimentally validated for consistent expression in BLCA patients. Following risk scoring, the low-risk group of BLCA patients had noticeably higher overall survival (OS) than the high-risk group (p &lt; 0.001). Risk score was still an independent prognostic factor (HR = 1.728, 95% CI =1.289–2.315, p &lt; 0.001). In addition, we found remarkable correlations among the expression of pyroptosis-related prognostic genes and the immune infiltration of CD4+ T cells, CD8+ T cells, B cells, dendritic cells, macrophages, and neutrophils.

Conclusion

Genes (CRTAC1, GSDMB, AIM2, and FOXO3) associated with pyroptosis are potential BLCA prognostic biomarkers that act as an essential part in the predictive prognosis of survival and immunotherapy of BLCA.

Gardnerella Exposures Alter Bladder Gene Expression and Augment Uropathogenic Escherichia coli Urinary Tract Infection in Mice

The anaerobic actinobacterium Gardnerella was first isolated from the bladder by suprapubic aspiration more than 50 years ago. Since then, Gardnerella has been increasingly recognized as a common and often abundant member of the female urinary microbiome (urobiome). Some studies even suggest that the presence of Gardnerella is associated with urological disorders in women. We recently reported that inoculation of Gardnerella into the bladders of mice results in urothelial exfoliation. Here, we performed whole bladder RNA-seq in our mouse model to identify additional host pathways involved in the response to Gardnerella bladder exposure. The transcriptional response to Gardnerella reflected the urothelial turnover that is a consequence of exfoliation while also illustrating the activation of pathways involved in inflammation and immunity. Additional timed exposure experiments in mice provided further evidence of a potentially clinically relevant consequence of bladder exposure to Gardnerella-increased susceptibility to subsequent UTI caused by uropathogenic Escherichia coli. Together, these data provide a broader picture of the bladder's response to Gardnerella and lay the groundwork for future studies examining the impact of Gardnerella on bladder health.

N6-methyladenosine modifications of mRNAs and long noncoding RNAs in oxygen-induced retinopathy in mice

Retinal neovascular diseases are major causes of blindness worldwide. As a common epitranscriptomic modification of eukaryotic RNAs, N⁶-methyladenosine (m⁶A) is associated with the pathogenesis of many diseases, including angiogenesis, through the regulation of RNA metabolism and functions. The aim of this study was to identify m⁶A modifications of mRNAs and long noncoding RNAs (lncRNAs) and determine their potential roles in retinal neovascularization. The transcriptome-wide m⁶A profiles of mRNAs and lncRNAs in the retinal tissues of mice with oxygen-induced retinopathy (OIR) and controls were identified by microarray analysis of immunoprecipitated methylated RNAs. The m⁶A methylation levels of mRNAs and lncRNAs identified in the microarray data were validated by MeRIP-qPCR. A total of 1321 mRNAs (151 hypermethylated and 1170 hypomethylated) and 192 lncRNAs (15 hypermethylated and 177 hypomethylated) were differentially methylated with the m⁶A modification in OIR and control mice. Gene ontology analysis showed that hypermethylated mRNAs were enriched in the regulation of multicellular organismal process, intracellular organelle, and protein binding, while hypomethylated mRNAs were enriched in cellular metabolic process, intracellular process, and binding. Kyoto Encyclopedia of Genes and Genomes pathway analysis indicated that hypermethylated mRNAs were involved in dopaminergic synapses, glutamatergic synapse, and PI3K-Akt signaling pathway, while hypomethylated mRNAs were involved in autophagy, ubiquitin-mediated proteolysis, and spliceosome. Moreover, the altered levels of m⁶A methylation of ANGPT2, GNG12, ROBO4, and ENSMUST00000153785 were validated by MeRIP-qPCR. The results revealed an altered m⁶A epitranscriptome in OIR retinas. These methylated RNAs may act as novel modulators and targets in retinal neovascularization.

Prognosis Risk Model Based on Pyroptosis-Related lncRNAs for Bladder Cancer

Objective. Bladder cancer (BC) is the most common malignancy in the urinary system and is prone to recurrence and metastasis. Pyroptosis is a kind of cell necrosis that is triggered by the gasdermin protein family. lncRNAs are noncoding RNAs that are more than 200 nucleotides long. Both pyroptosis and lncRNAs are associated with tumor development and progression. This study is aimed at exploring and establishing a prognostic signature of BC based on pyroptosis-related lncRNAs. Methods. In this study, The Cancer Genome Atlas (TCGA) database provided us with the RNA sequencing transcriptome data of bladder cancer patients, and we identified differentially expressed pyroptosis-related lncRNAs in bladder cancer. Then, the prognostic significance of these lncRNAs was assessed using univariate Cox regression analysis and LASSO regression analysis. Subsequently, 4 pyroptosis-related lncRNAs, namely, AL121652.1, AL161729.4, AC007128.1, and AC124312.3, were identified by multivariate Cox regression analysis, thus constructing the prognostic risk model. Then, we compared the levels of immune infiltration, differences in cell function, immune checkpoints, and m6A-related gene expression levels between the high- and low-risk groups. Result. Patients were divided into low-risk or high-risk groups based on the median risk score. Kaplan–Meier survival analysis indicated that the overall survival of bladder cancer patients in the low-risk group was substantially superior to that in the high-risk group ( $p<0.001$ ). The receiver operating characteristic (ROC) curve further confirmed the credibility of our model. Moreover, gene set enrichment analysis (GSEA) indicated that these were different signal pathways significantly enriched between the two groups. Immune infiltration, immune checkpoint, and N6-methyladenosine-related gene analysis also reflected that there were notable differences between the two groups. Conclusion. Therefore, this prognostic risk model is based on the level of pyroptotic lncRNAs, which is conducive to individualized assessment of the risk of patients and provides a reference for clinical treatment. This will also help provide insights into the prognosis and treatment of bladder cancer.