Comprehensive evaluation of coding region point mutations in microsatellite-unstable colorectal cancer.

Microsatellite instability (MSI) leads to accumulation of an excessive number of mutations in the genome, mostly small insertions and deletions. MSI colorectal cancers (CRCs), however, also contain more point mutations than microsatellite-stable (MSS) tumors, yet they have not been as comprehensively studied. To identify candidate driver genes affected by point mutations in MSI CRC, we ranked genes based on mutation significance while correcting for replication timing and gene expression utilizing an algorithm, MutSigCV Somatic point mutation data from the exome kit-targeted area from 24 exome-sequenced sporadic MSI CRCs and respective normals, and 12 whole-genome-sequenced sporadic MSI CRCs and respective normals were utilized. The top 73 genes were validated in 93 additional MSI CRCs. The MutSigCV ranking identified several well-established MSI CRC driver genes and provided additional evidence for previously proposed CRC candidate genes as well as shortlisted genes that have to our knowledge not been linked to CRC before. Two genes, SMARCB1 and STK38L, were also functionally scrutinized, providing evidence of a tumorigenic role, for SMARCB1 mutations in particular.

Glycomic Profiling Highlights Increased Fucosylation in Pseudomyxoma Peritonei.

Pseudomyxoma peritonei (PMP) is a subtype of mucinous adenocarcinoma that most often originates from the appendix, and grows in the peritoneal cavity filling it with mucinous ascites. KRAS and GNAS mutations are frequently found in PMP, but other common driver mutations are infrequent. As altered glycosylation can promote carcinogenesis, we compared N-linked glycan profiles of PMP tissues to those of normal appendix. Glycan profiles of eight normal appendix samples and eight low-grade and eight high-grade PMP specimens were analyzed by mass spectrometry. Our results show differences in glycan profiles between PMP and the controls, especially in those of neutral glycans, and the most prominent alteration was increased fucosylation. We further demonstrate up-regulated mRNA expression of four fucosylation-related enzymes, the core fucosylation performing fucosyltransferase 8 and three GDP-fucose biosynthetic enzymes in PMP tissues when compared with the controls. Up-regulated protein expression of the latter three enzymes was further observed in PMP cells by immunohistochemistry. We also demonstrate that restoration of fucosylation either by salvage pathway or by introduction of an expression of intact GDP-mannose 4,6-dehydratase enhance expression of MUC2, which is the predominant mucin molecule secreted by the PMP cells, in an intestinal-derived adenocarcinoma cell line with defective fucosylation because of deletion in the GDP-mannose 4,6-dehydratase gene. Thus, altered glycosylation especially in the form of fucosylation is linked to the characteristic mucin production of PMP. Glycomic data are available via ProteomeXchange with identifier PXD010086.

Upregulation of Early and Downregulation of Terminal Pathway Complement Genes in Subcutaneous Adipose Tissue and Adipocytes in Acquired Obesity.

Inflammation is an important mediator of obesity-related complications such as the metabolic syndrome but its causes and mechanisms are unknown. As the complement system is a key mediator of inflammation, we studied whether it is activated in acquired obesity in subcutaneous adipose tissue (AT) and isolated adipocytes. We used a special study design of genetically matched controls of lean and heavy groups, rare monozygotic twin pairs discordant for body mass index (BMI) [n = 26, within-pair difference (Δ) in body mass index, BMI >3 kg/m2] with as much as 18 kg mean Δweight. Additionally, 14 BMI-concordant (BMI <3 kg/m2) served as a reference group. The detailed measurements included body composition (DEXA), fat distribution (MRI), glucose, insulin, adipokines, C3a and SC5b-9 levels, and the expression of complement and insulin signaling pathway-related genes in AT and adipocytes. In both AT and isolated adipocytes, the classical and alternative pathway genes were upregulated, and the terminal pathway genes downregulated in the heavier co-twins of the BMI-discordant pairs. The upregulated genes included C1q, C1s, C2, ficolin-1, factor H, receptors for C3a and C5a (C5aR1), and the iC3b receptor (CR3). While the terminal pathway components C5 and C6 were downregulated, its inhibitor clusterin was upregulated. Complement gene upregulation in AT and adipocytes correlated positively with adiposity and hyperinsulinemia and negatively with the expression of insulin signaling-related genes. Plasma C3a, but not SC5b-9, levels were elevated in the heavier co-twins. There were no differences between the co-twins in BMI-concordant pairs. Obesity is associated with increased expression of the early, but not late, complement pathway components and of key receptors. The twins with acquired obesity have therefore an inflated inflammatory activity in the AT. The results suggest that complement is likely involved in orchestrating clearance of apoptotic debris and inflammation in the AT.

