Pyroptosis‐related molecular classification and immune microenvironment infiltration in breast cancer: A novel therapeutic target

The underlying role of pyroptosis in breast cancer (BC) remains unknown. Herein, we investigated the correlations of 33 pyroptosis‐related genes (PRGs) with immune checkpoints and immune cell infiltrations in BC patients based on The Cancer Genome Atlas cohort (n = 996) and Gene Expression Omnibus cohort (n = 3,262). Enrichment analysis revealed that these PRGs mainly functioned in pyroptosis, inflammasomes and regulation of autophagy pathway. Four prognostic independent PRGs (CASP9, TIRAP, GSDMC and IL18) were identified. Then, cluster 1/2 was recognized using consensus clustering for these four PRGs. Patients from cluster 1 had a favourable prognosis and diverse immune cell infiltrations. A nomogram was developed based on age, TNM stage, tumour subtype and pyroptosis score. Patients with the high‐risk group exhibited worse 5‐year OS, and the result was consistent in the external cohort. Additionally, high‐risk group patients were associated with downregulated immune checkpoint expression. Further analysis suggested that the high‐risk group patients were associated with a higher IC50 of paclitaxel, doxorubicin, cisplatin, methotrexate and vinorelbine. In summarizing, the pyroptosis score‐based nomogram might serve as an independent prognostic predictor and could guide medication for chemotherapy. Additionally, it may bring novel insight into the regulation of tumour immune microenvironment in BC and help to achieve precision immunotherapy.     

Effects of dietary macronutrients on the hepatic transcriptome and serum metabolome in mice

Dietary macronutrient composition influences both hepatic function and aging. Previous work suggested that longevity and hepatic gene expression levels were highly responsive to dietary protein, but almost unaffected by other macronutrients. In contrast, we found expression of 4005, 4232, and 4292 genes in the livers of mice were significantly associated with changes in dietary protein (5%–30%), fat (20%–60%), and carbohydrate (10%–75%), respectively. More genes in aging‐related pathways (notably mTOR, IGF‐1, and NF‐kappaB) had significant correlations with dietary fat intake than protein and carbohydrate intake, and the pattern of gene expression changes in relation to dietary fat intake was in the opposite direction to the effect of graded levels of caloric restriction consistent with dietary fat having a negative impact on aging. We found 732, 808, and 995 serum metabolites were significantly correlated with dietary protein (5%–30%), fat (8.3%–80%), and carbohydrate (10%–80%) contents, respectively. Metabolomics pathway analysis revealed sphingosine‐1‐phosphate signaling was the significantly affected pathway by dietary fat content which has also been identified as significant changed metabolic pathway in the previous caloric restriction study. Our results suggest dietary fat has major impact on aging‐related gene and metabolic pathways compared with other macronutrients.     

Long noncoding RNA CCAT2 reduces chemosensitivity to 5‐fluorouracil in breast cancer cells by activating the mTOR axis

Breast cancer (BC) is the most prevalent cancer in women and the second leading cause for cancer‐related death in women. LncRNA CCAT2 is involved in BC cell drug sensitivity. Drug resistance of BC cells after chemotherapy is the main obstacle to therapeutic effects. This study explored whether BC cell drug sensitivity to 5‐Fu was related to lncRNA CCAT2‐regulated mTOR pathway. Normal breast tissues and BC tissues before/after neoadjuvant chemotherapy were collected, and CCAT2 expression was detected by RT‐qPCR. Correlation between CCATA2 expression and neoadjuvant chemotherapy efficacy was analysed using the Kendall''s tau‐b correlation analysis. Normal breast epithelial cells and BC cell lines were cultured. BC cell lines were treated with 5‐Fu, and CCAT2 mRNA level in cells was detected. The 5‐Fu‐resistant MCF‐7/5‐Fu and MDA‐MB‐231/5‐Fu cells were treated with CCAT2 overexpression/knockdown or CCI‐779 (the mTOR pathway inhibitor). The mTOR pathway levels were detected. Expression of apoptosis‐related factors was identified. A subcutaneous xenograft model was carried out. High CCAT2 expression was detected in BC tissues and BC drug‐resistant cells after neoadjuvant chemotherapy, and a negative link was revealed between CCAT2 expression and efficacy of neoadjuvant chemotherapy. p‐mTOR/mTOR in 5‐Fu‐resistant BC cells with inhibited CCAT2 was decreased, while CCAT2 overexpression activated the mTOR pathway. IC50 value, proliferation, cells in S phase increased and apoptosis reduced after CCAT2 overexpression. After si‐CCAT2 or CCI‐779 treatment, the growth rate of transplanted tumours was inhibited, while promoted after CCAT2 overexpression. CCAT2 may reduce BC cell chemosensitivity to 5‐Fu by activating the mTOR pathway.     

