Role of vascular endothelial growth factor D in lung adenocarcinoma immunotherapy response.

OBJECTIVE: To identify key genes associated with tumor-associated macrophages (TAMs), tumor immunotherapy, in the prognosis of lung adenocarcinoma (LUAD). METHODS: The mRNA expression profiles of LUAD samples were obtained from The Cancer Genome Atlas (TCGA) database. The "CIBERSORT" R package was employed to calculate the proportion of innate immune cell infiltration in both tumor and adjacent normal tissues. TAM-associated genes in LUAD were identified to construct a prognostic risk model using weighted gene correlation network analysis (WGCNA), Least Absolute Shrinkage and Selection Operator (LASSO), and multivariate Cox regression analyses (COX). The IMvigor210 cohort was utilized to validate the roles of these genes as predictors of immunotherapy response. Tissue microarrays, immunofluorescence staining, and mRNA level detection methods were used to determine the correlation of risk factors in LUAD tissues. RESULTS: CIBERSORT analysis revealed significant differences in innate immune cells between tumor and adjacent tissues. Seventy-four differential genes linked to these cells were identified from WGCNA. Four hub genes (endothelin receptor type B, vascular endothelial growth factor D (VEGFD), latent transforming growth factor beta binding protein 4 (LTBP4), and fibroblast growth factor receptor 4 (FGFR4)) in the TAM prognostic model were identified as independent prognostic risk factors (P < 0.05). VEGFD expression was identified as a low-risk factor for LUAD prognosis prediction (P < 0.05). Moreover, low-risk patients exhibited higher sensitivity to anti-PD-L1 therapy compared to high-risk patients (P < 0.05). VEGFD levels were negatively correlated with programmed cell death 1 (PD-1) levels (r = -0.363; P < 0.05), suggesting that VEGFD may serve as a predictor for anti-PD-1 treatment. CONCLUSIONS: VEGFD is associated with innate immunity in LUAD, it can predict LUAD prognosis, and therefor may be a potential predictor for anti-PD-1 treatment in patients with LUAD.

Integrative multi-omics analysis reveals the underlying toxicological mechanisms of enrofloxacin on the growth of wheat seedling roots.

The continuous release of antibiotics into agroecosystems has raised concerns about the potential negative effects of antibiotic residues on crops. In this study, the toxicological effects of enrofloxacin (ENR) on wheat seedlings were analyzed using a combination of morpho-physiological, transcriptomic, proteomic, and metabolomic approaches. ENR inhibited the growth of wheat (Triticum aestivum L.) roots and induced oxidative stress. In particular, ENR downregulated the oxidative phosphorylation pathway, while it enhanced glycolysis and the tricarboxylic acid cycle, thereby regulating the balance of intracellular energy metabolism. In addition, sustained exposure to excessive reactive oxygen species (ROS) resulted in an increase in reduced glutathione (GSH), a slight decrease in ascorbic acid (AsA), and a significant decrease in the ratio of GSH to oxidized glutathione (GSSG), which imbalanced the AsA-GSH cycle. In addition, the resulting increase in abnormal proteins triggered ubiquitin-independent proteasomal degradation pathways. Further, an increase in abscisic acid (ABA) and a decrease in jasmonic acid (JA) and its derivatives alleviated the inhibitory effect of ENR on the growth of wheat roots. In conclusion, direct damage and signaling by ROS, hormonal regulation, a decrease in the GSH to GSSG ratio, and insufficient energy supply were identified as key factors for the significant inhibition of wheat root growth under ENR stress.

Unraveling Immune Heterogeneity across Pan-Cancer and Deep Insights in Lung Adenocarcinoma based on Alternative Splicing.

