Search

Background Mosaic loss of chromosome Y (LOY) in leukocytes is the most prevalent somatic aneuploidy in aging humans. Men with LOY have increased risks of all-cause mortality and the major causes of death, including many forms of cancer. It has been suggested that the association between LOY and disease risk depends on what type of leukocyte is affected with Y loss, with prostate cancer patients showing higher levels of LOY in CD4 + T lymphocytes. In previous studies, Y loss has however been observed at relatively low levels in this cell type. This motivated us to investigate whether specific subsets of CD4 + T lymphocytes are particularly affected by LOY. Publicly available, T lymphocyte enriched, single-cell RNA sequencing datasets from patients with liver, lung or colorectal cancer were used to study how LOY affects different subtypes of T lymphocyte. To validate the observations from the public data, we also generated a single-cell RNA sequencing dataset comprised of 23 PBMC samples and 32 CD4 + T lymphocytes enriched samples. Results Regulatory T cells had significantly more LOY than any other studied T lymphocytes subtype. Furthermore, LOY in regulatory T cells increased the ratio of regulatory T cells compared with other T lymphocyte subtypes, indicating an effect of Y loss on lymphocyte differentiation. This was supported by developmental trajectory analysis of CD4 + T lymphocytes culminating in the regulatory T cells cluster most heavily affected by LOY. Finally, we identify dysregulation of 465 genes in regulatory T cells with Y loss, many involved in the immunosuppressive functions and development of regulatory T cells. Conclusions Here, we show that regulatory T cells are particularly affected by Y loss, resulting in an increased fraction of regulatory T cells and dysregulated immune functions. Considering that regulatory T cells plays a critical role in the process of immunosuppression; this enrichment for regulatory T cells with LOY might contribute to, De två första författarna delar förstaförfattarskapet

Background Mosaic loss of chromosome Y (LOY) in leukocytes is the most prevalent somatic aneuploidy in aging humans. Men with LOY have increased risks of all-cause mortality and the major causes of death, including many forms of cancer. It has been suggested that the association between LOY and disease risk depends on what type of leukocyte is affected with Y loss, with prostate cancer patients showing higher levels of LOY in CD4 + T lymphocytes. In previous studies, Y loss has however been observed at relatively low levels in this cell type. This motivated us to investigate whether specific subsets of CD4 + T lymphocytes are particularly affected by LOY. Publicly available, T lymphocyte enriched, single-cell RNA sequencing datasets from patients with liver, lung or colorectal cancer were used to study how LOY affects different subtypes of T lymphocyte. To validate the observations from the public data, we also generated a single-cell RNA sequencing dataset comprised of 23 PBMC samples and 32 CD4 + T lymphocytes enriched samples. Results Regulatory T cells had significantly more LOY than any other studied T lymphocytes subtype. Furthermore, LOY in regulatory T cells increased the ratio of regulatory T cells compared with other T lymphocyte subtypes, indicating an effect of Y loss on lymphocyte differentiation. This was supported by developmental trajectory analysis of CD4 + T lymphocytes culminating in the regulatory T cells cluster most heavily affected by LOY. Finally, we identify dysregulation of 465 genes in regulatory T cells with Y loss, many involved in the immunosuppressive functions and development of regulatory T cells. Conclusions Here, we show that regulatory T cells are particularly affected by Y loss, resulting in an increased fraction of regulatory T cells and dysregulated immune functions. Considering that regulatory T cells plays a critical role in the process of immunosuppression; this enrichment for regulatory T cells with LOY might contribute to, De två första författarna delar förstaförfattarskapet

Background Mosaic loss of chromosome Y (LOY) in leukocytes is the most prevalent somatic aneuploidy in aging humans. Men with LOY have increased risks of all-cause mortality and the major causes of death, including many forms of cancer. It has been suggested that the association between LOY and disease risk depends on what type of leukocyte is affected with Y loss, with prostate cancer patients showing higher levels of LOY in CD4 + T lymphocytes. In previous studies, Y loss has however been observed at relatively low levels in this cell type. This motivated us to investigate whether specific subsets of CD4 + T lymphocytes are particularly affected by LOY. Publicly available, T lymphocyte enriched, single-cell RNA sequencing datasets from patients with liver, lung or colorectal cancer were used to study how LOY affects different subtypes of T lymphocyte. To validate the observations from the public data, we also generated a single-cell RNA sequencing dataset comprised of 23 PBMC samples and 32 CD4 + T lymphocytes enriched samples. Results Regulatory T cells had significantly more LOY than any other studied T lymphocytes subtype. Furthermore, LOY in regulatory T cells increased the ratio of regulatory T cells compared with other T lymphocyte subtypes, indicating an effect of Y loss on lymphocyte differentiation. This was supported by developmental trajectory analysis of CD4 + T lymphocytes culminating in the regulatory T cells cluster most heavily affected by LOY. Finally, we identify dysregulation of 465 genes in regulatory T cells with Y loss, many involved in the immunosuppressive functions and development of regulatory T cells. Conclusions Here, we show that regulatory T cells are particularly affected by Y loss, resulting in an increased fraction of regulatory T cells and dysregulated immune functions. Considering that regulatory T cells plays a critical role in the process of immunosuppression; this enrichment for regulatory T cells with LOY might contribute to, De två första författarna delar förstaförfattarskapet

Although sows do not directly enter the market, they play an important role in piglet breeding on farms. They consume large amounts of feed, resulting in a significant environmental burden. Pig farms can increase their income and reduce environmental pollution by increasing the litter size (LS) of swine. PCR-RFLP/SSCP and GWAS are common methods to evaluate single-nucleotide polymorphisms (SNPs) in candidate genes. We conducted a systematic meta-analysis of the effect of SNPs on pig LS. We collected and analysed data published over the past 30 years using traditional and network meta-analyses. Trial sequential analysis (TSA) was used to analyse population data. Gene set enrichment analysis and protein–protein interaction network analysis were used to analyse the GWAS dataset. The results showed that the candidate genes were positively correlated with LS, and defects in PCR-RFLP/SSCP affected the reliability of candidate gene results. However, the genotypes with high and low LSs did not have a significant advantage. Current breeding and management practices for sows should consider increasing the LS while reducing lactation length and minimizing the sows’ non-pregnancy period as much as possible.

Background Mosaic loss of chromosome Y (LOY) in leukocytes is the most prevalent somatic aneuploidy in aging humans. Men with LOY have increased risks of all-cause mortality and the major causes of death, including many forms of cancer. It has been suggested that the association between LOY and disease risk depends on what type of leukocyte is affected with Y loss, with prostate cancer patients showing higher levels of LOY in CD4 + T lymphocytes. In previous studies, Y loss has however been observed at relatively low levels in this cell type. This motivated us to investigate whether specific subsets of CD4 + T lymphocytes are particularly affected by LOY. Publicly available, T lymphocyte enriched, single-cell RNA sequencing datasets from patients with liver, lung or colorectal cancer were used to study how LOY affects different subtypes of T lymphocyte. To validate the observations from the public data, we also generated a single-cell RNA sequencing dataset comprised of 23 PBMC samples and 32 CD4 + T lymphocytes enriched samples. Results Regulatory T cells had significantly more LOY than any other studied T lymphocytes subtype. Furthermore, LOY in regulatory T cells increased the ratio of regulatory T cells compared with other T lymphocyte subtypes, indicating an effect of Y loss on lymphocyte differentiation. This was supported by developmental trajectory analysis of CD4 + T lymphocytes culminating in the regulatory T cells cluster most heavily affected by LOY. Finally, we identify dysregulation of 465 genes in regulatory T cells with Y loss, many involved in the immunosuppressive functions and development of regulatory T cells. Conclusions Here, we show that regulatory T cells are particularly affected by Y loss, resulting in an increased fraction of regulatory T cells and dysregulated immune functions. Considering that regulatory T cells plays a critical role in the process of immunosuppression; this enrichment for regulatory T cells with LOY might contribute to, De två första författarna delar förstaförfattarskapet

