Germline variants associated with leukocyte genes predict tumor recurrence in breast cancer patients

Breast Tumor Metastasis and Its Microenvironment: It Takes Both Seed and Soil to Grow a Tumor and Target It for Treatment

Metastasis remains a major challenge in treating breast cancer. Breast tumors metastasize to organ-specific locations such as the brain, lungs, and bone, but why some organs are favored over others remains unclear. Breast tumors also show heterogeneity, plasticity, and distinct microenvironments. This contributes to treatment failure and relapse. The interaction of breast cancer cells with their metastatic microenvironment has led to the concept that primary breast cancer cells act as seeds, whereas the metastatic tissue microenvironment (TME) is the soil. Improving our understanding of this interaction could lead to better treatment strategies for metastatic breast cancer. Targeted treatments for different subtypes of breast cancers have improved overall patient survival, even with metastasis. However, these targeted treatments are based upon the biology of the primary tumor and often these patients' relapse, after therapy, with metastatic tumors. The advent of immunotherapy allowed the immune system to target metastatic tumors. Unfortunately, immunotherapy has not been as effective in metastatic breast cancer relative to other cancers with metastases, such as melanoma. This review will describe the heterogeneic nature of breast cancer cells and their microenvironments. The distinct properties of metastatic breast cancer cells and their microenvironments that allow interactions, especially in bone and brain metastasis, will also be described. Finally, we will review immunotherapy approaches to treat metastatic breast tumors and discuss future therapeutic approaches to improve treatments for metastatic breast cancer.

Germline modifiers of the tumor immune microenvironment implicate drivers of cancer risk and immunotherapy response

With the continued promise of immunotherapy for treating cancer, understanding how host genetics contributes to the tumor immune microenvironment (TIME) is essential to tailoring cancer screening and treatment strategies. Here, we study 1084 eQTLs affecting the TIME found through analysis of The Cancer Genome Atlas and literature curation. These TIME eQTLs are enriched in areas of active transcription, and associate with gene expression in specific immune cell subsets, such as macrophages and dendritic cells. Polygenic score models built with TIME eQTLs reproducibly stratify cancer risk, survival and immune checkpoint blockade (ICB) response across independent cohorts. To assess whether an eQTL-informed approach could reveal potential cancer immunotherapy targets, we inhibit CTSS, a gene implicated by cancer risk and ICB response-associated polygenic models; CTSS inhibition results in slowed tumor growth and extended survival in vivo. These results validate the potential of integrating germline variation and TIME characteristics for uncovering potential targets for immunotherapy.

Roles of Skeletal Muscle in Development: A Bioinformatics and Systems Biology Overview

The ability to assess various cellular events consequent to perturbations, such as genetic mutations, disease states and therapies, has been recently revolutionized by technological advances in multiple "omics" fields. The resulting deluge of information has enabled and necessitated the development of tools required to both process and interpret the data. While of tremendous value to basic researchers, the amount and complexity of the data has made it extremely difficult to manually draw inference and identify factors key to the study objectives. The challenges of data reduction and interpretation are being met by the development of increasingly complex tools that integrate disparate knowledge bases and synthesize coherent models based on current biological understanding. This chapter presents an example of how genomics data can be integrated with biological network analyses to gain further insight into the developmental consequences of genetic perturbations. State of the art methods for conducting similar studies are discussed along with modern methods used to analyze and interpret the data.

TGFBR1*6A as a modifier of breast cancer risk and progression: advances and future prospects

There is growing evidence that germline mutations in certain genes influence cancer susceptibility, tumor evolution, as well as clinical outcomes. Identification of a disease-causing genetic variant enables testing and diagnosis of at-risk individuals. For breast cancer, several genes such as BRCA1, BRCA2, PALB2, ATM, and CHEK2 act as high- to moderate-penetrance cancer susceptibility genes. Genotyping of these genes informs genetic risk assessment and counseling, as well as treatment and management decisions in the case of high-penetrance genes. TGFBR1*6A (rs11466445) is a common variant of the TGF-β receptor type I (TGFBR1) that has a global minor allelic frequency (MAF) of 0.051 according to the 1000 Genomes Project Consortium. It is emerging as a high frequency, low penetrance tumor susceptibility allele associated with increased cancer risk among several cancer types. The TGFBR1*6A allele has been associated with increased breast cancer risk in women, OR 1.15 (95% CI 1.01–1.31). Functionally, TGFBR1*6A promotes breast cancer cell proliferation, migration, and invasion through the regulation of the ERK pathway and Rho-GTP activation. This review discusses current findings on the genetic, functional, and mechanistic associations between TGFBR1*6A and breast cancer risk and proposes future directions as it relates to genetic association studies and mechanisms of action for tumor growth, metastasis, and immune suppression.