Aging‐related cell type‐specific pathophysiologic immune responses that exacerbate disease severity in aged COVID‐19 patients

Coronavirus disease 2019 (COVID‐19) is especially severe in aged patients, defined as 65 years or older, for reasons that are currently unknown. To investigate the underlying basis for this vulnerability, we performed multimodal data analyses on immunity, inflammation, and COVID‐19 incidence and severity as a function of age. Our analysis leveraged age‐specific COVID‐19 mortality and laboratory testing from a large COVID‐19 registry, along with epidemiological data of ~3.4 million individuals, large‐scale deep immune cell profiling data, and single‐cell RNA‐sequencing data from aged COVID‐19 patients across diverse populations. We found that decreased lymphocyte count and elevated inflammatory markers (C‐reactive protein, D‐dimer, and neutrophil–lymphocyte ratio) are significantly associated with age‐specific COVID‐19 severities. We identified the reduced abundance of naïve CD8 T cells with decreased expression of antiviral defense genes (i.e., IFITM3 and TRIM22) in aged severe COVID‐19 patients. Older individuals with severe COVID‐19 displayed type I and II interferon deficiencies, which is correlated with SARS‐CoV‐2 viral load. Elevated expression of SARS‐CoV‐2 entry factors and reduced expression of antiviral defense genes (LY6E and IFNAR1) in the secretory cells are associated with critical COVID‐19 in aged individuals. Mechanistically, we identified strong TGF‐beta‐mediated immune–epithelial cell interactions (i.e., secretory‐non‐resident macrophages) in aged individuals with critical COVID‐19. Taken together, our findings point to immuno‐inflammatory factors that could be targeted therapeutically to reduce morbidity and mortality in aged COVID‐19 patients.     

Association between IGF‐1 levels ranges and all‐cause mortality: A meta‐analysis

The association between IGF‐1 levels and mortality in humans is complex with low levels being associated with both low and high mortality. The present meta‐analysis investigates this complex relationship between IGF‐1 and all‐cause mortality in prospective cohort studies. A systematic literature search was conducted in PubMed/MEDLINE, Scopus, and Cochrane Library up to September 2019. Published studies were eligible for the meta‐analysis if they had a prospective cohort design, a hazard ratio (HR) and 95% confidence interval (CI) for two or more categories of IGF‐1 and were conducted among adults. A random‐effects model with a restricted maximum likelihood heterogeneity variance estimator was used to find combined HRs for all‐cause mortality. Nineteen studies involving 30,876 participants were included. Meta‐analysis of the 19 eligible studies showed that with respect to the low IGF‐1 category, higher IGF‐1 was not associated with increased risk of all‐cause mortality (HR = 0.84, 95% CI = 0.68–1.05). Dose–response analysis revealed a U‐shaped relation between IGF‐1 and mortality HR. Pooled results comparing low vs. middle IGF‐1 showed a significant increase of all‐cause mortality (HR = 1.33, 95% CI = 1.14–1.57), as well as comparing high vs. middle IGF‐1 categories (HR = 1.23, 95% CI = 1.06–1.44). Finally, we provide data on the association between IGF‐1 levels and the intake of proteins, carbohydrates, certain vitamins/minerals, and specific foods. Both high and low levels of IGF‐1 increase mortality risk, with a specific 120–160 ng/ml range being associated with the lowest mortality. These findings can explain the apparent controversy related to the association between IGF‐1 levels and mortality.     

Combining bioinformatics techniques to explore the molecular mechanisms involved in pancreatic cancer metastasis and prognosis

