RNA-Seq based genome-wide analysis reveals loss of inter-chromosomal regulation in breast cancer

Coordinated inflammation and immune response transcriptional regulation in breast cancer molecular subtypes

Breast cancer, characterized by its complexity and diversity, presents significant challenges in understanding its underlying biology. In this study, we employed gene co-expression network analysis to investigate the gene composition and functional patterns in breast cancer subtypes and normal breast tissue. Our objective was to elucidate the detailed immunological features distinguishing these tumors at the transcriptional level and to explore their implications for diagnosis and treatment. The analysis identified nine distinct gene module clusters, each representing unique transcriptional signatures within breast cancer subtypes and normal tissue. Interestingly, while some clusters exhibited high similarity in gene composition between normal tissue and certain subtypes, others showed lower similarity and shared traits. These clusters provided insights into the immune responses within breast cancer subtypes, revealing diverse immunological functions, including innate and adaptive immune responses. Our findings contribute to a deeper understanding of the molecular mechanisms underlying breast cancer subtypes and highlight their unique characteristics. The immunological signatures identified in this study hold potential implications for diagnostic and therapeutic strategies. Additionally, the network-based approach introduced herein presents a valuable framework for understanding the complexities of other diseases and elucidating their underlying biology.

The network structure of hematopoietic cancers

Hematopoietic cancers (HCs) are a heterogeneous group of malignancies that affect blood, bone marrow and lymphatic system. Here, by analyzing 1960 RNA-Seq samples from three independent datasets, we explored the co-expression landscape in HCs, by inferring gene co-expression networks (GCNs) with four cancer phenotypes (B and T-cell acute leukemia -BALL, TALL-, acute myeloid leukemia -AML-, and multiple myeloma -MM-) as well as non-cancer bone marrow. We characterized their structure (topological features) and function (enrichment analyses). We found that, as in other types of cancer, the highest co-expression interactions are intra-chromosomal, which is not the case for control GCNs. We also detected a highly co-expressed group of overexpressed pseudogenes in HC networks. The four GCNs present only a small fraction of common interactions, related to canonical functions, like immune response or erythrocyte differentiation. With this approach, we were able to reveal cancer-specific features useful for detection of disease manifestations.

Methylation-related genes involved in renal carcinoma progression

Renal carcinomas are a group of malignant tumors often originating in the cells lining the small tubes in the kidney responsible for filtering waste from the blood and urine production. Kidney tumors arise from the uncontrolled growth of cells in the kidneys and are responsible for a large share of global cancer-related morbidity and mortality. Understanding the molecular mechanisms driving renal carcinoma progression results crucial for the development of targeted therapies leading to an improvement of patient outcomes. Epigenetic mechanisms such as DNA methylation are known factors underlying the development of several cancer types. There is solid experimental evidence of relevant biological functions modulated by methylation-related genes, associated with the progression of different carcinomas. Those mechanisms can often be associated to different epigenetic marks, such as DNA methylation sites or chromatin conformation patterns. Currently, there is no definitive method to establish clear relations between genetic and epigenetic factors that influence the progression of cancer. Here, we developed a data-driven method to find methylation-related genes, so we could find relevant bonds between gene co-expression and methylation-wide-genome regulation patterns able to drive biological processes during the progression of clear cell renal carcinoma (ccRC). With this approach, we found out genes such as ITK oncogene that appear hypomethylated during all four stages of ccRC progression and are strongly involved in immune response functions. Also, we found out relevant tumor suppressor genes such as RAB25 hypermethylated, thus potentially avoiding repressed functions in the AKT signaling pathway during the evolution of ccRC. Our results have relevant implications to further understand some epigenetic-genetic-affected roles underlying the progression of renal cancer.

Drug repurposing for Basal breast cancer subpopulations using modular network signatures

Breast cancer is characterized as being a heterogeneous pathology with a broad phenotype variability. Breast cancer subtypes have been developed in order to capture some of this heterogeneity. Each of these breast cancer subtypes, in turns retains varied characteristic features impacting diagnostic, prognostic and therapeutics. Basal breast tumors, in particular have been challenging in these regards. Basal breast cancer is often more aggressive, of rapid evolution and no tailor-made targeted therapies are available yet to treat it. Arguably, epigenetic variability is behind some of these intricacies. It is possible to further classify basal breast tumor in groups based on their non-coding transcriptome and methylome profiles. It is expected that these groups will have differences in survival as well as in sensitivity to certain classes of drugs. With this in mind, we implemented a computational learning approach to infer different subpopulations of basal breast cancer (from TCGA multi-omic data) based on their epigenetic signatures. Such epigenomic signatures were associated with different survival profiles; we then identified their associated gene co-expression network structure, extracted a signature based on modules within these networks, and use these signatures to find and prioritize drugs (in the LINCS dataset) that may be used to target these types of cancer. In this way we are introducing the analytical workflow for an epigenomic signature-based drug repurposing structure.

CNVs in 8q24.3 do not influence gene co-expression in breast cancer subtypes

Gene co-expression networks are a useful tool in the study of interactions that have allowed the visualization and quantification of diverse phenomena, including the loss of co-expression over long distances in cancerous samples. This characteristic, which could be considered fundamental to cancer, has been widely reported in various types of tumors. Since copy number variations (CNVs) have previously been identified as causing multiple genetic diseases, and gene expression is linked to them, they have often been mentioned as a probable cause of loss of co-expression in cancerous networks. In order to carry out a comparative study of the validity of this statement, we took 477 protein-coding genes from chromosome 8, and the CNVs of 101 genes, also protein-coding, belonging to the 8q24.3 region, a cytoband that is particularly active in the appearance of breast cancer. We created CNVS-conditioned co-expression networks of each of the 101 genes in the 8q24.3 region using conditional mutual information. The study was carried out using the four molecular subtypes of breast cancer (Luminal A, Luminal B, Her2, and Basal), as well as a case corresponding to healthy samples. We observed that in all cancer cases, the measurement of the Kolmogorov-Smirnov statistic shows that there are no significant differences between one and other values of the CNVs for any case. Furthermore, the co-expression interactions are stronger in all cancer subtypes than in the control networks. However, the control network presents a homogeneously distributed set of co-expression interactions, while for cancer networks, the highest interactions are more confined to specific cytobands, in particular 8q24.3 and 8p21.3. With this approach, we demonstrate that despite copy number alterations in the 8q24 region being a common trait in breast cancer, the loss of long-distance co-expression in breast cancer is not determined by CNVs.