Alternative splicing (AS) participates in tumor development and tumor microenvironment formation. However, the landscape of immune infiltrating AS events (IIASE) in pan-cancer and mechanisms of AS in lung adenocarcinoma (LUAD) have not been comprehensively characterized. We systematically profiled the IIASE landscape of pan-cancer using data from The Cancer Genome Atlas (TCGA), analyzing both commonalities and specific characteristics among different cancer types. We found that AS events tend to occur specifically in one cancer type rather than in multiple cancer types. AS events were used to classify 512 LUAD samples into two subtypes by unsupervised clustering: aberrant splicing subtype (ABS) and immune infiltrating subtype (IIS). The two subtypes showed significant differences in clinicopathology, prognosis, transcriptomics, genomics and immune microenvironment. We constructed a classification signature comprising 10 genes involved in 14 AS events using Logistic regression. The robustness of the signature was validated in three independent datasets using survival analysis. To explore AS mechanisms in LUAD, we constructed subtype-specific co-expression networks using Pearson correlation analysis. AS event of AKT3 regulated by splicing factor ENOX1 was associated with poor prognosis in LUAD. Overall, we outline AS events associated with immune infiltration in pan-cancer and this study provides insights into AS mechanisms in LUAD patient classification.

Individualized detection of TMPRSS2-ERG fusion status in prostate cancer: a rank-based qualitative transcriptome signature.

BACKGROUND: TMPRSS2-ERG (T2E) fusion is highly related to aggressive clinical features in prostate cancer (PC), which guides individual therapy. However, current fusion prediction tools lacked enough accuracy and biomarkers were unable to be applied to individuals across different platforms due to their quantitative nature. This study aims to identify a transcriptome signature to detect the T2E fusion status of PC at the individual level. METHODS: Based on 272 high-throughput mRNA expression profiles from the Sboner dataset, we developed a rank-based algorithm to identify a qualitative signature to detect T2E fusion in PC. The signature was validated in 1223 samples from three external datasets (Setlur, Clarissa, and TCGA). RESULTS: A signature, composed of five mRNAs coupled to ERG (five ERG-mRNA pairs, 5-ERG-mRPs), was developed to distinguish T2E fusion status in PC. 5-ERG-mRPs reached 84.56% accuracy in Sboner dataset, which was verified in Setlur dataset (n = 455, accuracy = 82.20%) and Clarissa dataset (n = 118, accuracy = 81.36%). Besides, for 495 samples from TCGA, two subtypes classified by 5-ERG-mRPs showed a higher level of significance in various T2E fusion features than subtypes obtained through current fusion prediction tools, such as STAR-Fusion. CONCLUSIONS: Overall, 5-ERG-mRPs can robustly detect T2E fusion in PC at the individual level, which can be used on any gene measurement platform without specific normalization procedures. Hence, 5-ERG-mRPs may serve as an auxiliary tool for PC patient management.

Up-regulated SPP1 increases the risk from IPF to lung cancer via activating the pro-tumor macrophages.

The incidence of lung cancer (LC) in Idiopathic Pulmonary Fibrosis (IPF) patients is more than twice that in non-IPF. This study aims to investigate IPF-to-LC pathogenesis and to develop a predictor for detecting IPF predisposing patients to LC. We conducted unsupervised clustering to detect high-risk subtypes from IPF to LC. Subsequently, we performed single-cell RNA-seq analysis to characterize high-risk IPF by examining the immune microenvironment. We identified 42 common immune function-related pathogenic genes between IPF and LC. We developed an LC risk classifier for IPF patients, comprising five genes: SPP1, MMP9, MMP12, FABP4, and IL1B. The five-gene classifier can successfully distinguish the high-risk population from IPF patients. High-risk IPF patients exhibited an immunosuppressive microenvironment with higher oncogene expression than low-risk patients. Single-cell analysis revealed that SPP1+ macrophages at the terminal of macrophages' developmental trajectory may promote the progression from IPF to LC. The strong crosstalk between SPP1+ macrophages and inflammation-related cancer-associated fibroblasts promoted the tumorigenic process in IPF. In vitro, assays showed that co-culturing macrophages overexpressing SPP1 with MRC-5 cells induced the transition of fibroblasts into cancer-associated fibroblasts. SPP1 produced by macrophages promoted epithelial-mesenchymal transition in alveolar epithelial cells via stimulating the upregulation of N-cadherin and Vimentin in MLE-12 cells. This study provided a novel method to identify the LC risk population from IPF, revealing the cellular interactions involved in the transition from IPF to LC. Our findings highlighted SPP1 as a critical driver in IPF progression, offering a potential target for therapy in fibrosis.