Comparative analysis of protein-protein interaction networks in metastatic breast cancer

Metastatic lesions leading causes of the majority of deaths in patients with the breast cancer. The present study aimed to provide a comprehensive analysis of the differentially expressed genes (DEGs) in the brain (MDA-MB-231 BrM2) and lung (MDA-MB-231 LM2) metastatic cell lines obtained from breast cancer patients compared with those who have primary breast cancer. We identified 981 and 662 DEGs for brain and lung metastasis, respectively. Protein-protein interaction (PPI) analysis revealed seven shared (PLCB1, FPR1, FPR2, CX3CL1, GABBR2, GPR37, and CXCR4) hub genes between brain and lung metastasis in breast cancer. Moreover, GNG2 and CXCL8, C3, and PTPN6 in the brain and SAA1 and CCR5 in lung metastasis were found as unique hub genes. Besides, five co-regulation of clusters via seven important co-expression genes (COL1A2, LUM, SPARC, THBS2, IL1B, CXCL8, THY1) were identified in the brain PPI network. Clusters screening followed by biological process (BP) function and pathway enrichment analysis for both metastatic cell lines showed that complement receptor signalling, acetylcholine receptor signalling, and gastric acid secretion pathways were common between these metastases, whereas other pathways were site-specific. According to our findings, there are a set of genes and functional pathways that mark and mediate breast cancer metastasis to the brain and lungs, which may enable us understand the molecular basis of breast cancer development in a deeper levele to the brain and lungs, which may help us gain a more complete understanding of the molecular underpinnings of breast cancer development.

Identification and analysis of genes associated with epithelial ovarian cancer by integrated bioinformatics methods

Background

Though considerable efforts have been made to improve the treatment of epithelial ovarian cancer (EOC), the prognosis of patients has remained poor. Identifying differentially expressed genes (DEGs) involved in EOC progression and exploiting them as novel biomarkers or therapeutic targets is of great value.

Methods

Overlapping DEGs were screened out from three independent gene expression omnibus (GEO) datasets and were subjected to Gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) pathway enrichment analyses. The protein-protein interactions (PPI) network of DEGs was constructed based on the STRING database. The expression of hub genes was validated in GEPIA and GEO. The relationship of hub genes expression with tumor stage and overall survival and progression-free survival of EOC patients was investigated using the cancer genome atlas data.

Results

A total of 306 DEGs were identified, including 265 up-regulated and 41 down-regulated. Through PPI network analysis, the top 20 genes were screened out, among which 4 hub genes, which were not researched in depth so far, were selected after literature retrieval, including CDC45, CDCA5, KIF4A, ESPL1. The four genes were up-regulated in EOC tissues compared with normal tissues, but their expression decreased gradually with the continuous progression of EOC. Survival curves illustrated that patients with a lower level of CDCA5 and ESPL1 had better overall survival and progression-free survival statistically.

Conclusion

Two hub genes, CDCA5 and ESPL1, identified as probably playing tumor-promotive roles, have great potential to be utilized as novel therapeutic targets for EOC treatment.

The Role of BRCA1/2-Mutated Tumor Microenvironment in Breast Cancer

