Computational Models for Clinical Applications in Personalized Medicine-Guidelines and Recommendations for Data Integration and Model Validation.

The future development of personalized medicine depends on a vast exchange of data from different sources, as well as harmonized integrative analysis of large-scale clinical health and sample data. Computational-modelling approaches play a key role in the analysis of the underlying molecular processes and pathways that characterize human biology, but they also lead to a more profound understanding of the mechanisms and factors that drive diseases; hence, they allow personalized treatment strategies that are guided by central clinical questions. However, despite the growing popularity of computational-modelling approaches in different stakeholder communities, there are still many hurdles to overcome for their clinical routine implementation in the future. Especially the integration of heterogeneous data from multiple sources and types are challenging tasks that require clear guidelines that also have to comply with high ethical and legal standards. Here, we discuss the most relevant computational models for personalized medicine in detail that can be considered as best-practice guidelines for application in clinical care. We define specific challenges and provide applicable guidelines and recommendations for study design, data acquisition, and operation as well as for model validation and clinical translation and other research areas.

The citrullinated/native index of autoantibodies against hnRNP-DL predicts an individual "window of treatment success" in RA patients.

BACKGROUND: There is a need for biomarker to identify patients "at risk" for rheumatoid arthritis (risk-RA) and to better predict the therapeutic response and in this study we tested the hypothesis that novel native and citrullinated heterogeneous nuclear ribonucleoprotein (hnRNP)-DL autoantibodies could be possible biomarkers. METHODS: Using protein macroarray and ELISA, epitope recognition against hnRNP-DL was analysed in sera from different developed RA disease and diagnosed SLE patients. Toll-like receptor (TLR) 7/9 and myeloid differentiation primary response gene 88 (MyD88)-dependency were studied in sera from murine disease models. HnRNP-DL expression in cultivated cells and synovial tissue was analysed by indirect immunofluorescence, immunoblot and immunohistochemistry. RESULTS: HnRNP-DL was highly expressed in stress granules, citrullinated in the rheumatoid joint and targeted by autoantibodies either as native or citrullinated proteins in patient subsets with different developed RA disease. Structural citrullination dependent epitopes (SCEs) of hnRNP-DL were detected in 58% of the SLE patients although 98% of these sera were α-CCP-2-negative. To obtain a specific citrullinated signal value, we subtracted the native antibody value from the citrullinated signal. The citrullinated/native index of autoantibodies against hnRNP-DL (CNDL-Index) was identified as a new value for an "individual window of treatment success" in early RA and for the detection of RF IgM/α-CCP-2 seronegative RA patients (24-46%). Negative CNDL-index was found in SLE patients, risk-RA and early RA cohorts such as EIRA where the majority of these patients are DAS28-responders to methotrexate (MTX) treatment (87%). High positive CNDL-values were associated with more severe RA, shared epitope and parenchymal changes in the lung. Specifically, native α-hnRNP-DL is TLR7/9-dependent, associated with pain and ROC analysis revealed an association to initial MTX or etanercept treatment response, especially in seronegative RA patients. CONCLUSION: CNDL-index defines people at risk to develop RA and the "window of treatment success" thereby closing the sensitivity gap in RA.

Yeast XRS2 and human NBN gene: Experimental evidence for homology using codon optimized cDNA.

The genes, XRS2 in Saccharomyces cerevisiae and NBN in mammals, have little sequence identity at the amino acid level. Nevertheless, they are both found together with MRE11 and RAD50 in a highly conserved protein complex which functions in the repair of DNA double-strand breaks. Here, we have examined the evolutionary and functional relationship of these two genes by cross-complementation experiments. These experiments necessitated sequence correction for specific codon usage before they could be successfully conducted. We present evidence that despite extreme sequence divergence nibrin can, at least partially, replace Xrs2 in the cellular DNA damage response, and Xrs2 is able to promote nuclear localization of MRE11 in NBS cells. We discuss that the extreme sequence divergence reflects a unique adaptive pressure during evolution related to the specific eukaryotic role for both Xrs2 and nibrin in the subcellular localisation of the DNA repair complex. This, we suggest, is of particular relevance when cells are infected by viruses. The conflict hypothesis of co-evolution of DNA repair genes and DNA viruses may thus explain the very low sequence identity of these two homologous genes.