The study of gene expression, the abundance of RNA transcripts in cells, has played a critical role in expanding knowledge of the molecular underpinnings of cancer and the ways in which treatments affect it. Because genes operate in interlinking pathways with nuanced and context-dependent behavior, deriving the phenotypic implication of changes in abundances of individual genes can be challenging. For this reason, Gene Set Enrichment Analysis (GSEA) was developed. GSEA is a statistical method that quantifies the activation of pathways or processes as represented by a priori annotated gene sets and plays a critical role in constructing pathway-level understandings of cell states in diseases such as cancer. In the use of GSEA and the expansion of its companion database of gene sets, the Molecular Signatures Database (MSigDB), two challenges have emerged. First, some gene sets lack the properties of context-specificity and coordinate regulation, that is, not all their members be collectively more expressed in samples that have a specific phenotype to which the gene set should correspond. Second, the new technologies for measuring gene expression at single cell resolution yield data with different properties than the “bulk” gene expression data for which GSEA was initially developed. I will address the first challenge in Chapter 1, in which I propose a data-driven method for gene set refinement and show that this method yields new gene sets that are both more coordinately regulated and context-specific. I address the second challenge in Chapter 2, in which I present a characterization of the performance of the single sample version of GSEA (ssGSEA) in multiple datasets and propose adaptations which I will show produce enrichment scores that are more stable and certain in the single cell context. Finally, Chapter 3 is an application of gene set based pathway characterization to the problem of understanding the development of chemoresistance in ovarian cancer. Using a cel

Background Mosaic loss of chromosome Y (LOY) in leukocytes is the most prevalent somatic aneuploidy in aging humans. Men with LOY have increased risks of all-cause mortality and the major causes of death, including many forms of cancer. It has been suggested that the association between LOY and disease risk depends on what type of leukocyte is affected with Y loss, with prostate cancer patients showing higher levels of LOY in CD4 + T lymphocytes. In previous studies, Y loss has however been observed at relatively low levels in this cell type. This motivated us to investigate whether specific subsets of CD4 + T lymphocytes are particularly affected by LOY. Publicly available, T lymphocyte enriched, single-cell RNA sequencing datasets from patients with liver, lung or colorectal cancer were used to study how LOY affects different subtypes of T lymphocyte. To validate the observations from the public data, we also generated a single-cell RNA sequencing dataset comprised of 23 PBMC samples and 32 CD4 + T lymphocytes enriched samples. Results Regulatory T cells had significantly more LOY than any other studied T lymphocytes subtype. Furthermore, LOY in regulatory T cells increased the ratio of regulatory T cells compared with other T lymphocyte subtypes, indicating an effect of Y loss on lymphocyte differentiation. This was supported by developmental trajectory analysis of CD4 + T lymphocytes culminating in the regulatory T cells cluster most heavily affected by LOY. Finally, we identify dysregulation of 465 genes in regulatory T cells with Y loss, many involved in the immunosuppressive functions and development of regulatory T cells. Conclusions Here, we show that regulatory T cells are particularly affected by Y loss, resulting in an increased fraction of regulatory T cells and dysregulated immune functions. Considering that regulatory T cells plays a critical role in the process of immunosuppression; this enrichment for regulatory T cells with LOY might contribute to, De två första författarna delar förstaförfattarskapet

Breast cancer stands as a prominent global health concern, consistently ranking among the most prevalent types of cancer affecting women. In this work, the focus is put on invasive ductal carcinoma and invasive lobular carcinoma, which are two most common types of breast cancer. Transfer learning techniques are applied to implement a Convolutional Neural Network (CNN) architecture that accurately classifies histological images accurately based on the subtype of cancer. Next, an integrative analysis is done with the artificial descriptors obtained from the CNN and gene expression RNA-Seq data to see the associations between the two datasets. Finally, a gene set enrichment analysis is performed with the aim of providing a biological interpretation to the artificial descriptors obtained from the CNN, by correlating them with biological pathways relevant to these two subtypes of breast cancer.

Sphingosine 1-phosphate lyase (SGPL1) insufficiency (SPLIS) is a syndrome which presents with adrenal insufficiency, steroid-resistant nephrotic syndrome, hypothyroidism, neurological disease, and ichthyosis. Where a skin phenotype is reported, 94% had abnormalities such as ichthyosis, acanthosis, and hyperpigmentation. To elucidate the disease mechanism and the role SGPL1 plays in the skin barrier we established clustered regularly interspaced short palindromic repeats-Cas9 SGPL1 KO and a lentiviral-induced SGPL1 overexpression (OE) in telomerase reverse-transcriptase immortalised human keratinocytes (N/TERT-1) and thereafter organotypic skin equivalents. Loss of SGPL1 caused an accumulation of S1P, sphingosine, and ceramides, while its overexpression caused a reduction of these species. RNAseq analysis showed perturbations in sphingolipid pathway genes, particularly in SGPL1_KO, and our gene set enrichment analysis revealed polar opposite differential gene expression between SGPL1_KO and _OE in keratinocyte differentiation and Ca2+ signaling genesets. SGPL1_KO upregulated differentiation markers, while SGPL1_OE upregulated basal and proliferative markers. The advanced differentiation of SGPL1_KO was confirmed by 3D organotypic models that also presented with a thickened and retained stratum corneum and a breakdown of E-cadherin junctions. We conclude that SPLIS associated ichthyosis is a multifaceted disease caused possibly by sphingolipid imbalance and excessive S1P signaling, leading to increased differentiation and an imbalance of the lipid lamellae throughout the epidermis.

Group-aggregated responses to tuberculosis (TB) have been well characterized on a molecular level. However, human beings differ and individual responses to infection vary. We have combined a novel approach to individual gene set analysis (GSA) with the clustering of transcriptomic profiles of TB patients from seven datasets in order to identify individual molecular endotypes of transcriptomic responses to TB. We found that TB patients differ with respect to the intensity of their hallmark interferon (IFN) responses, but they also show variability in their complement system, metabolic responses and multiple other pathways. This variability cannot be sufficiently explained with covariates such as gender or age, and the molecular endotypes are found across studies and populations. Using datasets from a Cynomolgus macaque model of TB, we revealed that transcriptional signatures of different molecular TB endotypes did not depend on TB progression post-infection. Moreover, we provide evidence that patients with molecular endotypes characterized by high levels of IFN responses (IFN-rich), suffered from more severe lung pathology than those with lower levels of IFN responses (IFN-low). Harnessing machine learning (ML) models, we derived gene signatures classifying IFN-rich and IFN-low TB endotypes and revealed that the IFN-low signature allowed slightly more reliable overall classification of TB patients from non-TB patients than the IFN-rich one. Using the paradigm of molecular endotypes and the ML-based predictions allows more precisely tailored treatment regimens, predicting treatment-outcome with higher accuracy and therefore bridging the gap between conventional treatment and precision medicine., Peer Reviewed

