Search

Triple-negative breast cancer (TNBC), accounting for 15-20% of all breast cancers, has one of the poorest prognoses and survival rates. Metastasis, a critical process in cancer progression, causes most cancer-related deaths, underscoring the need for alternative therapeutic approaches. This study explores the anti-migratory, anti-invasive, anti-tumoral, and antimetastatic effects of copper coordination compounds Casiopeína IIIia (CasIIIia) and Casiopeína IIgly (CasIIgly) on MDA-MB-231 and 4T1 breast carcinoma cell lines in vitro and in vivo. These emerging anticancer agents, mixed chelate copper(II) compounds, induce apoptosis by generating reactive oxygen species (ROS) and causing DNA damage. Whole-transcriptome analysis via gene expression arrays indicated that subtoxic concentrations of CasIIIia upregulate genes involved in metal response mechanisms. Casiopeínas ® reduced TNBC cell viability dose-dependently and more efficiently than Cisplatin. At subtoxic concentrations (IC 20 ), they inhibited random and chemotactic migration of MDA-MB-231 and 4T1 cells by 50-60%, similar to Cisplatin, as confirmed by transcriptome analysis. In vivo, CasIIIia and Cisplatin significantly reduced tumor growth, volume, and weight in a syngeneic breast cancer model with 4T1 cells. Furthermore, both compounds significantly decreased metastatic foci in treated mice compared to controls. Thus, CasIIIia and CasIIgly are promising chemotherapeutic candidates against TNBC.

Dave Garbers' work significantly contributed to our understanding of sperm's regulated motility, capacitation, and the acrosome reaction. These key sperm functions involve complex multistep signaling pathways engaging numerous finely orchestrated elements. Despite significant progress, many parameters and interactions among these elements remain elusive. Mathematical modeling emerges as a potent tool to study sperm physiology, providing a framework to integrate experimental results and capture functional dynamics considering biochemical, biophysical, and cellular elements. Depending on research objectives, different modeling strategies, broadly categorized into continuous and discrete approaches, reveal valuable insights into cell function. These models allow the exploration of hypotheses regarding molecules, conditions, and pathways, whenever they become challenging to evaluate experimentally. This review presents an overview of current theoretical and experimental efforts to understand sperm motility regulation, capacitation, and the acrosome reaction. We discuss the strengths and weaknesses of different modeling strategies and highlight key findings and unresolved questions. Notable discoveries include the importance of specific ion channels, the role of intracellular molecular heterogeneity in capacitation and the acrosome reaction, and the impact of pH changes on acrosomal exocytosis. Ultimately, this review underscores the crucial importance of mathematical frameworks in advancing our understanding of sperm physiology and guiding future experimental investigations., (© 2024 The Author(s). Molecular Reproduction and Development published by Wiley Periodicals LLC.)

Background: Metformin, a widely prescribed antidiabetic drug, has shown several promising effects for cancer treatment. These effects have been shown to be mediated by dual modulation of the AMPK-mTORC1 axis, where AMPK acts upstream of mTORC1 to decrease its activity. Nevertheless, alternative pathways have been recently discovered suggesting that metformin can act through of different targets regulation., Methods: We performed a transcriptome screening analysis using HeLa xenograft tumors generated in NOD-SCID mice treated with or without metformin to examine genes regulated by metformin. Western Blot analysis, Immunohistochemical staining, and RT-qPCR were used to confirm alterations in gene expression. The TNMplot and GEPIA2 platform were used for in silico analysis of genes found up-regulated by metformin, in cervical cancer patients. We performed an AMPK knock-down using AMPK-targeted siRNAs and mTOR inhibition with rapamycin to investigate the molecular mechanisms underlying the effect of metformin in cervical cancer cell lines., Results: We shown that metformin decreases tumor growth and increased the expression of a group of antitumoral genes involved in DNA-binding transcription activator activity, hormonal response, and Dcp1-Dcp2 mRNA-decapping complex. We demonstrated that ZFP36 could act as a new molecular target increased by metformin. mTORC1 inhibition using rapamycin induces ZFP36 expression, which could suggest that metformin increases ZFP36 expression and requires mTORC1 inhibition for such effect. Surprisingly, in HeLa cells AMPK inhibition did not affect ZFP36 expression, suggesting that additional signal transducers related to suppressing mTORC1 activity, could be involved., Conclusions: These results highlight the importance of ZFP36 activation in response to metformin treatment involving mTORC1 inhibition., (© 2024. The Author(s).)

