Disparity in Tumor Immune Microenvironment of Breast Cancer and Prognostic Impact: Asian Versus Western Populations

Clinical outcomes and prognostic factors of concurrent chemoradiotherapy for anal squamous cell carcinoma in Japan

BACKGROUND

Concurrent chemoradiotherapy (CCRT) is the standard treatment for locoregional anal squamous cell carcinoma (ASCC) in western countries. However, there have been few reports on the clinical outcomes of CCRT in Japan. This study aimed to evaluate the clinical outcomes of CCRT, prognostic factors, and the clinical impact of programmed cell death-ligand 1 (PD-L1) expression of ASCC in Japan.

METHODS

Patients with locoregional ASCC were enrolled between 2007 and 2017. All patients received CCRT consisting of ≥ 45 Gy of radiation, 5-fluorouracil, and mitomycin C. Disease-free survival (DFS), overall survival (OS), and adverse events (AEs) were estimated. Expression of p16 and PD-L1 were assessed by immunohistochemical staining (IHC).

RESULTS

This study included 36 patients, of whom 30 (83.3%) were female. Among the participants, 32 (88.9%) achieved complete clinical remission, while six (16.7%) experienced recurrence. The five-year DFS and five-year OS were 72.2% and 84.7%, respectively. Grades ≥ 3 serious AEs included neutropenia in 10 (27.7%) and perianal dermatitis in eight (22.2%). In a univariate analysis, male sex, lymph node metastasis, and large tumor size were significantly associated with worse outcome. In a multivariate analysis, tumor size was an independent factor associated with short DFS. Of the 30 patients whose biopsy specimens were available for IHC, 29 (96.7%) were positive for p16, and 13 (43.3%) were positive for PD-L1. However, PD-L1 expression did not show any clinical impact.

CONCLUSIONS

The comparative etiology, clinical outcomes, and prognostic factors of CCRT observed in Japanese patients with locoregional ASCC were consistent with western data.

Genomic and transcriptomic profiling of inflammatory breast cancer reveals distinct molecular characteristics to non-inflammatory breast cancers

PURPOSE

Inflammatory breast cancer (IBC), a rare and highly aggressive form of breast cancer, accounts for 10% of breast cancer-related deaths. Previous omics studies of IBC have focused solely on one of genomics or transcriptomics and did not discover common differences that could distinguish IBC from non-IBC.

METHODS

Seventeen IBC patients and five non-IBC patients as well as additional thirty-three Asian breast cancer samples from TCGA-BRCA were included for the study. We performed whole-exon sequencing (WES) to investigate different somatic genomic alterations, copy number variants, and large structural variants between IBC and non-IBC. Bulk RNA sequencing (RNA-seq) was performed to examine the differentially expressed genes, pathway enrichment, and gene fusions. WES and RNA-seq data were further investigated in combination to discover genes that were dysregulated in both genomics and transcriptomics.

RESULTS

Copy number variation analysis identified 10 cytobands that showed higher frequency in IBC. Structural variation analysis showed more frequent deletions in IBC. Pathway enrichment and immune infiltration analysis indicated increased immune activation in IBC samples. Gene fusions including CTSC-RAB38 were found to be more common in IBC. We demonstrated more commonly dysregulated RAS pathway in IBC according to both WES and RNA-seq. Inhibitors targeting RAS signaling and its downstream pathways were predicted to possess promising effects in IBC treatment.

CONCLUSION

We discovered differences unique in Asian women that could potentially explain IBC etiology and presented RAS signaling pathway as a potential therapeutic target in IBC treatment.

Interaction of base excision repair gene polymorphism and estrogen-DNA adducts in breast cancer risk among East Asian women

PURPOSE

In East Asia, the incidence of breast cancer has been increasing rapidly, particularly among premenopausal women. An elevated ratio of estrogen-DNA adducts was linked to a higher risk of breast cancer. The present study explored the influence of the interaction between base excision repair (BER) gene polymorphisms and estrogen-DNA adducts on breast cancer risk.

METHODS

We conducted a case-control study comprising healthy volunteers and individuals with benign breast disease (control arm, n = 176) and patients with invasive carcinoma or carcinoma in situ (case arm, n = 177). Genotyping for BER-related genes, including SMUG1, OGG1, ERCC5, and APEX1, was performed. A logistic regression model, incorporating interactions between gene polymorphisms, estrogen-DNA adduct ratio, and clinical variables, was used to identify the risk factors for breast cancer.

RESULTS

Univariate analysis indicated marginal associations between breast cancer risk and APEX1 rs1130409 T > G (P = 0.057) and APEX1 rs1760944 T > G (P = 0.065). Multivariate regression analysis revealed significant associations with increased breast cancer risk for APEX1_rs1130409 (GT/GG versus TT) combined with a natural logarithmic value of the estrogen-DNA adduct ratio (estimated OR 1.164, P = 0.023) and premenopausal status with an estrogen-DNA adduct ratio > 2.93 (estimated OR 2.433, P = 0.001).

CONCLUSION

APEX1_rs1130409 (GT/GG versus TT) polymorphisms, which are related to decreased BER activity, combined with an increased ratio of estrogen-DNA adducts, increase the risk of breast cancer in East Asian women.

Identification and Validation of JAM-A as a Novel Prognostic and Immune Factor in Human Tumors

Junctional adhesion molecule-A (JAM-A), also known as F11 receptor (F11R), is a transmembrane glycoprotein that is involved in various biological processes, including cancer initiation and progression. However, the functional characteristics and significance of JAM-A in pan-cancer remain unexplored. In this study, we used multiple databases to gain a comprehensive understanding of JAM-A in human cancers. JAM-A was widely expressed in various tissues, mainly located on the microtubules and cell junctions. Aberrant expression of JAM-A was detected in multiple cancers at both mRNA and protein levels, which can be correlated with poorer prognosis and may be attributed to genetic alterations and down-regulated DNA methylation. JAM-A expression was also associated with immune infiltration and may affect immunotherapy responses in several cancers. Functional enrichment analysis indicated that JAM-A participated in tight junction and cancer-related pathways. In vitro experiments verified that JAM-A knockdown suppressed the proliferation and migration abilities of breast cancer cells and liver cancer cells. Overall, our study suggests that JAM-A is a pan-cancer regulator and a potential biomarker for predicting prognosis and immune-therapeutic responses for different tumors.

Clustering of HR + /HER2- breast cancer in an Asian cohort is driven by immune phenotypes