Taking into account the factors of high incidence rate, prevalence and mortality, breast cancer represents a crucial social and economic burden. Most cases of breast cancer develop as a consequence of somatic mutations accumulating in mammary epithelial cells throughout lifetime and approximately 5-10% can be ascribed to monogenic predispositions. Even though the role of genetic predispositions in breast cancer is well described in the context of genetics, very little is known about the role of the microenvironment carrying the same aberrant cells impaired by the germline mutation in the breast cancer development and progression. Based on the clinical observations, carcinomas carrying mutations in hereditary tumor-suppressor genes involved in maintaining genome integrity such as BRCA1/2 have worse prognosis and aggressive behavior. One of the mechanisms clarifying the aggressive nature of BRCA-associated tumors implies alterations within the surrounding adipose tissue itself. The objective of this review is to look at the role of BRCA1/2 mutations in the context of breast tumor microenvironment and plausible mechanisms by which it contributes to the aggressive behavior of the tumor cells.

Plasma exosome-derived microRNAs expression profiling and bioinformatics analysis under cross-talk between increased low-density lipoprotein cholesterol level and ATP-sensitive potassium channels variant rs1799858

Background

Exosome-derived microRNAs (exo-miRs) as messengers play important roles, in the cross-talk between genetic [ATP-sensitive potassium channels (KATP) genetic variant rs1799858] and environmental [elevated serum low-density lipoprotein cholesterol (LDL-C) level] factors, but the plasma exo-miRs expression profile and its role in biological processes from genotype to phenotype remain unclear.

Methods

A total of 14 subjects with increased LDL-C serum levels (≥ 1.8 mmol/L) were enrolled in the study. The KATP rs1799858 was genotyped by the Sequenom MassARRAY system. The plasma exo-miRs expression profile was identified by next-generation sequencing.

Results

64 exo-miRs were significantly differentially expressed (DE), among which 44 exo-miRs were up-regulated and 20 exo-miRs were down-regulated in those subjects carrying T-allele (TT + CT) of rs1799858 compared to those carrying CC genotype. The top 20 up-regulated DE-exo-miRs were miR-378 family, miR-320 family, miR-208 family, miR-483-5p, miR-22-3p, miR-490-3p, miR-6515-5p, miR-31-5p, miR-210-3p, miR-17-3p, miR-6807-5p, miR-497-5p, miR-33a-5p, miR-3611 and miR-126-5p. The top 20 down-regulated DE-exo-miRs were let-7 family, miR-221/222 family, miR-619-5p, miR-6780a-5p, miR-641, miR-200a-5p, miR-581, miR-605-3p, miR-548ar-3p, miR-135a-3p, miR-451b, miR-509-3-5p, miR-4664-3p and miR-224-5p. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were subsequently implemented to identify the top 10 DE-exo-miRs related specific target genes and signaling pathways. Only 5 DE-exo-miRs were validated by qRT-PCR as follows: miR-31-5p, miR-378d, miR-619-5p, miR-320a-3p and let-7a-5p (all P < 0.05).

Conclusion

These results firstly indicated the plasma exo-miRs expression profile bridging the link between genotype (KATP rs1799858) and phenotype (higher LDL-C serum level), these 5 DE-exo-miRs may be potential target intermediates for molecular intervention points.

Development of cancer prognostic signature based on pan-cancer proteomics

Utilizing genomic data to predict cancer prognosis was insufficient. Proteomics can improve our understanding of the etiology and progression of cancer and improve the assessment of cancer prognosis. And the Clinical Proteomic Tumor Analysis Consortium (CPTAC) has generated extensive proteomics data of the vast majority of tumors. Based on CPTAC, we can perform a proteomic pan-carcinoma analysis. We collected the proteomics data and clinical features of cancer patients from CPTAC. Then, we screened 69 differentially expressed proteins (DEPs) with R software in five cancers: hepatocellular carcinoma (HCC), children’s brain tumor tissue consortium (CBTTC), clear cell renal cell carcinoma (CCRC), lung adenocarcinoma (LUAD) and uterine corpus endometrial carcinoma (UCEC). GO and KEGG analysis were performed to clarify the function of these proteins. We also identified their interactions. The DEPs-based prognostic model for predicting over survival was identified by least absolute shrinkage and selection operator (LASSO)-Cox regression model in training cohort. Then, we used the time-dependent receiver operating characteristics analysis to evaluate the ability of the prognostic model to predict overall survival and validated it in validation cohort. The results showed that the DEPs-based prognostic model could accurately and effectively predict the survival rate of most cancers.

A Multi-Gene Model Effectively Predicts the Overall Prognosis of Stomach Adenocarcinomas With Large Genetic Heterogeneity Using Somatic Mutation Features