Breast cancer encompasses a diverse array of subtypes, each exhibiting distinct clinical characteristics and treatment responses. Unraveling the underlying regulatory mechanisms that govern gene expression patterns in these subtypes is essential for advancing our understanding of breast cancer biology. Gene co-expression networks (GCNs) help us identify groups of genes that work in coordination. Previous research has revealed a marked reduction in the interaction of genes located on different chromosomes within GCNs for breast cancer, as well as for lung, kidney, and hematopoietic cancers. However, the reasons behind why genes on the same chromosome often co-express remain unclear. In this study, we investigate the role of transcription factors in shaping gene co-expression networks within the four main breast cancer subtypes: Luminal A, Luminal B, HER2+, and Basal, along with normal breast tissue. We identify communities within each GCN and calculate the transcription factors that may regulate these communities, comparing the results across different phenotypes. Our findings indicate that, in general, regulatory behavior is to a large extent similar among breast cancer molecular subtypes and even in healthy networks. This suggests that transcription factor motif usage does not fully determine long-range co-expression patterns. Specific transcription factor motifs, such as CCGGAAG, appear frequently across all phenotypes, even involving multiple highly connected transcription factors. Additionally, certain transcription factors exhibit unique actions in specific subtypes but with limited influence. Our research demonstrates that the loss of inter-chromosomal co-expression is not solely attributable to transcription factor regulation. Although the exact mechanism responsible for this phenomenon remains elusive, this work contributes to a better understanding of gene expression regulatory programs in breast cancer.

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., (© 2023. The Author(s).)

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., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Zamora-Fuentes, Hernández-Lemus and Espinal-Enríquez.)

Pseudogenes are duplicates of protein-coding genes that have accumulated multiple detrimental alterations, rendering them unable to produce the protein they encode. Initially disregarded as "junk DNA" due to their perceived lack of functionality, research on their biological roles has been hindered by this assumption. Nevertheless, recent focus has shifted towards these molecules due to their abnormal expression in cancer phenotypes. In this review, our objective is to provide a thorough overview of the current understanding of pseudogene formation, the mechanisms governing their expression, and the roles they may play in promoting tumorigenesis.

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., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Hernández-Gómez, Hernández-Lemus and Espinal-Enríquez.)

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., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Zamora-Fuentes, Hernández-Lemus and Espinal-Enríquez.)

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.

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., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Hernández-Gómez, Hernández-Lemus and Espinal-Enríquez.)

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., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 González-Espinoza, Zamora-Fuentes, Hernández-Lemus and Espinal-Enríquez.)

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., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Alcalá-Corona, Sandoval-Motta, Espinal-Enríquez and Hernández-Lemus.)

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., (© 2021. The Author(s).)

Breast cancer is a complex, heterogeneous disease at the phenotypic and molecular level. In particular, the transcriptional regulatory programs are known to be significantly affected and such transcriptional alterations are able to capture some of the heterogeneity of the disease, leading to the emergence of breast cancer molecular subtypes. Recently, it has been found that network biology approaches to decipher such abnormal gene regulation programs, for instance by means of gene co-expression networks, have been able to recapitulate the differences between breast cancer subtypes providing elements to further understand their functional origins and consequences. Network biology approaches may be extended to include other co-expression patterns, like those found between genes and non-coding transcripts such as microRNAs (miRs). As is known, miRs play relevant roles in the establishment of normal and anomalous transcription processes. Commodore miRs (cdre-miRs) have been defined as miRs that, based on their connectivity and redundancy in co-expression networks, are potential control elements of biological functions. In this work, we reconstructed miR-gene co-expression networks for each breast cancer molecular subtype, from high throughput data in 424 samples from the Cancer Genome Atlas consortium. We identified cdre-miRs in three out of four molecular subtypes. We found that in each subtype, each cdre-miR was linked to a different set of associated genes, as well as a different set of associated biological functions. We used a systematic literature validation strategy, and identified that the associated biological functions to these cdre-miRs are hallmarks of cancer such as angiogenesis, cell adhesion, cell cycle and regulation of apoptosis. The relevance of such cdre-miRs as actionable molecular targets in breast cancer is still to be determined from functional studies., (© 2021 The Authors.)

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., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 García-Cortés, Hernández-Lemus and Espinal-Enríquez.)

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., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Andonegui-Elguera, Zamora-Fuentes, Espinal-Enríquez and Hernández-Lemus.)