Oncogenic Role of miR-217 During Clear Cell Renal Carcinoma Progression

Clear cell renal carcinoma (ccRC) comprises a set of heterogeneous, fast-progressing pathologies with poor prognosis. Analyzing ccRC progression in terms of modifications at the molecular level may provide us with a broader understanding of the disease, paving the way for improved diagnostics and therapeutics. The role of micro-RNAs (miRs) in cancer by targeting both oncogenes and tumor suppressor genes is widely known. Despite this knowledge, the role of specific miRs and their targets in the progression of ccRC is still unknown. To evaluate the action of miRs and their target genes during ccRC progression, here we implemented a three-step method for constructing miR–gene co-expression networks for each progression stage of ccRC as well as for adjacent-normal renal tissue (NT). In the first step, we inferred all miR–gene co-expression interactions for each progression stage of ccRC and for NT. Afterwards, we filtered the whole miR–gene networks by differential gene and miR expression between successive stages: stage I with non-tumor, stage II with stage I, and so on. Finally, all miR–gene interactions whose relationships were inversely proportional (overexpressed miR and underexpressed genes and vice versa) were kept and removed otherwise. We found that miR-217 is differentially expressed in all contrasts; however, its targets were different depending on the ccRC stage. Furthermore, the target genes of miR-217 have a known role in cancer progression—for instance, in stage II network, GALNTL6 is overexpressed, and it is related to cell signaling, survival, and proliferation. In the stage III network, WNK2, a widely known tumor suppressor, is underexpressed. For the stage IV network, IGF2BP2, a post-transcriptional regulator of MYC and PTEN, is overexpressed. This data-driven network approach has allowed us to discover miRs that have different targets through ccRC progression, thus providing a method for searching possible stage-dependent therapeutic targets in this and other types of cancer.

The Breast Cancer Protein Co-Expression Landscape

Breast cancer is a complex phenotype (or better yet, several complex phenotypes) characterized by the interplay of a large number of cellular and biomolecular entities. Biological networks have been successfully used to capture some of the heterogeneity of intricate pathophenotypes, including cancer. Gene coexpression networks, in particular, have been used to study large-scale regulatory patterns. Ultimately, biological processes are carried out by proteins and their complexes. However, to date, most of the tumor profiling research has focused on the genomic and transcriptomic information. Here, we tried to expand this profiling through the analysis of open proteomic data via mutual information co-expression networks' analysis. We could observe that there are distinctive biological processes associated with communities of these networks and how some transcriptional co-expression phenomena are lost at the protein level. These kinds of data and network analyses are a broad resource to explore cellular behavior and cancer research.

Introduction to Genomic Network Reconstruction for Cancer Research

High-throughput genomic technologies have revolutionized the study of cancer. Current research in oncology is now limited more for the capacity of analyzing and interpreting data, rather than the availability of data itself. Integrative approaches to obtain functional information from data are at the core of the disciplines gathered under the systems biology banner. In this context, network models have been used to study cancer, from the identification of key molecules involved in the disease to the discovery of functional alterations associated with specific manifestations of the disease.In this chapter, we describe the state of the art of network reconstruction from genomic data, with an emphasis in gene expression experiments. We explore the strengths and limitations of correlation, Bayesian, and information theoretic approaches to network reconstruction. We present tools that leverage the flexibility of network science to gain a deeper understanding of cancer biology.

The Role of Copy Number Variants in Gene Co-Expression Patterns for Luminal B Breast Tumors

Gene co-expression networks have become a usual approach to integrate the vast amounts of information coming from gene expression studies in cancer cohorts. The reprogramming of the gene regulatory control and the molecular pathways depending on such control are central to the characterization of the disease, aiming to unveil the consequences for cancer prognosis and therapeutics. There is, however, a multitude of factors which have been associated with anomalous control of gene expression in cancer. In the particular case of co-expression patterns, we have previously documented a phenomenon of loss of long distance co-expression in several cancer types, including breast cancer. Of the many potential factors that may contribute to this phenomenology, copy number variants (CNVs) have been often discussed. However, no systematic assessment of the role that CNVs may play in shaping gene co-expression patterns in breast cancer has been performed to date. For this reason we have decided to develop such analysis. In this study, we focus on using probabilistic modeling techniques to evaluate to what extent CNVs affect the phenomenon of long/short range co-expression in Luminal B breast tumors. We analyzed the co-expression patterns in chromosome 8, since it is known to be affected by amplifications/deletions during cancer development. We found that the CNVs pattern in chromosome 8 of Luminal B network does not alter the co-expression patterns significantly, which means that the co-expression program in this cancer phenotype is not determined by CNV structure. Additionally, we found that region 8q24.3 is highly dense in interactions, as well as region p21.3. The most connected genes in this network belong to those cytobands and are associated with several manifestations of cancer in different tissues. Interestingly, among the most connected genes, we found MAF1 and POLR3D, which may constitute an axis of regulation of gene transcription, in particular for non-coding RNA species. We believe that by advancing on our knowledge of the molecular mechanisms behind gene regulation in cancer, we will be better equipped, not only to understand tumor biology, but also to broaden the scope of diagnostic, prognostic and therapeutic interventions to ultimately benefit oncologic patients.