A Clb/Cdk1-mediated regulation of Fkh2 synchronizes CLB expression in the budding yeast cell cycle.

Precise timing of cell division is achieved by coupling waves of cyclin-dependent kinase (Cdk) activity with a transcriptional oscillator throughout cell cycle progression. Although details of transcription of cyclin genes are known, it is unclear which is the transcriptional cascade that modulates their expression in a timely fashion. Here, we demonstrate that a Clb/Cdk1-mediated regulation of the Fkh2 transcription factor synchronizes the temporal mitotic CLB expression in budding yeast. A simplified kinetic model of the cyclin/Cdk network predicts a linear cascade where a Clb/Cdk1-mediated regulation of an activator molecule drives CLB3 and CLB2 expression. Experimental validation highlights Fkh2 as modulator of CLB3 transcript levels, besides its role in regulating CLB2 expression. A Boolean model based on the minimal number of interactions needed to capture the information flow of the Clb/Cdk1 network supports the role of an activator molecule in the sequential activation, and oscillatory behavior, of mitotic Clb cyclins. This work illustrates how transcription and phosphorylation networks can be coupled by a Clb/Cdk1-mediated regulation that synchronizes them.

The bromodomain protein BRD4 regulates splicing during heat shock.

The cellular response to heat stress is an ancient and evolutionarily highly conserved defence mechanism characterised by the transcriptional up-regulation of cyto-protective genes and a partial inhibition of splicing. These features closely resemble the proteotoxic stress response during tumor development. The bromodomain protein BRD4 has been identified as an integral member of the oxidative stress as well as of the inflammatory response, mainly due to its role in the transcriptional regulation process. In addition, there are also several lines of evidence implicating BRD4 in the splicing process. Using RNA-sequencing we found a significant increase in splicing inhibition, in particular intron retentions (IR), following heat treatment in BRD4-depleted cells. This leads to a decrease of mRNA abundancy of the affected transcripts, most likely due to premature termination codons. Subsequent experiments revealed that BRD4 interacts with the heat shock factor 1 (HSF1) such that under heat stress BRD4 is recruited to nuclear stress bodies and non-coding SatIII RNA transcripts are up-regulated. These findings implicate BRD4 as an important regulator of splicing during heat stress. Our data which links BRD4 to the stress induced splicing process may provide novel mechanisms of BRD4 inhibitors in regard to anti-cancer therapies.

Quantitative Global Proteomics of Yeast PBP1 Deletion Mutants and Their Stress Responses Identifies Glucose Metabolism, Mitochondrial, and Stress Granule Changes.

The yeast protein PBP1 is implicated in very diverse pathways. Intriguingly, its deletion mitigates the toxicity of human neurodegeneration factors. Here, we performed label-free quantitative global proteomics to identify crucial downstream factors, either without stress or under cell stress conditions (heat and NaN3). Compared to the wildtype BY4741 strain, PBP1 deletion always triggered downregulation of the key bioenergetics enzyme KGD2 and the prion protein RNQ1 as well as upregulation of the leucine biosynthesis enzyme LEU1. Without stress, enrichment of stress response factors was consistently detected for both deletion mutants; upon stress, these factors were more pronounced. The selective analysis of components of stress granules and P-bodies revealed a prominent downregulation of GIS2. Our yeast data are in good agreement with a global proteomics and metabolomics publication that the PBP1 ortholog ATAXIN-2 (ATXN2) knockout (KO) in mouse results in mitochondrial deficits in leucine/fatty acid catabolism and bioenergetics, with an obesity phenotype. Furthermore, our data provide the completely novel insight that PBP1 mutations in stress periods involve GIS2, a plausible scenario in view of previous data that both PBP1 and GIS2 relocalize from ribosomes to stress granules, interact with poly(A)-binding protein in translation regulation and prevent mitochondrial precursor overaccumulation stress (mPOS). This may be relevant for human diseases like spinocerebellar ataxias, amyotrophic lateral sclerosis, and the metabolic syndrome.

PRMT1-mediated arginine methylation controls ATXN2L localization.