This article aims to explore the underlying molecular mechanisms and prognosis‐related genes in pancreatic cancer metastasis. Pancreatic cancer metastasis‐related gene chip data were downloaded from GENE EXPRESSION OMNIBUS(GEO)database. Differentially expressed genes were screened after R‐package pre‐treatment. Functional annotations and related signalling pathways were analysed using DAVID software. GEPIA (Gene Expression Profiling Interactive Analysis) was used to perform prognostic analysis, and differential genes associated with prognosis were screened and validated using data from GEO. We screened 40 healthy patients, 40 primary pancreatic cancer and 40 metastatic pancreatic cancer patients, collected serum, designed primers and used qPCR to test the expression of prognosis‐related genes in each group. 109 differentially expressed genes related with pancreatic cancer metastasis were screened, of which 49 were up‐regulated and 60 were down‐regulated. Functional annotation and pathway analysis revealed differentially expressed genes were mainly concentrated in protein activation cascade, extracellular matrix construction, decomposition, etc In the biological process, it is mainly involved in signalling pathways such as PPAR, PI3K‐Akt and ECM receptor interaction. Prognostic analysis showed the expression levels of four genes were significantly correlated with the overall survival time of patients with pancreatic cancer, namely SCG5, CRYBA2, CPE and CHGB. qPCR experiments showed the expression of these four genes was decreased in both the primary pancreatic cancer group and the metastatic pancreatic cancer group, and the latter was more significantly reduced. Pancreatic cancer metastasis is closely related to the activation of PPAR pathway, PI3K‐Akt pathway and ECM receptor interaction. SCG5, CRYBA2, CPE and CHGB genes are associated with the prognosis of pancreatic cancer, and their low expression suggests a poor prognosis.     

Asprin‐loaded strontium‐containing α‐calcium sulphate hemihydrate/nano‐hydroxyapatite composite promotes regeneration of critical bone defects

Our laboratory originally synthesized strontium(Sr)‐containing α‐calcium sulphate hemihydrate/nano‐hydroxyapatite composite (Sr‐α‐CSH/n‐HA) and demonstrated its ability to repair critical bone defects. This study attempted to incorporate aspirin into it to produce a better bone graft material for critical bone defects. After 5% Sr‐α‐CSH was prepared by coprecipitation and hydrothermal methods, it was mixed with aspirin solution of different concentrations (50 μg/ml, 200 μg/ml, 800 μg/ml and 3200 μg/ml) at a fixed liquid‐solid ratio (0.54 v/w) to obtain aspirin‐loaded Sr‐α‐CSH/n‐HA composite. In vitro experiments were performed on the composite extracts. The tibial defects (3 mm*5 mm) in SD rat model were filled with the composite for 4 weeks and 12 weeks to evaluate its osteogenic capacity in vivo. Our results showed its capability of proliferation, migration and osteogenesis of BMSCs in vitro got improved. In vivo treatment with 800 μg/ml aspirin–loaded Sr‐α‐CSH/n‐HA composite led to significantly more new bone formation in the defects compared with Sr‐α‐CSH/n‐HA composite and significantly promoted the expression of osteogenic‐related genes and inhibited osteoclast activity. In general, our research suggests that aspirin‐loaded Sr‐α‐CSH/n‐HA composite may have a greater capacity of repairing tibial defects in SD rats than simple Sr‐α‐CSH/n‐HA composite.     

Comparative mRNA/micro-RNA co-expression network drives melanomagenesis by promoting epithelial–mesenchymal transition and vasculogenic mimicry signaling

This study aimed to identify a novel disease-associated differentially co-expressed mRNA-microRNA (miRNA) that is associated with vasculogenic mimicry (VM) and epithelial-to-mesenchymal transition (EMT) network at different stages of melanoma. By applying weighted gene co-expression network analysis, we constructed a VM+EMT biological network with the available microarray dataset downloaded from a public database. Quantitative real-time PCR, immunohistochemical staining, and CD31-periodic acid solution dual staining were performed to confirm the expression of genes associated with EMT and VM formation in subjects with malignant melanoma (n = 18) and primary melanoma (n = 13) and in healthy subjects (n = 10). Our findings suggested that phosphatidylserine-specific phospholipase A1-alpha (PLA1A) and dermokine (DMKN) genes function as oncogenes that trigger VM and EMT processes during melanomagenesis on interaction with miR-370, miR-563, and miR-770–5p. PLA1A and DMKN genes can be considered potential VM+EMT network-based diagnostic biomarkers for distinguishing between melanoma patients. We postulate that a network with altered PLA1A/miR-563 and DMNK/miR-770–5p/miR-370 may contribute to melanomagenesis by triggering the EMT signaling pathway and VM formation. This study provides a potentially valuable approach for the early diagnosis and prognosis of melanoma progression.     