Group-aggregated responses to tuberculosis (TB) have been well characterized on a molecular level. However, human beings differ and individual responses to infection vary. We have combined a novel approach to individual gene set analysis (GSA) with the clustering of transcriptomic profiles of TB patients from seven datasets in order to identify individual molecular endotypes of transcriptomic responses to TB. We found that TB patients differ with respect to the intensity of their hallmark interferon (IFN) responses, but they also show variability in their complement system, metabolic responses and multiple other pathways. This variability cannot be sufficiently explained with covariates such as gender or age, and the molecular endotypes are found across studies and populations. Using datasets from a Cynomolgus macaque model of TB, we revealed that transcriptional signatures of different molecular TB endotypes did not depend on TB progression post-infection. Moreover, we provide evidence that patients with molecular endotypes characterized by high levels of IFN responses (IFN-rich), suffered from more severe lung pathology than those with lower levels of IFN responses (IFN-low). Harnessing machine learning (ML) models, we derived gene signatures classifying IFN-rich and IFN-low TB endotypes and revealed that the IFN-low signature allowed slightly more reliable overall classification of TB patients from non-TB patients than the IFN-rich one. Using the paradigm of molecular endotypes and the ML-based predictions allows more precisely tailored treatment regimens, predicting treatment-outcome with higher accuracy and therefore bridging the gap between conventional treatment and precision medicine., Peer Reviewed

Транскриптоміка — це сукупність методів, що займаються аналізом даних експресії генів. Більшість ранніх транскриптомних методів цікавилася широкою групою генів. На додаток до проведення дослідницьких аналізів, заснованих на нормалізації рівня експресії, було спрямовано чимало зусиль на оптимізацію статистичних тестів, що вирішують, чи є різниця в експресії вибраного гену між двома (або більше) умовами (фенотипами). Такий аналіз називається аналізом диференціальної експресії (Differential Gene Expression analysis, DGE). Двома основними завданнями всіх інструментів DGE є: 1. Оцінити величину диференціальної експресії між двома або більше умовами на основі підрахунку читання від повторених зразків, тобто обчислити зміну кількості зчитувань у ході секвенування. 2. Оцінити статистичну значущість різниці. Оглянувши типові види візуалізацій, ми можемо встановити, що існуючі методи: складні в розумінні та вимагають значного поглиблення для повного розуміння. Зазвичай не дозволяють функціонально оцінити набори генів на генному рівні. Я створив нову візуалізацію, яка дозволяє це зробити. Крім того, було створено програмне забезпечення для відображення цієї та інших відомих візуалізацій., Transcriptomics consists of an array of methods that involve gene expression analysis. Most of the early transcriptomics methods were interested in a broad set of genes. In addition to research analyses based on expression level normalization, a lot of effort was given to the statistical tests optimization that decide whether there is a difference in expression of a given gene between two or more phenotypes. This analysis is called Differential Gene Expression (DGE) analysis. Two main goals of DGE are: 1. To estimate size of differential expression between two or more conditions based on calculations of used samples 2. To estimate the statistical significance of this change. Looking at the typical visualisations, we can establish that existing methods are: 1. Hard to understand and require significant dive to fully comprehend. 2. Don’t allow to functionally evaluate gene sets on a gene level. I have created a novel visualization that allows me to fix these problems. Besides that, new software was created to display this and other known visualizations.

BACKGROUND:Executive function (EF) is vital to human beings. It has been linked to many genes and family environmental factors in separate studies, but few studies have examined the potential interactions between gene(s) and environmental factor(s). The current study explored the whole genome to identify SNPs, genes, and pathways that interacted with parental warmth (PW) on EF. RESULTS:Nine EF tasks were used to measure its three components (common EF, updating, shifting) based on the model proposed by Miyake et al. (2000). We found that rs111605473, LAMP5, SLC4A7, and LRRK1 interacted significantly with PW to affect the updating component of EF, and the GSE43955 pathway interacted significantly with PW to affect the common EF component. CONCLUSIONS:The current study is the first to identify genes that interacted with PW to affect EF. Further studies are needed to reveal the underlying mechanism.

Background: Fetal growth restriction (FGR) is associated with an increased susceptibility for various noncommunicable diseases in adulthood, including cardiovascular and renal disease. During FGR, reduced uteroplacental blood flow, oxygen and nutrient supply to the fetus are hypothesized to detrimentally influence cardiovascular and renal programming. This study examined whether developmental programming profiles, especially related to the cardiovascular and renal system, differ in human umbilical vein endothelial cells (HUVECs) collected from pregnancies complicated by placental insufficiency-induced FGR compared to normal growth pregnancies. Our approach, involving transcriptomic profiling by RNA-sequencing and gene set enrichment analysis focused on cardiovascular and renal gene sets and targeted DNA methylation assays, contributes to the identification of targets underlying long-term cardiovascular and renal diseases. Results: Gene set enrichment analysis showed several downregulated gene sets, most of them involved in immune or inflammatory pathways or cell cycle pathways. seven of the 22 significantly upregulated gene sets related to kidney development and four gene sets involved with cardiovascular health and function were downregulated in FGR (n = 11) versus control (n = 8). Transcriptomic profiling by RNA-sequencing revealed downregulated expression of LGALS1, FPR3 and NRM and upregulation of lincRNA RP5-855F14.1 in FGR compared to controls. DNA methylation was similar for LGALS1 between study groups, but relative hypomethylation of FPR3 and hypermethylation of NRM were present in FGR, especially in male offspring. Absolute differences in methylation were, however, small. Conclusion: This study showed upregulation of gene sets related to renal development in HUVECs collected from pregnancies complicated by FGR compared to control donors. The differentially expressed gene sets related to cardiovascular function and health might be in line with the downregula

BACKGROUND:Executive function (EF) is vital to human beings. It has been linked to many genes and family environmental factors in separate studies, but few studies have examined the potential interactions between gene(s) and environmental factor(s). The current study explored the whole genome to identify SNPs, genes, and pathways that interacted with parental warmth (PW) on EF. RESULTS:Nine EF tasks were used to measure its three components (common EF, updating, shifting) based on the model proposed by Miyake et al. (2000). We found that rs111605473, LAMP5, SLC4A7, and LRRK1 interacted significantly with PW to affect the updating component of EF, and the GSE43955 pathway interacted significantly with PW to affect the common EF component. CONCLUSIONS:The current study is the first to identify genes that interacted with PW to affect EF. Further studies are needed to reveal the underlying mechanism.

Gall wasps (Hymenoptera: Cynipidae) induce complex galls on oaks, roses, and other plants, but the mechanism of gall induction is still unknown. Here, we take a comparative genomic approach to revealing the genetic basis of gall induction. We focus on Synergus itoensis, a species that induces galls inside oak acorns. Previous studies suggested that this species evolved the ability to initiate gall formation recently, as it is deeply nested within the genus Synergus, whose members are mostly inquilines that develop inside the galls of other species. We compared the genome of S. itoensis with that of three related Synergus inquilines to identify genomic changes associated with the origin of gall induction. We used a novel Bayesian selection analysis, which accounts for branch-specific and gene-specific selection effects, to search for signatures of selection in 7,600 single-copy orthologous genes shared by the four Synergus species. We found that the terminal branch leading to S. itoensis had more genes with a significantly elevated dN/dS ratio (positive signature genes) than the other terminal branches in the tree; the S. itoensis branch also had more genes with a significantly decreased dN/dS ratio. Gene set enrichment analysis showed that the positive signature gene set of S. itoensis, unlike those of the inquiline species, is enriched in several biological process Gene Ontology terms, the most prominent of which is “Ovarian Follicle Cell Development.” Our results indicate that the origin of gall induction is associated with distinct genomic changes, and provide a good starting point for further characterization of the genes involved.