Arginine methylation is a posttranslational modification that is of importance in diverse cellular processes. Recent proteomic mass spectrometry studies reported arginine methylation of ataxin-2-like (ATXN2L), the paralog of ataxin-2, a protein that is implicated in the neurodegenerative disorder spinocerebellar ataxia type 2. Here, we investigated the methylation state of ATXN2L and its significance for ATXN2L localization. We first confirmed that ATXN2L is asymmetrically dimethylated in vivo, and observed that the nuclear localization of ATXN2L is altered under methylation inhibition. We further discovered that ATXN2L associates with the protein arginine-N-methyltransferase 1 (PRMT1). Finally, we showed that neither mutation of the arginine-glycine-rich motifs of ATXN2L nor methylation inhibition alters ATXN2L localization to stress granules, suggesting that methylation of ATXN2L is probably not mandatory.

5-Fluorouracil affects assembly of stress granules based on RNA incorporation.

The antimetabolite 5-fluorouracil is a widely used chemotherapeutic for the treatment of several solid cancers. However, resistance to 5-fluorouracil remains a major drawback in its clinical use. In this study we report that treatment of HeLa cells with 5-fluorouracil resulted in de novo assembly of stress granules. Moreover, we revealed that stress granule assembly under stress conditions as well as disassembly is altered in cells treated with 5-fluorouracil. Notably, we discovered that RACK1, a protein mediating cell survival and apoptosis, is a component of 5-fluorouracil-induced stress granules. To explore the mode of action of 5-fluorouracil accountable for de novo stress granule assembly, we analyzed 5-fluorouracil metabolites and noticed that stress granule assembly is caused by RNA, not DNA incorporating 5-fluorouracil metabolites. Interestingly, we observed that other RNA incorporating drugs also cause assembly of stress granules. Thus, our results suggest that incorporation of chemotherapeutics into RNA may result in stress granule assembly with potential significance in chemoresistance.

High-throughput miRNA and mRNA sequencing of paired colorectal normal, tumor and metastasis tissues and bioinformatic modeling of miRNA-1 therapeutic applications.

MiRNAs are discussed as diagnostic and therapeutic molecules. However, effective miRNA drug treatments with miRNAs are, so far, hampered by the complexity of the miRNA networks. To identify potential miRNA drugs in colorectal cancer, we profiled miRNA and mRNA expression in matching normal, tumor and metastasis tissues of eight patients by Illumina sequencing. We validated six miRNAs in a large tissue screen containing 16 additional tumor entities and identified miRNA-1, miRNA-129, miRNA-497 and miRNA-215 as constantly de-regulated within the majority of cancers. Of these, we investigated miRNA-1 as representative in a systems-biology simulation of cellular cancer models implemented in PyBioS and assessed the effects of depletion as well as overexpression in terms of miRNA-1 as a potential treatment option. In this system, miRNA-1 treatment reverted the disease phenotype with different effectiveness among the patients. Scoring the gene expression changes obtained through mRNA-Seq from the same patients we show that the combination of deep sequencing and systems biological modeling can help to identify patient-specific responses to miRNA treatments. We present this data as guideline for future pre-clinical assessments of new and personalized therapeutic options.

Fkh1 and Fkh2 associate with Sir2 to control CLB2 transcription under normal and oxidative stress conditions.

The Forkhead (Fkh) box family of transcription factors is evolutionary conserved from yeast to higher eukaryotes and its members are involved in many physiological processes including metabolism, DNA repair, cell cycle, stress resistance, apoptosis, and aging. In budding yeast, four Fkh transcription factors were identified, namely Fkh1, Fkh2, Fhl1, and Hcm1, which are implicated in chromatin silencing, cell cycle regulation, and stress response. These factors impinge transcriptional regulation during cell cycle progression, and histone deacetylases (HDACs) play an essential role in this process, e.g., the nuclear localization of Hcm1 depends on Sir2 activity, whereas Sin3/Rpd3 silence cell cycle specific gene transcription in G2/M phase. However, a direct involvement of Sir2 in Fkh1/Fkh2-dependent regulation of target genes is at present unknown. Here, we show that Fkh1 and Fkh2 associate with Sir2 in G1 and M phase, and that Fkh1/Fkh2-mediated activation of reporter genes is antagonized by Sir2. We further report that Sir2 overexpression strongly affects cell growth in an Fkh1/Fkh2-dependent manner. In addition, Sir2 regulates the expression of the mitotic cyclin Clb2 through Fkh1/Fkh2-mediated binding to the CLB2 promoter in G1 and M phase. We finally demonstrate that Sir2 is also enriched at the CLB2 promoter under stress conditions, and that the nuclear localization of Sir2 is dependent on Fkh1 and Fkh2. Taken together, our results show a functional interplay between Fkh1/Fkh2 and Sir2 suggesting a novel mechanism of cell cycle repression. Thus, in budding yeast, not only the regulation of G2/M gene expression but also the protective response against stress could be directly coordinated by Fkh1 and Fkh2.