Gene Co-Expression in Breast Cancer: A Matter of Distance

Gene regulatory and signaling phenomena are known to be relevant players underlying the establishment of cellular phenotypes. It is also known that such regulatory programs are disrupted in cancer, leading to the onset and development of malignant phenotypes. Gene co-expression matrices have allowed us to compare and analyze complex phenotypes such as breast cancer (BrCa) and their control counterparts. Global co-expression patterns have revealed, for instance, that the highest gene-gene co-expression interactions often occur between genes from the same chromosome (cis-), meanwhile inter-chromosome (trans-) interactions are scarce and have lower correlation values. Furthermore, strength of cis- correlations have been shown to decay with the chromosome distance of gene couples. Despite this loss of long-distance co-expression has been clearly identified, it has been observed only in a small fraction of the whole co-expression landscape, namely the most significant interactions. For that reason, an approach that takes into account the whole interaction set results appealing. In this work, we developed a hybrid method to analyze whole-chromosome Pearson correlation matrices for the four BrCa subtypes (Luminal A, Luminal B, HER2+ and Basal), as well as adjacent normal breast tissue derived matrices. We implemented a systematic method for clustering gene couples, by using eigenvalue spectral decomposition and the k–medoids algorithm, allowing us to determine a number of clusters without removing any interaction. With this method we compared, for each chromosome in the five phenotypes: a) Whether or not the gene-gene co-expression decays with the distance in the breast cancer subtypes b) the chromosome location of cis- clusters of gene couples, and c) whether or not the loss of long-distance co-expression is observed in the whole range of interactions. We found that in the correlation matrix for the control phenotype, positive and negative Pearson correlations deviate from a random null model independently of the distance between couples. Conversely, for all BrCa subtypes, in all chromosomes, positive correlations decay with distance, and negative correlations do not differ from the null model. We also found that BrCa clusters are distance-dependent, meanwhile for the control phenotype, chromosome location does not determine the clustering. To our knowledge, this is the first time that a dependence on distance is reported for gene clusters in breast cancer. Since this method uses the whole cis- interaction geneset, combination with other -omics approaches may provide further evidence to understand in a more integrative fashion, the mechanisms that disrupt gene regulation in cancer.

Modularity in Biological Networks

Network modeling, from the ecological to the molecular scale has become an essential tool for studying the structure, dynamics and complex behavior of living systems. Graph representations of the relationships between biological components open up a wide variety of methods for discovering the mechanistic and functional properties of biological systems. Many biological networks are organized into a modular structure, so methods to discover such modules are essential if we are to understand the biological system as a whole. However, most of the methods used in biology to this end, have a limited applicability, as they are very specific to the system they were developed for. Conversely, from the statistical physics and network science perspective, graph modularity has been theoretically studied and several methods of a very general nature have been developed. It is our perspective that in particular for the modularity detection problem, biology and theoretical physics/network science are less connected than they should. The central goal of this review is to provide the necessary background and present the most applicable and pertinent methods for community detection in a way that motivates their further usage in biological research.

Spatial variation in gene expression of Tasmanian devil facial tumors despite minimal host transcriptomic response to infection

BACKGROUND

Transmissible cancers lie at the intersection of oncology and infectious disease, two traditionally divergent fields for which gene expression studies are particularly useful for identifying the molecular basis of phenotypic variation. In oncology, transcriptomics studies, which characterize the expression of thousands of genes, have identified processes leading to heterogeneity in cancer phenotypes and individual prognoses. More generally, transcriptomics studies of infectious diseases characterize interactions between host, pathogen, and environment to better predict population-level outcomes. Tasmanian devils have been impacted dramatically by a transmissible cancer (devil facial tumor disease; DFTD) that has led to widespread population declines. Despite initial predictions of extinction, populations have persisted at low levels, due in part to heterogeneity in host responses, particularly between sexes. However, the processes underlying this variation remain unknown.

RESULTS

We sequenced transcriptomes from healthy and DFTD-infected devils, as well as DFTD tumors, to characterize host responses to DFTD infection, identify differing host-tumor molecular interactions between sexes, and investigate the extent to which tumor gene expression varies among host populations. We found minimal variation in gene expression of devil lip tissues, either with respect to DFTD infection status or sex. However, 4088 genes were differentially expressed in tumors among our sampling localities. Pathways that were up- or downregulated in DFTD tumors relative to normal tissues exhibited the same patterns of expression with greater intensity in tumors from localities that experienced DFTD for longer. No mRNA sequence variants were associated with expression variation.

CONCLUSIONS

Expression variation among localities may reflect morphological differences in tumors that alter ratios of normal-to-tumor cells within biopsies. Phenotypic variation in tumors may arise from environmental variation or differences in host immune response that were undetectable in lip biopsies, potentially reflecting variation in host-tumor coevolutionary relationships among sites that differ in the time since DFTD arrival.

Spatial genome architecture and the emergence of malignancy

Human chromosomes are large spatially and hierarchically structured entities, the integrity of which needs to be preserved throughout the lifespan of the cell and in conjunction with cell cycle progression. Preservation of chromosomal structure is important for proper deployment of cell type-specific gene expression programs. Thus, aberrations in the integrity and structure of chromosomes will predictably lead to disease, including cancer. Here, we provide an updated standpoint with respect to chromatin misfolding and the emergence of various cancer types. We discuss recent studies implicating the disruption of topologically associating domains, switching between active and inactive compartments, rewiring of promoter-enhancer interactions in malignancy as well as the effects of single nucleotide polymorphisms in non-coding regions involved in long-range regulatory interactions. In light of these findings, we argue that chromosome conformation studies may now also be useful for patient diagnosis and drug target discovery.