Mutant RB1 enhances therapeutic efficacy of PARPis in lung adenocarcinoma by triggering the cGAS/STING pathway.

Poly (ADP-ribose) polymerase inhibitors (PARPis) are approved for cancer therapy according to their synthetic lethal interactions, and clinical trials have been applied in non-small cell lung cancer. However, the therapeutic efficacy of PARPis in lung adenocarcinoma (LUAD) is still unknown. We explored the effect of a mutated retinoblastoma gene (RB1) on PARPi sensitivity in LUAD. Bioinformatic screening was performed to identify PARPi-sensitive biomarkers. Here, we showed that viability of LUAD cell lines with mutated RB1 was significantly decreased by PARPis (niraparib, rucaparib, and olaparib). RB1 deficiency induced genomic instability, prompted cytosolic double-stranded DNA (dsDNA) formation, activated the cGAS/STING pathway, and upregulated downstream chemokines CCL5 and CXCL10, triggering immune cell infiltration. Xenograft experiments indicated that PARPi treatment reduced tumorigenesis in RB1-KO mice. Additionally, single-cell RNA sequencing analysis showed that malignant cells with downregulated expression of RB1 had more communications with other cell types, exhibiting activation of specific signaling such as GAS, IFN response, and antigen-presenting and cytokine activities. Our findings suggest that RB1 mutation mediates the sensitivity to PARPis through a synthetic lethal effect by triggering the cGAS/STING pathway and upregulation of immune infiltration in LUAD, which may be a potential therapeutic strategy.

Identification of a coagulation-related signature correlated with immune infiltration and their prognostic implications in lung adenocarcinoma.

BACKGROUND: Lung adenocarcinoma (LUAD) is a fatal form of lung cancer with a poor prognosis. Coagulation system had been confirmed closely related to tumor progression and the hypercoagulable state encouraged the immune infiltration and development of tumor cells, leading to a poor prognosis in cancer patients. However, the use of the coagulation-related genes (CRGs) for prognosis in LUAD has yet to be determined. In this study, we constructed an immune-related signature (CRRS) and identified a potential coagulation-related biomarker (P2RX1). METHODS: We obtained a total of 209 CRGs based on two coagulation-related KEGG pathways, then developed the CRRS signature by using the TCGA-LUAD RNA-seq data via the procedure of LASSO-Cox regression, stepwise-Cox regression, univariate and multivariate Cox regression. Grouped by the CRRS, Kaplan-Meier survival curves and receiver operating characteristic curves were drawn for the training and validation sets, respectively. In addition, single-sample gene set enrichment analysis was exploited to explore immune infiltration level. Moreover, immunophenotypes and immunotherapy grouped by CRRS were further analyzed. RESULTS: We developed an immune-related signature (CRRS) composed of COL1A2, F2, PLAUR, C4BPA, and P2RX1 in LUAD. CRRS was an independent risk factor for overall survival and displayed stable and powerful performance. Additionally, CRRS possessed distinctly superior accuracy than traditional clinical variables and molecular features. Functional analysis indicated that the differentially high expressed genes in the low-risk group significantly enriched in T cell and B cell receptor signaling pathways. The low-risk group was sensitive to anti-PD-1/PD-L1 immunotherapy and displayed abundant immune infiltration and immune checkpoint gene expression. Finally, we identified an independent prognostic gene P2RX1. Low expression of P2RX1 associated with poor overall survival and decreased immune infiltration. CONCLUSIONS: Our study revealed a significant correlation between CRRS and immune infiltration. CRRS could serve as a promising tool to improve the clinical outcomes for individual LUAD patients.