Evaluating otter reintroduction outcomes using genetic spatial capture–recapture modified for dendritic networks

Monitoring the demographics and genetics of reintroduced populations is critical to evaluating reintroduction success, but species ecology and the landscapes that they inhabit often present challenges for accurate assessments. If suitable habitats are restricted to hierarchical dendritic networks, such as river systems, animal movements are typically constrained and may violate assumptions of methods commonly used to estimate demographic parameters. Using genetic detection data collected via fecal sampling at latrines, we demonstrate applicability of the spatial capture–recapture (SCR) network distance function for estimating the size and density of a recently reintroduced North American river otter (Lontra canadensis) population in the Upper Rio Grande River dendritic network in the southwestern United States, and we also evaluated the genetic outcomes of using a small founder group (n = 33 otters) for reintroduction. Estimated population density was 0.23–0.28 otter/km, or 1 otter/3.57–4.35 km, with weak evidence of density increasing with northerly latitude (β = 0.33). Estimated population size was 83–104 total otters in 359 km of riverine dendritic network, which corresponded to average annual exponential population growth of 1.12–1.15/year since reintroduction. Growth was ≥40% lower than most reintroduced river otter populations and strong evidence of a founder effect existed 8–10 years post‐reintroduction, including 13–21% genetic diversity loss, 84%–87% genetic effective population size decline, and rapid divergence from the source population (F _ST accumulation = 0.06/generation). Consequently, genetic restoration via translocation of additional otters from other populations may be necessary to mitigate deleterious genetic effects in this small, isolated population. Combined with non‐invasive genetic sampling, the SCR network distance approach is likely widely applicable to demogenetic assessments of both reintroduced and established populations of multiple mustelid species that inhabit aquatic dendritic networks, many of which are regionally or globally imperiled and may warrant reintroduction or augmentation efforts.     

Genetic association study of intron variants in the forkhead box protein P3 gene in Chinese patients diagnosed with cervical cancer

The aim of this study was to investigate the effects of forkhead box protein P3 (FOXP3) intron single nucleotide variants (SNVs) in high‐risk human papilloma virus (HR‐HPV) infection and cervical cancer (CC) malignant lesions. We performed FOXP3 genotyping in 350 patients with CC and 350 healthy controls using the ImLDR multiple single nucleotide polymorphism genotyping technology. The heterozygous mutation TC in rs2294021 decreased the risk of HR‐HPV infection and CC malignant lesions (TC vs. TT: OR = 0.71, 95% CI = 0.51–0.99); the dominant model TC+CC and allele C in rs2294021 decreased the risk of CC malignant lesions (TC+CC vs. TT: OR = 0.69, 95% CI = 0.50–0.95; C vs. T: OR = 0.78, 95% CI = 0.63–0.97). The heterozygous mutation GA, dominant model GA+AA and allele A in rs3761549 also decreased the risk of HR‐HPV infection and CC malignant lesions (GA vs. GG: OR = 0.70, 95% CI = 0.51–0.96; GA+AA vs. GG: OR = 0.69, 95% CI = 0.51–0.94; A vs. G: OR = 0.75, 95% CI = 0.58–0.96). Patients with CC and HR‐HPV infection carrying rs2294021 TC and rs3761549 GA had lower expression of FOXP3 protein. Haplotype analysis revealed that T‐C‐A decreased the risk of HR‐HPV infection. Furthermore, we found a significant association between immune cells infiltration and prognosis in patients with CC. Our findings demonstrated that rs2294021 and rs3761549 variants may protect against HR‐HPV and CC malignant lesions by downregulating FOXP3 and that FOXP3 was associated with immune cells infiltration, which affected the prognosis of CC.     

Single-nucleotide polymorphisms in one-carbon metabolism genes,Mediterranean diet and breast cancer risk: a case–control study in the Greek-Cypriot female population