Understanding the spatial and temporal patterns of mortality rates in a highly heterogeneous metropolis, is a matter of public policy interest. In this context, there is no, to the best of our knowledge, previous studies that correlate both spatio-temporal and age-specific mortality rates in Mexico City. Spatio-temporal Kriging modeling was used over five age-specific mortality rates (from the years 2000 to 2016 in Mexico City), to gain both spatial (borough and neighborhood) and temporal (year and trimester) data level description. Mortality age-specific patterns have been modeled using multilevel modeling for longitudinal data. Posterior tests were carried out to compare mortality averages between geo-spatial locations. Mortality correlation extends in all study groups for as long as 12 years and as far as 13.27 km. The highest mortality rate takes place in the Cuauhtémoc borough, the commercial, touristic and cultural core downtown of Mexico City. On the contrary, Tlalpan borough is the one with the lowest mortality rates in all the study groups. Post-productive mortality is the first age-specific cause of death, followed by infant, productive, pre-school and scholar groups. The combinations of spatio-temporal Kriging estimation and time-evolution linear mixed-effect models, allowed us to unveil relevant time and location trends that may be useful for public policy planning in Mexico City., Competing Interests: The authors have declared that no competing interests exist.

Assessment of the air quality in metropolitan areas is a major challenge in environmental sciences. Issues related include the distribution of monitoring stations, their spatial range, or missing information. In Mexico City, stations have been located spanning the entire Metropolitan zone for pollutants, such as CO, NO 2 , O 3 , SO 2 , PM 2.5 , PM 10 , NO, NO x , and PM CO . A fundamental question is whether the number and location of such stations are adequate to optimally cover the city. By analyzing spatio-temporal correlations for pollutant measurements, we evaluated the distribution and performance of monitoring stations in Mexico City from 2009 to 2018. Based on our analysis, air quality evaluation of those contaminants is adequate to cover the 16 boroughs of Mexico City, with the exception of SO 2 , since its spatial range is shorter than the one needed to cover the whole surface of the city. We observed that NO and NO x concentrations must be taken into account since their long-range dispersion may have relevant consequences for public health. With this approach, we may be able to propose policy based on systematic criteria to locate new monitoring stations., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Baca-López, Fresno, Espinal-Enríquez, Martínez-García, Camacho-López, Flores-Merino and Hernández-Lemus.)