Breast cancer exhibits significant heterogeneity, manifesting in various subtypes that are critical in guiding treatment decisions. This study aimed to investigate the existence of distinct subtypes of breast cancer within the Asian population, by analysing the transcriptomic profiles of 934 breast cancer patients from a Malaysian cohort. Our findings reveal that the HR + /HER2- breast cancer samples display a distinct clustering pattern based on immune phenotypes, rather than conforming to the conventional luminal A-luminal B paradigm previously reported in breast cancers from women of European descent. This suggests that the activation of the immune system may play a more important role in Asian HR + /HER2- breast cancer than has been previously recognized. Analysis of somatic mutations by whole exome sequencing showed that counter-intuitively, the cluster of HR + /HER2- samples exhibiting higher immune scores was associated with lower tumour mutational burden, lower homologous recombination deficiency scores, and fewer copy number aberrations, implicating the involvement of non-canonical tumour immune pathways. Further investigations are warranted to determine the underlying mechanisms of these pathways, with the potential to develop innovative immunotherapeutic approaches tailored to this specific patient population.

Population-enriched innate immune variants may identify candidate gene targets at the intersection of cancer and cardio-metabolic disease

Both cancer and cardio-metabolic disease disparities exist among specific populations in the US. For example, African Americans experience the highest rates of breast and prostate cancer mortality and the highest incidence of obesity. Native and Hispanic Americans experience the highest rates of liver cancer mortality. At the same time, Pacific Islanders have the highest death rate attributed to type 2 diabetes (T2D), and Asian Americans experience the highest incidence of non-alcoholic fatty liver disease (NAFLD) and cancers induced by infectious agents. Notably, the pathologic progression of both cancer and cardio-metabolic diseases involves innate immunity and mechanisms of inflammation. Innate immunity in individuals is established through genetic inheritance and external stimuli to respond to environmental threats and stresses such as pathogen exposure. Further, individual genomes contain characteristic genetic markers associated with one or more geographic ancestries (ethnic groups), including protective innate immune genetic programming optimized for survival in their corresponding ancestral environment(s). This perspective explores evidence related to our working hypothesis that genetic variations in innate immune genes, particularly those that are commonly found but unevenly distributed between populations, are associated with disparities between populations in both cancer and cardio-metabolic diseases. Identifying conventional and unconventional innate immune genes that fit this profile may provide critical insights into the underlying mechanisms that connect these two families of complex diseases and offer novel targets for precision-based treatment of cancer and/or cardio-metabolic disease.

Pembrolizumab Plus Chemotherapy Followed by Pembrolizumab in Patients With Early Triple-Negative Breast Cancer: A Secondary Analysis of a Randomized Clinical Trial

Importance

In the phase 3 KEYNOTE-522 study, addition of pembrolizumab to neoadjuvant chemotherapy followed by adjuvant pembrolizumab significantly increased pathologic complete response (pCR) and event-free survival (EFS) vs neoadjuvant chemotherapy in patients with early triple-negative breast cancer.

Objective

To evaluate efficacy and safety outcomes for patients enrolled in East/Southeast Asia (Asia) in KEYNOTE-522.

Design, Setting, and Participants

KEYNOTE-522, a multicenter, double-blind, randomized clinical trial, enrolled 1174 patients between March 7, 2017, and September 13, 2018. For interim EFS and overall survival (OS) analyses (data cutoff, March 23, 2021), median follow-up was 39.8 months (range, 30.4-46.9 months) for pembrolizumab plus chemotherapy and 40.8 months (range, 30.1-46.9 months) for placebo plus chemotherapy. Data cutoff for pCR analysis was September 24, 2018. This secondary analysis included adults enrolled in Asia with newly diagnosed, previously untreated, nonmetastatic triple-negative breast cancer (tumor stage T1c and nodal stage N1-2 or tumor stage T2-4 and nodal stage N0-2) and Eastern Cooperative Oncology Group performance status of 0 to 1, regardless of programmed cell death ligand 1 (PD-L1) status.

Intervention

Patients were randomized 2:1 to 4 cycles of pembrolizumab (200 mg every 3 weeks) or placebo plus carboplatin and paclitaxel and another 4 cycles of pembrolizumab or placebo plus doxorubicin or epirubicin and cyclophosphamide before surgery. After definitive surgery, patients received pembrolizumab or placebo every 3 weeks for 9 cycles or until recurrence or unacceptable toxic effects.

Main Outcomes and Measures

The main outcome was pCR (no evidence of primary tumor after neoadjuvant therapy or carcinoma in situ after neoadjuvant therapy and no regional lymph node involvement after neoadjuvant therapy) at the time of definitive surgery and EFS.

Results

A total of 216 of 1174 randomized patients (all female; median [range] age, 46.0 [24.0-71.0] years) were from Korea, Japan, Taiwan, and Singapore (136 in the pembrolizumab plus chemotherapy group and 80 in the placebo plus chemotherapy group). Of these patients, 104 (76.5%) in the pembrolizumab plus chemotherapy group and 60 (75.0%) in the placebo plus chemotherapy group had a tumor PD-L1 combined positive score of 1 or greater. Pathologic complete response was 58.7% (95% CI, 46.7%-69.9%) with pembrolizumab plus chemotherapy and 40.0% (95% CI, 26.4%-54.8%) with placebo plus chemotherapy; benefit was observed regardless of PD-L1 status. Thirteen patients (9.6%) in the pembrolizumab plus chemotherapy group and 20 patients (25.0%) in the placebo plus chemotherapy group had EFS events (hazard ratio, 0.35; 95% CI, 0.17-0.71). The 36-month EFS rate was 91.2% (95% CI, 85.0%-94.9%) with pembrolizumab plus chemotherapy and 77.2% (95% CI, 66.3%-85.0%) with placebo plus chemotherapy. Grade 3 to 4 treatment-related adverse events occurred in 109 patients (80.1%) receiving pembrolizumab plus chemotherapy and 64 patients (81.0%) receiving placebo plus chemotherapy.

Conclusions and Relevance

In this subgroup analysis of patients enrolled in Asia in KEYNOTE-522, neoadjuvant pembrolizumab plus chemotherapy followed by adjuvant pembrolizumab led to clinically meaningful improvements in pCR and EFS vs neoadjuvant chemotherapy alone. These findings support the use of neoadjuvant pembrolizumab plus chemotherapy followed by adjuvant pembrolizumab as a standard-of-care therapy for patients in Asian countries with early triple-negative breast cancer.

Trial Registration

ClinicalTrials.gov Identifier: NCT03036488.

Gut microbiome influences incidence and outcomes of breast cancer by regulating levels and activity of steroid hormones in women

BACKGROUND

Breast cancer, the leading cancer type in women worldwide, is affected by reproductive and nonreproductive factors. Estrogen and progesterone influence the incidence and progression of breast cancer. The microbiome of the gut, a complex organ that plays a vital role in digestion and homeostasis, enhances availability of estrogen and progesterone in the host. Thus, an altered gut microbiome may influence the hormone-induced breast cancer incidence. This review describes the current understanding of the roles of gut microbiome in influencing the incidence and progression of breast cancer, with an emphasis on the microbiome-induced metabolism of estrogen and progesterone.

RECENT FINDINGS