k-core genes underpin structural features of breast cancer

Gene co-expression networks (GCNs) have been developed as relevant analytical tools for the study of the gene expression patterns behind complex phenotypes. Determining the association between structure and function in GCNs is a current challenge in biomedical research. Several structural differences between GCNs of breast cancer and healthy phenotypes have been reported. In a previous study, using co-expression multilayer networks, we have shown that there are abrupt differences in the connectivity patterns of the GCN of basal-like breast cancer between top co-expressed gene-pairs and the remaining gene-pairs. Here, we compared the top-100,000 interactions networks for the four breast cancer phenotypes (Luminal-A, Luminal-B, Her2+ and Basal), in terms of structural properties. For this purpose, we used the graph-theoretical k-core of a network (maximal sub-network with nodes of degree at least k). We developed a comprehensive analysis of the network k-core ([Formula: see text]) structures in cancer, and its relationship with biological functions. We found that in the Top-100,000-edges networks, the majority of interactions in breast cancer networks are intra-chromosome, meanwhile inter-chromosome interactions serve as connecting bridges between clusters. Moreover, core genes in the healthy network are strongly associated with processes such as metabolism and cell cycle. In breast cancer, only the core of Luminal A is related to those processes, and genes in its core are over-expressed. The intersection of the core nodes in all subtypes of cancer is composed only by genes in the chr8q24.3 region. This region has been observed to be highly amplified in several cancers before, and its appearance in the intersection of the four breast cancer k-cores, may suggest that local co-expression is a conserved phenomenon in cancer. Considering the many intricacies associated with these phenomena and the vast amount of research in epigenomic regulation which is currently undergoing, there is a need for further research on the epigenomic effects on the structure and function of gene co-expression networks in cancer.

Luminal A Breast Cancer Co-expression Network: Structural and Functional Alterations

Luminal A is the most common breast cancer molecular subtype in women worldwide. These tumors have characteristic yet heterogeneous alterations at the genomic and transcriptomic level. Gene co-expression networks (GCNs) have contributed to better characterize the cancerous phenotype. We have previously shown an imbalance in the proportion of intra-chromosomal (cis-) over inter-chromosomal (trans-) interactions when comparing cancer and healthy tissue GCNs. In particular, for breast cancer molecular subtypes (Luminal A included), the majority of high co-expression interactions connect gene-pairs in the same chromosome, a phenomenon that we have called loss of trans- co-expression. Despite this phenomenon has been described, the functional implication of this specific network topology has not been studied yet. To understand the biological role that communities of co-expressed genes may have, we constructed GCNs for healthy and Luminal A phenotypes. Network modules were obtained based on their connectivity patterns and they were classified according to their chromosomal homophily (proportion of cis-/trans- interactions). A functional overrepresentation analysis was performed on communities in both networks to observe the significantly enriched processes for each community. We also investigated possible mechanisms for which the loss of trans- co-expression emerges in cancer GCN. To this end we evaluated transcription factor binding sites, CTCF binding sites, differential gene expression and copy number alterations (CNAs) in the cancer GCN. We found that trans- communities in Luminal A present more significantly enriched categories than cis- ones. Processes, such as angiogenesis, cell proliferation, or cell adhesion were found in trans- modules. The differential expression analysis showed that FOXM1, CENPA, and CIITA transcription factors, exert a major regulatory role on their communities by regulating expression of their target genes in other chromosomes. Finally, identification of CNAs, displayed a high enrichment of deletion peaks in cis- communities. With this approach, we demonstrate that network topology determine, to at certain extent, the function in Luminal A breast cancer network. Furthermore, several mechanisms seem to be acting together to avoid trans- co-expression. Since this phenomenon has been observed in other cancer tissues, a remaining question is whether the loss of long distance co-expression is a novel hallmark of cancer.

An Information Theoretical Multilayer Network Approach to Breast Cancer Transcriptional Regulation

Breast cancer is a complex, highly heterogeneous disease at multiple levels ranging from its genetic origins and molecular processes to clinical manifestations. This heterogeneity has given rise to the so-called intrinsic or molecular breast cancer subtypes. Aside from classification, these subtypes have set a basis for differential prognosis and treatment. Multiple regulatory mechanisms-involving a variety of biomolecular entities-suffer from alterations leading to the diseased phenotypes. Information theoretical approaches have been found to be useful in the description of these complex regulatory programs. In this work, we identified the interactions occurring between three main mechanisms of regulation of the gene expression program: transcription factor regulation, regulation via noncoding RNA, and epigenetic regulation through DNA methylation. Using data from The Cancer Genome Atlas, we inferred probabilistic multilayer networks, identifying key regulatory circuits able to (partially) explain the alterations that lead from a healthy phenotype to different manifestations of breast cancer, as captured by its molecular subtype classification. We also found some general trends in the topology of the multi-omic regulatory networks: Tumor subtype networks present longer shortest paths than their normal tissue counterpart; epigenomic regulation has frequently focused on genes enriched for certain biological processes; CpG methylation and miRNA interactions are often part of a regulatory core of conserved interactions. The use of probabilistic measures to infer information regarding theoretical-derived multilayer networks based on multi-omic high-throughput data is hence presented as a useful methodological approach to capture some of the molecular heterogeneity behind regulatory phenomena in breast cancer, and potentially other diseases.

Loss of Long Distance Co-Expression in Lung Cancer