Multiple components of PKA and TGF-ß pathways are mutated in pseudomyxoma peritonei.

Pseudomyxoma peritonei (PMP) is a subtype of mucinous adenocarcinoma mainly restricted to the peritoneal cavity and most commonly originating from the appendix. The genetic background of PMP is poorly understood and no targeted treatments are currently available for this fatal disease. While RAS signaling pathway is affected in most if not all PMP cases and over half of them also have a mutation in the GNAS gene, other genetic alterations and affected pathways are, to a large degree, poorly known. In this study, we sequenced whole coding genome of nine PMP tumors and paired normal tissues in order to identify additional, commonly mutated genes and signaling pathways affected in PMP. These exome sequencing results were validated with an ultra-deep amplicon sequencing method, leading to 14 validated variants. The validated results contain seven genes that contribute to the protein kinase A (PKA) pathway. PKA pathway, which also contains GNAS, is a major player of overproduction of mucin, which is the characteristic feature of PMP. In addition to PKA pathway, we identified mutations in six genes that belong to the transforming growth factor beta (TGF-ß) pathway, which is a key regulator of cell proliferation. Since either GNAS mutation or an alternative mutation in the PKA pathway was identified in 8/9 patients, inhibition of the PKA pathway might reduce mucin production in most of the PMP patients and potentially suppress disease progression.

Complete androgen insensitivity syndrome caused by a deep intronic pseudoexon-activating mutation in the androgen receptor gene.

Mutations in the X-linked androgen receptor (AR) gene underlie complete androgen insensitivity syndrome (CAIS), the most common cause of 46,XY sex reversal. Molecular genetic diagnosis of CAIS, however, remains uncertain in patients who show normal coding region of AR. Here, we describe a novel mechanism of AR disruption leading to CAIS in two 46,XY sisters. We analyzed whole-genome sequencing data of the patients for pathogenic variants outside the AR coding region. Patient fibroblasts from the genital area were used for AR cDNA analysis and protein quantification. Analysis of the cDNA revealed aberrant splicing of the mRNA caused by a deep intronic mutation (c.2450-118A>G) in the intron 6 of AR. The mutation creates a de novo 5' splice site and a putative exonic splicing enhancer motif, which leads to the preferential formation of two aberrantly spliced mRNAs (predicted to include a premature stop codon). Patient fibroblasts contained no detectable AR protein. Our results show that patients with CAIS and normal AR coding region need to be examined for deep intronic mutations that can lead to pseudoexon activation.

Genomic profile of pseudomyxoma peritonei analyzed using next-generation sequencing and immunohistochemistry.

Pseudomyxoma peritonei (PMP) is a relatively rare clinical syndrome characterized by neoplastic epithelial cells growing in the peritoneal cavity and secreting mucinous ascites. Our aim was to explore the molecular events behind this fatal but under-investigated disease. We extracted DNA from 19 appendix-derived PMP tumors and nine corresponding normal tissues, and analyzed the mutational hotspot areas of 48 cancer-related genes by amplicon-based next-generation sequencing (NGS). Further, we analyzed the protein expression of V600E mutated BRAF, MLH1, MSH2, MSH6 and p53 from a larger set of PMP tumors (n = 74) using immunohistochemistry. With NGS, we detected activating somatic KRAS mutations in all of the tumors studied. GNAS was mutated in 63% of the tumors with no marked difference between low-grade and high-grade tumors. Only one (5.3%) tumor showed oncogenic PIK3CA mutation, one showed oncogenic AKT1 mutation, three (15.8%) showed SMAD4 mutations and none showed an APC mutation. P53 protein was aberrantly expressed in higher proportion of high-grade tumors as compared with low-grade ones (31.3 vs. 7.1%, respectively; p = 0.012) and aberrant expression was an independent factor for reduced overall survival (p = 0.002). BRAF V600E mutation was only found in one (1.4%) high-grade tumor by immunohistochemistry (n = 74). All the studied tumors expressed mismatch repair proteins MLH1, MSH2 and MSH6. Our results indicate that KRAS mutations are evident in all and GNAS mutations in most of the PMPs, but BRAF V600E, PIK3CA and APC mutations are rare. Aberrantly expressed p53 is associated with high-grade histology and reduced survival.