Erratum: [Corrigendum] MicroRNA­23a inhibits endometrial cancer cell development by targeting SIX1.

[This corrects the article DOI: 10.3892/ol.2019.10694.].

Synthetic viability induces resistance to immune checkpoint inhibitors in cancer cells.

BACKGROUND: Immune checkpoint inhibitors (ICI) have revolutionized the treatment for multiple cancers. However, most of patients encounter resistance. Synthetic viability (SV) between genes could induce resistance. In this study, we established SV signature to predict the efficacy of ICI treatment for melanoma. METHODS: We collected features and predicted SV gene pairs by random forest classifier. This work prioritized SV gene pairs based on CRISPR/Cas9 screens. SV gene pairs signature were constructed to predict the response to ICI for melanoma patients. RESULTS: This study predicted robust SV gene pairs based on 14 features. Filtered by CRISPR/Cas9 screens, we identified 1,861 SV gene pairs, which were also related with prognosis across multiple cancer types. Next, we constructed the six SV pairs signature to predict resistance to ICI for melanoma patients. This study applied the six SV pairs signature to divide melanoma patients into high-risk and low-risk. High-risk melanoma patients were associated with worse response after ICI treatment. Immune landscape analysis revealed that high-risk melanoma patients had lower natural killer cells and CD8+ T cells infiltration. CONCLUSIONS: In summary, the 14 features classifier accurately predicted robust SV gene pairs for cancer. The six SV pairs signature could predict resistance to ICI.

A Combination of Genome Mining with an OSMAC Approach Facilitates the Discovery of and Contributions to the Biosynthesis of Melleolides from the Basidiomycete Armillaria tabescens.

Genome mining revealed that the genomes of basidiomycetes may include a considerable number of biosynthetic gene clusters (BGCs), yet numerous clusters remain unidentified. Herein, we report a combination of genome mining with an OSMAC (one strain, many compounds) approach to characterize the spectrum of melleolides produced by Armillaria tabescens CPCC 401429. Using F1 fermentation medium, the metabolic pathway of the gene cluster mel was successfully upregulated. From the extracts of the wild-type strain, two new melleolides (1 and 2), along with five new orsellinic acid-derived lactams (10-14), were isolated, and their structures were elucidated by LC-HR-ESIMS/MS and 2D-NMR. Several melleolides exhibited moderate anti-carcinoma (A549, NCI-H520, and H1299) effects with IC50 values of 4.0-48.8 µM. RNA-sequencing based transcriptomic profiling broadened our knowledge of the genetic background, regulation, and mechanisms of melleolide biosynthesis. These results may promote downstream metabolic engineering studies of melleolides. Our study demonstrates the approach is effective for discovering new secondary metabolites from Armillaria sp. and will facilitate the mining of the unexploited biosynthetic potential in other basidiomycetes.

Polyglycerol Functionalized 10 B Enriched Boron Carbide Nanoparticle as an Effective Bimodal Anticancer Nanosensitizer for Boron Neutron Capture and Photothermal Therapies.