Homicide is without doubt one of Mexico's most important security problems, with data showing that this dismal kind of violence sky-rocketed shortly after the war on drugs was declared in 2007. Since then, violent war-like zones have appeared and disappeared throughout Mexico, causing unfathomable human, social and economic losses. One of the most emblematic of these zones is the Monterrey metropolitan area (MMA), a central scenario in the narco-war. Being an important metropolitan area in Mexico and a business hub, MMA has counted hundreds to thousands of casualties. In spite of several approaches being developed to understand and analyze crime in general, and homicide in particular, the lack of accurate spatio-temporal homicide data results in incomplete descriptions. In order to describe the manner in which violence has evolved and spread in time and space through the city, here we propose a network-based approach. For this purpose, we define a homicide network where nodes are geographical entities that are connected through spatial and temporal relationships. We analyzed the time series of homicides in different municipalities and neighborhoods of the MMA, to observe whether or not a global correlation appeared. We studied the spatial correlation between neighborhoods where homicides took place, to observe whether distance is a factor of influence in the frequency of homicides. We constructed yearly co-occurrence networks, by correlating neighborhoods with homicides happening within a same week, and counting the co-occurrences of these neighborhood pairs in 1 year. We also constructed a crime network by aggregating all data of homicides, eliminating the temporal correlation, in order to observe whether homicide clusters appeared, and what those clusters were distributed geographically. Finally, we correlated the location and frequency of homicides with roads, freeways and highways, to observe if a trend in the homicidal violence appeared. Our network approach in the homicide evolution of MMA allows us to identify that (1) analyzing the whole 86-month period, we observed a correlation between close cities, which decreases in distant places. (2) at neighborhood level, correlations are not distance-dependent, on the contrary, highest co-occurrences appeared between distant neighborhoods and a polygon formed by close neighborhoods in downtown Monterrey. Moreover, (3) An elevated number of homicides occur close to the 85th freeway, which connects MMA with the US border. (4) Some socioeconomic barriers determine the presence of homicide violence. Finally, (5) we show a relation between homicidal crime and the urban landscape by studying the distance of safe and violent neighborhoods to the closest highway and by studying the evolution of highway and crime distance over the cartel-related years and the following period. With this approach, we are able to describe the spatial and temporal evolution of homicidal crime in a metropolitan area.

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., (Copyright © 2020 Zamora-Fuentes, Hernández-Lemus and Espinal-Enríquez.)

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., (Copyright © 2020 García-Cortés, de Anda-Jáuregui, Fresno, Hernández-Lemus and Espinal-Enríquez.)

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., (Copyright © 2020 Serrano-Carbajal, Espinal-Enríquez and Hernández-Lemus.)

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.

Casiopeinas are a group of copper-based compounds designed to be used as less toxic, more efficient chemotherapeutic agents. In this study, we analyzed the in vitro effects of Casiopeina II-gly on the expression of canonical biological pathways. Using microarray data from HeLa cell lines treated with Casiopeina II-gly, we identified biological pathways that are perturbed after treatment. We present a novel approach integrating pathway analysis and network theory: The Pathway Crosstalk Network. We constructed a network with deregulated pathways, featuring links between those pathways that crosstalk with each other. We identified modules grouping deregulated pathways that are functionally related. Through this approach, we were able to identify three features of Casiopeina treatment: (a) Perturbation of signaling pathways, related to induction of apoptosis; (b) perturbation of metabolic pathways, and (c) activation of immune responses. These findings can be useful to drive new experimental exploration on their role in adverse effects and efficacy of Casiopeinas., (Copyright © 2018 Elsevier Ltd. All rights reserved.)

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.

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.

Breast cancer is a major public health problem which treatment needs new pharmacological options. In the last decades, during the postgenomic era new theoretical and technological tools that give us novel and promising ways to address these problems have emerged. In this work, we integrate several tools that exploit disease-specific experimental transcriptomic results in addition to information from biological and pharmacological data bases obtaining a contextual prioritization of pathways and drugs in breast cancer subtypes. The usefulness of these results should be evaluated in terms of drug repurposing in each breast cancer molecular subtype therapy. In favor of breast cancer patients, this methodology could be further developed to provide personalized treatment schemes. The latter are particularly needed in those breast cancer subtypes with limited therapeutic options or those who have developed resistance to the current pharmacological schemes.

Alterations to transcriptional regulation are an important factor in breast cancer. Noncoding RNA, such as microRNA (miR), have very influential roles in the transcriptional regulation of genes. Transcriptional regulation can be successfully modeled and analyzed using complex network theory. Particularly, interactions between two distinct classes of biological elements, such as miR and genes, can be approached through the bipartite network formalism. Based on bipartite network properties, it is possible to identify highly influential miRs in the network, such as those that have a large number of connections indicating regulation of a large set of genes. Some miRs in a network are nonredundant, which indicates that they are solely responsible of the regulation of a particular set of genes, which in turn may be associated to a particular biological process. We hypothesize that highly influential, nonredundant miRs, which we call Commodore miRs (Cdre-miRs), have an important role on the control of biological functions through transcriptional networks. In this work, we analyze the regulation of gene expression by miRs in healthy and cancerous breast tissue using bipartite miR-gene networks inferred from the Cancer Genome Atlas (TCGA) expression data. We observe differences in the degree, clustering coefficient and redundancy distributions for miRs and genes in the network, indicating differences in the way that these elements interact with each other. Furthermore, we identify a small set of five Cdre-miRs in the breast cancer network: miR-190b, miR-let7i, miR-292-b, miR-511, and miR-141. The neighborhood of genes controlled by each of these miRs is involved in particular biological functions such as dynein structure-associated processes, immune response, angiogenesis, cytokine activity, and cell motility. We propose that these Cdre-miRs are important control elements of biological functions deregulated in breast cancer.

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.