Single-nucleotide polymorphisms (SNPs) within genes of the one-carbon metabolism pathway have been shown to interact with dietary folate intake to modify breast cancer (BC) risk. Our group has previously demonstrated that the Mediterranean dietary pattern, rich in beneficial one-carbon metabolism micronutrients, protects against BC in Greek-Cypriot women. We aimed to investigate whether SNPs in the MTHFR (rs1801133 and rs1801131) and MTR (rs1805087) genes modify the effect of the Mediterranean dietary pattern on BC risk. Dietary intake data were obtained using a 32-item food-frequency questionnaire. A dietary pattern specific to the Greek-Cypriot population, which closely resembles the Mediterranean diet, was derived using principal component analysis (PCA) and used as our dietary variable. Genotyping was performed on subjects from the MASTOS study, a case–control study of BC in Cyprus, using TaqMan assays. Adjusted odds ratios (ORs) were estimated using logistic regression analyses. High adherence to the PCA-derived Mediterranean dietary pattern further reduced BC risk with increasing number of variant MTHFR 677T alleles (OR_Q4vs.Q1 for 677TT = 0.37, 95 % CI 0.20–0.69, for 677 CT = 0.60, 95 % CI 0.42–0.86). Additionally, high adherence to the Mediterranean dietary pattern decreased BC risk in subjects with at least one MTR 2756A allele (OR_Q4vs.Q1 for 2756AA = 0.59, 95 % CI 0.43–0.81, for 2756AG = 0.59, 95 % CI 0.39–0.91) and in subjects with the MTHFR 1298CC genotype (OR_Q4vs.Q1 0.44, 95 % CI 0.30–0.65). Overall P-interaction values, however, were not statistically significant. Our study suggests that these MTHFR and MTR SNPs may act as effect modifiers, highlighting their biological significance in the association between Mediterranean diet, the one-carbon metabolism pathway and BC.

Electronic supplementary material

The online version of this article (doi:10.1007/s12263-015-0453-7) contains supplementary material, which is available to authorized users.     

CTNNAL1 participates in the regulation of mucus overproduction in HDM‐induced asthma mouse model through the YAP‐ROCK2 pathway

Our previous study indicated that adhesion molecule catenin alpha‐like 1(CTNNAL1) is downregulated in airway epithelial cells of asthma patients and asthma animal model but little is known about how the CTNNAL1 affects asthma pathogenesis. To reveal the direct relationship between asthma and CTNNAL1, CTNNAL1‐deficient mouse model in bronchopulmonary tissue was constructed by introducing CTNNAL1‐siRNA sequence using adeno‐associated virus (AAV) as vector. The mouse model of asthma was established by stimulation of house dust mite (HDM). After HDM‐challenged, there was marked airway inflammation, especially mucus hypersecretion in the CTNNAL1‐deficient mice. In addition, the CTNNAL1‐deficient mice exhibited an increase of lung IL‐4 and IL‐13 levels, as well as a significant increase of goblet cell hyperplasia and MUC5AC after HDM exposure. The expression of Yes‐associated protein (YAP), protein that interacted with α‐catenin, was downregulated after CTNNAL1 silencing and was upregulated due to its overexpression. In addition, the interaction between CTNNAL1 and YAP was confirmed by CO‐IP. Besides, inhibition of YAP could decrease the secretion of MUC5AC, IL‐4 and IL‐13 in CTNNAL1‐deficient 16HBE14o‐cells. Above results indicated us that CTNNAL1 regulated mucus hypersecretion through YAP pathway. In addition, the expression of ROCK2 increased when CTNNAL1 was silenced and decreased after YAP silencing, and inhibition of YAP decreased the expression of ROCK2 in CTNNAL1‐deficient HBE cells. Inhibition of ROCK2 decreased MUC5AC expression and IL‐13 secretion. In all, our study demonstrates that CTNNAL1 plays an important role in HDM‐induced asthma, mediating mucus secretion through the YAP‐ROCK2 pathway.     

Robust prognostic model based on immune infiltration‐related genes and clinical information in ovarian cancer

