Breast cancer is the leading cause of cancer related death among women. Breast cancers are generally diagnosed and treated based on clinical and histopathological features, along with subtype ...classification determined by the Prosigna Breast Cancer Prognostic Gene Signature Assay (also known as PAM50). Currently the copy number alteration (CNA) landscape of the tumour is not considered. We set out to examine the role of genomic instability (GI) in breast cancer survival since CNAs reflect GI and correlate with survival in other cancers. We focused on the 70% of breast cancers classified as luminal and carried out a comprehensive survival and association analysis using Molecular Taxonomy of Breast Cancer International Consortium (METABRIC) data to determine whether CNA Score Quartiles derived from absolute CNA counts are associated with survival. Analysis revealed that patients diagnosed with luminal A breast cancer have a CNA landscape associated with disease specific survival, suggesting that CNA Score can provide a statistically robust prognostic factor. Furthermore, stratification of patients into subtypes based on gene expression has shown that luminal A and B cases overlap, and it is in this region we largely observe luminal A cases with reduced survival outlook. Therefore, luminal A breast cancer patients with quantitatively elevated CNA counts may benefit from more aggressive therapy. This demonstrates how individual genomic landscapes can facilitate personalisation of therapeutic interventions to optimise survival outcomes.
WT1, a critical regulator of kidney development, is a tumor suppressor for nephroblastoma but in some contexts functions as an oncogene. A limited number of direct transcriptional targets of WT1 have ...been identified to explain its complex roles in tumorigenesis and organogenesis. In this study we performed genome-wide screening for direct WT1 targets, using a combination of ChIP-ChIP and expression arrays. Promoter regions bound by WT1 were highly G-rich and resembled the sites for a number of other widely expressed transcription factors such as SP1, MAZ, and ZNF219. Genes directly regulated by WT1 were implicated in MAPK signaling, axon guidance, and Wnt pathways. Among directly bound and regulated genes by WT1, nine were identified in the Wnt signaling pathway, suggesting that WT1 modulates a subset of Wnt components and responsive genes by direct binding. To prove the biological importance of the interplay between WT1 and Wnt signaling, we showed that WT1 blocked the ability of Wnt8 to induce a secondary body axis during Xenopus embryonic development. WT1 inhibited TCF-mediated transcription activated by Wnt ligand, wild type and mutant, stabilized β-catenin by preventing TCF4 loading onto a promoter. This was neither due to direct binding of WT1 to the TCF binding site nor to interaction between WT1 and TCF4, but by competition of WT1 and TCF4 for CBP. WT1 interference with Wnt signaling represents an important mode of its action relevant to the suppression of tumor growth and guidance of development.
DNA methylation is an epigenetic mark that can stably repress gene expression. Because of its biological and clinical significance, several methods have been developed to compare genome-wide patterns ...of methylation between groups of samples. The application of gene set analysis to identify relevant groups of genes that are enriched for differentially methylated genes is often a major component of the analysis of these data. This can be used, for example, to identify processes or pathways that are perturbed in disease development. We show that gene-set analysis, as it is typically applied to genome-wide methylation assays, is severely biased as a result of differences in the numbers of CpG sites associated with different classes of genes and gene promoters.
We demonstrate this bias using published data from a study of differential CpG island methylation in lung cancer and a dataset we generated to study methylation changes in patients with long-standing ulcerative colitis. We show that several of the gene sets that seem enriched would also be identified with randomized data. We suggest two existing approaches that can be adapted to correct the bias. Accounting for the bias in the lung cancer and ulcerative colitis datasets provides novel biological insights into the role of methylation in cancer development and chronic inflammation, respectively. Our results have significant implications for many previous genome-wide methylation studies that have drawn conclusions on the basis of such strongly biased analysis.
cathal.seoighe@nuigalway.ie
Supplementary data are available at Bioinformatics online.
Abstract
We simulate possible stellar coronal mass ejection (CME) scenarios over the magnetic cycle of
ϵ
Eridani (18 Eridani; HD 22049). We use three separate epochs from 2008, 2011, and 2013, and ...estimate the radio emission frequencies associated with these events. These stellar eruptions have proven to be elusive, although a promising approach to detect and characterize these phenomena are low-frequency radio observations of potential type II bursts as CME-induced shocks propagate through the stellar corona. Stellar type II radio bursts are expected to emit below 450 MHz, similarly to their solar counterparts. We show that the length of time these events remain above the ionospheric cutoff is not necessarily dependent on the stellar magnetic cycle, but more on the eruption location relative to the stellar magnetic field. We find that these type II bursts would remain within the frequency range of LOFAR for a maximum of 20–30 minutes post-eruption for the polar CMEs (50 minutes for second harmonics). We find evidence of slower equatorial CMEs, which result in slightly longer observable windows for the 2008 and 2013 simulations. Stellar magnetic geometry and strength have a significant effect on the detectability of these events. We place the CMEs in the context of the stellar mass-loss rate (27–48× solar mass-loss rate), showing that they can amount to 3%–50% of the stellar wind mass-loss rate for
ϵ
Eridani. Continuous monitoring of likely stellar CME candidates with low-frequency radio telescopes will be required to detect these transient events.
Familial binding profiles (FBPs) represent the average binding specificity for a group of structurally related DNA-binding proteins. The construction of such profiles allows the classification of ...novel motifs based on similarity to known families, can help to reduce redundancy in motif databases and de novo prediction algorithms, and can provide valuable insights into the evolution of binding sites. Many current approaches to automated motif clustering rely on progressive tree-based techniques, and can suffer from so-called frozen sub-alignments, where motifs which are clustered early on in the process remain 'locked' in place despite the potential for better placement at a later stage. In order to avoid this scenario, we have developed a genetic-k-medoids approach which allows motifs to move freely between clusters at any point in the clustering process.
We demonstrate the performance of our algorithm, GMACS, on multiple benchmark motif datasets, comparing results obtained with current leading approaches. The first dataset includes 355 position weight matrices from the TRANSFAC database and indicates that the k-mer frequency vector approach used in GMACS outperforms other motif comparison techniques. We then cluster a set of 79 motifs from the JASPAR database previously used in several motif clustering studies and demonstrate that GMACS can produce a higher number of structurally homogeneous clusters than other methods without the need for a large number of singletons. Finally, we show the robustness of our algorithm to noise on multiple synthetic datasets consisting of known motifs convolved with varying degrees of noise.
Our proposed algorithm is generally applicable to any DNA or protein motifs, can produce highly stable and biologically meaningful clusters, and, by avoiding the problem of frozen sub-alignments, can provide improved results when compared with existing techniques on benchmark datasets.
Assessing whole-body radiation injury and absorbed dose is essential for remediation efforts following accidental or deliberate exposure in medical, industrial, military, or terrorist incidents. We ...hypothesize that variations in specific metabolite concentrations extracted from blood plasma would correlate with whole-body radiation injury and dose.
Groups of C57BL/6 mice (n=12 per group) were exposed to 0, 2, 4, 8, and 10.4 Gy of whole-body gamma radiation. At 24 hours after treatment, all animals were euthanized, and both plasma and liver biopsy samples were obtained, the latter being used to identify a distinct hepatic radiation injury response within plasma. A semiquantitative, untargeted metabolite/lipid profile was developed using gas chromatography-mass spectrometry and liquid chromatography-tandem mass spectrometry, which identified 354 biochemical compounds. A second set of C57BL/6 mice (n=6 per group) were used to assess a subset of identified plasma markers beyond 24 hours.
We identified a cohort of 37 biochemical compounds in plasma that yielded the optimal separation of the irradiated sample groups, with the most correlated metabolites associated with pyrimidine (positively correlated) and tryptophan (negatively correlated) metabolism. The latter were predominantly associated with indole compounds, and there was evidence that these were also correlated between liver and plasma. No evidence of saturation as a function of dose was observed, as has been noted for studies involving metabolite analysis of urine.
Plasma profiling of specific metabolites related to pyrimidine and tryptophan pathways can be used to differentiate whole-body radiation injury and dose response. As the tryptophan-associated indole compounds have their origin in the intestinal microbiome and subsequently the liver, these metabolites particularly represent an attractive marker for radiation injury within blood plasma.
Unrestrained receptor tyrosine kinase (RTK) signaling and epigenetic deregulation are root causes of tumorigenesis. We establish linkage between these processes by demonstrating that aberrant RTK ...signaling unleashed by oncogenic HRasG12V or loss of negative feedback through Sprouty gene deletion remodels histone modifications associated with active typical and super-enhancers. However, although both lesions disrupt the Ras-Erk axis, the expression programs, enhancer signatures, and transcription factor networks modulated upon HRasG12V transformation or Sprouty deletion are largely distinct. Oncogenic HRasG12V elevates histone 3 lysine 27 acetylation (H3K27ac) levels at enhancers near the transcription factor Gata4 and the kinase Prkcb, as well as their expression levels. We show that Gata4 is necessary for the aberrant gene expression and H3K27ac marking at enhancers, and Prkcb is required for the oncogenic effects of HRasG12V-driven cells. Taken together, our findings demonstrate that dynamic reprogramming of the cellular enhancer landscape is a major effect of oncogenic RTK signaling.
Display omitted
•Unrestrained Ras-Erk signaling deregulates histone marking at active enhancers•Sprouty loss and oncogenic HRasG12V alter distinct enhancers and target genes•Gata4 is necessary for HRasG12V-driven enhancer marking and target gene expression•Prkcb and BET bromodomains are necessary for the oncogenic effects of HRasG12V
Aberrant receptor tyrosine kinase signaling mediated by oncogenic Ras or loss of Sprouty promotes tumorigenesis. Nabet et al. find that unrestrained receptor tyrosine signaling driven by these lesions alters distinct super-enhancers, transcription factors, and target genes. Gata4 and Prkcb are identified as mediators of the oncogenic program upon Ras transformation.
•RNA-seq data taken from patient tumor and normal samples are used to derive GARD (Genome Adjusted Radiation Dose) radiosensitivity indices for a luminal B breast cancer patient (n = 1), and a cohort ...of prostate cancer patients (n = 14)•Clear differences in radiosensitivity were evident between tumor and normal tissues for the luminal B breast cancer patient•A treatment of 50 Gy/25fx would yield a GARD value passing the threshold for an optimal therapeutic index for this breast cancer patient•The majority of prostate cancer patients show tumour radiosensitivity differences and so potential suitability for hypofractionation•A treatment of 72 Gy/36fx results in variable GARD values for the prostate cancer patients, highlighting the benefit of higher optimal dosages in some cases
The use of a 10 gene transcriptional signature as part of the GARD model has been shown to be predictive of radiotherapy benefit for a range of cancers, with the potential to determine an optimal overall dose per patient. We used publicly available RNA-seq transcriptomics data from a luminal B breast cancer patient and from 14 prostate cancer patients to explore the radiosensitivity indices (RSI) and so GARD estimates of both tumour and proximal normal biopsies from each individual. Clear differences of clinical relevance in derived radiobiological properties between tumour and proximal normal tissues were evident for the breast cancer patient, whilst such differences across the prostate cancer cohort were more equivocal. Using the prostate cancer cohort’s median tumour predicted GARD value as a threshold for high therapeutic effect for radiotherapy, we found evidence that a higher overall prescribed dose than the widely used 72 Gy/36fx could benefit half of these patients. This exploratory study demonstrates the potential combining the GARD model with sequencing based transcriptomics could have in informing personalised radiotherapeutic practise for both breast and prostate cancer patients.
The diagnostic criteria for antibody-mediated rejection (AMR) are continuously evolving. Here we investigated the clinical and molecular significance of different Banff microvascular inflammation ...(MVI) scores in transplant kidney biopsies. A total of 356 patients with clinically indicated kidney transplant biopsies were classified into three groups based on MVI scores of 0, 1, 2, or more for Groups 1–3, respectively. Gene expression profiles were assessed using arrays on a representative subset of 93 patients. The incidence of donor-specific anti-HLA antibodies was increased from 25% in Group 1 to 36% in Group 2 and to 54% in Group 3. Acute and chronic AMR were significantly more frequent in Group 3 (15% and 35%) compared with the Group 2 (3% and 15%) and Group 1 (0% and 5%), respectively. Gene expression profiles showed increased interferon-γ and rejection-induced, cytotoxic and regulatory T-cell, natural killer cell–associated and donor-specific antibody (DSA)–selective transcripts in Group 3 compared with Groups 1 and 2. There was no significant difference in gene expression profiles between the Groups 1 and 2. Increased intragraft expression of DSA–selective transcripts was found in the biopsies of C4d− Group 3 patients. Thus, an MVI score of 2 or more was significantly associated with a histological diagnosis of acute and chronic antibody-mediated rejection. Hence, increased intragraft DSA–selective gene transcripts may be used as molecular markers for AMR, especially in C4d− biopsies.
RNA:DNA hybrids represent a non-canonical nucleic acid structure that has been associated with a range of human diseases and potential transcriptional regulatory functions. Mapping of RNA:DNA hybrids ...in human cells reveals them to have a number of characteristics that give insights into their functions.
We find RNA:DNA hybrids to occupy millions of base pairs in the human genome. A directional sequencing approach shows the RNA component of the RNA:DNA hybrid to be purine-rich, indicating a thermodynamic contribution to their in vivo stability. The RNA:DNA hybrids are enriched at loci with decreased DNA methylation and increased DNase hypersensitivity, and within larger domains with characteristics of heterochromatin formation, indicating potential transcriptional regulatory properties. Mass spectrometry studies of chromatin at RNA:DNA hybrids shows the presence of the ILF2 and ILF3 transcription factors, supporting a model of certain transcription factors binding preferentially to the RNA:DNA conformation.
Overall, there is little to indicate a dependence for RNA:DNA hybrids forming co-transcriptionally, with results from the ribosomal DNA repeat unit instead supporting the intriguing model of RNA generating these structures in trans. The results of the study indicate heterogeneous functions of these genomic elements and new insights into their formation and stability in vivo.