Lung cancer is one of the deadliest, most aggressive cancers. Abrupt changes in gene expression represent an important challenge to understand and fight the disease. Gene co-expression networks (GCNs) have been widely used to study the genomic regulatory landscape of human cancer. Here, based on 1,143 RNA-Seq experiments from the TCGA collaboration, we constructed GCN for the most common types of lung tumors: adenocarcinoma (TAD) and squamous cells (TSCs) as well as their respective control networks (NAD and NSC). We compared the number of intra-chromosome (cis-) and inter-chromosome (trans-) co-expression interactions in normal and cancer GCNs. We compared the number of shared interactions between TAD and TSC, as well as in NAD and NSC, to observe which phenotypes were more alike. By means of an over-representation analysis, we associated network topology features with biological functions. We found that TAD and TSC present mostly cis- small disconnected components, whereas in control GCNs, both types have a giant trans- component. In both cancer networks, we observed cis- components in which genes not only belong to the same chromosome but to the same cytoband or to neighboring cytobands. This supports the hypothesis that in lung cancer, gene co-expression is constrained to small neighboring regions. Despite this loss of distant co-expression observed in TAD and TSC, there are some remaining trans- clusters. These clusters seem to play relevant roles in the carcinogenic processes. For instance, some clusters in TAD and TSC are associated with the immune system, response to virus, or control of gene expression. Additionally, other non-enriched trans- clusters are composed of one gene and several associated pseudo-genes, as in the case of the FTH1 gene. The appearance of those common trans- clusters reflects that the gene co-expression program in lung cancer conserves some aspects for cell maintenance. Unexpectedly, 0.48% of the edges are shared between control networks; conversely, 35% is shared between lung cancer GCNs, a 73-fold larger intersection. This suggests that in lung cancer a process of de-differentiation may be occurring. To further investigate the implications of the loss of distant co-expression, it will become necessary to broaden the investigation with other omic-based approaches. However, the present approach provides a basis for future work toward an integrative perspective of abnormal transcriptional regulatory programs in lung cancer.

Pathway-Based Drug-Repurposing Schemes in Cancer: The Role of Translational Bioinformatics

Cancer is a set of complex pathologies that has been recognized as a major public health problem worldwide for decades. A myriad of therapeutic strategies is indeed available. However, the wide variability in tumor physiology, response to therapy, added to multi-drug resistance poses enormous challenges in clinical oncology. The last years have witnessed a fast-paced development of novel experimental and translational approaches to therapeutics, that supplemented with computational and theoretical advances are opening promising avenues to cope with cancer defiances. At the core of these advances, there is a strong conceptual shift from gene-centric emphasis on driver mutations in specific oncogenes and tumor suppressors-let us call that the silver bullet approach to cancer therapeutics-to a systemic, semi-mechanistic approach based on pathway perturbations and global molecular and physiological regulatory patterns-we will call this the shrapnel approach. The silver bullet approach is still the best one to follow when clonal mutations in driver genes are present in the patient, and when there are targeted therapies to tackle those. Unfortunately, due to the heterogeneous nature of tumors this is not the common case. The wide molecular variability in the mutational level often is reduced to a much smaller set of pathway-based dysfunctions as evidenced by the well-known hallmarks of cancer. In such cases "shrapnel gunshots" may become more effective than "silver bullets". Here, we will briefly present both approaches and will abound on the discussion on the state of the art of pathway-based therapeutic designs from a translational bioinformatics and computational oncology perspective. Further development of these approaches depends on building collaborative, multidisciplinary teams to resort to the expertise of clinical oncologists, oncological surgeons, and molecular oncologists, but also of cancer cell biologists and pharmacologists, as well as bioinformaticians, computational biologists and data scientists. These teams will be capable of engaging on a cycle of analyzing high-throughput experiments, mining databases, researching on clinical data, validating the findings, and improving clinical outcomes for the benefits of the oncological patients.

Gene Expression and Co-expression Networks Are Strongly Altered Through Stages in Clear Cell Renal Carcinoma

Clear cell renal carcinoma (ccRC) is a highly heterogeneous and progressively malignant disease. Analyzing ccRC progression in terms of modifications at the molecular and genetic level may help us to develop a broader understanding of its patho-physiology and may give us a glimpse toward improved therapeutics. In this work, by using TCGA data, we studied the molecular progression of the four main ccRC stages (i, ii, iii, iv) in two different yet complementary approaches: (a) gene expression and (b) gene co-expression. For (a) we analyzed the differential gene expression between each stage and the control non-cancer group. We compared the progression molecular signature between stages, and observed those genes that change their expression patterns through progression stages. For (b) we constructed and analyzed co-expression networks for the four ccRC progression stages, as well as for the control phenotype, to observe whether and how the co-expression landscape changes with progression. We separated genomic interactions into intra-chromosome (cis-) and inter-chromosome (trans-). Finally, we intersected those networks and performed functional enrichment analysis. All calculations were made over different network sizes, from the top 100 edges to top 1,000,000. We show that differential expression is quite similar between ccRC progression stages. However, interestingly, two genes, namely SLC6A19 and PLG show a significant progressive decrease in their expression according to ccRC stage, meanwhile two other genes, SAA2-SAA4 and CXCL13 show progressive increase. Despite the high similarity between gene expression profiles, all networks are substantially different between them in terms of their topological features. Control network has a larger proportion of trans- interactions, meanwhile for any stage, the amount of cis- interactions is higher, independent of the network cut-off. The majority of interactions in any network are phenotype-specific. Only 189 interactions are shared between the five networks, and 533 edges are ccRC-specific, independent of the stage. The small resulting connected components in both cases are formed by genes with the same differential expression trend, and are associated with important biological processes, such as cell cycle or immune system, suggesting that activity of these categories follows the differential expression trend. With this approach we have shown that, even if the expression program is similar during ccRC progression, the co-expression programs strongly differ. More research is needed to understand the delicate interplay between expression and co-expression, but this is a first approach to enclose both approaches in an integrative view aimed at a deeper understanding in gene regulation in tumor evolution.

Gene Co-expression Is Distance-Dependent in Breast Cancer

