RESUMEN

Inferring and characterizing gene co-expression networks has led to important insights on the molecular mechanisms of complex diseases. Most co-expression analyses to date have been performed on gene expression data collected from bulk tissues with different cell type compositions across samples. As a result, the co-expression estimates only offer an aggregated view of the underlying gene regulations and can be confounded by heterogeneity in cell type compositions, failing to reveal gene coordination that may be distinct across different cell types. In this paper, we introduce a flexible framework for estimating cell-type-specific gene co-expression networks from bulk sample data, without making specific assumptions on the distributions of gene expression profiles in different cell types. We develop a novel sparse least squares estimator, referred to as CSNet, that is efficient to implement and has good theoretical properties. Using CSNet, we analyzed the bulk gene expression data from a cohort study on Alzheimer's disease and identified previously unknown cell-type-specific co-expressions among Alzheimer's disease risk genes, suggesting cell-type-specific disease mechanisms.

RESUMEN

MOTIVATION: Systems biology analyses often use correlations in gene expression profiles to infer co-expression networks that are then used as input for gene regulatory network inference or to identify functional modules of co-expressed or putatively co-regulated genes. While systematic biases, including batch effects, are known to induce spurious associations and confound differential gene expression analyses (DE), the impact of batch effects on gene co-expression has not been fully explored. Methods have been developed to adjust expression values, ensuring conditional independence of mean and variance from batch or other covariates for each gene, resulting in improved fidelity of DE analysis. However, such adjustments do not address the potential for spurious differential co-expression (DC) between groups. Consequently, uncorrected, artifactual DC can skew the correlation structure, leading to the identification of false, non-biological associations, even when the input data is corrected using standard batch correction. RESULTS: In this work, we demonstrate the persistence of confounders in covariance after standard batch correction using synthetic and real-world gene expression data examples. We then introduce Co-expression Batch Reduction Adjustment (COBRA), a method for computing a batch-corrected gene co-expression matrix based on estimating a conditional covariance matrix. COBRA estimates a reduced set of parameters expressing the co-expression matrix as a function of the sample covariates, allowing control for continuous and categorical covariates. COBRA is computationally efficient, leveraging the inherently modular structure of genomic data to estimate accurate gene regulatory associations and facilitate functional analysis for high-dimensional genomic data. AVAILABILITY AND IMPLEMENTATION: COBRA is available under the GLP3 open source license in R and Python in netZoo (https://netzoo.github.io). SUPPLEMENTARY INFORMATION: Supplementary information is available at Bioinformatics online.

RESUMEN

The omnigenic model posits that genetic risk for traits with complex heritability involves cumulative effects of peripheral genes on mechanistic "core genes," suggesting that in a network of genes, those closer to clusters including core genes should have higher GWAS signals. In gene co-expression networks, we confirmed that GWAS signals accumulate in genes more connected to risk-enriched gene clusters, highlighting across-network risk convergence. This was strongest in adult psychiatric disorders, especially schizophrenia (SCZ), spanning 70% of network genes, suggestive of super-polygenic architecture. In snRNA-seq cell type networks, SCZ risk convergence was strongest in L2/L3 excitatory neurons. We prioritized genes most connected to SCZ-GWAS genes, which showed robust association to a CRISPRa measure of PGC3 regulation and were consistently identified across several brain regions. Several genes, including dopamine-associated ones, were prioritized specifically in the striatum. This strategy thus retrieves current drug targets and can be used to prioritize other potential drug targets.

RESUMEN

A neurological dogma is that the contralateral effects of brain injury are set through crossed descending neural tracts. We have recently identified a novel topographic neuroendocrine system (T-NES) that operates via a humoral pathway and mediates the left-right side-specific effects of unilateral brain lesions. In rats with completely transected thoracic spinal cords, unilateral injury to the sensorimotor cortex produced contralateral hindlimb flexion, a proxy for neurological deficit. Here, we investigated in acute experiments whether T-NES consists of left and right counterparts and whether they differ in neural and molecular mechanisms. We demonstrated that left- and right-sided hormonal signaling is differentially blocked by the δ-, κ- and µ-opioid antagonists. Left and right neurohormonal signaling differed in targeting the afferent spinal mechanisms. Bilateral deafferentation of the lumbar spinal cord abolished the hormone-mediated effects of the left-brain injury but not the right-sided lesion. The sympathetic nervous system was ruled out as a brain-to-spinal cord-signaling pathway since hindlimb responses were induced in rats with cervical spinal cord transections that were rostral to the preganglionic sympathetic neurons. Analysis of gene-gene co-expression patterns identified the left- and right-side-specific gene co-expression networks that were coordinated via the humoral pathway across the hypothalamus and lumbar spinal cord. The coordination was ipsilateral and disrupted by brain injury. These findings suggest that T-NES is bipartite and that its left and right counterparts contribute to contralateral neurological deficits through distinct neural mechanisms, and may enable ipsilateral regulation of molecular and neural processes across distant neural areas along the neuraxis.

Asunto(s)

RESUMEN

BACKGROUND: A widely used approach for extracting information from gene expression data employs the construction of a gene co-expression network and the subsequent computational detection of gene clusters, called modules. WGCNA and related methods are the de facto standard for module detection. The purpose of this work is to investigate the applicability of more sophisticated algorithms toward the design of an alternative method with enhanced potential for extracting biologically meaningful modules. RESULTS: We present self-learning gene clustering pipeline (SGCP), a spectral method for detecting modules in gene co-expression networks. SGCP incorporates multiple features that differentiate it from previous work, including a novel step that leverages gene ontology (GO) information in a self-leaning step. Compared with widely used existing frameworks on 12 real gene expression datasets, we show that SGCP yields modules with higher GO enrichment. Moreover, SGCP assigns highest statistical importance to GO terms that are mostly different from those reported by the baselines. CONCLUSION: Existing frameworks for discovering clusters of genes in gene co-expression networks are based on relatively simple algorithmic components. SGCP relies on newer algorithmic techniques that enable the computation of highly enriched modules with distinctive characteristics, thus contributing a novel alternative tool for gene co-expression analysis.

Asunto(s)

RESUMEN

Objectives: To develop a multi-omics prognostic model integrating transcriptomics and radiomics for predicting overall survival in patients with glioblastoma multiforme (GBM), and investigate the biological pathways of radiomics patterns. Materials and methods: Transcription profiles of GBM patients and normal controls were used to obtain differentially expressed mRNAs and long non-coding RNAs (lncRNAs). Radiomics features were extracted from magnetic resonance imaging (MRI). Least absolute shrinkage and selection operator (LASSO) Cox regression was employed to select survival-associated features for the construction of transcriptomics and radiomics signatures. Genes associated with GBM prognosis were identified through the analysis of lncRNA-mRNA co-expression networks and Weighted Gene Co-expression Network Analysis (WGCNA), and their biological pathways were investigated using Genomes enrichment analysis. Transcriptomics, radiomics, and clinical data were integrated to evaluate the multi-omics prognostic model's performance. Results: LASSO Cox regression yielded 21 survival-related features, including 19 transcriptomics features and 2 radiomics features. Based on transcriptomics and radiomics signature, GBM patients were classified as high-risk or low-risk. The genes obtained from the co-expression network screen were associated with microtubule binding, while those from the WGCNA screen were associated with growth factor receptor binding. In the training set, the AUC values for the multi-omics model and clinical model were 0.964 and 0.830, respectively, while in the validation set, they were 0.907 and 0.787. The multi-omics prognostic model outperformed the clinical prognostic model. Conclusions: The co-expression network and WGCNA methods revealed genes associated with multiple biological pathways in GBM. The multi-omics prognostic model demonstrated excellent performance and indicated significant potential for clinical application.

RESUMEN

In reproductive technologies, uncovering the molecular aspects of oocyte and embryo competence under different conditions is crucial for refining protocols and enhancing efficiency. RNA-seq generates high-throughput data and provides transcriptomes that can undergo additional computational analyses. This study presented the transcriptomic profiles of in vitro matured oocytes and blastocysts produced in vitro from buffalo crossbred (Bubalus bubalis), coupled with gene co-expression and module preservation analysis. Cumulus Oophorus Complexes, obtained from slaughterhouse-derived ovaries, were subjected to in vitro maturation to yield metaphase II oocytes (616) or followed in vitro fertilization and culture to yield blastocysts for sequencing (526). Oocyte maturation (72%, ±3.34 sd) and embryo development (21.3%, ±4.18 sd) rates were obtained from three in vitro embryo production routines following standard protocols. Sequencing of 410 metaphase II oocytes and 70 hatched blastocysts (grade 1 and 2) identified a total of 13,976 genes, with 62% being ubiquitously expressed (8,649). Among them, the differentially expressed genes (4,153) and the strongly variable genes with the higher expression (fold-change above 11) were highlighted in oocytes (BMP15, UCHL1, WEE1, NLRPs, KPNA7, ZP2, and ZP4) and blastocysts (APOA1, KRT18, ANXA2, S100A14, SLC34A2, PRSS8 and ANXA2) as representative indicators of molecular quality. Additionally, genes exclusively found in oocytes (224) and blastocysts (2,200) with specific biological functions were identified. Gene co-expression network and module preservation analysis revealed strong preservation of functional modules related to exosome components, steroid metabolism, cell proliferation, and morphogenesis. However, cell cycle and amino acid transport modules exhibited weak preservation, which may reflect differences in embryo development kinetics and the activation of cell signaling pathways between buffalo and bovine. This comprehensive transcriptomic profile serves as a valuable resource for assessing the molecular quality of buffalo oocytes and embryos in future in vitro embryo production assays.

RESUMEN

BACKGROUND: Malaria results in more than 550,000 deaths each year due to drug resistance in the most lethal Plasmodium (P.) species P. falciparum. A full P. falciparum genome was published in 2002, yet 44.6% of its genes have unknown functions. Improving the functional annotation of genes is important for identifying drug targets and understanding the evolution of drug resistance. RESULTS: Genes function by interacting with one another. So, analyzing gene co-expression networks can enhance functional annotations and prioritize genes for wet lab validation. Earlier efforts to build gene co-expression networks in P. falciparum have been limited to a single network inference method or gaining biological understanding for only a single gene and its interacting partners. Here, we explore multiple inference methods and aim to systematically predict functional annotations for all P. falciparum genes. We evaluate each inferred network based on how well it predicts existing gene-Gene Ontology (GO) term annotations using network clustering and leave-one-out crossvalidation. We assess overlaps of the different networks' edges (gene co-expression relationships), as well as predicted functional knowledge. The networks' edges are overall complementary: 47-85% of all edges are unique to each network. In terms of the accuracy of predicting gene functional annotations, all networks yielded relatively high precision (as high as 87% for the network inferred using mutual information), but the highest recall reached was below 15%. All networks having low recall means that none of them capture a large amount of all existing gene-GO term annotations. In fact, their annotation predictions are highly complementary, with the largest pairwise overlap of only 27%. We provide ranked lists of inferred gene-gene interactions and predicted gene-GO term annotations for future use and wet lab validation by the malaria community. CONCLUSIONS: The different networks seem to capture different aspects of the P. falciparum biology in terms of both inferred interactions and predicted gene functional annotations. Thus, relying on a single network inference method should be avoided when possible. SUPPLEMENTARY DATA: Attached.

Asunto(s)

RESUMEN

Pediatric T-cell Acute Lymphoblastic Leukemia (T-ALL) relapses are still associated with a dismal outcome, justifying the search for new therapeutic targets and relapse biomarkers. Using single-cell RNA sequencing (scRNAseq) data from three paired samples of pediatric T-ALL at diagnosis and relapse, we first conducted a high-dimensional weighted gene co-expression network analysis (hdWGCNA). This analysis highlighted several gene co-expression networks (GCNs) and identified relapse-associated hub genes, which are considered potential driver genes. Shared relapse-expressed genes were found to be related to antigen presentation (HLA, B2M), cytoskeleton remodeling (TUBB, TUBA1B), translation (ribosomal proteins, EIF1, EEF1B2), immune responses (MIF, EMP3), stress responses (UBC, HSP90AB1/AA1), metabolism (FTH1, NME1/2, ARCL4C), and transcriptional remodeling (NF-κB family genes, FOS-JUN, KLF2, or KLF6). We then utilized sparse partial least squares discriminant analysis to select from a pool of 481 unique leukemic hub genes, which are the genes most discriminant between diagnosis and relapse states (comprising 44, 35, and 31 genes, respectively, for each patient). Applying a Cox regression method to these patient-specific genes, along with transcriptomic and clinical data from the TARGET-ALL AALL0434 cohort, we generated three model gene signatures that efficiently identified relapsed patients within the cohort. Overall, our approach identified new potential relapse-associated genes and proposed three model gene signatures associated with lower survival rates for high-score patients.

RESUMEN

Bud Rot, caused by Phytophthora palmivora, is considered one of the main diseases affecting African oil palm (Elaeis guineensis). In this study, we investigated the in vitro molecular dynamics of the pathogen-host interaction by analyzing gene expression profiles from oil palm genotypes that were either susceptible or resistant to the disease. We observed distinct interactions of P. palmivora with resistant and susceptible oil palms through co-expression network analysis. When interacting with susceptible genotypes, P. palmivora exhibited upregulation of carbohydrate and sulfate transport genes. These genes demonstrated co-expression with apoplastic and cytoplasmic effectors, including cell wall degrading enzymes, elicitins, and RxLR motif effectors. The pathogen manipulated susceptible oil palm materials, exacerbating the response and compromising the phenylpropanoid pathway, ultimately leading to susceptibility. In contrast, resistant materials exhibited control over their response through putative Heat Shock Proteins (HSP) that maintained homeostasis between primary metabolism and biotic defense. Co-expressed genes related to flavonoids, WRKY transcripts, lectin-type receptors, and LRR receptors may play important roles in pathogen control. Overall, the study provides new knowledge of the molecular mechanisms underlying the interaction between E. guineensis and P. palmivora, which can contribute to controlling Bud Rot in oil palms and gives new insights into the interactions of P. palmivora with their hosts.

RESUMEN

Salt stress is a significant challenge that severely hampers rice growth, resulting in decreased yield and productivity. Over the years, researchers have identified biomarkers associated with salt stress to enhance rice tolerance. However, the understanding of the mechanism underlying salt tolerance in rice remains incomplete due to the involvement of multiple genes. Given the vast amount of genomics and transcriptomics data available today, it is crucial to integrate diverse datasets to identify key genes that play essential roles during salt stress in rice. In this study, we propose an integration of multiple datasets to identify potential key transcription factors. This involves utilizing network analysis based on weighted co-expression networks, focusing on gene-centric measurement and differential co-expression relationships among genes. Consequently, our analysis reveals 86 genes located in markers from previous meta-QTL analysis. Moreover, six transcription factors, namely LOC_Os03g45410 (OsTBP2), LOC_Os07g42400 (OsGATA23), LOC_Os01g13030 (OsIAA3), LOC_Os05g34050 (OsbZIP39), LOC_Os09g29930 (OsBIM1), and LOC_Os10g10990 (transcription initiation factor IIF), exhibited significantly altered co-expression relationships between salt-sensitive and salt-tolerant rice networks. These identified genes hold potential as crucial references for further investigation into the functions of salt stress response in rice plants and could be utilized in the development of salt-resistant rice cultivars. Overall, our findings shed light on the complex genetic regulation underlying salt tolerance in rice and contribute to the broader understanding of rice's response to salt stress.

Asunto(s)

RESUMEN

Background: Gene co-expression networks (GCNs) describe relationships among expressed genes key to maintaining cellular identity and homeostasis. However, the small sample size of typical RNA-seq experiments which is several orders of magnitude fewer than the number of genes is too low to infer GCNs reliably. recount3, a publicly available dataset comprised of 316,443 uniformly processed human RNA-seq samples, provides an opportunity to improve power for accurate network reconstruction and obtain biological insight from the resulting networks. Results: We compared alternate aggregation strategies to identify an optimal workflow for GCN inference by data aggregation and inferred three consensus networks: a universal network, a non-cancer network, and a cancer network in addition to 27 tissue context-specific networks. Central network genes from our consensus networks were enriched for evolutionarily constrained genes and ubiquitous biological pathways, whereas central context-specific network genes included tissue-specific transcription factors and factorization based on the hubs led to clustering of related tissue contexts. We discovered that annotations corresponding to context-specific networks inferred from aggregated data were enriched for trait heritability beyond known functional genomic annotations and were significantly more enriched when we aggregated over a larger number of samples. Conclusion: This study outlines best practices for network GCN inference and evaluation by data aggregation. We recommend estimating and regressing confounders in each data set before aggregation and prioritizing large sample size studies for GCN reconstruction. Increased statistical power in inferring context-specific networks enabled the derivation of variant annotations that were enriched for concordant trait heritability independent of functional genomic annotations that are context-agnostic. While we observed strictly increasing held-out log-likelihood with data aggregation, we noted diminishing marginal improvements. Future directions aimed at alternate methods for estimating confounders and integrating orthogonal information from modalities such as Hi-C and ChIP-seq can further improve GCN inference.

RESUMEN

Next generation risk assessment of chemicals revolves around the use of mechanistic information without animal experimentation. In this regard, toxicogenomics has proven to be a useful tool to elucidate the underlying mechanisms of adverse effects of xenobiotics. In the present study, two widely used human in vitro hepatocyte culture systems, namely primary human hepatocytes (PHH) and human hepatoma HepaRG cells, were exposed to liver toxicants known to induce liver cholestasis, steatosis or necrosis. Benchmark concentration-response modelling was applied to transcriptomics gene co-expression networks (modules) to derive benchmark concentrations (BMCs) and to gain mechanistic insight into the hepatotoxic effects. BMCs derived by concentration-response modelling of gene co-expression modules recapitulated concentration-response modelling of individual genes. Although PHH and HepaRG cells showed overlap in deregulated genes and modules by the liver toxicants, PHH demonstrated a higher responsiveness, based on the lower BMCs of co-regulated gene modules. Such BMCs can be used as transcriptomics point of departure (tPOD) for assessing module-associated cellular (stress) pathways/processes. This approach identified clear tPODs of around maximum systemic concentration (Cmax) levels for the tested drugs, while for cosmetics ingredients the BMCs were 10-100-fold higher than the estimated plasma concentrations. This approach could serve next generation risk assessment practice to identify early responsive modules at low BMCs, that could be linked to key events in liver adverse outcome pathways. In turn, this can assist in delineating potential hazards of new test chemicals using in vitro systems and used in a risk assessment when BMCs are paired with chemical exposure assessment.

Risk assessment of chemicals has traditionally been focused on animal experiments. In contrast, next generation risk assessment uses biological information obtained from experiments in cell culture models without animals to identify potential hazards. Since the liver is the main target organ of toxicity, many liver cell (hepatocyte) models have been developed and applied for hazard assessment. In this study, two widely used human hepatocyte cell models, PHH and HepaRG, were exposed to liver toxic chemicals. Biological changes in gene expression were measured in a concentration range to identify the concentration at which a biological response was perturbed using concentration response modelling. Genes belonging to the same biological process were joined based on co-expression to derive an average concentration of this process. This animal-free approach could be applied for risk assessment when biological response concentrations were related to the expected human exposure to identify potential hazard of the test chemicals.

Asunto(s)

RESUMEN

The prefrontal cortex is a crucial regulator of alcohol drinking, and dependence, and other behavioral phenotypes associated with AUD. Comprehensive identification of cell-type specific transcriptomic changes in alcohol dependence will improve our understanding of mechanisms underlying the excessive alcohol use associated with alcohol dependence and will refine targets for therapeutic development. We performed single nucleus RNA sequencing (snRNA-seq) and Visium spatial gene expression profiling on the medial prefrontal cortex (mPFC) obtained from C57BL/6 J mice exposed to the two-bottle choice-chronic intermittent ethanol (CIE) vapor exposure (2BC-CIE, defined as dependent group) paradigm which models phenotypes of alcohol dependence including escalation of alcohol drinking. Gene co-expression network analysis and differential expression analysis identified highly dysregulated co-expression networks in multiple cell types. Dysregulated modules and their hub genes suggest novel understudied targets for studying molecular mechanisms contributing to the alcohol dependence state. A subtype of inhibitory neurons was the most alcohol-sensitive cell type and contained a downregulated gene co-expression module; the hub gene for this module is Cpa6, a gene previously identified by GWAS to be associated with excessive alcohol consumption. We identified an astrocytic Gpc5 module significantly upregulated in the alcohol-dependent group. To our knowledge, there are no studies linking Cpa6 and Gpc5 to the alcohol-dependent phenotype. We also identified neuroinflammation related gene expression changes in multiple cell types, specifically enriched in microglia, further implicating neuroinflammation in the escalation of alcohol drinking. Here, we present a comprehensive atlas of cell-type specific alcohol dependence mediated gene expression changes in the mPFC and identify novel cell type-specific targets implicated in alcohol dependence.

Asunto(s)

RESUMEN

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.

Asunto(s)

RESUMEN

Gene networks allow researchers to understand the underlying mechanisms between diseases and genes while reducing the need for wet lab experiments. Numerous gene network inference (GNI) algorithms have been presented in the literature to infer accurate gene networks. We proposed a hybrid GNI algorithm, k-Strong Inference Algorithm (ksia), to infer more reliable and robust gene networks from omics datasets. To increase reliability, ksia integrates Pearson correlation coefficient (PCC) and Spearman rank correlation coefficient (SCC) scores to determine mutual information scores between molecules to increase diversity of relation predictions. To infer a more robust gene network, ksia applies three different elimination steps to remove redundant and spurious relations between genes. The performance of ksia was evaluated on microbe microarrays database in the overlap analysis with other GNI algorithms, namely ARACNE, C3NET, CLR, and MRNET. Ksia inferred less number of relations due to its strict elimination steps. However, ksia generally performed better on Escherichia coli (E.coli) and Saccharomyces cerevisiae (yeast) gene expression datasets due to F- measure and precision values. The integration of association estimator scores and three elimination stages slightly increases the performance of ksia based gene networks. Users can access ksia R package and user manual of package via https://github.com/ozgurcingiz/ksia .

RESUMEN

Gene regulatory networks (GRNs) and gene co-expression networks (GCNs) allow genome-wide exploration of molecular regulation patterns in health and disease. The standard approach for obtaining GRNs and GCNs is to infer them from gene expression data, using computational network inference methods. However, since network inference methods are usually applied on aggregate data, distortion of the networks by demographic confounders might remain undetected, especially because gene expression patterns are known to vary between different demographic groups. In this paper, we present a computational framework to systematically evaluate the influence of demographic confounders on network inference from gene expression data. Our framework compares similarities between networks inferred for different demographic groups with similarity distributions obtained for random splits of the expression data. Moreover, it allows to quantify to which extent demographic groups are represented by networks inferred from the aggregate data in a confounder-agnostic way. We apply our framework to test four widely used GRN and GCN inference methods as to their robustness w. r. t. confounding by age, ethnicity and sex in cancer. Our findings based on more than $ {44000}$ inferred networks indicate that age and sex confounders play an important role in network inference for certain cancer types, emphasizing the importance of incorporating an assessment of the effect of demographic confounders into network inference workflows. Our framework is available as a Python package on GitHub: https://github.com/bionetslab/grn-confounders.

Asunto(s)

RESUMEN

Glioblastoma multiforme (GBM) is the most aggressive form of brain tumor characterized by inter and intra-tumor heterogeneity and complex tumor microenvironment. To uncover the molecular targets in this milieu, we systematically identified immune and stromal interactions at the glial cell type level that leverages on RNA-sequencing data of GBM patients from The Cancer Genome Atlas. The perturbed genes between the high vs low immune and stromal scored patients were subjected to weighted gene co-expression network analysis to identify the glial cell type specific networks in immune and stromal infiltrated patients. The intramodular connectivity analysis identified the highly connected genes in each module. Combining it with univariable and multivariable prognostic analysis revealed common vital gene ITGB2, between the immune and stromal infiltrated patients enriched in microglia and newly formed oligodendrocytes. We found following unique hub genes in immune infiltrated patients; COL6A3 (microglia), ITGAM (oligodendrocyte precursor cells), TNFSF9 (microglia), and in stromal infiltrated patients, SERPINE1 (microglia) and THBS1 (newly formed oligodendrocytes, oligodendrocyte precursor cells). To validate these hub genes, we used external GBM patient single cell RNA-sequencing dataset and this identified ITGB2 to be significantly enriched in microglia, newly formed oligodendrocytes, T-cells, macrophages and adipocyte cell types in both immune and stromal datasets. The tumor infiltration analysis of ITGB2 showed that it is correlated with myeloid dendritic cells, macrophages, monocytes, neutrophils, B-cells, fibroblasts and adipocytes. Overall, the systematic screening of tumor microenvironment components at glial cell types uncovered ITGB2 as a potential target in primary GBM.

Asunto(s)

RESUMEN

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.

RESUMEN

Neuroblastoma represents a neoplastic expansion of neural crest cells in the developing sympathetic nervous system and is childhood's most common extracranial solid tumor. The heterogeneity of gene expression in different types of cancer is well-documented, and genetic features of neuroblastoma have been described by classification, development stage, malignancy, and progression of tumors. Here, we aim to analyze RNA sequencing datasets, publicly available in the GDC data portal, of neuroblastoma tumor samples from various patients and compare them with normal adrenal gland tissue from the GTEx data portal to elucidate the gene expression profile and regulation networks they share. Our results from the differential expression, weighted correlation network, and functional enrichment analyses that we performed with the count data from neuroblastoma and standard normal gland samples indicate that the analysis of transcriptome data from 58 patients diagnosed with high-risk neuroblastoma shares the expression pattern of 104 genes. More importantly, our analyses identify the co-expression relationship and the role of these genes in multiple biological processes and signaling pathways strongly associated with this disease phenotype. Our approach proposes a group of genes and their biological functions to be further investigated as essential molecules and possible therapeutic targets of neuroblastoma regardless of the etiology of individual tumors.