Bone metastasis is the lethal end-stage of prostate cancer (PC), but the biology of bone metastases is poorly understood. The overall aim of this study was therefore to explore molecular variability in PC bone metastases of potential importance for therapy. Specifically, genome-wide expression profiles of bone metastases from untreated patients (n = 12) and patients treated with androgen-deprivation therapy (ADT, n = 60) were analyzed in relation to patient outcome and to morphological characteristics in metastases and paired primary tumors. Principal component analysis and unsupervised classification were used to identify sample clusters based on mRNA profiles. Clusters were characterized by gene set enrichment analysis and related to histological and clinical parameters using univariate and multivariate statistics. Selected proteins were analyzed by immunohistochemistry in metastases and matched primary tumors (n = 52) and in transurethral resected prostate (TUR-P) tissue of a separate cohort (n = 59). Three molecular subtypes of bone metastases (MetA-C) characterized by differences in gene expression pattern, morphology, and clinical behavior were identified. MetA (71% of the cases) showed increased expression of androgen receptor-regulated genes, including prostate-specific antigen (PSA), and glandular structures indicating a luminal cell phenotype. MetB (17%) showed expression profiles related to cell cycle activity and DNA damage, and a pronounced cellular atypia. MetC (12%) exhibited enriched stroma-epithelial cell interactions. MetB patients had the lowest serum PSA levels and the poorest prognosis after ADT. Combined analysis of PSA and Ki67 immunoreactivity (proliferation) in bone metastases, paired primary tumors, and TUR-P samples was able to differentiate MetA-like (high PSA, low Ki67) from MetB-like (low PSA, high Ki67) tumors and demonstrate their different prognosis. In conclusion, bone metastases from PC patients are separated based on gene expressi

Bone metastasis is the lethal end-stage of prostate cancer (PC), but the biology of bone metastases is poorly understood. The overall aim of this study was therefore to explore molecular variability in PC bone metastases of potential importance for therapy. Specifically, genome-wide expression profiles of bone metastases from untreated patients (n = 12) and patients treated with androgen-deprivation therapy (ADT, n = 60) were analyzed in relation to patient outcome and to morphological characteristics in metastases and paired primary tumors. Principal component analysis and unsupervised classification were used to identify sample clusters based on mRNA profiles. Clusters were characterized by gene set enrichment analysis and related to histological and clinical parameters using univariate and multivariate statistics. Selected proteins were analyzed by immunohistochemistry in metastases and matched primary tumors (n = 52) and in transurethral resected prostate (TUR-P) tissue of a separate cohort (n = 59). Three molecular subtypes of bone metastases (MetA-C) characterized by differences in gene expression pattern, morphology, and clinical behavior were identified. MetA (71% of the cases) showed increased expression of androgen receptor-regulated genes, including prostate-specific antigen (PSA), and glandular structures indicating a luminal cell phenotype. MetB (17%) showed expression profiles related to cell cycle activity and DNA damage, and a pronounced cellular atypia. MetC (12%) exhibited enriched stroma-epithelial cell interactions. MetB patients had the lowest serum PSA levels and the poorest prognosis after ADT. Combined analysis of PSA and Ki67 immunoreactivity (proliferation) in bone metastases, paired primary tumors, and TUR-P samples was able to differentiate MetA-like (high PSA, low Ki67) from MetB-like (low PSA, high Ki67) tumors and demonstrate their different prognosis. In conclusion, bone metastases from PC patients are separated based on gene expressi

Bone metastasis is the lethal end-stage of prostate cancer (PC), but the biology of bone metastases is poorly understood. The overall aim of this study was therefore to explore molecular variability in PC bone metastases of potential importance for therapy. Specifically, genome-wide expression profiles of bone metastases from untreated patients (n = 12) and patients treated with androgen-deprivation therapy (ADT, n = 60) were analyzed in relation to patient outcome and to morphological characteristics in metastases and paired primary tumors. Principal component analysis and unsupervised classification were used to identify sample clusters based on mRNA profiles. Clusters were characterized by gene set enrichment analysis and related to histological and clinical parameters using univariate and multivariate statistics. Selected proteins were analyzed by immunohistochemistry in metastases and matched primary tumors (n = 52) and in transurethral resected prostate (TUR-P) tissue of a separate cohort (n = 59). Three molecular subtypes of bone metastases (MetA-C) characterized by differences in gene expression pattern, morphology, and clinical behavior were identified. MetA (71% of the cases) showed increased expression of androgen receptor-regulated genes, including prostate-specific antigen (PSA), and glandular structures indicating a luminal cell phenotype. MetB (17%) showed expression profiles related to cell cycle activity and DNA damage, and a pronounced cellular atypia. MetC (12%) exhibited enriched stroma-epithelial cell interactions. MetB patients had the lowest serum PSA levels and the poorest prognosis after ADT. Combined analysis of PSA and Ki67 immunoreactivity (proliferation) in bone metastases, paired primary tumors, and TUR-P samples was able to differentiate MetA-like (high PSA, low Ki67) from MetB-like (low PSA, high Ki67) tumors and demonstrate their different prognosis. In conclusion, bone metastases from PC patients are separated based on gene expressi

Bone metastasis is the lethal end-stage of prostate cancer (PC), but the biology of bone metastases is poorly understood. The overall aim of this study was therefore to explore molecular variability in PC bone metastases of potential importance for therapy. Specifically, genome-wide expression profiles of bone metastases from untreated patients (n = 12) and patients treated with androgen-deprivation therapy (ADT, n = 60) were analyzed in relation to patient outcome and to morphological characteristics in metastases and paired primary tumors. Principal component analysis and unsupervised classification were used to identify sample clusters based on mRNA profiles. Clusters were characterized by gene set enrichment analysis and related to histological and clinical parameters using univariate and multivariate statistics. Selected proteins were analyzed by immunohistochemistry in metastases and matched primary tumors (n = 52) and in transurethral resected prostate (TUR-P) tissue of a separate cohort (n = 59). Three molecular subtypes of bone metastases (MetA-C) characterized by differences in gene expression pattern, morphology, and clinical behavior were identified. MetA (71% of the cases) showed increased expression of androgen receptor-regulated genes, including prostate-specific antigen (PSA), and glandular structures indicating a luminal cell phenotype. MetB (17%) showed expression profiles related to cell cycle activity and DNA damage, and a pronounced cellular atypia. MetC (12%) exhibited enriched stroma-epithelial cell interactions. MetB patients had the lowest serum PSA levels and the poorest prognosis after ADT. Combined analysis of PSA and Ki67 immunoreactivity (proliferation) in bone metastases, paired primary tumors, and TUR-P samples was able to differentiate MetA-like (high PSA, low Ki67) from MetB-like (low PSA, high Ki67) tumors and demonstrate their different prognosis. In conclusion, bone metastases from PC patients are separated based on gene expressi