Breast carcinomas are characterized by anomalous gene regulatory programs. As is well-known, gene expression programs are able to shape phenotypes. Hence, the understanding of gene co-expression may shed light on the underlying mechanisms behind the transcriptional regulatory programs affecting tumor development and evolution. For instance, in breast cancer, there is a clear loss of inter-chromosomal (trans-) co-expression, compared with healthy tissue. At the same time cis- (intra-chromosomal) interactions are favored in breast tumors. In order to have a deeper understanding of regulatory phenomena in cancer, here, we constructed Gene Co-expression Networks by using TCGA-derived RNA-seq whole-genome samples corresponding to the four breast cancer molecular subtypes, as well as healthy tissue. We quantify the cis-/trans- co-expression imbalance in all phenotypes. Additionally, we measured the association between co-expression and physical distance between genes, and characterized the ratio of intra/inter-cytoband interactions per phenotype. We confirmed loss of trans- co-expression in all molecular subtypes. We also observed that gene cis- co-expression decays abruptly with distance in all tumors in contrast with healthy tissue. We observed co-expressed gene hotspots, that tend to be connected at cytoband regions, and coincide accurately with already known copy number altered regions, such as Chr17q12, or Chr8q24.3 for all subtypes. Our methodology recovered different alterations already reported for specific breast cancer subtypes, showing how co-expression network approaches might help to capture distinct events that modify the cell regulatory program.

Transcriptomic insight into salinomycin mechanisms in breast cancer cell lines: synergistic effects with dasatinib and induction of estrogen receptor β

BACKGROUND

Tumors are heterogeneous in nature, composed of different cell populations with various mutations and/or phenotypes. Using a single drug to encounter cancer progression is generally ineffective. To improve the treatment outcome, multiple drugs of distinctive mechanisms but complementary anticancer activities (combination therapy) are often used to enhance antitumor efficacy and minimize the risk of acquiring drug resistance. We report here the synergistic effects of salinomycin (a polyether antibiotic) and dasatinib (a Src kinase inhibitor).

METHODS

Functionally, both drugs induce cell cycle arrest, intracellular reactive oxygen species (iROS) production, and apoptosis. We rationalized that an overlapping of the drug activities should offer an enhanced anticancer effect, either through vertical inhibition of the Src-STAT3 axis or horizontal suppression of multiple pathways. We determined the toxicity induced by the drug combination and studied the kinetics of iROS production by fluorescence imaging and flow cytometry. Using genomic and proteomic techniques, including RNA-sequencing (RNA-seq), reverse transcription-quantitative polymerase chain reaction (RT-qPCR), and Western Blot, we subsequently identified the responsible pathways that contributed to the synergistic effects of the drug combination.

RESULTS

Compared to either drug alone, the drug combination showed enhanced potency against MDA-MB-468, MDA-MB-231, and MCF-7 human breast cancer (BC) cell lines and tumor spheroids. The drug combination induces both iROS generation and apoptosis in a time-dependent manner, following a 2-step kinetic profile. RNA-seq data revealed that the drug combination exhibited synergism through horizontal suppression of multiple pathways, possibly through a promotion of cell cycle arrest at the G1/S phase via the estrogen-mediated S-phase entry pathway, and partially via the BRCA1 and DNA damage response pathway.

CONCLUSION

Transcriptomic analyses revealed for the first time, that the estrogen-mediated S-phase entry pathway partially contributed to the synergistic effect of the drug combination. More importantly, our studies led to the discoveries of new potential therapeutic targets, such as E2F2, as well as a novel drug-induced targeting of estrogen receptor β (ESR2) approach for triple-negative breast cancer treatment, currently lacking of targeted therapies.

OncoOmics approaches to reveal essential genes in breast cancer: a panoramic view from pathogenesis to precision medicine

Breast cancer (BC) is the leading cause of cancer-related death among women and the most commonly diagnosed cancer worldwide. Although in recent years large-scale efforts have focused on identifying new therapeutic targets, a better understanding of BC molecular processes is required. Here we focused on elucidating the molecular hallmarks of BC heterogeneity and the oncogenic mutations involved in precision medicine that remains poorly defined. To fill this gap, we established an OncoOmics strategy that consists of analyzing genomic alterations, signaling pathways, protein-protein interactome network, protein expression, dependency maps in cell lines and patient-derived xenografts in 230 previously prioritized genes to reveal essential genes in breast cancer. As results, the OncoOmics BC essential genes were rationally filtered to 140. mRNA up-regulation was the most prevalent genomic alteration. The most altered signaling pathways were associated with basal-like and Her2-enriched molecular subtypes. RAC1, AKT1, CCND1, PIK3CA, ERBB2, CDH1, MAPK14, TP53, MAPK1, SRC, RAC3, BCL2, CTNNB1, EGFR, CDK2, GRB2, MED1 and GATA3 were essential genes in at least three OncoOmics approaches. Drugs with the highest amount of clinical trials in phases 3 and 4 were paclitaxel, docetaxel, trastuzumab, tamoxifen and doxorubicin. Lastly, we collected ~3,500 somatic and germline oncogenic variants associated with 50 essential genes, which in turn had therapeutic connectivity with 73 drugs. In conclusion, the OncoOmics strategy reveals essential genes capable of accelerating the development of targeted therapies for precision oncology.

Targeting Metabolic Deregulation Landscapes in Breast Cancer Subtypes