Boron neutron capture therapy (BNCT) is a non-invasive cancer treatment with little adverse effect utilizing nuclear fission of 10 B upon neutron irradiation. While neutron source has been developed from a nuclear reactor to a compact accelerator, only two kinds of drugs, boronophenylalanine and sodium borocaptate, have been clinically used for decades despite their low tumor specificity and/or retentivity. To overcome these challenges, various boron-containing nanomaterials, or "nanosensitizers", have been designed based on micelles, (bio)polymers and inorganic nanoparticles. Among them, inorganic nanoparticles such as boron carbide can include a much higher 10 B content, but successful in vivo applications are very limited. Additionally, recent reports on the photothermal effect of boron carbide are motivating for the addition of another modality of photothermal therapy. In this study, 10 B enriched boron carbide (10 B4 C) nanoparticle is functionalized with polyglycerol (PG), giving 10 B4 C-PG with enough dispersibility in a physiological environment. Pharmacokinetic experiments show that 10 B4 C-PG fulfills the following three requirements for BNCT; 1) low intrinsic toxicity, 2) 10 B in tumor/tumor tissue (wt/wt) ≥ 20 ppm, and 3) 10 B concentrations in tumor/blood ≥ 3. In vivo study reveals that neutron irradiation after intravenous administration of 10 B4 C-PG suppresses cancer growth significantly and eradicates cancer with the help of near-infrared light irradiation.

A transcriptional signature detects homologous recombination deficiency in pancreatic cancer at the individual level.

Pancreatic cancer (PC) with homologous recombination deficiency (HRD) has been reported to benefit from poly ADP-ribose polymerase (PARP) inhibitors. However, accurate identification of HRD status for PC patients from the transcriptional level is still a great challenge. Here, based on a relative expression ordering (REO)-based algorithm, we developed an HRD signature including 24 gene pairs (24-GPS) using PC transcriptional profiles from The Cancer Genome Atlas (TCGA). HRD samples classified by 24-GPS showed worse overall survival (p = 4.4E-3 for TCGA; p = 1.2E-3 for International Cancer Genome Consortium-Australia cohort; p = 6.4E-2 for GSE17891; p = 7.5E-2 for GSE57495) and higher HRD scores than non-HRD samples (p = 1.4E-4). HRD samples showed highly unstable genomic characteristics and also displayed HRD-related alterations at the epigenomic and proteomic levels. Moreover, HRD cell lines identified by 24-GPS tended to be sensitive to PARP inhibitors (p = 6.6E-2 for olaparib; p = 2.6E-3 for niraparib). Compared with the non-HRD group, the HRD group presented lower immune scores and CD4/CD8 T cell infiltration proportion. Interestingly, PC tumor cells with co-inhibition of PARP-related genes and ATR showed reduced survival ability. In conclusion, 24-GPS can robustly identify PC patients with HRD status at the individualized level.

Bioinformatics Analysis Reveals MCM3 as an Important Prognostic Marker in Cervical Cancer.

The minichromosome maintenance complex 3 (MCM3) is essential for the regulation of DNA replication and cell cycle progression. However, the expression and prognostic values of MCM3 in cervical cancer (CC) have not been well-studied. Herein, we investigated the expression patterns and survival data of MCM3 in cervical cancer patients from the ONCOMINE, GEPIA, Human Protein Atlas, UALCAN, Kaplan-Meier Plotter, and LinkedOmics databases. The expression level of MCM3 is negatively correlated with advanced tumor stage and metastatic status. Specifically, MCM3 is significantly differentially expressed between patients in stage 1 and stage 3 cervical cancer with p value 0.0138. Similarly, the p values between stage 1 and stage 4 cervical cancer, between stage 2 and stage 3, and between stage 2 and stage 4 are 0.00089, 0.0244, and 0.00197, respectively. Not only that, cervical cancer patients with high mRNA expression of MCM3 may indicate longer overall survival but indicate shorter relapse-free survival. PRIM2 and MCM6 are positively correlated genes of MCM3. Bioinformatics analysis revealed that MCM3 might be considered a biological indicator for prognostic evaluation of cervical cancer. However, it is currently limited to bioinformatics analysis, and more clinical tissue specimens and cell experiments are needed to further explore the role of MCM3 in the occurrence and progression of cervical cancer.

Nickel-Catalyzed Cross-Electrophile C(sp3)-Si Coupling of Unactivated Alkyl Bromides with Vinyl Chlorosilanes.