Bone metastasis is the lethal end-stage of prostate cancer (PC), but the biology of bone metastases is poorly understood. The overall aim of this study was therefore to explore molecular variability in PC bone metastases of potential importance for therapy. Specifically, genome-wide expression profiles of bone metastases from untreated patients (n = 12) and patients treated with androgen-deprivation therapy (ADT, n = 60) were analyzed in relation to patient outcome and to morphological characteristics in metastases and paired primary tumors. Principal component analysis and unsupervised classification were used to identify sample clusters based on mRNA profiles. Clusters were characterized by gene set enrichment analysis and related to histological and clinical parameters using univariate and multivariate statistics. Selected proteins were analyzed by immunohistochemistry in metastases and matched primary tumors (n = 52) and in transurethral resected prostate (TUR-P) tissue of a separate cohort (n = 59). Three molecular subtypes of bone metastases (MetA-C) characterized by differences in gene expression pattern, morphology, and clinical behavior were identified. MetA (71% of the cases) showed increased expression of androgen receptor-regulated genes, including prostate-specific antigen (PSA), and glandular structures indicating a luminal cell phenotype. MetB (17%) showed expression profiles related to cell cycle activity and DNA damage, and a pronounced cellular atypia. MetC (12%) exhibited enriched stroma-epithelial cell interactions. MetB patients had the lowest serum PSA levels and the poorest prognosis after ADT. Combined analysis of PSA and Ki67 immunoreactivity (proliferation) in bone metastases, paired primary tumors, and TUR-P samples was able to differentiate MetA-like (high PSA, low Ki67) from MetB-like (low PSA, high Ki67) tumors and demonstrate their different prognosis. In conclusion, bone metastases from PC patients are separated based on gene expressi

Bone metastasis is the lethal end-stage of prostate cancer (PC), but the biology of bone metastases is poorly understood. The overall aim of this study was therefore to explore molecular variability in PC bone metastases of potential importance for therapy. Specifically, genome-wide expression profiles of bone metastases from untreated patients (n = 12) and patients treated with androgen-deprivation therapy (ADT, n = 60) were analyzed in relation to patient outcome and to morphological characteristics in metastases and paired primary tumors. Principal component analysis and unsupervised classification were used to identify sample clusters based on mRNA profiles. Clusters were characterized by gene set enrichment analysis and related to histological and clinical parameters using univariate and multivariate statistics. Selected proteins were analyzed by immunohistochemistry in metastases and matched primary tumors (n = 52) and in transurethral resected prostate (TUR-P) tissue of a separate cohort (n = 59). Three molecular subtypes of bone metastases (MetA-C) characterized by differences in gene expression pattern, morphology, and clinical behavior were identified. MetA (71% of the cases) showed increased expression of androgen receptor-regulated genes, including prostate-specific antigen (PSA), and glandular structures indicating a luminal cell phenotype. MetB (17%) showed expression profiles related to cell cycle activity and DNA damage, and a pronounced cellular atypia. MetC (12%) exhibited enriched stroma-epithelial cell interactions. MetB patients had the lowest serum PSA levels and the poorest prognosis after ADT. Combined analysis of PSA and Ki67 immunoreactivity (proliferation) in bone metastases, paired primary tumors, and TUR-P samples was able to differentiate MetA-like (high PSA, low Ki67) from MetB-like (low PSA, high Ki67) tumors and demonstrate their different prognosis. In conclusion, bone metastases from PC patients are separated based on gene expressi

Bone metastasis is the lethal end-stage of prostate cancer (PC), but the biology of bone metastases is poorly understood. The overall aim of this study was therefore to explore molecular variability in PC bone metastases of potential importance for therapy. Specifically, genome-wide expression profiles of bone metastases from untreated patients (n = 12) and patients treated with androgen-deprivation therapy (ADT, n = 60) were analyzed in relation to patient outcome and to morphological characteristics in metastases and paired primary tumors. Principal component analysis and unsupervised classification were used to identify sample clusters based on mRNA profiles. Clusters were characterized by gene set enrichment analysis and related to histological and clinical parameters using univariate and multivariate statistics. Selected proteins were analyzed by immunohistochemistry in metastases and matched primary tumors (n = 52) and in transurethral resected prostate (TUR-P) tissue of a separate cohort (n = 59). Three molecular subtypes of bone metastases (MetA-C) characterized by differences in gene expression pattern, morphology, and clinical behavior were identified. MetA (71% of the cases) showed increased expression of androgen receptor-regulated genes, including prostate-specific antigen (PSA), and glandular structures indicating a luminal cell phenotype. MetB (17%) showed expression profiles related to cell cycle activity and DNA damage, and a pronounced cellular atypia. MetC (12%) exhibited enriched stroma-epithelial cell interactions. MetB patients had the lowest serum PSA levels and the poorest prognosis after ADT. Combined analysis of PSA and Ki67 immunoreactivity (proliferation) in bone metastases, paired primary tumors, and TUR-P samples was able to differentiate MetA-like (high PSA, low Ki67) from MetB-like (low PSA, high Ki67) tumors and demonstrate their different prognosis. In conclusion, bone metastases from PC patients are separated based on gene expressi

Bone metastasis is the lethal end-stage of prostate cancer (PC), but the biology of bone metastases is poorly understood. The overall aim of this study was therefore to explore molecular variability in PC bone metastases of potential importance for therapy. Specifically, genome-wide expression profiles of bone metastases from untreated patients (n = 12) and patients treated with androgen-deprivation therapy (ADT, n = 60) were analyzed in relation to patient outcome and to morphological characteristics in metastases and paired primary tumors. Principal component analysis and unsupervised classification were used to identify sample clusters based on mRNA profiles. Clusters were characterized by gene set enrichment analysis and related to histological and clinical parameters using univariate and multivariate statistics. Selected proteins were analyzed by immunohistochemistry in metastases and matched primary tumors (n = 52) and in transurethral resected prostate (TUR-P) tissue of a separate cohort (n = 59). Three molecular subtypes of bone metastases (MetA-C) characterized by differences in gene expression pattern, morphology, and clinical behavior were identified. MetA (71% of the cases) showed increased expression of androgen receptor-regulated genes, including prostate-specific antigen (PSA), and glandular structures indicating a luminal cell phenotype. MetB (17%) showed expression profiles related to cell cycle activity and DNA damage, and a pronounced cellular atypia. MetC (12%) exhibited enriched stroma-epithelial cell interactions. MetB patients had the lowest serum PSA levels and the poorest prognosis after ADT. Combined analysis of PSA and Ki67 immunoreactivity (proliferation) in bone metastases, paired primary tumors, and TUR-P samples was able to differentiate MetA-like (high PSA, low Ki67) from MetB-like (low PSA, high Ki67) tumors and demonstrate their different prognosis. In conclusion, bone metastases from PC patients are separated based on gene expressi