Metabolic deregulation is an emergent hallmark of cancer. Altered patterns of metabolic pathways result in exacerbated synthesis of macromolecules, increased proliferation, and resistance to treatment via alteration of drug processing. In addition, molecular heterogeneity creates a barrier to therapeutic options. In breast cancer, this broad variation in molecular metabolism constitutes, simultaneously, a source of prognostic and therapeutic challenges and a doorway to novel interventions. In this work, we investigated the metabolic deregulation landscapes in breast cancer molecular subtypes. Such landscapes are the regulatory signatures behind subtype-specific metabolic features. n = 735 breast cancer samples of the Luminal A, Luminal B, Her2+, and Basal subtypes, as well as n = 113 healthy breast tissue samples were analyzed. By means of a single-sample-based algorithm, deregulation for all metabolic pathways in every sample was determined. Deregulation levels match almost perfectly with the molecular classification, indicating that metabolic anomalies are closely associated with gene-expression signatures. Luminal B tumors are the most deregulated but are also the ones with higher within-subtype variance. We argued that this variation may underlie the fact that Luminal B tumors usually present the worst prognosis, a high rate of recurrence, and the lowest response to treatment in the long term. Finally, we designed a therapeutic scheme to regulate purine metabolism in breast cancer, independently of the molecular subtype. This scheme is founded on a computational tool that provides a set of FDA-approved drugs to target pathway-specific differentially expressed genes. By providing metabolic deregulation patterns at the single-sample level in breast cancer subtypes, we have been able to further characterize tumor behavior. This approach, together with targeted therapy, may open novel avenues for the design of personalized diagnostic, prognostic, and therapeutic strategies.

Master Regulators of Signaling Pathways: An Application to the Analysis of Gene Regulation in Breast Cancer

Analysis of gene regulatory networks allows the identification of master transcriptional factors that control specific groups of genes. In this work, we inferred a gene regulatory network from a large dataset of breast cancer samples to identify the master transcriptional factors that control the genes within signal transduction pathways. The focus in a particular subset of relevant genes constitutes an extension of the original Master Regulator Inference Algorithm (MARINa) analysis. This modified version of MARINa utilizes a restricted molecular signature containing genes from the 25 human pathways in KEGG's signal transduction category. Our breast cancer RNAseq expression dataset consists of 881 samples comprising tumors and normal mammary gland tissue. The top 10 master transcriptional factors found to regulate signal transduction pathways in breast cancer we identified are: TSHZ2, HOXA2, MEIS2, HOXA3, HAND2, HOXA5, TBX18, PEG3, GLI2, and CLOCK. The functional enrichment of the regulons of these master transcriptional factors showed an important proportion of processes related to morphogenesis. Our results suggest that, as part of the aberrant regulation of signaling pathways in breast cancer, pathways similar to the regulation of cell differentiation, cardiovascular system development, and vasculature development may be dysregulated and co-opted in favor of tumor development through the action of these transcription factors.

The Many Faces of Gene Regulation in Cancer: A Computational Oncogenomics Outlook

Cancer is a complex disease at many different levels. The molecular phenomenology of cancer is also quite rich. The mutational and genomic origins of cancer and their downstream effects on processes such as the reprogramming of the gene regulatory control and the molecular pathways depending on such control have been recognized as central to the characterization of the disease. More important though is the understanding of their causes, prognosis, and therapeutics. There is a multitude of factors associated with anomalous control of gene expression in cancer. Many of these factors are now amenable to be studied comprehensively by means of experiments based on diverse omic technologies. However, characterizing each dimension of the phenomenon individually has proven to fall short in presenting a clear picture of expression regulation as a whole. In this review article, we discuss some of the more relevant factors affecting gene expression control both, under normal conditions and in tumor settings. We describe the different omic approaches that we can use as well as the computational genomic analysis needed to track down these factors. Then we present theoretical and computational frameworks developed to integrate the amount of diverse information provided by such single-omic analyses. We contextualize this within a systems biology-based multi-omic regulation setting, aimed at better understanding the complex interplay of gene expression deregulation in cancer.

Spatial Organization of the Gene Regulatory Program: An Information Theoretical Approach to Breast Cancer Transcriptomics

Gene regulation may be studied from an information-theoretic perspective. Gene regulatory programs are representations of the complete regulatory phenomenon associated to each biological state. In diseases such as cancer, these programs exhibit major alterations, which have been associated with the spatial organization of the genome into chromosomes. In this work, we analyze intrachromosomal, or cis-, and interchromosomal, or trans-gene regulatory programs in order to assess the differences that arise in the context of breast cancer. We find that using information theoretic approaches, it is possible to differentiate cis-and trans-regulatory programs in terms of the changes that they exhibit in the breast cancer context, indicating that in breast cancer there is a loss of trans-regulation. Finally, we use these programs to reconstruct a possible spatial relationship between chromosomes.

Unveiling the Link Between Inflammation and Adaptive Immunity in Breast Cancer

Inflammation has been recognized as an important driver in the development and growth of malignancies. Inflammatory signaling in cancer emerges from the combinatorial interaction of several deregulated pathways. Pathway deregulation is often driven by changes in the underlying gene regulatory networks. Confronted with such complex scenario, it can be argued that a closer analysis of the structure of such regulatory networks will shed some light on how gene deregulation led to sustained inflammation in cancer. Here, we inferred an inflammation-associated gene regulatory network from 641 breast cancer and 78 healthy samples. A modular structure analysis of the regulatory network was carried out, revealing a hierarchical modular structure. Modules show significant overrepresentation score p-values for biological processes unveiling a definite association between inflammatory processes and adaptive immunity. Other modules are enriched for T-cell activation, differentiation of CD8⁺ lymphocytes and immune cell migration, thus reinforcing the aforementioned association. These analyses suggest that in breast cancer tumors, the balance between antitumor response and immune tolerance involving CD8⁺ T cells is tipped in favor of the tumor. One possible mechanism is the induction of tolerance and anergization of these cells by persistent antigen exposure.

Gene prioritization, communality analysis, networking and metabolic integrated pathway to better understand breast cancer pathogenesis