Microbiome has been recognized as a promising hallmark of cancer. Next-generation sequencing technologies have aided in rapid identification of components of the gut microbiome that are capable of metabolizing estrogen and progesterone. Moreover, studies have indicated a wider role of the gut microbiome in metabolizing chemotherapeutic and hormonal therapy agents and reducing their efficacy in patients with breast cancer, with a predominant effect in postmenopausal women.

CONCLUSION

The gut microbiome and variations in its composition significantly alter the incidence and therapy outcomes of patients with breast cancer. Thus, a healthy and diverse microbiome is required for better response to anticancer therapies. Finally, the review emphasizes the requirement of studies to elucidate mechanisms that may aid in improving the gut microbiome composition, and hence, survival outcomes of patients with breast cancer.

Disparity in survival benefits of pembrolizumab between Asian and non-Asian patients with advanced cancers: A systematic review and meta-regression analysis

BACKGROUND

Immune checkpoint inhibitors have revolutionized the treatment of malignancies. However, disproportionate enrollment among races and ethnicities places the generalizability of global trial results in doubt.

METHODS

In this systematic review, phase 3 randomized controlled trials investigating pembrolizumab in advanced cancers and providing subgroup analyses of Asian and non-Asian participants were included. The primary and secondary effect measures were the mean differences (MDs) in the natural logarithms of the hazard ratios (HRs) for overall survival (OS) and progression-free survival (PFS) between these two subgroups, respectively. We used random-effects meta-analysis to calculate the pooled ratios of HRs (i.e., exp(MD)) and implemented a meta-regression analysis to identify significant covariates.

RESULTS

A total of 17 and 11 trials were included in the meta-analyses of OS and PFS, respectively. These trials included 2732 (25.49%) Asian and 7000 (65.32%) non-Asian participants in the OS analysis and 1438 (22.5%) Asian and 4129 (64.61%) non-Asian participants in the PFS analysis. The pooled ratio of HRs for OS was 0.87 (95% CI: 0.76-0.99; p = 0.0391), favoring Asian participants, but no significant difference was found in PFS (pooled ratio of HRs: 0.93; 95% CI: 0.82-1.07; p = 0.2391). Both linear meta-regression analyses revealed an open-label design as a crucial covariate, which indicated more benefits for non-Asian participants.

CONCLUSIONS

Compared with non-Asian patients, Asian patients with advanced cancers may derive superior OS benefits from pembrolizumab. Although the results warrant further exploration, this meta-analysis provides insight into clinical research design.

Landmark Series: The Cancer Genome Atlas and the Study of Breast Cancer Disparities

Race-related variation in breast cancer incidence and mortality are well-documented in the United States. The effect of genetic ancestry on disparities in tumor genomics, risk factors, treatment, and outcomes of breast cancer is less understood. The Cancer Genome Atlas (TCGA) is a publicly available resource that has allowed for the recent emergence of genome analysis research seeking to characterize tumor DNA and protein expression by ancestry as well as the social construction of race and ethnicity. Results from TCGA based studies support previous clinical evidence that demonstrates that American women with African ancestry are more likely to be afflicted with breast cancers featuring aggressive biology and poorer outcomes compared with women with other backgrounds. Data from TCGA based studies suggest that Asian women have tumors with favorable immune microenvironments and may experience better disease-free survival compared with white Americans. TCGA contains limited data on Hispanic/Latinx patients due to small sample size. Overall, TCGA provides important opportunities to define the molecular, biologic, and germline genetic factors that contribute to breast cancer disparities.

Psychological distress influences lung cancer: Advances and perspectives on the immune system and immunotherapy

Lung cancer has the highest incidence rate and mortality worldwide. Moreover, multiple factors may cause heterogeneity in the efficacy of immunotherapy for lung cancer, and preclinical studies have gradually uncovered the promotive effects of psychological distress (PD) on tumor hallmarks. Therefore, treatment targeted at PD may be a vital factor in adjusting and improving immunotherapy for lung cancer. Here, by focusing on the central nervous system, as well as stress-related crucial neurotransmitters and hormones, we highlight the effects of PD on the lung immune system, the lung tumor microenvironment (TME) and immunotherapy, which brings a practicable means and psychosocial perspective to lung cancer treatment.

Outcomes in breast cancer-does ethnicity matter?

Ethnic or racial differences in breast cancer (BC) survival outcomes have been reported, but current data are largely restricted to comparisons between African Americans and non-Hispanic whites. Most analyses have traditionally been based on self-reported race which may not always be accurate, or are oversimplified in their classification. With increasing globalization, quantification of the genetic ancestry from genomic data may offer a solution to infer the complex makeup from admixture of races. Focusing on the larger and the latest studies, we will discuss recent findings on the differing host and tumor biology that may be driving these disparities, in addition to the extrinsic environmental or lifestyle factors. Socioeconomic disparities with lower cancer literacy may lead to late presentation, poorer adherence to treatment, and other lifestyle factors such as unhealthy diet, obesity, and inadequate physical activity. These hardships may also result in greater allostatic load, which is in turn associated with aggressive BC features in disadvantaged populations. Epigenetic reprogramming may mediate the effects of the environment or lifestyle factors on gene expression, with ensuing differences in BC characteristics and outcome. There is increasing evidence that germline genetics can influence somatic gene alterations or expression, as well as modulate the tumor or immune microenvironment. Although the precise mechanisms remain elusive, this may account for the varying distribution of different BC subtypes across ethnicities. These gaps in our knowledge highlight the need to interrogate the multiomics landscape of BC in diverse populations, ideally in large-scale collaborative settings with standardized methodology for the comparisons to be statistically robust. Together with improving BC awareness and access to good quality health care, a holistic approach with insights of the biological underpinnings is much needed to eradicate ethnic disparities in BC outcomes.

Concomitant analyses of intratumoral microbiota and genomic features reveal distinct racial differences in breast cancer

Racial disparities are most accentuated among Black women as their lifetime risk of breast cancer incidence is lower than white and Asian women but their breast cancer related mortality is the highest among all races. Black women are more likely to develop triple-negative breast cancer at a younger age and harbor more aggressive tumors. In addition to tumor-centric alterations, tumor growth is also influenced by multiple other tumor microenvironment-related features, including resident immune cells and microbiota. Hence, in this study, we conduct concurrent genomic and metagenomic analyses, and uncover distinctive intratumoral microbial community compositions and tumor immune microenvironment-related traits in breast tumors from Asian, Black and white women. Interestingly, unique racially associated genomic nodes are found in the breast tumors from Asian, Black and white women. Examination of the cellular heterogeneity show differential enrichment of 11 out of 64 immune and stroma cell types in the breast tumors from different racial groups. In terms of microbial diversity, significant differences are revealed in alpha and beta-diversity measures. Intriguingly, potential race-specific microbial biomarkers of breast cancer are identified which significantly correlate with genes involved with tumor aggressiveness, angiogenesis, tumor cell migration and metastasis as well as oncogenic pathways-GLI and Notch. Investigating the metabolic features of intratumoral microbes, we find a significant differential enrichment of environmental information processing pathways, oncogenic pathways, and lipid metabolism pathways. Concomitantly investigating tumor-centric, tumor immune microenvironment-related and microbial alterations, our study provides a comprehensive understanding of racial disparities in breast cancer and warrants further exploration.