Bone metastasis is the lethal end-stage of prostate cancer (PC), but the biology of bone metastases is poorly understood. The overall aim of this study was therefore to explore molecular variability in PC bone metastases of potential importance for therapy. Specifically, genome-wide expression profiles of bone metastases from untreated patients (n = 12) and patients treated with androgen-deprivation therapy (ADT, n = 60) were analyzed in relation to patient outcome and to morphological characteristics in metastases and paired primary tumors. Principal component analysis and unsupervised classification were used to identify sample clusters based on mRNA profiles. Clusters were characterized by gene set enrichment analysis and related to histological and clinical parameters using univariate and multivariate statistics. Selected proteins were analyzed by immunohistochemistry in metastases and matched primary tumors (n = 52) and in transurethral resected prostate (TUR-P) tissue of a separate cohort (n = 59). Three molecular subtypes of bone metastases (MetA-C) characterized by differences in gene expression pattern, morphology, and clinical behavior were identified. MetA (71% of the cases) showed increased expression of androgen receptor-regulated genes, including prostate-specific antigen (PSA), and glandular structures indicating a luminal cell phenotype. MetB (17%) showed expression profiles related to cell cycle activity and DNA damage, and a pronounced cellular atypia. MetC (12%) exhibited enriched stroma-epithelial cell interactions. MetB patients had the lowest serum PSA levels and the poorest prognosis after ADT. Combined analysis of PSA and Ki67 immunoreactivity (proliferation) in bone metastases, paired primary tumors, and TUR-P samples was able to differentiate MetA-like (high PSA, low Ki67) from MetB-like (low PSA, high Ki67) tumors and demonstrate their different prognosis. In conclusion, bone metastases from PC patients are separated based on gene expressi

Background: Unsupervised machine learning methods (deep learning) have shown their usefulness with noisy single cell mRNA-sequencing data (scRNA-seq), where the models generalize well, despite the zero-inflation of the data. A class of neural networks, namely autoencoders, has been useful for denoising of single cell data, imputation of missing values and dimensionality reduction. Results: Here, we present a striking feature with the potential to greatly increase the usability of autoencoders: With specialized training, the autoencoder is not only able to generalize over the data, but also to tease apart biologically meaningful modules, which we found encoded in the representation layer of the network. Our model can, from scRNA-seq data, delineate biological meaningful modules that govern a dataset, as well as give information as to which modules are active in each single cell. Importantly, most of these modules can be explained by known biological functions, as provided by the Hallmark gene sets. Conclusions: We discover that tailored training of an autoencoder makes it possible to deconvolute biological modules inherent in the data, without any assumptions. By comparisons with gene signatures of canonical pathways we see that the modules are directly interpretable. The scope of this discovery has important implications, as it makes it possible to outline the drivers behind a given effect of a cell. In comparison with other dimensionality reduction methods, or supervised models for classification, our approach has the benefit of both handling well the zero-inflated nature of scRNA-seq, and validating that the model captures relevant information, by establishing a link between input and decoded data. In perspective, our model in combination with clustering methods is able to provide information about which subtype a given single cell belongs to, as well as which biological functions determine that membership.

Bone metastasis is the lethal end-stage of prostate cancer (PC), but the biology of bone metastases is poorly understood. The overall aim of this study was therefore to explore molecular variability in PC bone metastases of potential importance for therapy. Specifically, genome-wide expression profiles of bone metastases from untreated patients (n = 12) and patients treated with androgen-deprivation therapy (ADT, n = 60) were analyzed in relation to patient outcome and to morphological characteristics in metastases and paired primary tumors. Principal component analysis and unsupervised classification were used to identify sample clusters based on mRNA profiles. Clusters were characterized by gene set enrichment analysis and related to histological and clinical parameters using univariate and multivariate statistics. Selected proteins were analyzed by immunohistochemistry in metastases and matched primary tumors (n = 52) and in transurethral resected prostate (TUR-P) tissue of a separate cohort (n = 59). Three molecular subtypes of bone metastases (MetA-C) characterized by differences in gene expression pattern, morphology, and clinical behavior were identified. MetA (71% of the cases) showed increased expression of androgen receptor-regulated genes, including prostate-specific antigen (PSA), and glandular structures indicating a luminal cell phenotype. MetB (17%) showed expression profiles related to cell cycle activity and DNA damage, and a pronounced cellular atypia. MetC (12%) exhibited enriched stroma-epithelial cell interactions. MetB patients had the lowest serum PSA levels and the poorest prognosis after ADT. Combined analysis of PSA and Ki67 immunoreactivity (proliferation) in bone metastases, paired primary tumors, and TUR-P samples was able to differentiate MetA-like (high PSA, low Ki67) from MetB-like (low PSA, high Ki67) tumors and demonstrate their different prognosis. In conclusion, bone metastases from PC patients are separated based on gene expressi

Pathway analysis is an important task for gaining novel insights into the molecular architecture of many complex diseases. With the advancement of new sequencing technologies, a large amount of quantitative gene expression data have been continuously acquired. The springing up omics data sets such as proteomics has facilitated the investigation on disease relevant pathways. Although much work has previously been done to explore the single omics data, little work has been reported using multi-omics data integration, mainly due to methodological and technological limitations. While a single omic data can provide useful information about the underlying biological processes, multi-omics data integration would be much more comprehensive about the cause-effect processes responsible for diseases and their subtypes. This project investigates the combination of miRNAseq, proteomics, and RNAseq data on seven types of muscular dystrophies and control group. These unique multi-omics data sets provide us with the opportunity to identify disease-specific and most relevant biological pathways. We first perform t-test and OVEPUG test separately to define the differential expressed genes in protein and mRNA data sets. In multi-omics data sets, miRNA also plays a significant role in muscle development by regulating their target genes in mRNA dataset. To exploit the relationship between miRNA and gene expression, we consult with the commonly used gene library - Targetscan to collect all paired miRNA-mRNA and miRNA-protein co-expression pairs. Next, by conducting statistical analysis such as Pearson's correlation coefficient or t-test, we measured the biologically expected correlation of each gene with its upstream miRNAs and identify those showing negative correlation between the aforementioned miRNA-mRNA and miRNA-protein pairs. Furthermore, we identify and assess the most relevant disease-specific pathways by inputting the differential expressed genes and negative correlated genes i

Integrative analysis based on quantitative representation of whole slide images (WSIs) in a large histology cohort may provide predictive models of clinical outcome. On one hand, the efficiency and effectiveness of such representation is hindered as a result of large technical variations (e.g., fixation, staining) and biological heterogeneities (e.g., cell type, cell state) that are always present in a large cohort. On the other hand, perceptual interpretation/validation of important multivariate phenotypic signatures are often difficult due to the loss of visual information during feature transformation in hyperspace. To address these issues, we propose a novel approach for integrative analysis based on cellular morphometric context, which is a robust representation of WSI, with the emphasis on tumor architecture and tumor heterogeneity, built upon cellular level morphometric features within the spatial pyramid matching (SPM) framework. The proposed approach is applied to The Cancer Genome Atlas (TCGA) lower grade glioma (LGG) cohort, where experimental results (i) reveal several clinically relevant cellular morphometric types, which enables both perceptual interpretation/validation and further investigation through gene set enrichment analysis; and (ii) indicate the significantly increased survival rates in one of the cellular morphometric context subtypes derived from the cellular morphometric context.

Gene set enrichment analysis aims to discover sets of genes, such as biological pathways or protein complexes, which may show moderate but coordinated differentiation across experimental conditions. The existing gene set enrichment approaches utilize single gene statistic as a measure of differentiation for individual genes. These approaches do not utilize any inter-gene correlations, but it has been known that genes in a pathway often interact with each other. Motivated by the need for taking gene dependence into account, we propose a novel gene set enrichment algorithm, where the gene-gene correlation is addressed via a gene-pair representation strategy. Relying on an appropriately defined gene pair statistic, the gene set statistic is formulated using a competitive null hypothesis. Extensive simulation studies show that our proposed approach can correctly control the type I error (false positive rate), and retain good statistical power for detecting true differential expression. The new method is also applied to analyze several gene expression datasets.