Cross-electrophile C-Si coupling has emerged as a promising tool for the construction of organosilanes, but the potential of this method remains largely unexplored. Herein, we report a C(sp3)-Si coupling of unactivated alkyl bromides with vinyl chlorosilanes. The reaction proceeds under mild conditions, and it offers a new approach to alkylsilanes. Functionalities such as Grignard-sensitive groups (e.g., acid, amide, alcohol, ketone, and ester), acid-sensitive groups (e.g., ketal and THP protection), alkyl fluoride and chloride, aryl bromide, alkyl tosylate and mesylate, silyl ether, and amine were tolerated. Incorporation of the -Si(vinyl)R2 moiety into complex molecules and the immobilization of a glass surface by formed organosilanes were demonstrated.

Discovering a qualitative transcriptional signature of homologous recombination defectiveness for prostate cancer.

The discovery of homologous recombination deficiency (HRD) biomarkers in prostate cancer is important for patients who will benefit from poly ADP-ribose polymerase inhibitor (PARPi). Here, we developed a transcriptional homologous recombination defectiveness (HRDness) signature, comprising 16 gene pairs (16-GPS), for prostate cancer by a relative expression ordering (REO)-based discovery procedure. Subsequently, two newly subtypes classified by 16-GPS showed a higher significance level in various clinicopathological and HRD features than subtypes obtained by other methods, such as HRDetect. HRDness subtype also displayed more aggressive features and higher genomics scores than non-HRDness in three independent datasets. HRDness prostate cancer cells were more sensitive to PARPi than non-HRDness. Moreover, the HRDness samples showed distinct multi-omics characteristics related to homologous recombination repair function loss. Overall, the newly proposed qualitative signature can robustly determine the HRD status for prostate cancer at the personalized level, and especially be an auxiliary tool for PARPi treatment strategy.

Revealing biomarkers associated with PARP inhibitors based on genetic interactions in cancer genome.

Poly (ADPribose) polymerase inhibitors (PARPis) are clinically approved drugs designed according to the concept of synthetic lethality (SL) interaction. It is crucial to expand the scale of patients who can benefit from PARPis, and overcome drug resistance associated with it. Genetic interactions (GIs) include SL and synthetic viability (SV) that participate in drug response in cancer cells. Based on the hypothesis that mutated genes with SL or SV interactions with PARP1/2/3 are potential sensitive or resistant PARPis biomarkers, respectively, we developed a novel computational method to identify them. We analyzed fitness variation of cell lines to identify PARP1/2/3-related GIs according to CRISPR/Cas9 and RNA interference functional screens. Potential resistant/sensitive mutated genes were identified using pharmacogenomic datasets. We identified 41 candidate resistant and 130 candidate sensitive PARPi-response related genes, and observed that EGFR with gain-of-function mutation induced PARPi resistance, and predicted a combination therapy with PARP inhibitor (veliparib) and EGFR inhibitor (erlotinib) for lung cancer. We also revealed that a resistant gene set (TNN, PLEC, and TRIP12) in lower grade glioma and a sensitive gene set (BRCA2, TOP3A, and ASCC3) in ovarian cancer, which were associated with prognosis. Thus, cancer genome-derived GIs provide new insights for identifying PARPi biomarkers and a new avenue for precision therapeutics.

Virtual or real: lifelike cinematic rendering of adrenal tumors.

The adrenal gland is small in size and hidden in location. Adrenal tumors are relatively difficult to diagnose due to the wide variety of tumors and partial overlap of image features. Cinematic rendering (CR) is a novel, three-dimensional post-processing technology that simulates how light propagates in the real world, providing high-resolution visualizations that truly present subtle anatomical details. We retrospectively collected a series of pathologically confirmed adrenal tumor cases, raw data was introduced into the post-processing workstation, and different tools and templates of CR software were used for reconstruction and rendering. Compared with traditional black and white two-dimensional images and three-dimensional volume rendering (VR) images, CR images were more colorful, layered, and closer to the truth. CR has potential in diagnosing and preoperative planning of adrenal tumors, allowing vivid and realistic visualization of tumor location, morphology, different components (solid, cystic, fat, calcification, etc.), the pattern of enhancement, and the relationship with surrounding tissues and organs.