Ataxin-2-like is a regulator of stress granules and processing bodies.

Paralogs for several proteins implicated in neurodegenerative disorders have been identified and explored to further facilitate the identification of molecular mechanisms contributing to disease pathogenesis. For the disease-causing protein in spinocerebellar ataxia type 2, ataxin-2, a paralog of unknown function, termed ataxin-2-like, has been described. We discovered that ataxin-2-like associates with known interaction partners of ataxin-2, the RNA helicase DDX6 and the poly(A)-binding protein, and with ataxin-2 itself. Furthermore, we found that ataxin-2-like is a component of stress granules. Interestingly, sole ataxin-2-like overexpression led to the induction of stress granules, while a reduction of stress granules was detected in case of a low ataxin-2-like level. Finally, we observed that overexpression of ataxin-2-like as well as its reduction has an impact on the presence of microscopically visible processing bodies. Thus, our results imply a functional overlap between ataxin-2-like and ataxin-2, and further indicate a role for ataxin-2-like in the regulation of stress granules and processing bodies.

FOX-2 dependent splicing of ataxin-2 transcript is affected by ataxin-1 overexpression.

Alternative splicing is a fundamental posttranscriptional mechanism for controlling gene expression, and splicing defects have been linked to various human disorders. The splicing factor FOX-2 is part of a main protein interaction hub in a network related to human inherited ataxias, however, its impact remains to be elucidated. Here, we focused on the reported interaction between FOX-2 and ataxin-1, the disease-causing protein in spinocerebellar ataxia type 1. In this line, we further evaluated this interaction by yeast-2-hybrid analyses and co-immunoprecipitation experiments in mammalian cells. Interestingly, we discovered that FOX-2 localization and splicing activity is affected in the presence of nuclear ataxin-1 inclusions. Moreover, we observed that FOX-2 directly interacts with ataxin-2, a protein modulating spinocerebellar ataxia type 1 pathogenesis. Finally, we provide evidence that splicing of pre-mRNA of ataxin-2 depends on FOX-2 activity, since reduction of FOX-2 levels led to increased skipping of exon 18 in ataxin-2 transcripts. Most striking, we observed that ataxin-1 overexpression has an effect on this splicing event as well. Thus, our results demonstrate that FOX-2 is involved in splicing of ataxin-2 transcripts and that this splicing event is altered by overexpression of ataxin-1.

Sic1 plays a role in timing and oscillatory behaviour of B-type cyclins.

Budding yeast cell cycle oscillates between states of low and high cyclin-dependent kinase activity, driven by association of Cdk1 with B-type (Clb) cyclins. Various Cdk1-Clb complexes are activated and inactivated in a fixed, temporally regulated sequence, inducing the behaviour known as "waves of cyclins". The transition from low to high Clb activity is triggered by degradation of Sic1, the inhibitor of Cdk1-Clb complexes, at the entry to S phase. The G(1) phase is characterized by low Clb activity and high Sic1 levels. High Clb activity and Sic1 proteolysis are found from the beginning of the S phase until the end of mitosis. The mechanism regulating the appearance on schedule of Cdk1-Clb complexes is currently unknown. Here, we analyse oscillations of Clbs, focusing on the role of their inhibitor Sic1. We compare mathematical networks differing in interactions that Sic1 may establish with Cdk1-Clb complexes. Our analysis suggests that the wave-like cyclins pattern derives from the binding of Sic1 to all Clb pairs rather than from Clb degradation. These predictions are experimentally validated, showing that Sic1 indeed interacts and coexists in time with Clbs. Intriguingly, a sic1Δ strain looses cell cycle-regulated periodicity of Clbs, which is observed in the wild type, whether a SIC1-0P strain delays the formation of Clb waves. Our results highlight an additional role for Sic1 in regulating Cdk1-Clb complexes, coordinating their appearance.