Over the last years, microRNAs (miRs) have shown to be crucial for breast tumour establishment and progression. To understand the influence that miRs have over transcriptional regulation in breast cancer, we constructed mutual information networks from 86 TCGA matched breast invasive carcinoma and control tissue RNA-Seq and miRNA-Seq sequencing data. We show that miRs are determinant for tumour and control data network structure. In tumour data network, miR-200, miR-199 and neighbour miRs seem to cooperate on the regulation of the acquisition of epithelial and mesenchymal traits by the biological processes: Epithelial-Mesenchymal Transition (EMT) and Mesenchymal to Epithelial Transition (MET). Despite structural differences between tumour and control networks, we found a conserved set of associations between miR-200 family members and genes such as VIM, ZEB-1/2 and TWIST-1/2. Further, a large number of miRs observed in tumour network mapped to a specific chromosomal location in DLK1-DIO3 (Chr14q32); some of those miRs have also been associated with EMT and MET regulation. Pathways related to EMT and TGF-beta reinforce the relevance of miR-200, miR-199 and DLK1-DIO3 cluster in breast cancer. With this approach, we stress that miR inclusion in gene regulatory network construction improves our understanding of the regulatory mechanisms underlying breast cancer biology.

Spermatozoa sea urchin swimming behaviour is regulated by small peptides from the egg outer envelope. Speract, such a peptide, after binding to its receptor in Strongylocentrotus purpuratus sperm flagella, triggers a signaling pathway that culminates with a train of intracellular calcium oscillations, correlated with changes in sperm swimming pattern. This pathway has been widely studied but not fully characterized. Recent work on Arbacia punctulata sea urchin spermatozoa has documented the presence of the Ca 2+ CatSper channel in their flagella and its involvement in chemotaxis. However, if other calcium channels participate in chemotaxis remains unclear. Here, based on an experimentally-backed logical network model, we conclude that CatSper is fundamental in the S. purpuratus speract-activated sea urchin sperm signaling cascade, although other Ca 2+ channels could still be relevant. We also present for the first time experimental corroboration of its active presence in S. purpuratus sperm flagella. We argue, prompted by in silico knock-out calculations, that CatSper is the main generator of calcium oscillations in the signaling pathway and that other calcium channels, if present, have a complementary role. The approach adopted here allows us to unveil processes, which are hard to detect exclusively by experimental procedures.

Breast cancer is a complex heterogeneous disease. Common hallmark features of cancer can be found. Their origin may be traced back to their intricate relationships governing regulatory programs during the development of this disease. To unveil distinctive features of the transcriptional regulation program in breast cancer, a pipeline for RNA-seq analysis in 780 breast cancer and 101 healthy breast samples, at gene expression and network level, was implemented. Inter-chromosomal relationships between genes resulted strikingly scarce in a cancer network, in comparison to its healthy counterpart. We suggest that inter-chromosomal regulation loss may be a novel feature in breast cancer. Additional evidence was obtained by independent validation in microarray and Hi-C data as well as supplementary computational analyses. Functional analysis showed upregulation in processes related to cell cycle and division; while migration, adhesion and cell-to-cell communication, were downregulated. Both the BRCA1 DNA repairing signalling and the Estrogen-mediated G1/S phase entry pathways were found upregulated. In addition, a synergistic underexpression of the γ-protocadherin complex, located at Chr5q31 is also shown. This region has previously been reported to be hypermethylated in breast cancer. These findings altogether provide further evidence for the central role of transcriptional regulatory programs in shaping malignant phenotypes.