Background

Stomach adenocarcinoma (STAD) is one of the most common malignancies worldwide with poor prognosis. It remains unclear whether the prognosis is associated with somatic gene mutations.

Methods

In this research, we collected two independent STAD cohorts with both genetic profiling and clinical follow-up data, systematically investigated the association between the prognosis and somatic mutations, and analyzed the influence of heterogeneity on the prognosis-genetics association.

Results

Typical association was identified between somatic mutations and overall prognosis for individual cohorts. In The Cancer Genome Atlas (TCGA) cohort, a list of 24 genes was also identified that tended to mutate within cases of the poorest prognosis. The association showed apparent heterogeneity between different cohorts, although common signatures could be identified. A machine-learning model was trained with 20 common genes that showed a similar mutation rate difference between prognostic groups in the two cohorts, and it classified the cases in each cohort into two groups with significantly different prognosis. The model outperformed both single-gene models and TNM-based staging system significantly.

Conclusion

The study made a systematic analysis on the association between STAD prognosis and somatic mutations, identified signature genes that showed mutation preference in different prognostic groups, and developed an effective multi-gene model that can effectively predict the overall prognosis of STAD in different cohorts.

Patient-specific multi-omics models and the application in personalized combination therapy

The rapid advancement of high-throughput technologies and sharp decrease in cost have opened up the possibility to generate large amount of multi-omics data on an individual basis. The development of high-throughput -omics, including genomics, epigenomics, transcriptomics, proteomics, metabolomics and microbiomics, enables the application of multi-omics technologies in the clinical settings. Combination therapy, defined as disease treatment with two or more drugs to achieve efficacy with lower doses or lower drug toxicity, is the basis for the care of diseases like cancer. Patient-specific multi-omics data integration can help the identification and development of combination therapies. In this review, we provide an overview of different -omics platforms, and discuss the methods for multi-omics, high-throughput, data integration, personalized combination therapy.

Construction and Analysis of a Long Non-Coding RNA-Associated Competing Endogenous RNA Network Identified Potential Prognostic Biomarkers in Luminal Breast Cancer

Purpose

To construct a competing endogenous RNA (ceRNA) topology network of RNA-seq data and micro RNA-seq (miRNA-seq) data to identify key prognostic long non-coding RNA (lncRNAs) in luminal breast cancer, and validate the results by human luminal breast cancer samples.

Materials and Methods

The RNA-seq data and miRNA-seq data of luminal A breast cancer in the The Cancer Genome Atlas (TCGA) database were downloaded and compared with those in the miRcode database to obtain lncRNA–miRNA relationship pairs. Final target genes were predicted by all three databases (miRDB, miRTarBase, and TargetScan), thereby obtaining the miRNA-messenger RNA (miRNA-mRNA) relationship pairs and a ceRNA topology network was constructed, then mRNA enrichment analysis, ceRNA topological and stability analysis, univariate and multivariate Cox regression analysis were performed. Overall survival (OS) was evaluated and the key prognostic RNAs were identified. The expression difference between normal and tumor, as well as the correlation of high expression in tumor with pathological parameters (Ki-67, Grade, tumor diameter) were validated by human breast cancer specimens.

Results

A ceRNA topology network was constructed and six lncRNAs were finally identified (The higher expression of PART1, IGF2.AS, WT1.AS, OIP5.AS1, and SLC25A5.AS1 was associated with poor prognosis while AL035706.1 was adverse) and the poor prognostic ones were higher expressed in tumor tissue and correlated with a higher Ki-67 (>10%), tumor grades (II, III) and tumor diameters (>1.5 cm). Using six lncRNAs, we constructed a prognostic model, which performed well for the classification of prognosis in the module.

Conclusion

We identified and verified six biomarkers (OS-predicting) in luminal breast cancer, which significantly enriched the prediction and potential targets of this subtype.

Methylation-Based Signatures for Gastroesophageal Tumor Classification