UnlabelledGene set enrichment analysis (GSEA) approaches are widely used to identify coordinately regulated genes associated with phenotypes of interest. Here, we present Constellation Map, a tool to visualize and interpret the results when enrichment analyses yield a long list of significantly enriched gene sets. Constellation Map identifies commonalities that explain the enrichment of multiple top-scoring gene sets and maps the relationships between them. Constellation Map can help investigators take full advantage of GSEA and facilitates the biological interpretation of enrichment results.AvailabilityConstellation Map is freely available as a GenePattern module at http://www.genepattern.org.

UnlabelledGene set enrichment analysis (GSEA) approaches are widely used to identify coordinately regulated genes associated with phenotypes of interest. Here, we present Constellation Map, a tool to visualize and interpret the results when enrichment analyses yield a long list of significantly enriched gene sets. Constellation Map identifies commonalities that explain the enrichment of multiple top-scoring gene sets and maps the relationships between them. Constellation Map can help investigators take full advantage of GSEA and facilitates the biological interpretation of enrichment results.AvailabilityConstellation Map is freely available as a GenePattern module at http://www.genepattern.org.

The Molecular Signatures Database (MSigDB) is one of the most widely used and comprehensive databases of gene sets for performing gene set enrichment analysis. Since its creation, MSigDB has grown beyond its roots in metabolic disease and cancer to include >10,000 gene sets. These better represent a wider range of biological processes and diseases, but the utility of the database is reduced by increased redundancy across, and heterogeneity within, gene sets. To address this challenge, here we use a combination of automated approaches and expert curation to develop a collection of "hallmark" gene sets as part of MSigDB. Each hallmark in this collection consists of a "refined" gene set, derived from multiple "founder" sets, that conveys a specific biological state or process and displays coherent expression. The hallmarks effectively summarize most of the relevant information of the original founder sets and, by reducing both variation and redundancy, provide more refined and concise inputs for gene set enrichment analysis.

Circadian rhythms control metabolism and energy homeostasis, but the role of the skeletal muscle clock has never been explored. We generated conditional and inducible mouse lines with muscle-specific ablation of the core clock gene Bmal1. Skeletal muscles from these mice showed impaired insulin-stimulated glucose uptake with reduced protein levels of GLUT4, the insulin-dependent glucose transporter, and TBC1D1, a Rab-GTPase involved in GLUT4 translocation. Pyruvate dehydrogenase (PDH) activity was also reduced due to altered expression of circadian genes Pdk4 and Pdp1, coding for PDH kinase and phosphatase, respectively. PDH inhibition leads to reduced glucose oxidation and diversion of glycolytic intermediates to alternative metabolic pathways, as revealed by metabolome analysis. The impaired glucose metabolism induced by muscle-specific Bmal1 knockout suggests that a major physiological role of the muscle clock is to prepare for the transition from the rest/fasting phase to the active/feeding phase, when glucose becomes the predominant fuel for skeletal muscle.

Circadian rhythms control metabolism and energy homeostasis, but the role of the skeletal muscle clock has never been explored. We generated conditional and inducible mouse lines with muscle-specific ablation of the core clock gene Bmal1. Skeletal muscles from these mice showed impaired insulin-stimulated glucose uptake with reduced protein levels of GLUT4, the insulin-dependent glucose transporter, and TBC1D1, a Rab-GTPase involved in GLUT4 translocation. Pyruvate dehydrogenase (PDH) activity was also reduced due to altered expression of circadian genes Pdk4 and Pdp1, coding for PDH kinase and phosphatase, respectively. PDH inhibition leads to reduced glucose oxidation and diversion of glycolytic intermediates to alternative metabolic pathways, as revealed by metabolome analysis. The impaired glucose metabolism induced by muscle-specific Bmal1 knockout suggests that a major physiological role of the muscle clock is to prepare for the transition from the rest/fasting phase to the active/feeding phase, when glucose becomes the predominant fuel for skeletal muscle.

The mechanisms generating epileptic neuronal networks following insults such as severe seizures are unknown. We have previously shown that interfering with the function of the neuron-restrictive silencer factor (NRSF/REST), an important transcription factor that influences neuronal phenotype, attenuated development of this disorder. In this study, we found that epilepsy-provoking seizures increased the low NRSF levels in mature hippocampus several fold yet surprisingly, provoked repression of only a subset (∼10%) of potential NRSF target genes. Accordingly, the repressed gene-set was rescued when NRSF binding to chromatin was blocked. Unexpectedly, genes selectively repressed by NRSF had mid-range binding frequencies to the repressor, a property that rendered them sensitive to moderate fluctuations of NRSF levels. Genes selectively regulated by NRSF during epileptogenesis coded for ion channels, receptors, and other crucial contributors to neuronal function. Thus, dynamic, selective regulation of NRSF target genes may play a role in influencing neuronal properties in pathological and physiological contexts.

Circadian rhythms control metabolism and energy homeostasis, but the role of the skeletal muscle clock has never been explored. We generated conditional and inducible mouse lines with muscle-specific ablation of the core clock gene Bmal1. Skeletal muscles from these mice showed impaired insulin-stimulated glucose uptake with reduced protein levels of GLUT4, the insulin-dependent glucose transporter, and TBC1D1, a Rab-GTPase involved in GLUT4 translocation. Pyruvate dehydrogenase (PDH) activity was also reduced due to altered expression of circadian genes Pdk4 and Pdp1, coding for PDH kinase and phosphatase, respectively. PDH inhibition leads to reduced glucose oxidation and diversion of glycolytic intermediates to alternative metabolic pathways, as revealed by metabolome analysis. The impaired glucose metabolism induced by muscle-specific Bmal1 knockout suggests that a major physiological role of the muscle clock is to prepare for the transition from the rest/fasting phase to the active/feeding phase, when glucose becomes the predominant fuel for skeletal muscle.

The mechanisms generating epileptic neuronal networks following insults such as severe seizures are unknown. We have previously shown that interfering with the function of the neuron-restrictive silencer factor (NRSF/REST), an important transcription factor that influences neuronal phenotype, attenuated development of this disorder. In this study, we found that epilepsy-provoking seizures increased the low NRSF levels in mature hippocampus several fold yet surprisingly, provoked repression of only a subset (∼10%) of potential NRSF target genes. Accordingly, the repressed gene-set was rescued when NRSF binding to chromatin was blocked. Unexpectedly, genes selectively repressed by NRSF had mid-range binding frequencies to the repressor, a property that rendered them sensitive to moderate fluctuations of NRSF levels. Genes selectively regulated by NRSF during epileptogenesis coded for ion channels, receptors, and other crucial contributors to neuronal function. Thus, dynamic, selective regulation of NRSF target genes may play a role in influencing neuronal properties in pathological and physiological contexts.