Prostatic acid phosphatase is the main acid phosphatase with 5'-ectonucleotidase activity in the male mouse saliva and regulates salivation.

We have previously shown that in addition to the well-known secreted isoform of prostatic acid phosphatase (sPAP), a transmembrane isoform exists (TMPAP) that interacts with snapin (a SNARE-associated protein) and regulates the endo-/exocytic pathways. We have also shown that PAP has 5'-ectonucleotidase and thiamine monophosphatase activity and elicits antinociceptive effects in mouse models of chronic inflammatory and neuropathic pain. Therefore, to determine the physiological role of PAP in a typical exocrine organ, we studied the submandibular salivary gland (SMG) of PAP(-/-) and wild-type C57BL/6J mice by microarray analyses, microRNA sequencing, activity tests, immunohistochemistry, and biochemical and physiological analyses of saliva. We show that PAP is the main acid phosphatase in the wild-type male mouse saliva, accounting for 50% of the total acid phosphatase activity, and that it is expressed only in the granular convoluted tubules of the SMGs, where it is the only 5'-ectonucleotidase. The lack of PAP in male PAP(-/-) mice was associated with a significant increase in the salivation volume under secretagogue stimulation, overexpression of genes related to cell proliferation (Mki67, Aurkb, Birc5) and immune response (Irf7, Cxcl9, Ccl3, Fpr2), and upregulation of miR-146a in SMGs. An increased and sustained acinar cell proliferation was detected without signs of glandular hyperplasia. Our results indicate that in PAP(-/-) mice, SMG homeostasis is maintained by an innate immune response. Additionally, we suggest that in male mice, PAP via its 5'-ectonucleotidase activity and production of adenosine can elicit analgesic effects when animals lick their wounds.

Analysis of BMP4 and BMP7 signaling in breast cancer cells unveils time-dependent transcription patterns and highlights a common synexpression group of genes.

BACKGROUND: Bone morphogenetic proteins (BMPs) are members of the TGF-beta superfamily of growth factors. They are known for their roles in regulation of osteogenesis and developmental processes and, in recent years, evidence has accumulated of their crucial functions in tumor biology. BMP4 and BMP7, in particular, have been implicated in breast cancer. However, little is known about BMP target genes in the context of tumor. We explored the effects of BMP4 and BMP7 treatment on global gene transcription in seven breast cancer cell lines during a 6-point time series, using a whole-genome oligo microarray. Data analysis included hierarchical clustering of differentially expressed genes, gene ontology enrichment analyses and model based clustering of temporal data. RESULTS: Both ligands had a strong effect on gene expression, although the response to BMP4 treatment was more pronounced. The cellular functions most strongly affected by BMP signaling were regulation of transcription and development. The observed transcriptional response, as well as its functional outcome, followed a temporal sequence, with regulation of gene expression and signal transduction leading to changes in metabolism and cell proliferation. Hierarchical clustering revealed distinct differences in the response of individual cell lines to BMPs, but also highlighted a synexpression group of genes for both ligands. Interestingly, the majority of the genes within these synexpression groups were shared by the two ligands, probably representing the core molecular responses common to BMP4 and BMP7 signaling pathways. CONCLUSIONS: All in all, we show that BMP signaling has a remarkable effect on gene transcription in breast cancer cells and that the functions affected follow a logical temporal pattern. Our results also uncover components of the common cellular transcriptional response to BMP4 and BMP7. Most importantly, this study provides a list of potential novel BMP target genes relevant in breast cancer.

The phytoestrogen genistein is a tissue-specific androgen receptor modulator.