Targeted high throughput sequencing in clinical cancer settings: formaldehyde fixed-paraffin embedded (FFPE) tumor tissues, input amount and tumor heterogeneity.

BACKGROUND: Massively parallel sequencing technologies have brought an enormous increase in sequencing throughput. However, these technologies need to be further improved with regard to reproducibility and applicability to clinical samples and settings. METHODS: Using identification of genetic variations in prostate cancer as an example we address three crucial challenges in the field of targeted re-sequencing: Small nucleotide variation (SNV) detection in samples of formalin-fixed paraffin embedded (FFPE) tissue material, minimal amount of input sample and sampling in view of tissue heterogeneity. RESULTS: We show that FFPE tissue material can supplement for fresh frozen tissues for the detection of SNVs and that solution-based enrichment experiments can be accomplished with small amounts of DNA with only minimal effects on enrichment uniformity and data variance.Finally, we address the question whether the heterogeneity of a tumor is reflected by different genetic alterations, e.g. different foci of a tumor display different genomic patterns. We show that the tumor heterogeneity plays an important role for the detection of copy number variations. CONCLUSIONS: The application of high throughput sequencing technologies in cancer genomics opens up a new dimension for the identification of disease mechanisms. In particular the ability to use small amounts of FFPE samples available from surgical tumor resections and histopathological examinations facilitates the collection of precious tissue materials. However, care needs to be taken in regard to the locations of the biopsies, which can have an influence on the prediction of copy number variations. Bearing these technological challenges in mind will significantly improve many large-scale sequencing studies and will - in the long term - result in a more reliable prediction of individual cancer therapies.

Ortho2ExpressMatrix--a web server that interprets cross-species gene expression data by gene family information.

BACKGROUND: The study of gene families is pivotal for the understanding of gene evolution across different organisms and such phylogenetic background is often used to infer biochemical functions of genes. Modern high-throughput experiments offer the possibility to analyze the entire transcriptome of an organism; however, it is often difficult to deduct functional information from that data. RESULTS: To improve functional interpretation of gene expression we introduce Ortho2ExpressMatrix, a novel tool that integrates complex gene family information, computed from sequence similarity, with comparative gene expression profiles of two pre-selected biological objects: gene families are displayed with two-dimensional matrices. Parameters of the tool are object type (two organisms, two individuals, two tissues, etc.), type of computational gene family inference, experimental meta-data, microarray platform, gene annotation level and genome build. Family information in Ortho2ExpressMatrix bases on computationally different protein family approaches such as EnsemblCompara, InParanoid, SYSTERS and Ensembl Family. Currently, respective all-against-all associations are available for five species: human, mouse, worm, fruit fly and yeast. Additionally, microRNA expression can be examined with respect to miRBase or TargetScan families. The visualization, which is typical for Ortho2ExpressMatrix, is performed as matrix view that displays functional traits of genes (differential expression) as well as sequence similarity of protein family members (BLAST e-values) in colour codes. Such translations are intended to facilitate the user's perception of the research object. CONCLUSIONS: Ortho2ExpressMatrix integrates gene family information with genome-wide expression data in order to enhance functional interpretation of high-throughput analyses on diseases, environmental factors, or genetic modification or compound treatment experiments. The tool explores differential gene expression in the light of orthology, paralogy and structure of gene families up to the point of ambiguity analyses. Results can be used for filtering and prioritization in functional genomic, biomedical and systems biology applications. The web server is freely accessible at http://bioinf-data.charite.de/o2em/cgi-bin/o2em.pl.

Effects of EpCAM overexpression on human breast cancer cell lines.