Understanding the general principles underlying genetic regulation in eukaryotes is an incomplete and challenging endeavor. The lack of experimental information regarding the regulation of the whole set of transcription factors and their targets in different cell types is one of the main reasons to this incompleteness. So far, there is a small set of curated known interactions between transcription factors and their downstream genes. Here, we built a transcription factor network for human monocytic THP-1 myeloid cells based on the experimentally curated FANTOM4 database where nodes are genes and the experimental interactions correspond to links. We present the topological parameters which define the network as well as some global structural features and introduce a relative inuence parameter to quantify the relevance of a transcription factor in the context of induction of a phenotype. Genes like ZHX2, ADNP, or SMAD6 seem to be highly regulated to avoid an avalanche transcription event. We compare these results with those of RegulonDB , a highly curated transcriptional network for the prokaryotic organism E. coli , finding similarities between general hallmarks on both transcriptional programs. We believe that an approach, such as the one shown here, could help to understand the one regulation of transcription in eukaryotic cells.

Breast cancer heterogeneity is evident at the clinical, histological and molecular level. High throughput technologies allowed the identification of intrinsic subtypes that capture transcriptional differences among tumors. A remaining question is whether said differences are associated to a particular transcriptional program which involves different connections between the same molecules. In other words, whether particular transcriptional network architectures can be linked to specific phenotypes. In this work we infer, construct and analyze transcriptional networks from whole-genome gene expression microarrays, by using an information theory approach. We use 493 samples of primary breast cancer tissue classified in four molecular subtypes: Luminal A, Luminal B, Basal and HER2-enriched. For comparison, a network for non-tumoral mammary tissue (61 samples) is also inferred and analyzed. Transcriptional networks present particular architectures in each breast cancer subtype as well as in the non-tumor breast tissue. We find substantial differences between the non-tumor network and those networks inferred from cancer samples, in both structure and gene composition. More importantly, we find specific network architectural features associated to each breast cancer subtype. Based on breast cancer networks' centrality, we identify genes previously associated to the disease, either, generally (i.e., CNR2) or to a particular subtype (such as LCK). Similarly, we identify LUZP4, a gene barely explored in breast cancer, playing a role in transcriptional networks with subtype-specific relevance. With this approach we observe architectural differences between cancer and non-cancer at network level, as well as differences between cancer subtype networks which might be associated with breast cancer heterogeneity. The centrality measures of these networks allow us to identify genes with potential biomedical implications to breast cancer.

Gene regulatory networks are useful to understand the activity behind the complex mechanisms in transcriptional regulation. A main goal in contemporary biology is using such networks to understand the systemic regulation of gene expression. In this work, we carried out a systematic study of a transcriptional regulatory network derived from a comprehensive selection of all potential transcription factor interactions downstream from MEF2C, a human transcription factor master regulator. By analyzing the connectivity structure of such network, we were able to find different biologically functional processes and specific biochemical pathways statistically enriched in communities of genes into the network, such processes are related to cell signaling, cell cycle and metabolism. In this way we further support the hypothesis that structural properties of biological networks encode an important part of their functional behavior in eukaryotic cells.

The understanding of the mechanisms associated with the action of chemotherapeutic agents is fundamental to assess and account for possible side-effects of such treatments. Casiopeínas have demonstrated a cytotoxic effect by activation of pro-apoptotic processes in malignant cells. Such processes have been proved to activate the apoptotic intrinsic route, as well as cell cycle arrest. Despite this knowledge, the whole mechanism of action of Casiopeínas is yet to be completely understood. In this work we implement a systems biology approach based on two pathway analysis tools (Over-Representation Analysis and Causal Network Analysis) to observe changes in some hallmarks of cancer, induced by this copper-based chemotherapeutic agent in HeLa cell lines. We find that the metabolism of metal ions is exacerbated, as well as cell division processes being globally diminished. We also show that cellular migration and proliferation events are decreased. Moreover, the molecular mechanisms of liver protection are increased in the cell cultures under the actions of Casiopeínas, unlike the case in many other cytotoxic drugs. We argue that this chemotherapeutic agent may be promising, given its protective hepatic function, concomitant with its cytotoxic participation in the onset of apoptotic processes in malignant cells.