Racial/Ethnic Differences Among Tumor-Infiltrating Lymphocytes in Breast Cancer Tumors

PURPOSE

Tumor-infiltrating lymphocytes (TILs) have emerged as a predictor of breast cancer treatment response and patient outcomes. Current studies investigating racial/ethnic differences in TILs and immune profiles in breast cancer offer varying results. Our study provides some preliminary data in the breast cancer tumor microenvironment where there is a paucity of information, from Asian and Native Hawaiian/Pacific Islander (NHPI) racial/ethnic groups, not well represented in the literature.

METHODS

We reviewed 183 cases of women diagnosed with early stage breast cancer who received neoadjuvant treatment at 2 large health systems in Hawaii between 2008 and 2020. We evaluated clinical and demographic information including: age at diagnosis, self-reported race/ethnicity, tumor stage, tumor subtype according to ER, PR, and HER2 receptor status, the amount of TILs and pathologic complete response (pCR).

RESULTS

We found a significantly greater amount of TILs in Asians (37.7%, P = .01) and NHPI (37.2%, P = .02) patients compared to White patients on multivariate analysis. We found no significant differences in pCR among the different racial/ethnic groups.

CONCLUSIONS

Racial/ethnic differences in the amount of TILs in breast cancer tumors may suggest differences in the breast tumor microenvironment. This may in part contribute to known outcome disparities in these populations and should be further evaluated.

Prognostic value of tumor-infiltrating lymphocytes in DCIS: a meta-analysis

Background

Tumor infiltrating lymphocytes (TILs) have been shown to be associated with the prognosis of breast ductal carcinoma in situ (DCIS). In this systematic review and meta-analysis, we investigated the role of TILs and TIL subsets in predicting the recurrence risk of DCIS.

Method

PubMed, Medline, Web of Science, Embase and Cochrane were searched to identify publications investigating the prognostic role of TILs in DCIS. After study screening, data extraction and risk of bias assessment, a meta-analysis was performed to assess the association between TILs (total TILs, CD4+, CD8+, FOXP3+, PD-L1+ TILs) and the risk of DCIS recurrence.

Results

A pooled analysis indicated that dense stromal TILs in DCIS were associated with a higher recurrence risk (HR 2.11 (95% CI 1.35–3.28)). Subgroup analysis showed that touching TILs (HR 4.73 (95% CI 2.28–9.80)) was more precise than the TIL ratio (HR 1.49 (95% CI 1.11–1.99)) in estimating DCIS recurrence risk. Moreover, the prognostic value of TILs seemed more suitable for patients who are diagnosed with DCIS and then undergo surgery (HR 2.77, (95% CI 1.26–6.07)) or surgery accompanied by radiotherapy (HR 2.26, (95% CI 1.29–3.95)), than for patients who receive comprehensive adjuvant therapies (HR 1.16, (95% CI 1.35–3.28)). Among subsets of TILs, dense stromal PD-L1+ TILs were valuable in predicting higher recurrence risk of DCIS.

Conclusion

This systematic review and meta-analysis suggested a non-favorable prognosis of TILs and stromal PD-L1+ TILs in DCIS and indicated an appropriate assessment method for TILs and an eligible population.

The impact of race and ethnicity in breast cancer-disparities and implications for precision oncology

Breast cancer is the most commonly diagnosed cancer worldwide and is one of the leading causes of cancer death. The incidence, pathological features, and clinical outcomes in breast cancer differ by geographical distribution and across racial and ethnic populations. Importantly, racial and ethnic diversity in breast cancer clinical trials is lacking, with both Blacks and Hispanics underrepresented. In this forum article, breast cancer researchers from across the globe discuss the factors contributing to racial and ethnic breast cancer disparities and highlight specific implications of precision oncology approaches for equitable provision of breast cancer care to improve outcomes and address disparities.

Clinicopathological characteristics and treatment outcomes of luminal B1 breast cancer in Taiwan

BACKGROUND

Hormone receptor-positive, human epidermal growth factor receptor II (HER2)-negative luminal B1 breast cancer is associated with a higher risk of disease relapse than luminal A breast cancer. Therefore, we assessed and compared the distant metastasis pattern and clinical outcomes associated with luminal B1 and luminal A breast cancer in an Asian population.

METHODS

In this observational study, we assessed patients with estrogen receptor-positive, HER2-negative breast cancer who underwent surgery from 2009 to 2016. Patients were classified into luminal A or luminal B1 subsets via immunohistochemical analysis. Disease-free survival, post-metastasis survival, and overall survival were estimated; time to disease relapse and patterns of distant metastasis were compared. Risk of relapse and mortality were assessed using Cox proportional hazards model.

RESULTS

Patients with luminal B1 breast cancer (n = 677) were significantly younger and had larger tumors and a higher degree of affected axillary lymph nodes, lymphovascular invasion, and tumor necrosis than those with luminal A breast cancer (n = 630). Higher rates of local recurrence and distant metastasis were observed for luminal B1 (both p < 0.05); however, no difference was observed in the specific distant metastatic sites. We observed a significant increase in disease relapse risk in luminal B1 patients compared with that in luminal A (hazard ratio: 2.157, 95% CI: 1.340-3.473, p < 0.05). Patient age, tumor size, stage, lymphovascular invasion, and receiving chemotherapy and hormone therapy were independent risk factors for metastasis and recurrence. Only the luminal B1 subtype (hazard ratio: 5.653, 95% CI: 1.166-27.409, p < 0.05) and stage (hazard ratio: 3.400, 95% CI: 1.512-7.649, p < 0.05) were identified as independent risk factors for post metastatic mortality.

CONCLUSION

Luminal B1 breast cancer has aggressive tumor biology compared with luminal A breast cancer in the follow-up period. However, there was no significant difference in the disease relapse pattern between the groups.

International Validation of the SORG Machine-learning Algorithm for Predicting the Survival of Patients with Extremity Metastases Undergoing Surgical Treatment

BACKGROUND

The Skeletal Oncology Research Group machine-learning algorithms (SORG-MLAs) estimate 90-day and 1-year survival in patients with long-bone metastases undergoing surgical treatment and have demonstrated good discriminatory ability on internal validation. However, the performance of a prediction model could potentially vary by race or region, and the SORG-MLA must be externally validated in an Asian cohort. Furthermore, the authors of the original developmental study did not consider the Eastern Cooperative Oncology Group (ECOG) performance status, a survival prognosticator repeatedly validated in other studies, in their algorithms because of missing data.

QUESTIONS/PURPOSES