Immune infiltration of ovarian cancer (OV) is a critical factor in determining patient''s prognosis. Using data from TCGA and GTEx database combined with WGCNA and ESTIMATE methods, 46 genes related to OV occurrence and immune infiltration were identified. Lasso and multivariate Cox regression were applied to define a prognostic score (IGCI score) based on 3 immune genes and 3 types of clinical information. The IGCI score has been verified by K‐M curves, ROC curves and C‐index on test set. In test set, IGCI score (C‐index = 0.630) is significantly better than AJCC stage (C‐index = 0.541, p < 0.05) and CIN25 (C‐index = 0.571, p < 0.05). In addition, we identified key mutations to analyse prognosis of patients and the process related to immunity. Chi‐squared tests revealed that 6 mutations are significantly (p < 0.05) related to immune infiltration: BRCA1, ZNF462, VWF, RBAK, RB1 and ADGRV1. According to mutation survival analysis, we found 5 key mutations significantly related to patient prognosis (p < 0.05): CSMD3, FLG2, HMCN1, TOP2A and TRRAP. RB1 and CSMD3 mutations had small p‐value (p < 0.1) in both chi‐squared tests and survival analysis. The drug sensitivity analysis of key mutation showed when RB1 mutation occurs, the efficacy of six anti‐tumour drugs has changed significantly (p < 0.05).     

Comparison of two treatment strategies for irradiation of regional lymph nodes in patients with breast cancer: Lymph flow guided portals versus standard radiation fields

Aim and Background

Radiotherapy being an essential part of breast cancer treatment, we evaluate various radiotherapy strategies in patients with breast cancer.

Materials and methods

Lymph node (LN) scintigraphy was performed in 172 primary patients with BC. LN visualization started 30–360 min after intratumoral injection of 75–150 MBq of 99mTc-nanocolloids.Our standard recommendation for postoperative radiotherapy in patients with LN invasion by BC were as follows: for patients with external localization of tumour – breast + axillary (Ax) + sub-supraclavicular (SSCL) regions; with internal localization – all above + internal mammary nodes (IM). Proposed strategy of lymph flow guided radiotherapy is based on the assumption that only regions that contain ‘hot’ LNs must be included in a treatment volume.

Results

Among 110 patients with external localization of BC, Ax LNs were visualized in all cases and in 62 patients it was the only region with ‘hot’ LN. Twenty-three patients (20.9%) had drainage to Ax + SSCL, 12 (10.9%) – Ax + IM, 13 (11.8%) – Ax + SSCL + IM regions. After the visualization of lymph flow patterns, standard treatment volume was changed in 87/110 cases (79.1%): in 56.4%, reduced, in 22.7%, enlarged or changed.In 62 patients with tumours in internal quadrants, we revealed the following patterns of lymph-flow: only to the Ax region in 23 (37.1%); Ax + IM, 13 (21%); Ax + SSCL, 15 (24.2%); Ax + IM + ISSCL, 11 (17.7%) cases. After lymph-flow visualization, the standard irradiation volume was reduced in 53/62 (85.5%) cases.

Conclusion

Visualization of an individual lymph flow pattern from BC can be used for the optimization of standard fields used for irradiation of regional LNs.     

The deacetylation of Foxk2 by Sirt1 reduces chemosensitivity to cisplatin

In multiple types of cancer, decreased tumour cell apoptosis during chemotherapy is indicative of decreased chemosensitivity. Forkhead box K2 (FOXK2), which is essential for cell fate, regulates cancer cell apoptosis through several post‐translational modifications. However, FOXK2 acetylation has not been extensively studied. Here, we evaluated the effects of sirtiun 1 (SIRT1) on FOXK2 deacetylation. Our findings demonstrated that SIRT1 inhibition increased FOXK2‐induced chemosensitivity to cisplatin and that K223 in FOXK2 was acetylated. Furthermore, FOXK2 K223 deacetylation reduced chemosensitivity to cisplatin in vitro and in vivo. Mechanistically, FOXK2 was acetylated by the acetyltransferase cAMP response element binding protein and deacetylated by SIRT1. Furthermore, cisplatin attenuated the interaction between FOXK2 and SIRT1. Cisplatin or SIRT1 inhibition enhanced FOXK2 acetylation, thereby reducing the nuclear distribution of FOXK2. Additionally, FOXK2 K223 acetylation significantly affected the expression of cell cycle–related and apoptosis‐related genes in cisplatin‐stimulated cancer cells, and FOXK2 K223 hyperacetylation promoted mitotic catastrophe, which enhanced chemosensitivity to cisplatin. Overall, our results provided insights into the mechanisms of SIRT1‐mediated FOXK2 deacetylation, which was involved in chemosensitivity to cisplatin.