BACKGROUND: Recently, EpCAM has attracted major interest as a target for antibody- and vaccine-based cancer immunotherapies. In breast cancer, the EpCAM antigen is overexpressed in 30-40% of all cases and this increased expression correlates with poor prognosis. The use of EpCAM-specific monoclonal antibodies is a promising treatment approach in these patients. METHODS: In order to explore molecular changes following EpCAM overexpression, we investigated changes of the transcriptome upon EpCAM gene expression in commercially available human breast cancer cells lines Hs578T and MDA-MB-231. To assess cell proliferation, a tetrazolium salt based assay was performed. A TCF/LEF Reporter Kit was used to measure the transcriptional activity of the Wnt/ß-catenin pathway. To evaluate the accumulation of ß-catenin in the nucleus, a subcellular fractionation assay was performed. RESULTS: For the first time we could show that expression profiling data of EpCAM transfected cell lines Hs578TEpCAM and MDA-MB-231EpCAM indicate an association of EpCAM overexpression with the downregulation of the Wnt signaling inhibitors SFRP1 and TCF7L2. Confirmation of increased Wnt signaling was provided by a TCF/LEF reporter kit and by the finding of the nuclear accumulation of ß-catenin for MDA-MB-231 EpCAM but not Hs578T EpCAM cells. In Hs578T cells, an increase of proliferation and chemosensitivity to Docetaxel was associated with EpCAM overexpression. CONCLUSIONS: These data show a cell type dependent modification of Wnt signaling components after EpCAM overexpression in breast cancer cell lines, which results in marginal functional changes. Further investigations on the interaction of EpCAM with SFRP1 and TCF7L2 and on additional factors, which may be causal for changes upon EpCAM overexpression, will help to characterize unique molecular properties of EpCAM-positive breast cancer cells.

The KRAB-containing zinc-finger transcriptional regulator ZBRK1 activates SCA2 gene transcription through direct interaction with its gene product, ataxin-2.

Gene transcription is controlled by transcriptional regulators acting with specific co-regulators to allow gene activation and repression. Here, we report the identification of the KRAB-containing zinc-finger transcriptional regulator, ZBRK1, as an interaction partner of the SCA2 gene product ataxin-2. Furthermore, we discovered that an elevated ZBRK1 level resulted in increased ataxin-2 levels, whereas interference on transcriptional and protein levels of ZBRK1 yielded reduced ataxin-2 levels, suggesting that a complex comprising ZBRK1 and ataxin-2 regulates SCA2 gene transcription. A bioinformatic analysis utilizing the known ZBRK1 consensus DNA-binding motif revealed ZBRK1-binding sites in the SCA2 promoter. These predicted sites were experimentally validated by chromatin-immunoprecipitation experiments along with luciferase-based promoter analyses corroborating that SCA2 gene transcription is controlled by a ZBRK1/ataxin-2 complex. Finally, we demonstrate that SCA2 gene transcription is significantly reduced in colon tumors possessing low ZBRK1 transcripts. Thus, our results provide first evidence that ataxin-2 acts as a co-regulator of ZBRK1 activating its own transcription, thereby representing the first identified ZBRK1 co-activator.

Huntington's disease: From molecular basis to therapeutic advances.

Huntington's disease is an autosomal dominant genetic neurodegenerative disorder, which is characterized by progressive motor dysfunction, emotional disturbances, dementia, and weight loss. The disease is caused by pathological CAG-triplet repeat extension(s), encoding polyglutamines, within the gene product, huntingtin. Huntingtin is ubiquitously expressed through the body and is a protein of uncertain molecular function(s). Mutant huntingtin, containing pathologically extended polyglutamines causes the earliest and most dramatic neuropathologic changes in the neostriatum and cerebral cortex. Extended polyglutamines confer structural conformational changes to huntingtin, which gains novel properties, resulting in aberrant interactions with multiple cellular components. The diverse and variable aberrations mediated by mutant huntingtin perturb many cellular functions essential for neuronal homeostasis and underlie pleiotropic mechanisms of Huntington's disease pathogenesis. The only approved drug for Huntington's disease is a symptomatic treatment, tetrabenazine; thus, novel neuroprotective strategies, slowing, blocking and possibly reversing disease progression, are vital for developing effective therapies.