(1) Is the SORG-MLA generalizable to Taiwanese patients for predicting 90-day and 1-year mortality? (2) Is the ECOG score an independent factor associated with 90-day and 1-year mortality while controlling for SORG-MLA predictions?

METHODS

All 356 patients who underwent surgery for long-bone metastases between 2014 and 2019 at one tertiary care center in Taiwan were included. Ninety-eight percent (349 of 356) of patients were of Han Chinese descent. The median (range) patient age was 61 years (25 to 95), 52% (184 of 356) were women, and the median BMI was 23 kg/m2 (13 to 39 kg/m2). The most common primary tumors were lung cancer (33% [116 of 356]) and breast cancer (16% [58 of 356]). Fifty-five percent (195 of 356) of patients presented with a complete pathologic fracture. Intramedullary nailing was the most commonly performed type of surgery (59% [210 of 356]), followed by plate screw fixation (23% [81 of 356]) and endoprosthetic reconstruction (18% [65 of 356]). Six patients were lost to follow-up within 90 days; 30 were lost to follow-up within 1 year. Eighty-five percent (301 of 356) of patients were followed until death or for at least 2 years. Survival was 82% (287 of 350) at 90 days and 49% (159 of 326) at 1 year. The model's performance metrics included discrimination (concordance index [c-index]), calibration (intercept and slope), and Brier score. In general, a c-index of 0.5 indicates random guess and a c-index of 0.8 denotes excellent discrimination. Calibration refers to the agreement between the predicted outcomes and the actual outcomes, with a perfect calibration having an intercept of 0 and a slope of 1. The Brier score of a prediction model must be compared with and ideally should be smaller than the score of the null model. A decision curve analysis was then performed for the 90-day and 1-year prediction models to evaluate their net benefit across a range of different threshold probabilities. A multivariate logistic regression analysis was used to evaluate whether the ECOG score was an independent prognosticator while controlling for the SORG-MLA's predictions. We did not perform retraining/recalibration because we were not trying to update the SORG-MLA algorithm in this study.

RESULTS

The SORG-MLA had good discriminatory ability at both timepoints, with a c-index of 0.80 (95% confidence interval 0.74 to 0.86) for 90-day survival prediction and a c-index of 0.84 (95% CI 0.80 to 0.89) for 1-year survival prediction. However, the calibration analysis showed that the SORG-MLAs tended to underestimate Taiwanese patients' survival (90-day survival prediction: calibration intercept 0.78 [95% CI 0.46 to 1.10], calibration slope 0.74 [95% CI 0.53 to 0.96]; 1-year survival prediction: calibration intercept 0.75 [95% CI 0.49 to 1.00], calibration slope 1.22 [95% CI 0.95 to 1.49]). The Brier score of the 90-day and 1-year SORG-MLA prediction models was lower than their respective null model (0.12 versus 0.16 for 90-day prediction; 0.16 versus 0.25 for 1-year prediction), indicating good overall performance of SORG-MLAs at these two timepoints. Decision curve analysis showed SORG-MLAs provided net benefits when threshold probabilities ranged from 0.40 to 0.95 for 90-day survival prediction and from 0.15 to 1.0 for 1-year prediction. The ECOG score was an independent factor associated with 90-day mortality (odds ratio 1.94 [95% CI 1.01 to 3.73]) but not 1-year mortality (OR 1.07 [95% CI 0.53 to 2.17]) after controlling for SORG-MLA predictions for 90-day and 1-year survival, respectively.

CONCLUSION

SORG-MLAs retained good discriminatory ability in Taiwanese patients with long-bone metastases, although their actual survival time was slightly underestimated. More international validation and incremental value studies that address factors such as the ECOG score are warranted to refine the algorithms, which can be freely accessed online at https://sorg-apps.shinyapps.io/extremitymetssurvival/.

LEVEL OF EVIDENCE

Level III, therapeutic study.

CX3CR1 Acts as a Protective Biomarker in the Tumor Microenvironment of Colorectal Cancer

The tumor microenvironment (TME) plays an important role in the pathogenesis of many cancers. We aimed to screen the TME-related hub genes of colorectal adenoma (CRAD) and identify possible prognostic biomarkers. The gene expression profiles and clinical data of 464 CRAD patients in The Cancer Genome Atlas (TCGA) database were downloaded. The Estimation of STromal and Immune cells in MAlignant Tumours using Expression data (ESTIMATE) algorithm was performed to calculate the ImmuneScore, StromalScore, and EstimateScore. Thereafter, differentially expressed genes (DEGs) were screened. Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway, and protein–protein interaction (PPI) analysis were performed to explore the roles of DEGs. Furthermore, univariate and multivariate Cox analyses were accomplished to identify independent prognostic factors of CRAD. CX3CR1 was selected as a hub gene, and the expression was confirmed in colorectal cancer (CRC) patients and cell lines. The correlations between CX3CR1 and tumor-infiltrating immune cells were estimated by Tumor IMmune Estimation Resource database (TIMER) and CIBERSORT analysis. Besides, we investigated the effects of coculture with THP-1-derived macrophages with HCT8 cells with low CX3CR1 expression on immune marker expression, cell viability, and migration. There were significant differences in the ImmuneScore and EstimateScore among different stages. Patients with low scores presented significantly lower lifetimes than those in the high-score group. Moreover, we recognized 1,578 intersection genes in ImmuneScore and StromalScore, and these genes were mainly enriched in numerous immune-related biological processes. CX3CR1 was found to be associated with immune cell infiltration levels, immune marker expression, and macrophage polarization. Simultaneous silencing of CX3CR1 and coculture with THP-1 cells further regulated macrophage polarization and promoted the cell proliferation and migration of CRC cells. CX3CR1 was decreased in CRAD tissues and cell lines and was related to T and N stages, tumor differentiation, and prognosis. Our results suggest that CX3CR1 contributes to the recruitment and regulation of immune-infiltrating cells and macrophage polarization in CRC and TAM-induced CRC progression. CX3CR1 may act as a prognostic biomarker in CRC.

Development of a 15-Gene Signature Model as a Prognostic Tool in Sex Hormone-Dependent Cancers