Pathway analysis is a set of widely used tools for research in life sciences intended to give meaning to high-throughput biological data. The methodology of these tools settles in the gathering and usage of knowledge that comprise biomolecular functioning, coupled with statistical testing and other algorithms. Despite their wide employment, pathway analysis foundations and overall background may not be fully understood, leading to misinterpretation of analysis results. This review attempts to comprise the fundamental knowledge to take into consideration when using pathway analysis as a hypothesis generation tool. We discuss the key elements that are part of these methodologies, their capabilities and current deficiencies. We also present an overview of current and all-time popular methods, highlighting different classes across them. In doing so, we show the exploding diversity of methods that pathway analysis encompasses, point out commonly overlooked caveats, and direct attention to a potential new class of methods that attempt to zoom the analysis scope to the sample scale.

Steroid hormones are involved on cell growth, development and differentiation. Such effects are often mediated by steroid receptors. One paradigmatic example of this coupling is the estrogen signaling pathway. Its dysregulation is involved in most tumors of the mammary gland. It is thus an important pharmacological target in breast cancer. This pathway, however, crosstalks with several other molecular pathways, a fact that may have consequences for the effectiveness of hormone modulating drug therapies, such as tamoxifen. For this work, we performed a systematic analysis of the major routes involved in crosstalk phenomena with the estrogen pathway - based on gene expression experiments (819 samples) and pathway analysis (493 samples) - for biopsy-captured tissue and contrasted in two independent datasets with in vivo and in vitro pharmacological stimulation. Our results confirm the presence of a number of crosstalk events across the estrogen signaling pathway with others that are dysregulated in different molecular subtypes of breast cancer. These may be involved in proliferation, invasiveness and apoptosis-evasion in patients. The results presented may open the way to new designs of adjuvant and neoadjuvant therapies for breast cancer treatment., (Copyright © 2015 Elsevier Ltd. All rights reserved.)

Gene regulatory networks account for the delicate mechanisms that control gene expression. Under certain circumstances, gene regulatory programs may give rise to amplification cascades. Such transcriptional cascades are events in which activation of key-responsive transcription factors called master regulators trigger a series of gene expression events. The action of transcriptional master regulators is then important for the establishment of certain programs like cell development and differentiation. However, such cascades have also been related with the onset and maintenance of cancer phenotypes. Here we present a systematic implementation of a series of algorithms aimed at the inference of a gene regulatory network and analysis of transcriptional master regulators in the context of primary breast cancer cells. Such studies were performed in a highly curated database of 880 microarray gene expression experiments on biopsy-captured tissue corresponding to primary breast cancer and healthy controls. Biological function and biochemical pathway enrichment analyses were also performed to study the role that the processes controlled - at the transcriptional level - by such master regulators may have in relation to primary breast cancer. We found that transcription factors such as AGTR2, ZNF132, TFDP3 and others are master regulators in this gene regulatory network. Sets of genes controlled by these regulators are involved in processes that are well-known hallmarks of cancer. This kind of analyses may help to understand the most upstream events in the development of phenotypes, in particular, those regarding cancer biology., (Copyright © 2015 Elsevier Ltd. All rights reserved.)

Genome-wide association studies and other genetic analyses have identified a large number of genes and variants implicating a variety of disease etiological mechanisms. It is imperative for the study of human diseases to put these genetic findings into a coherent functional context. Here we use system biology tools to examine disease connections of five master genes for CD4+ T cell subtypes (TBX21, GATA3, RORC, BCL6, and FOXP3). We compiled a list of genes functionally interacting (protein-protein interaction, or by acting in the same pathway) with the master genes, then we surveyed the disease connections, either by experimental evidence or by genetic association. Embryonic lethal genes (also known as essential genes) are over-represented in master genes and their interacting genes (55% versus 40% in other genes). Transcription factors are significantly enriched among genes interacting with the master genes (63% versus 10% in other genes). Predicted haploinsufficiency is a feature of most these genes. Disease-connected genes are enriched in this list of genes: 42% of these genes have a disease connection according to Online Mendelian Inheritance in Man (OMIM) (versus 23% in other genes), and 74% are associated with some diseases or phenotype in a Genome Wide Association Study (GWAS) (versus 43% in other genes). Seemingly, not all of the diseases connected to genes surveyed were immune related, which may indicate pleiotropic functions of the master regulator genes and associated genes., (Copyright © 2015 Elsevier Ltd. All rights reserved.)