Contention exists within the field of oncology with regards to gastroesophageal junction (GEJ) tumors, as in the past, they have been classified as gastric cancer, esophageal cancer, or a combination of both. Misclassifications of GEJ tumors ultimately influence treatment options, which may be rendered ineffective if treating for the wrong cancer attributes. It has been suggested that misclassification rates were as high as 45%, which is greater than reported for junctional cancer occurrences. Here, we aimed to use the methylation profiles of GEJ tumors to improve classifications of GEJ tumors. Four cohorts of DNA methylation profiles, containing ~27,000 (27k) methylation sites per sample, were collected from the Gene Expression Omnibus and The Cancer Genome Atlas. Tumor samples were assigned into discovery (nEC = 185, nGC = 395; EC, esophageal cancer; GC gastric cancer) and validation (nEC = 179, nGC = 369) sets. The optimized Multi-Survival Screening (MSS) algorithm was used to identify methylation biomarkers capable of distinguishing GEJ tumors. Three methylation signatures were identified: They were associated with protein binding, gene expression, and cellular component organization cellular processes, and achieved precision and recall rates of 94.7% and 99.2%, 97.6% and 96.8%, and 96.8% and 97.6%, respectively, in the validation dataset. Interestingly, the methylation sites of the signatures were very close (i.e., 170-270 base pairs) to their downstream transcription start sites (TSSs), suggesting that the methylations near TSSs play much more important roles in tumorigenesis. Here we presented the first set of methylation signatures with a higher predictive power for characterizing gastroesophageal tumors. Thus, they could improve the diagnosis and treatment of gastroesophageal tumors.

MeSHHeading2vec: a new method for representing MeSH headings as vectors based on graph embedding algorithm

Effectively representing Medical Subject Headings (MeSH) headings (terms) such as disease and drug as discriminative vectors could greatly improve the performance of downstream computational prediction models. However, these terms are often abstract and difficult to quantify. In this paper, we converted the MeSH tree structure into a relationship network and applied several graph embedding algorithms on it to represent these terms. Specifically, the relationship network consisting of nodes (MeSH headings) and edges (relationships), which can be constructed by the tree num. Then, five graph embedding algorithms including DeepWalk, LINE, SDNE, LAP and HOPE were implemented on the relationship network to represent MeSH headings as vectors. In order to evaluate the performance of the proposed methods, we carried out the node classification and relationship prediction tasks. The results show that the MeSH headings characterized by graph embedding algorithms can not only be treated as an independent carrier for representation, but also can be utilized as additional information to enhance the representation ability of vectors. Thus, it can serve as an input and continue to play a significant role in any computational models related to disease, drug, microbe, etc. Besides, our method holds great hope to inspire relevant researchers to study the representation of terms in this network perspective.

A Computational Study of Potential miRNA-Disease Association Inference Based on Ensemble Learning and Kernel Ridge Regression

As increasing experimental studies have shown that microRNAs (miRNAs) are closely related to multiple biological processes and the prevention, diagnosis and treatment of human diseases, a growing number of researchers are focusing on the identification of associations between miRNAs and diseases. Identifying such associations purely via experiments is costly and demanding, which prompts researchers to develop computational methods to complement the experiments. In this paper, a novel prediction model named Ensemble of Kernel Ridge Regression based MiRNA-Disease Association prediction (EKRRMDA) was developed. EKRRMDA obtained features of miRNAs and diseases by integrating the disease semantic similarity, the miRNA functional similarity and the Gaussian interaction profile kernel similarity for diseases and miRNAs. Under the computational framework that utilized ensemble learning and feature dimensionality reduction, multiple base classifiers that combined two Kernel Ridge Regression classifiers from the miRNA side and disease side, respectively, were obtained based on random selection of features. Then average strategy for these base classifiers was adopted to obtain final association scores of miRNA-disease pairs. In the global and local leave-one-out cross validation, EKRRMDA attained the AUCs of 0.9314 and 0.8618, respectively. Moreover, the model’s average AUC with standard deviation in 5-fold cross validation was 0.9275 ± 0.0008. In addition, we implemented three different types of case studies on predicting miRNAs associated with five important diseases. As a result, there were 90% (Esophageal Neoplasms), 86% (Kidney Neoplasms), 86% (Lymphoma), 98% (Lung Neoplasms), and 96% (Breast Neoplasms) of the top 50 predicted miRNAs verified to have associations with these diseases.