Sex hormone dependence is associated with tumor progression and prognosis. Here, we explored the molecular basis of luminal A-like phenotype in sex hormone-dependent cancers. RNA-sequencing data from 8 cancer types were obtained from The Cancer Genome Atlas (TCGA). We investigated the enrichment function of differentially expressed genes (DEGs) in luminal A breast cancer (BRCA). Weighted coexpression network analysis (WGCNA) was used to identify gene modules associated with the luminal A-like phenotype, and we calculated the module's preservation in 8 cancer types. Module hub genes screened using least absolute shrinkage and selection operator (LASSO) were used to construct a gene signature model for the luminal A-like phenotype, and we assessed the model's relationship with prognosis, enriched pathways, and immune infiltration using bioinformatics approaches. Compared to other BRCA subtypes, the enrichment functions of upregulated genes in luminal A BRCA were related to hormone biological processes and receptor activity, and the downregulated genes were associated with the cell cycle and nuclear division. A gene module significantly associated with luminal A BRCA was shared by uterine corpus endometrial carcinoma (UCEC), leading to a similar phenotype. Fifteen hub genes were used to construct a gene signature model for the assessment of the luminal A-like phenotype, and the corrected C-statistics and Brier scores were 0.986 and 0.023, respectively. Calibration plots showed good performance, and decision curve analysis indicated a high net benefit of the model. The 15-gene signature model was associated with better overall survival in BRCA and UCEC and was characterized by downregulation of DNA replication, cell cycle and activated CD4 T cells. In conclusion, our study elucidated that BRCA and UCEC share a similar sex hormone-dependent phenotype and constructed a 15-gene signature model for use as a prognostic tool to quantify the probability of the phenotype.

Breast and gut microbiome in health and cancer

The microbiota plays essential roles in health and disease, in both the intestine and the extra-intestine. Dysbiosis of the gut microbiota causes dysfunction in the intestine, which leads to inflammatory, immune, and infectious diseases. Dysbiosis is also associated with diseases beyond the intestine via microbial translocation or metabolisms. The in situ breast microbiome, which may be sourced from the gut through lactation and sexual contact, could be altered and cause breast diseases. In this review, we summarize the recent progress in understanding the interactions among the gut microbiome, breast microbiome, and breast diseases. We discuss the intestinal microbiota, microbial metabolites, and roles of microbiota in immune system. We emphasize the novel roles and mechanisms of the microbiome (both in situ and gastrointestinal sourced) and bacterial products in the development and progression of breast cancer. The intestinal microbial translocation suggests that the gut microbiome is translocated to the skin and subsequently to the breast tissue. The gut bacterial translocation is also due to the increased intestinal permeability. The breast and intestinal microbiota are important factors in maintaining healthy breasts. Micronutrition queuine (Q) is derived from a de novo synthesized metabolite in bacteria. All human cells use queuine and incorporate it into the wobble anticodon position of specific transfer RNAs. We have demonstrated that Q modification regulates genes critical in tight junctions and migration in human breast cancer cells and a breast tumor model. We further discuss the challenges and future perspectives that can move the field forward for prevention, diagnosis, and treatment of breast diseases.

Body mass index and type 2 diabetes and breast cancer survival: a Mendelian randomization study

The causal relationship between body mass index (BMI) and type 2 diabetes (T2D) and breast cancer prognosis is still ambiguous. The aim of this study was to investigate the prognostic effect of BMI and T2D on breast cancer disease-free survival (DFS) among Asian individuals. In this two-sample Mendelian randomization (MR) study, the instrumental variables (IVs) were identified using a genome-wide association study (GWAS) among 24,000 participants in the Taiwan Biobank. Importantly, the validity of these IVs was confirmed with a previous large-scale GWAS (Biobank Japan Project, BBJ). In this study, we found that a genetic predisposition toward higher BMI (as indicated by BMI IVs, F = 86.88) was associated with poor breast cancer DFS (hazard ratio [HR] = 6.11; P < 0.001). Furthermore, higher level of genetically predicted T2D (as indicated by T2D IVs) was associated with an increased risk of recurrence of and mortality from breast cancer (HR = 1.43; P < 0.001). Sensitivity analyses, including the weighted-median approach, MR-Egger regression, Radial regression and Mendelian randomization pleiotropy residual sum and outlier (MR-PRESSO) supported the consistency of our findings. Finally, the causal relationship between BMI and poor breast cancer prognosis was confirmed in a prospective cohort study. Our MR analyses demonstrated the causal relationship between the genetic prediction of elevated BMI and a greater risk of T2D with poor breast cancer prognosis. BMI and T2D have important clinical implications and may be used as prognostic indicators of breast cancer.

High prevalence of APOA1/C3/A4/A5 alterations in luminal breast cancers among young women in East Asia

In East Asia, the breast cancer incidence rate among women aged <50 years has rapidly increased. Emerging tumors are distinctly characterized by a high prevalence of estrogen receptor (ER)-positive/human epidermal growth factor receptor (HER2)-negative cancer. In the present study, we identified unique genetic alterations in these emerging tumors. We analyzed gene copy number variations (CNVs) in breast tumors from 120 Taiwanese patients, and obtained public datasets of CNV and gene expression (GE). The data regarding CNV and GE were separately compared between East Asian and Western patients, and the overlapping genes identified in the comparisons were explored to identify the gene-gene interaction networks. In the age <50 years/ER + /HER2- subgroup, tumors of East Asian patients exhibited a higher frequency of copy number loss in APOA1/C3/A4/A5, a lipid-metabolizing gene cluster (33 vs. 10%, P < .001) and lower APOA1/C3/A4/A5 expressions than tumors of Western patients. These copy number loss related- and GE-related results were validated in another Taiwanese cohort and in two GE datasets, respectively. The copy number loss was significantly associated with poor survival among Western patients, but not among East Asian patients. Lower APOA1, APOC3, and APOA5 expressions were associated with higher ESTIMATE immune scores, indicating an abundance of tumor-infiltrating immune cells. In conclusion, APOA1/C3/A4/A5 copy number loss was more prevalent in luminal breast tumors among East Asian women aged <50 years, and its immunomodulatory effect on the tumor microenvironment possibly plays various roles in the tumor biology of East Asian patients.

Bioinformatics-Based Identification of Tumor Microenvironment-Related Prognostic Genes in Pancreatic Cancer

Objective

Growing evidence has highlighted that the immune and stromal cells that infiltrate in pancreatic cancer microenvironment significantly influence tumor progression. However, reliable microenvironment-related prognostic gene signatures are yet to be established. The present study aimed to elucidate tumor microenvironment-related prognostic genes in pancreatic cancer.

Methods

We applied the ESTIMATE algorithm to categorize patients with pancreatic cancer from TCGA dataset into high and low immune/stromal score groups and determined their differentially expressed genes. Then, univariate and LASSO Cox regression was performed to identify overall survival-related differentially expressed genes (DEGs). And multivariate Cox regression analysis was used to screen independent prognostic genes and construct a risk score model. Finally, the performance of the risk score model was evaluated by Kaplan-Meier curve, time-dependent receiver operating characteristic and Harrell’s concordance index.

Results

The overall survival analysis demonstrated that high immune/stromal score groups were closely associated with poor prognosis. The multivariate Cox regression analysis indicated that the signatures of four genes, including TRPC7, CXCL10, CUX2, and COL2A1, were independent prognostic factors. Subsequently, the risk prediction model constructed by those genes was superior to AJCC staging as evaluated by time-dependent receiver operating characteristic and Harrell’s concordance index, and both KRAS and TP53 mutations were closely associated with high risk scores. In addition, CXCL10 was predominantly expressed by tumor associated macrophages and its receptor CXCR3 was highly expressed in T cells at the single-cell level.