Since December 2006, more than a thousand cities in México have suffered the effects of the war between several drug cartels, amongst themselves, as well as with Mexican armed forces. Sources are not in agreement about the number of casualties of this war, with reports varying from 30 to 100 thousand dead; the economic and social ravages are impossible to quantify. In this work we analyze the official report of casualties in terms of the location and the date of occurrence of the homicides. We show how the violence, as reflected by the number of casualties, has increased over time and spread across the country. Next, based on the correlations between cities in the changes of the monthly number of casualties attributed to organized crime, we construct a narco-war network where nodes are the affected cities and links represent correlations between them. We find that close geographical distance between violent cities does not imply a strong correlation amongst them. We observe that the dynamics of the conflict has evolved in short-term periods where a small core of violent cities determines the main theatre of the war at each stage. This kind of analysis may also help to describe the emergence and propagation of gang-related violence waves.

Background: Thyroid cancer (TC) is the most common malignant cancer of the Endocrine System. Histologically, there are three main subtypes of TC: follicular, papillary and anaplastic. Diagnosing a thyroid tumor subtype with a high level of accuracy and confidence is still a difficult task because genetic, molecular and cellular mechanisms underlying the transition from differentiated to undifferentiated thyroid tumors are not well understood. A genome-wide analysis of these three subtypes of thyroid carcinoma was carried out in order to identify significant differences in expression levels as well as enriched pathways for non-shared molecular and cellular features between subtypes., Results: Inhibition of matrix metalloproteinases pathway is a major event involved in thyroid cancer progression and its dysregulation may result crucial for invasiveness, migration and metastasis. This pathway is drastically altered in ATC while in FTC and PTC, the most important pathways are related to DNA-repair activation or cell to cell signaling events., Conclusion: A progression from FTC to PTC and then to ATC was detected and validated on two independent datasets. Moreover, PTX3, COLEC12 and PDGFRA genes were found as possible candidates for biomarkers of ATC while GPR110 could be tested to distinguish PTC over other tumor subtypes. The genome-wide analysis emphasizes the preponderance of pathway-dysregulation mechanisms over simple gene-malfunction as the main mechanism involved in the development of a cancer phenotype.

The motility of spermatozoa of both Lytechinus pictus and Strongylocentrotus purpuratus sea urchin species is modulated by the egg-derived decapeptide speract via an oscillatory [Ca2+]-dependent signaling pathway. Comprehension of this pathway is hence directly related to the understanding of regulated sperm swimming. Niflumic acid (NFA), a nonsteroidal anti-inflammatory drug alters several ion channels. Though unspecific, NFA profoundly affects how sea urchin sperm respond to speract, increasing the [Ca2+]i oscillation period, amplitude, peak and average level values of the responses in immobilized and swimming cells. A previous logical network model we developed for the [Ca2+] dynamics of speract signaling cascade in sea urchin sperm allows integrated dissection of individual and multiple actions of NFA. Among the channels affected by NFA are: hyperpolarization-activated and cyclic nucleotide gated Na+ channels (HCN), [Ca2+]-dependent Cl- channels (CaCC) and [Ca2+]-dependent K+ channels (CaKC), all present in the sea urchin genome. Here, using our model we investigated the effect of blocking in silico HCN and CaCC channels suggested by experiments. Regarding CaKC channels, arguments can be provided for either their blockage or activation by NFA. Our study yielded two scenarios compliant with experimental observations: i) under CaKC inhibition, this [Ca2+]-dependent K+ channel should be different from the Slo1 channel and ii) under activation of the CaKC channel, another [Ca2+] channel not considered previously in the network is required, such as the pH-dependent CatSper channel. Additionally, our findings predict cause-effect relations resulting from a selective inhibition of those channels. Knowledge of these relations may be of consequence for a variety of electrophysiological studies and have an impact on drug related investigations. Our study contributes to a better grasp of the network dynamics and suggests further experimental work.