Consensus strategy was proved to be highly efficient in the recognition of gene-disease association. Therefore, the main objective of this study was to apply theoretical approaches to explore genes and communities directly involved in breast cancer (BC) pathogenesis. We evaluated the consensus between 8 prioritization strategies for the early recognition of pathogenic genes. A communality analysis in the protein-protein interaction (PPi) network of previously selected genes was enriched with gene ontology, metabolic pathways, as well as oncogenomics validation with the OncoPPi and DRIVE projects. The consensus genes were rationally filtered to 1842 genes. The communality analysis showed an enrichment of 14 communities specially connected with ERBB, PI3K-AKT, mTOR, FOXO, p53, HIF-1, VEGF, MAPK and prolactin signaling pathways. Genes with highest ranking were TP53, ESR1, BRCA2, BRCA1 and ERBB2. Genes with highest connectivity degree were TP53, AKT1, SRC, CREBBP and EP300. The connectivity degree allowed to establish a significant correlation between the OncoPPi network and our BC integrated network conformed by 51 genes and 62 PPi. In addition, CCND1, RAD51, CDC42, YAP1 and RPA1 were functional genes with significant sensitivity score in BC cell lines. In conclusion, the consensus strategy identifies both well-known pathogenic genes and prioritized genes that need to be further explored.

The Hierarchical Modular Structure of HER2+ Breast Cancer Network

HER2-enriched breast cancer is a complex disease characterized by the overexpression of the ERBB2 amplicon. While the effects of this genomic aberration on the pathology have been studied, genome-wide deregulation patterns in this subtype of cancer are also observed. A novel approach to the study of this malignant neoplasy is the use of transcriptional networks. These networks generally exhibit modular structures, which in turn may be associated to biological processes. This modular regulation of biological functions may also exhibit a hierarchical structure, with deeper levels of modular organization accounting for more specific functional regulation. In this work, we identified the most probable (maximum likelihood) model of the hierarchical modular structure of the HER2-enriched transcriptional network as reconstructed from gene expression data, and analyzed the statistical associations of modules and submodules to biological functions. We found modular structures, independent from direct ERBB2 amplicon regulation, involved in different biological functions such as signaling, immunity, and cellular morphology. Higher resolution submodules were identified in more specific functions, such as micro-RNA regulation and the activation of viral-like immune response. We propose the approach presented here as one that may help to unveil mechanisms involved in the development of the pathology.

Transcriptome profiling of the interconnection of pathways involved in malignant transformation and response to hypoxia

In tumor tissues, hypoxia is a commonly observed feature resulting from rapidly proliferating cancer cells outgrowing their surrounding vasculature network. Transformed cancer cells are known to exhibit phenotypic alterations, enabling continuous proliferation despite a limited oxygen supply. The four-step isogenic BJ cell model enables studies of defined steps of tumorigenesis: the normal, immortalized, transformed, and metastasizing stages. By transcriptome profiling under atmospheric and moderate hypoxic (3% O₂) conditions, we observed that despite being highly similar, the four cell lines of the BJ model responded strikingly different to hypoxia. Besides corroborating many of the known responses to hypoxia, we demonstrate that the transcriptome adaptation to moderate hypoxia resembles the process of malignant transformation. The transformed cells displayed a distinct capability of metabolic switching, reflected in reversed gene expression patterns for several genes involved in oxidative phosphorylation and glycolytic pathways. By profiling the stage-specific responses to hypoxia, we identified ASS1 as a potential prognostic marker in hypoxic tumors. This study demonstrates the usefulness of the BJ cell model for highlighting the interconnection of pathways involved in malignant transformation and hypoxic response.

Novel insights into chromosomal conformations in cancer

Exploring gene function is critical for understanding the complexity of life. DNA sequences and the three-dimensional organization of chromatin (chromosomal interactions) are considered enigmatic factors underlying gene function, and interactions between two distant fragments can regulate transactivation activity via mediator proteins. Thus, a series of chromosome conformation capture techniques have been developed, including chromosome conformation capture (3C), circular chromosome conformation capture (4C), chromosome conformation capture carbon copy (5C), and high-resolution chromosome conformation capture (Hi-C). The application of these techniques has expanded to various fields, but cancer remains one of the major topics. Interactions mediated by proteins or long noncoding RNAs (lncRNAs) are typically found using 4C-sequencing and chromatin interaction analysis by paired-end tag sequencing (ChIA-PET). Currently, Hi-C is used to identify chromatin loops between cancer risk-associated single-nucleotide polymorphisms (SNPs) found by genome-wide association studies (GWAS) and their target genes. Chromosomal conformations are responsible for altered gene regulation through several typical mechanisms and contribute to the biological behavior and malignancy of different tumors, particularly prostate cancer, breast cancer and hematologic neoplasms. Moreover, different subtypes may exhibit different 3D-chromosomal conformations. Thus, C-tech can be used to help diagnose cancer subtypes and alleviate cancer progression by destroying specific chromosomal conformations. Here, we review the fundamentals and improvements in chromosome conformation capture techniques and their clinical applications in cancer to provide insight for future research.

Network Modularity in Breast Cancer Molecular Subtypes

Breast cancer is a heterogeneous and complex disease, a clear manifestation of this is its classification into different molecular subtypes. On the other hand, gene transcriptional networks may exhibit different modular structures that can be related to known biological processes. Thus, modular structures in transcriptional networks may be seen as manifestations of regulatory structures that tightly controls biological processes. In this work, we identify modular structures on gene transcriptional networks previously inferred from microarray data of molecular subtypes of breast cancer: luminal A, luminal B, basal, and HER2-enriched. We analyzed the modules (communities) found in each network to identify particular biological functions (described in the Gene Ontology database) associated to them. We further explored these modules and their associated functions to identify common and unique features that could allow a better level of description of breast cancer, particularly in the basal-like subtype, the most aggressive and poor prognosis manifestation. Our findings related to the immune system and a decrease in cell death-related processes in basal subtype could help to understand it and design strategies for its treatment.