Effect of Protein Corona on Mitochondrial Targeting Ability and Cytotoxicity of Triphenylphosphonium Conjugated with Polyglycerol-Functionalized Nanodiamond.

Functionalization of nanoparticles (NPs) with targeting moieties has a high potential to advance precision nanomedicine. However, the targeting moieties on a NP surface are known to be masked by a protein corona in biofluids, lowering the targeting efficiency. Although it has been demonstrated at the cellular level, little is known about the influence of the protein corona on the subcellular targeting. Herein, we adopted triphenylphosphonium (TPP) as a mitochondrial targeting moiety and investigated the effects of protein coronas from fetal bovine serum and human plasma on its targeting ability and cytotoxicity. Specifically, we introduced TPP in low (l) and high (h) densities on the surface of nanodiamond (ND) functionalized with polyglycerol (PG). Despite the "corona-free" PG interface, we found that the TPP moiety attracted proteins to form a corona layer with clear linearity between the TPP density and the protein amount. By performing investigations on human cervix epithelium (HeLa) and human lung epithelial carcinoma (A549) cells, we further demonstrated that (1) the protein corona alleviated the cytotoxicity of both ND-PG-TPP-l and -h, (2) a smaller amount of proteins on the surface of ND-PG-TPP-l did not affect its mitochondrial targeting ability, and (3) a larger amount of proteins on the surface of ND-PG-TPP-h diminished its targeting specificity by restricting the NDs inside the endosome and lysosome compartments. Our findings will provide in-depth insights into the design of NPs with active targeting moiety for more precise and safer delivery at the subcellular level.

Individualized lncRNA differential expression profile reveals heterogeneity of breast cancer.

Long non-coding RNAs (lncRNAs) play key regulatory roles in breast cancer. However, population-level differential expression analysis methods disregard the heterogeneous expression of lncRNAs in individual patients. Therefore, we individualized lncRNA expression profiles for breast invasive carcinoma (BRCA) using the method of LncRNA Individualization (LncRIndiv). After evaluating the robustness of LncRIndiv, we constructed an individualized differentially expressed lncRNA (IDElncRNA) profile for BRCA and investigated the subtype-specific IDElncRNAs. The breast cancer subtype-specific IDElncRNA showed frequent co-occurrence with alterations of protein-coding genes, including mutations, copy number variation and differential methylation. We performed hierarchical clustering to subdivide TNBC and revealed mesenchymal subtype and immune subtype for TNBC. The TNBC immune subtype showed a better prognosis than the TNBC mesenchymal subtype. LncRNA PTOV1-AS1 was the top differentially expressed lncRNA in the mesenchymal subtype. And biological experiments validated that the upregulation of PTOV1-AS1 could downregulate TJP1 (ZO-1) and E-Cadherin, and upregulate Vimentin, which suggests PTOV1-AS1 may promote epithelial-mesenchymal transition and lead to migration and invasion of TNBC cells. The mesenchymal subtype showed a higher fraction of M2 macrophages, whereas the immune subtype was more associated with CD4 + T cells. The immune subtype is characterized by genomic instability and upregulation of immune checkpoint genes, thereby suggesting a potential response to immunosuppressive drugs. Last, drug response analysis revealed lncRNA ENSG00000230082 (PRRT3-AS1) is a potential resistance biomarker for paclitaxel in BRCA treatment. Our analysis highlights that IDElncRNAs can characterize inter-tumor heterogeneity in BRCA and the new TNBC subtypes indicate novel insights into TNBC immunotherapy.