To enable studies of androgen signaling in different tissues in vivo, we generated an androgen receptor (AR) reporter mouse line by inserting a luciferase gene construct into the murine genome. The construct is driven by four copies of androgen-responsive elements from the mouse sex-limited protein gene (slp-HRE2) and a minimal thymidine kinase promoter. Luciferase activity was readily measurable in a number of murine tissues, including prostate, lung, testis, brain, and skeletal muscle, and testosterone administration elicited a significant increase in reporter gene activity in these tissues. Consumption of isoflavonoid genistein is linked to reduced risk of prostate cancer, but direct effects of genistein on the AR pathway are not well understood. To examine androgen-modulating activity of genistein in vivo, male mice received daily doses of genistein (10 mg/kg) for 5 d. In intact males, genistein was antiandrogenic in testis, prostate, and brain, and it attenuated reporter gene activity by 50-80%. In castrated males, genistein exhibited significant androgen agonistic activity in prostate and brain by increasing reporter gene activity over 2-fold in both tissues. No antiandrogenic action was seen in lung or skeletal muscle of intact males. Gene expression profiling of the murine prostate under the same experimental conditions revealed that genistein modulates androgen-dependent transcription program in prostate in a fashion similar to that observed in reporter mice by luciferase expression. In conclusion, genistein is a partial androgen agonist/antagonist in some but not in all mouse tissues and should be considered as a tissue-specific AR modulator.

Large-scale data integration framework provides a comprehensive view on glioblastoma multiforme.

BACKGROUND: Coordinated efforts to collect large-scale data sets provide a basis for systems level understanding of complex diseases. In order to translate these fragmented and heterogeneous data sets into knowledge and medical benefits, advanced computational methods for data analysis, integration and visualization are needed. METHODS: We introduce a novel data integration framework, Anduril, for translating fragmented large-scale data into testable predictions. The Anduril framework allows rapid integration of heterogeneous data with state-of-the-art computational methods and existing knowledge in bio-databases. Anduril automatically generates thorough summary reports and a website that shows the most relevant features of each gene at a glance, allows sorting of data based on different parameters, and provides direct links to more detailed data on genes, transcripts or genomic regions. Anduril is open-source; all methods and documentation are freely available. RESULTS: We have integrated multidimensional molecular and clinical data from 338 subjects having glioblastoma multiforme, one of the deadliest and most poorly understood cancers, using Anduril. The central objective of our approach is to identify genetic loci and genes that have significant survival effect. Our results suggest several novel genetic alterations linked to glioblastoma multiforme progression and, more specifically, reveal Moesin as a novel glioblastoma multiforme-associated gene that has a strong survival effect and whose depletion in vitro significantly inhibited cell proliferation. All analysis results are available as a comprehensive website. CONCLUSIONS: Our results demonstrate that integrated analysis and visualization of multidimensional and heterogeneous data by Anduril enables drawing conclusions on functional consequences of large-scale molecular data. Many of the identified genetic loci and genes having significant survival effect have not been reported earlier in the context of glioblastoma multiforme. Thus, in addition to generally applicable novel methodology, our results provide several glioblastoma multiforme candidate genes for further studies.Anduril is available at http://csbi.ltdk.helsinki.fi/anduril/The glioblastoma multiforme analysis results are available at http://csbi.ltdk.helsinki.fi/anduril/tcga-gbm/

Gene set enrichment analysis reveals several globally affected pathways due to SKI-1/S1P inhibition in HepG2 cells.

Sterol regulatory element-binding proteins (SREBPs) are transcription factors governing transcription of genes related to cholesterol and fatty acid metabolism. To become active, SREBPs must undergo a proteolytic cleavage to allow an active NH(2)-terminal segment to translocate into the nucleus. SKI-1/S1P is the first protease in the proteolytic activation cascade of SREBPs. SREBP inhibition may be useful, for example, in the treatment of liver steatosis caused by homocysteine-induced lipid synthesis. Accordingly, we overexpressed inhibitory prodomains (proSKI) of SKI-1/S1P in HepG2 cells to block SREBP activation to evaluate the potential of SKI-1/S1P in controlling cellular cholesterol synthesis. SKI-1/S1P inhibition resulted in reduced cholesterol synthesis and mRNA levels of the rate-limiting enzymes, HMG-CoA reductase and squalene epoxidase, in the cholesterol synthetic pathway. The inhibitory effect was maintained in the presence of homocysteine-induced endoplasmic reticulum stress. A gene set enrichment analysis was performed to elucidate other metabolic effects caused by SKI-1/S1P inhibition. SKI-1/S1P inhibition was observed to affect a number of other metabolic pathways, including glycolysis and citric acid cycle. These results demonstrate that inhibition of SREBPs decreases cholesterol synthesis in HepG2 cells both in the absence and presence of homocysteine. SKI-1/S1P inhibition may cause widespread changes in other key metabolic pathways.