Conclusions

This study comprehensively investigated the tumor microenvironment and verified immune/stromal-related biomarkers for pancreatic cancer.

Immune marker signature helps to predict survival in uveal melanoma

The detailed molecular function of tumor microenvironment (TEM) in uveal melanoma (UVM) remains unclear. This study generated the immune index and the stromal index scores by ESTIMATE algorithm based on RNA-sequencing data with 80 UVM patients. There was no correlation between the immune stromal index and clinical parameters. The differentially expressed genes related to the immune stromal index were calculated and were described by functional annotations and protein-protein interaction network diagrams. After univariate and multivariate Cox regression analyses, there were four genes (HLA-J, MMP12, HES6, and ADAMDEC1) with significant prognostic significance. The prognostic model was constructed using these four characteristic genes, and the KM curve and tROC curve were described to show that the model had a better ability to predict survival outcomes and prognosis. The verification results in GSE62075 showed that HLA-J and HES6 were expressed differently in the cancer group than in the non-cancer group. This study indicates that the risk signature based on the immune index can be used as an indicator to evaluate the prognosis of patients with UVM.

Mining immune-related genes with prognostic value in the tumor microenvironment of breast invasive ductal carcinoma

ABSTRACT

The tumor microenvironment (TME) plays an important role in the development of breast cancer. Due to limitations in experimental conditions, the molecular mechanism of TME in breast cancer has not yet been elucidated. With the development of bioinformatics, the study of TME has become convenient and reliable.Gene expression and clinical feature data were downloaded from The Cancer Genome Atlas database and the Molecular Taxonomy of Breast Cancer International Consortium database. Immune scores and stromal scores were calculated using the Estimation of Stromal and Immune Cells in Malignant Tumor Tissues Using Expression Data algorithm. The interaction of genes was examined with protein-protein interaction and co-expression analysis. The function of genes was analyzed by gene ontology enrichment analysis, Kyoto Encyclopedia of Genes and Genomes analysis and gene set enrichment analysis. The clinical significance of genes was assessed with Kaplan-Meier analysis and univariate/multivariate Cox regression analysis.Our results showed that the immune scores and stromal scores of breast invasive ductal carcinoma (IDC) were significantly lower than those of invasive lobular carcinoma. The immune scores were significantly related to overall survival of breast IDC patients and both the immune and stromal scores were significantly related to clinical features of these patients. According to the level of immune/stromal scores, 179 common differentially expressed genes and 5 hub genes with prognostic value were identified. In addition, the clinical significance of the hub genes was validated with data from the molecular taxonomy of breast cancer international consortium database, and gene set enrichment analysis analysis showed that these hub genes were mainly enriched in signaling pathways of the immune system and breast cancer.We identified five immune-related hub genes with prognostic value in the TME of breast IDC, which may partly determine the prognosis of breast cancer and provide some direction for development of targeted treatments in the future.

Development of a Prognostic Model for Ovarian Cancer Patients Based on Novel Immune Microenvironment Related Genes

Ovarian cancer (OV) has become the most lethal gynecological cancer. However, its treatment methods and staging system are far from ideal. In the present study, taking the advantage of large-scale public cohorts, we extracted a list of immune-related prognostic genes that differentially expressed in tumor and normal ovarian tissues. Importantly, an individualized immune-related gene based prognostic model (IPM) for OV patients were developed. Furthermore, we validated our IPM in Gene Expression Omnibus (GEO) repository and compared the immune landscape and pathways between high-risk and low-risk groups. The results of our study can serve as an important model to identify the immune subset of patients and has potential for use in immune therapeutic selection and patient management.

Immune Cell Infiltration-Based Characterization of Triple-Negative Breast Cancer Predicts Prognosis and Chemotherapy Response Markers

Breast cancer represents the number one cause of cancer-associated mortality globally. The most aggressive molecular subtype is triple negative breast cancer (TNBC), of which limited therapeutic options are available. It is well known that breast cancer prognosis and tumor sensitivity toward immunotherapy are dictated by the tumor microenvironment. Breast cancer gene expression profiles were extracted from the METABRIC dataset and two TNBC clusters displaying unique immune features were identified. Activated immune cells formed a large proportion of cells in the high infiltration cluster, which correlated to a good prognosis. Differentially expressed genes (DEGs) extracted between two heterogeneous subtypes were used to further explore the underlying immune mechanism and to identify prognostic biomarkers. Functional enrichment analysis revealed that the DEGs were predominately related to some processes involved in activation and regulation of innate immune signaling. Using network analysis, we identified two modules in which genes were selected for further prognostic investigation. Validation by independent datasets revealed that CXCL9 and CXCL13 were good prognostic biomarkers for TNBC. We also performed comparisons between the above two genes and immune markers (CYT, APM, TILs, and TIS), as well as cell checkpoint marker expressions, and found a statistically significant correlation between them in both METABRIC and TCGA datasets. The potential of CXCL9 and CXCL13 to predict chemotherapy sensitivity was also evaluated. We found that the CXCL9 and CXCL13 were good predictors for chemotherapy and their expressions were higher in chemotherapy-responsive patients in contrast to those who were not responsive. In brief, immune infiltrate characterization on TNBC revealed heterogeneous subtypes with unique immune features allowed for the identification of informative and reliable characteristics representative of the local immune tumor microenvironment and were potential candidates to guide the management of TNBC patients.

Prevalence of mammary Paget's disease in urban China in 2016

No national data have been available on descriptive epidemiology of mammary Paget's disease (MPD) in China. This population-based study aimed to estimate the prevalence of MPD and its pattens by sex, age and area in China. We conducted a population-based study using data in 2016 from China's Urban Employee Basic Medical Insurance and Urban Resident Basic Medical Insurance, covering approximately 0.43 billion residents. MPD cases were identified based on the diagnostic names and codes in claim data. A total of 825 patients of confirmed diagnosis of MPD were found during the study period. The prevalence of MPD in 2016 was 0.42 per 100,000 population (95% CI 0.19 to 0.73), with marked female predominance. The prevalence rates peaked at 40-59 years and ≥ 80 years in females and males, respectively. The prevalence rates varied among different regions, ranging from 0.06 (95% CI 0.00 to 0.23) in Northeast China to 1.21 (95% CI 0.07 to 3.72) in Northwest China. MPD showed marked female predominance in China. Chinese female patients were much younger, with lower prevalence than that in the United States. Obvious sex difference in the age pattern of MPD prevalence was also observed in China.

The molecular landscape of Asian breast cancers reveals clinically relevant population-specific differences