Pancreatic ductal adenocarcinoma (PDAC) remains one of the most lethal malignancies. It is typically detected at an advanced stage, at which the therapeutic options are very limited. One remarkable feature of PDAC that contributes to its resilience to treatment is the extreme stromal activation seen in these tumors. Often, the vast majority of tumor bulk consists of non-tumor cells that together provide a tumor-promoting environment. One of the signals that maintains and activates the stroma is the developmental protein Sonic Hedgehog (SHH). As the disease progresses, tumor cells produce increasing amounts of SHH, which activates the surrounding stroma to aid in tumor progression. To better understand this response and identify targets for inhibition, we aimed to elucidate the proteins that mediate the SHH-driven stromal response in PDAC. For this a novel mixed-species coculture model was set up in which the cancer cells are human, and the stroma is modeled by mouse fibroblasts. In conjunction with next-generation sequencing we were able to use the sequence difference between these species to genetically distinguish between the epithelial and stromal responses to SHH. The stromal SHH-dependent genes from this analysis were validated and their relevance for human disease was subsequently determined in two independent patient cohorts. In non-microdissected tissue from PDAC patients, in which a large amount of stroma is present, the targets were confirmed to associate with tumor stroma versus normal pancreatic tissue. Patient survival analysis and immunohistochemistry identified CDA, EDIL3, ITGB4, PLAUR and SPOCK1 as SHH-dependent stromal factors that are associated with poor prognosis in PDAC patients. Summarizing, the presented data provide insight into the role of the activated stroma in PDAC, and how SHH acts to mediate this response. In addition, the study has yielded several candidates that are interesting therapeutic targets for a disease for which treatment opt

Objectives: Primary osteoarthritis (OA) and Kashin–Beck disease (KBD) exhibit similar clinical manifestations and common articular cartilage lesions. Revealing the pathogenetic differences between OA and KBD is helpful for differential diagnosis and may provide new insights into the pathogenesis of OA and KBD. In this study, we compared the genome-wide gene ontology (GO) and pathway expression patterns of articular cartilage derived from both OA and KBD patients. Methods: Total RNA was isolated, amplified, labelled, and hybridized using Agilent whole genome microarray analysis. Gene set enrichment analysis (GSEA) was used to identify differentially expressed genes and pathways between OA and KBD. Nine differentially expressed GO categories and 85 differentially expressed pathways were identified by this study. Results: The reactive oxygen species (ROS)-related HOUSTIS_ROS pathway and the vascular endothelial growth factor (VEGF)-related ABE_VEGFA_TARGETS_2HR pathway were significantly up-regulated in OA compared to KBD. Higher expression levels of the collagen-related COLLAGEN GO, EXTRACELLULAR_MATRIX_PART GO, and nitric oxide (NO)-related BIOCARTA_NO1_PATHWAY pathways were detected in KBD than in OA. Conclusions: ROS-induced cartilage lesions seem to be more involved in the pathogenesis of OA whereas NO-mediated chondrocyte apoptosis contributes more to the development of KBD.

Pancreatic ductal adenocarcinoma (PDAC) remains one of the most lethal malignancies. It is typically detected at an advanced stage, at which the therapeutic options are very limited. One remarkable feature of PDAC that contributes to its resilience to treatment is the extreme stromal activation seen in these tumors. Often, the vast majority of tumor bulk consists of non-tumor cells that together provide a tumor-promoting environment. One of the signals that maintains and activates the stroma is the developmental protein Sonic Hedgehog (SHH). As the disease progresses, tumor cells produce increasing amounts of SHH, which activates the surrounding stroma to aid in tumor progression. To better understand this response and identify targets for inhibition, we aimed to elucidate the proteins that mediate the SHH-driven stromal response in PDAC. For this a novel mixed-species coculture model was set up in which the cancer cells are human, and the stroma is modeled by mouse fibroblasts. In conjunction with next-generation sequencing we were able to use the sequence difference between these species to genetically distinguish between the epithelial and stromal responses to SHH. The stromal SHH-dependent genes from this analysis were validated and their relevance for human disease was subsequently determined in two independent patient cohorts. In non-microdissected tissue from PDAC patients, in which a large amount of stroma is present, the targets were confirmed to associate with tumor stroma versus normal pancreatic tissue. Patient survival analysis and immunohistochemistry identified CDA, EDIL3, ITGB4, PLAUR and SPOCK1 as SHH-dependent stromal factors that are associated with poor prognosis in PDAC patients. Summarizing, the presented data provide insight into the role of the activated stroma in PDAC, and how SHH acts to mediate this response. In addition, the study has yielded several candidates that are interesting therapeutic targets for a disease for which treatment opt

Pancreatic ductal adenocarcinoma (PDAC) remains one of the most lethal malignancies. It is typically detected at an advanced stage, at which the therapeutic options are very limited. One remarkable feature of PDAC that contributes to its resilience to treatment is the extreme stromal activation seen in these tumors. Often, the vast majority of tumor bulk consists of non-tumor cells that together provide a tumor-promoting environment. One of the signals that maintains and activates the stroma is the developmental protein Sonic Hedgehog (SHH). As the disease progresses, tumor cells produce increasing amounts of SHH, which activates the surrounding stroma to aid in tumor progression. To better understand this response and identify targets for inhibition, we aimed to elucidate the proteins that mediate the SHH-driven stromal response in PDAC. For this a novel mixed-species coculture model was set up in which the cancer cells are human, and the stroma is modeled by mouse fibroblasts. In conjunction with next-generation sequencing we were able to use the sequence difference between these species to genetically distinguish between the epithelial and stromal responses to SHH. The stromal SHH-dependent genes from this analysis were validated and their relevance for human disease was subsequently determined in two independent patient cohorts. In non-microdissected tissue from PDAC patients, in which a large amount of stroma is present, the targets were confirmed to associate with tumor stroma versus normal pancreatic tissue. Patient survival analysis and immunohistochemistry identified CDA, EDIL3, ITGB4, PLAUR and SPOCK1 as SHH-dependent stromal factors that are associated with poor prognosis in PDAC patients. Summarizing, the presented data provide insight into the role of the activated stroma in PDAC, and how SHH acts to mediate this response. In addition, the study has yielded several candidates that are interesting therapeutic targets for a disease for which treatment opt

Background Self organizing maps (SOM) enable the straightforward portraying of high-dimensional data of large sample collections in terms of sample-specific images. The analysis of their texture provides so-called spot-clusters of co-expressed genes which require subsequent significance filtering and functional interpretation. We address feature selection in terms of the gene ranking problem and the interpretation of the obtained spot-related lists using concepts of molecular function. Results Different expression scores based either on simple fold change-measures or on regularized Student's t-statistics are applied to spot-related gene lists and compared with special emphasis on the error characteristics of microarray expression data. The spot-clusters are analyzed using different methods of gene set enrichment analysis with the focus on overexpression and/or overrepresentation of predefined sets of genes. Metagene-related overrepresentation of selected gene sets was mapped into the SOM images to assign gene function to different regions. Alternatively we estimated set-related overexpression profiles over all samples studied using a gene set enrichment score. It was also applied to the spot-clusters to generate lists of enriched gene sets. We used the tissue body index data set, a collection of expression data of human tissues as an illustrative example. We found that tissue related spots typically contain enriched populations of gene sets well corresponding to molecular processes in the respective tissues. In addition, we display special sets of housekeeping and of consistently weak and high expressed genes using SOM data filtering. Conclusions The presented methods allow the comprehensive downstream analysis of SOM-transformed expression data in terms of cluster-related gene lists and enriched gene sets for functional interpretation. SOM clustering implies the ability to define either new gene sets using selected SOM spots or to verify and/or to amend existing on