Comparative genomics of Bordetella pertussis reveals progressive gene loss in Finnish strains.

BACKGROUND: Bordetella pertussis is a gram-negative bacterium that infects the human respiratory tract and causes pertussis or whooping cough. The disease has resurged in many countries including Finland where the whole-cell pertussis vaccine has been used for more than 50 years. Antigenic divergence has been observed between vaccine strains and clinical isolates in Finland. To better understand genome evolution in B. pertussis circulating in the immunized population, we developed an oligonucleotide-based microarray for comparative genomic analysis of Finnish strains isolated during the period of 50 years. METHODOLOGY/PRINCIPAL FINDINGS: The microarray consisted of 3,582 oligonucleotides (70-mer) and covered 94% of 3,816 ORFs of Tohama I, the strain of which the genome has been sequenced. Twenty isolates from 1953 to 2004 were studied together with two Finnish vaccine strains and two international reference strains. The isolates were selected according to their characteristics, e.g. the year and place of isolation and pulsed-field gel electrophoresis profiles. Genomic DNA of the tested strains, along with reference DNA of Tohama I strain, was labelled and hybridized. The absence of genes as established with microarrays, was confirmed by PCR. Compared with the Tohama I strain, Finnish isolates lost 7 (8.6 kb) to 49 (55.3 kb) genes, clustered in one to four distinct loci. The number of lost genes increased with time, and one third of lost genes had functions related to inorganic ion transport and metabolism, or energy production and conversion. All four loci of lost genes were flanked by the insertion sequence element IS481. CONCLUSION/SIGNIFICANCE: Our results showed that the progressive gene loss occurred in Finnish B. pertussis strains isolated during a period of 50 years and confirmed that B. pertussis is dynamic and is continuously evolving, suggesting that the bacterium may use gene loss as one strategy to adapt to highly immunized populations.

A systems biology strategy reveals biological pathways and plasma biomarker candidates for potentially toxic statin-induced changes in muscle.

BACKGROUND: Aggressive lipid lowering with high doses of statins increases the risk of statin-induced myopathy. However, the cellular mechanisms leading to muscle damage are not known and sensitive biomarkers are needed to identify patients at risk of developing statin-induced serious side effects. METHODOLOGY: We performed bioinformatics analysis of whole genome expression profiling of muscle specimens and UPLC/MS based lipidomics analyses of plasma samples obtained in an earlier randomized trial from patients either on high dose simvastatin (80 mg), atorvastatin (40 mg), or placebo. PRINCIPAL FINDINGS: High dose simvastatin treatment resulted in 111 differentially expressed genes (1.5-fold change and p-value<0.05), while expression of only one and five genes was altered in the placebo and atorvastatin groups, respectively. The Gene Set Enrichment Analysis identified several affected pathways (23 gene lists with False Discovery Rate q-value<0.1) in muscle following high dose simvastatin, including eicosanoid synthesis and Phospholipase C pathways. Using lipidomic analysis we identified previously uncharacterized drug-specific changes in the plasma lipid profile despite similar statin-induced changes in plasma LDL-cholesterol. We also found that the plasma lipidomic changes following simvastatin treatment correlate with the muscle expression of the arachidonate 5-lipoxygenase-activating protein. CONCLUSIONS: High dose simvastatin affects multiple metabolic and signaling pathways in skeletal muscle, including the pro-inflammatory pathways. Thus, our results demonstrate that clinically used high statin dosages may lead to unexpected metabolic effects in non-hepatic tissues. The lipidomic profiles may serve as highly sensitive biomarkers of statin-induced metabolic alterations in muscle and may thus allow us to identify patients who should be treated with a lower dose to prevent a possible toxicity.