Molecular profiling of breast cancer has enabled the development of more robust molecular prognostic signatures and therapeutic options for breast cancer patients. However, non-Caucasian populations remain understudied. Here, we present the mutational, transcriptional, and copy number profiles of 560 Malaysian breast tumours and a comparative analysis of breast cancers arising in Asian and Caucasian women. Compared to breast tumours in Caucasian women, we show an increased prevalence of HER2-enriched molecular subtypes and higher prevalence of TP53 somatic mutations in ER+ Asian breast tumours. We also observe elevated immune scores in Asian breast tumours, suggesting potential clinical response to immune checkpoint inhibitors. Whilst HER2-subtype and enriched immune score are associated with improved survival, presence of TP53 somatic mutations is associated with poorer survival in ER+ tumours. Taken together, these population differences unveil opportunities to improve the understanding of this disease and lay the foundation for precision medicine in different populations.

M2 macrophage-induced lncRNA PCAT6 facilitates tumorigenesis and angiogenesis of triple-negative breast cancer through modulation of VEGFR2

As a common female malignancy, triple-negative breast cancer (TNBC) is the most malignant subtype of breast cancers (BC). This study further studied the role of long noncoding RNA (lncRNA) prostate cancer-associated transcript 6 (PCAT6) in TNBC. Functional assays, including EdU, wound healing, transwell, and immunofluorescence staining, revealed the effect of PCAT6 on cell proliferation, migration, and EMT process. The tube-formation assay disclosed the function of PCAT6 on angiogenesis. In vivo assays were also established to explore the impact of PCAT6 on tumor growth and microangiogenesis. The results revealed that PCAT6 boosted TNBC cell proliferation, migration, and angiogenesis both in vitro and in vivo. Then, this study unveiled that M2 macrophage secreted VEGF to stimulate the upregulation of PCAT6, thus promoting angiogenesis in TNBC. Next, through bioinformatics analysis and mechanism assays, we identified that PCAT6 positively regulated VEGFR2 expression via ceRNA pattern and then participated in VEGFR/AKT/mTOR signaling pathway to accelerate angiogenesis. Moreover, PCAT6 bound USP14, a deubiquitinase, to induce the deubiquitination of VEGFR2. On the whole, M2 macrophage-induced upregulation of PCAT6 facilitates TNBC tumorigenesis through modulation of VEGFR2 expression via ceRNA and deubiquitination patterns.

The Immune Landscape and Prognostic Immune Key Genes Potentially Involved in Modulating Synaptic Functions in Prostate Cancer

Background: Increasing evidence has indicated an association between differentially expressed genes (DEGs) in tumor-infiltrating immune cells (TIICs) and clinical outcome. The aim of this research is to investigate the influence of tumor microenvironment on the gene expression profile of TIICs and to identify their potential markers for modulating immune cell function in prostate cancer.

Methods: In our research, CIBERSORT algorithm was utilized to calculate the proportion of the TIICs in 164 tumor and 18 control samples from The Cancer Genome Atlas cohort. The differential expression analysis was conducted using R, and then the functional and the pathway enrichments of the DEGs were analyzed using Database for Annotation, Visualization, and Integrated Discovery, followed by integrated regulatory network analysis.

Results: As a result, nTreg, B cells, Th1, and DC cells were significantly increased, accompanied by largely decreased NK and NKT cells. The expressed immune-related gene correlation analysis showed that the signature gene expression extent of CD8 T cells was positively associated with CD4 memory activated T cells but negatively correlated with that of CD4 memory resting T cells. In addition, a total of 128 differentially expressed genes were identified. CytoHubba analysis obtained six hub genes, of which three prognostic-associated potential key molecules including CAV1, FLNA, and VCL were mainly involved in biological processes associated with the regulation of organic substance and synaptic connections.

Conclusions: This study provides a comprehensive understanding of the landscape of TIICs and the roles of the hub genes which may be valuable markers in prostate cancer diagnosis and immunotherapy.

Microenvironment-related gene TNFSF13B predicts poor prognosis in kidney renal clear cell carcinoma

Background

Kidney renal clear cell carcinoma (KIRC) affects the genitourinary system. Although treatment of KIRC in early stages can be highly successful, this type of cancer is difficult to detect until later stages of disease that are less easily treatable. Previous studies have focused on tumor cells, but have ignored the contributions of the tumor microenvironment.

Methods

We analyzed KIRC gene expression data from The Cancer Genome Atlas with the ESTIMATE algorithm to identify differentially expressed genes. Through protein-protein interaction network analysis, we identified clusters and picked genes from the clusters for further analysis. Differential expression, Kaplan-Meier, and univariate Cox analyses were used to select prognostic biomarkers. Gene Set Enrichment Analysis (GSEA) and Tumor Immune Estimation Resource (TIMER) analysis were used to explore the immune characteristics of these genes as biomarkers.

Results

Through the ESTIMATE algorithm and other system biology tools, TNFSF13B was identified as a prognostic biomarker. TNFSF13B is highly expressed in tumors, and high expression of TNFSF13B leads to poor prognosis. Further GSEA and TIMER analysis revealed that the expression of TNFSF13B was related to the immune signaling pathway and lymphocyte infiltration. Our findings strongly suggest that TNFSF13B may be a potential biomarker or target related to the tumor microenvironment for KIRC.

Immune system‑associated genes increase malignant progression and can be used to predict clinical outcome in patients with hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is one of the most malignant types of cancer, and is associated with high recurrence rates and a poor response to chemotherapy. Immune signatures in the microenvironment of HCC have not been well explored systematically. The aim of the present study was to identify prognostic immune signatures and build a nomogram for use in clinical evaluation. Using bioinformatics analysis, RNA‑seq data and overall survival (OS) information on 370 HCC cases from TCGA and 232 HCC cases from ICGC were analyzed. The differential expression of select immune genes, based on previously published studies, between HCC and adjacent tissue were analyzed using the limma package in R. Enrichment of pathways and gene ontology analysis was performed using clusterProfiler. Subsequently, univariate Cox regression analysis, Lasso penalty linear regression and multivariate Cox regression models were used to construct a model for immune risk score (IRS). The R packages, survival and survivalROC, were used to plot survival and the associated receiver operating characteristic curves. Infiltration of immune cells was calculated using Tumor IMmune Estimation Resource, with significance examined using a Pearson's correlation test. P<0.05 was considered significant. Based on the analysis, expression of 200 immune genes were upregulated and 47 immune genes were downregulated immune genes. In the multivariate Cox model, 5 genes (enhancer of zest homology 2, ferritin light chain, complement factor H related 3, isthmin 2, cyclin dependent kinase 5) were used to generate the IRS. By stratifying according to the median IRS, it was shown that patients with a high IRS had poor OS rates after 1, 2, 3 and 5 years, and this result was consistent across the testing, training and independent validation cohorts. Additionally, the IRS was correlated with the abundance of infiltrating immune cells. The nomogram built using IRS and clinical characteristics, was able to predict 1, 3 and 5 year OS with area under the curve values of >0.8. These results suggest that the model developed to calculate the IRS may be used to monitor the effectiveness of treatment strategies and for prognostic prediction.