Parkinson's disease and schizophrenia interactomes contain temporally distinct gene clusters underlying comorbid mechanisms and unique disease processes.

Genome-wide association studies suggest significant overlaps in Parkinson's disease (PD) and schizophrenia (SZ) risks, but the underlying mechanisms remain elusive. The protein-protein interaction network ('interactome') plays a crucial role in PD and SZ and can incorporate their spatiotemporal specificities. Therefore, to study the linked biology of PD and SZ, we compiled PD- and SZ-associated genes from the DisGeNET database, and constructed their interactomes using BioGRID and HPRD. We examined the interactomes using clustering and enrichment analyses, in conjunction with the transcriptomic data of 26 brain regions spanning foetal stages to adulthood available in the BrainSpan Atlas. PD and SZ interactomes formed four gene clusters with distinct temporal identities (Disease Gene Networks or 'DGNs'1-4). DGN1 had unique SZ interactome genes highly expressed across developmental stages, corresponding to a neurodevelopmental SZ subtype. DGN2, containing unique SZ interactome genes expressed from early infancy to adulthood, correlated with an inflammation-driven SZ subtype and adult SZ risk. DGN3 contained unique PD interactome genes expressed in late infancy, early and late childhood, and adulthood, and involved in mitochondrial pathways. DGN4, containing prenatally-expressed genes common to both the interactomes, involved in stem cell pluripotency and overlapping with the interactome of 22q11 deletion syndrome (comorbid psychosis and Parkinsonism), potentially regulates neurodevelopmental mechanisms in PD-SZ comorbidity. Our findings suggest that disrupted neurodevelopment (regulated by DGN4) could expose risk windows in PD and SZ, later elevating disease risk through inflammation (DGN2). Alternatively, variant clustering in DGNs may produce disease subtypes, e.g., PD-SZ comorbidity with DGN4, and early/late-onset SZ with DGN1/DGN2.

Flow blockage disrupts cilia-driven fluid transport in the epileptic brain.

A carpet of ependymal motile cilia lines the brain ventricular system, forming a network of flow channels and barriers that pattern cerebrospinal fluid (CSF) flow at the surface. This CSF transport system is evolutionary conserved, but its physiological function remains unknown. Here we investigated its potential role in epilepsy with studies focused on CDKL5 deficiency disorder (CDD), a neurodevelopmental disorder with early-onset epilepsy refractory to seizure medications and the most common cause of infant epilepsy. CDKL5 is a highly conserved X-linked gene suggesting its function in regulating cilia length and motion in the green alga Chlamydomonas might have implication in the etiology of CDD. Examination of the structure and function of airway motile cilia revealed both the CDD patients and the Cdkl5 knockout mice exhibit cilia lengthening and abnormal cilia motion. Similar defects were observed for brain ventricular cilia in the Cdkl5 knockout mice. Mapping ependymal cilia generated flow in the ventral third ventricle (v3V), a brain region with important physiological functions showed altered patterning of flow. Tracing of cilia-mediated inflow into v3V with fluorescent dye revealed the appearance of a flow barrier at the inlet of v3V in Cdkl5 knockout mice. Analysis of mice with a mutation in another epilepsy-associated kinase, Yes1, showed the same disturbance of cilia motion and flow patterning. The flow barrier was also observed in the Foxj1± and FOXJ1CreERT:Cdkl5y/fl mice, confirming the contribution of ventricular cilia to the flow disturbances. Importantly, mice exhibiting altered cilia-driven flow also showed increased susceptibility to anesthesia-induced seizure-like activity. The cilia-driven flow disturbance arises from altered cilia beating orientation with the disrupted polarity of the cilia anchoring rootlet meshwork. Together these findings indicate motile cilia disturbances have an essential role in CDD-associated seizures and beyond, suggesting cilia regulating kinases may be a therapeutic target for medication-resistant epilepsy.

Novel Protein-Protein Interactions Highlighting the Crosstalk between Hypoplastic Left Heart Syndrome, Ciliopathies and Neurodevelopmental Delays.

Hypoplastic left heart syndrome (HLHS) is a severe congenital heart disease (CHD) affecting 1 in 5000 newborns. We constructed the interactome of 74 HLHS-associated genes identified from a large-scale mouse mutagenesis screen, augmenting it with 408 novel protein-protein interactions (PPIs) using our High-Precision Protein-Protein Interaction Prediction (HiPPIP) model. The interactome is available on a webserver with advanced search capabilities. A total of 364 genes including 73 novel interactors were differentially regulated in tissue/iPSC-derived cardiomyocytes of HLHS patients. Novel PPIs facilitated the identification of TOR signaling and endoplasmic reticulum stress modules. We found that 60.5% of the interactome consisted of housekeeping genes that may harbor large-effect mutations and drive HLHS etiology but show limited transmission. Network proximity of diabetes, Alzheimer's disease, and liver carcinoma-associated genes to HLHS genes suggested a mechanistic basis for their comorbidity with HLHS. Interactome genes showed tissue-specificity for sites of extracardiac anomalies (placenta, liver and brain). The HLHS interactome shared significant overlaps with the interactomes of ciliopathy- and microcephaly-associated genes, with the shared genes enriched for genes involved in intellectual disability and/or developmental delay, and neuronal death pathways, respectively. This supported the increased burden of ciliopathy variants and prevalence of neurological abnormalities observed among HLHS patients with developmental delay and microcephaly, respectively.

Uncompensated mitochondrial oxidative stress underlies heart failure in an iPSC-derived model of congenital heart disease.

Hypoplastic left heart syndrome (HLHS) is a severe congenital heart disease with 30% mortality from heart failure (HF) in the first year of life, but the cause of early HF remains unknown. Induced pluripotent stem-cell-derived cardiomyocytes (iPSC-CM) from patients with HLHS showed that early HF is associated with increased apoptosis, mitochondrial respiration defects, and redox stress from abnormal mitochondrial permeability transition pore (mPTP) opening and failed antioxidant response. In contrast, iPSC-CM from patients without early HF showed normal respiration with elevated antioxidant response. Single-cell transcriptomics confirmed that early HF is associated with mitochondrial dysfunction accompanied with endoplasmic reticulum (ER) stress. These findings indicate that uncompensated oxidative stress underlies early HF in HLHS. Importantly, mitochondrial respiration defects, oxidative stress, and apoptosis were rescued by treatment with sildenafil to inhibit mPTP opening or TUDCA to suppress ER stress. Together these findings point to the potential use of patient iPSC-CM for modeling clinical heart failure and the development of therapeutics.

Circadian transcription factor NPAS2 and the NAD+ -dependent deacetylase SIRT1 interact in the mouse nucleus accumbens and regulate reward.

Substance use disorders are associated with disruptions to both circadian rhythms and cellular metabolic state. At the molecular level, the circadian molecular clock and cellular metabolic state may be interconnected through interactions with the nicotinamide adenine dinucleotide (NAD+ )-dependent deacetylase, sirtuin 1 (SIRT1). In the nucleus accumbens (NAc), a region important for reward, both SIRT1 and the circadian transcription factor neuronal PAS domain protein 2 (NPAS2) are highly enriched, and both are regulated by the metabolic cofactor NAD+ . Substances of abuse, like cocaine, greatly disrupt cellular metabolism and promote oxidative stress; however, their effects on NAD+ in the brain remain unclear. Interestingly, cocaine also induces NAc expression of both NPAS2 and SIRT1, and both have independently been shown to regulate cocaine reward in mice. However, whether NPAS2 and SIRT1 interact in the NAc and/or whether together they regulate reward is unknown. Here, we demonstrate diurnal expression of Npas2, Sirt1 and NAD+ in the NAc, which is altered by cocaine-induced upregulation. Additionally, co-immunoprecipitation reveals NPAS2 and SIRT1 interact in the NAc, and cross-analysis of NPAS2 and SIRT1 chromatin immunoprecipitation sequencing reveals several reward-relevant and metabolic-related pathways enriched among shared gene targets. Notably, NAc-specific Npas2 knock-down or a functional Npas2 mutation in mice attenuates SIRT1-mediated increases in cocaine preference. Together, our data reveal an interaction between NPAS2 and SIRT1 in the NAc, which may serve to integrate cocaine's effects on circadian and metabolic factors, leading to regulation of drug reward.

Generalized and social anxiety disorder interactomes show distinctive overlaps with striosome and matrix interactomes.

Mechanisms underlying anxiety disorders remain elusive despite the discovery of several associated genes. We constructed the protein-protein interaction networks (interactomes) of six anxiety disorders and noted enrichment for striatal expression among common genes in the interactomes. Five of these interactomes shared distinctive overlaps with the interactomes of genes that were differentially expressed in two striatal compartments (striosomes and matrix). Generalized anxiety disorder and social anxiety disorder interactomes showed exclusive and statistically significant overlaps with the striosome and matrix interactomes, respectively. Systematic gene expression analysis with the anxiety disorder interactomes constrained to contain only those genes that were shared with striatal compartment interactomes revealed a bifurcation among the disorders, which was influenced by the anterior cingulate cortex, nucleus accumbens, amygdala and hippocampus, and the dopaminergic signaling pathway. Our results indicate that the functionally distinct striatal pathways constituted by the striosome and the matrix may influence the etiological differentiation of various anxiety disorders.

Malignant Pleural Mesothelioma Interactome with 364 Novel Protein-Protein Interactions.

Malignant pleural mesothelioma (MPM) is an aggressive cancer affecting the outer lining of the lung, with a median survival of less than one year. We constructed an 'MPM interactome' with over 300 computationally predicted protein-protein interactions (PPIs) and over 2400 known PPIs of 62 literature-curated genes whose activity affects MPM. Known PPIs of the 62 MPM associated genes were derived from Biological General Repository for Interaction Datasets (BioGRID) and Human Protein Reference Database (HPRD). Novel PPIs were predicted by applying the HiPPIP algorithm, which computes features of protein pairs such as cellular localization, molecular function, biological process membership, genomic location of the gene, and gene expression in microarray experiments, and classifies the pairwise features as interacting or non-interacting based on a random forest model. We validated five novel predicted PPIs experimentally. The interactome is significantly enriched with genes differentially ex-pressed in MPM tumors compared with normal pleura and with other thoracic tumors, genes whose high expression has been correlated with unfavorable prognosis in lung cancer, genes differentially expressed on crocidolite exposure, and exosome-derived proteins identified from malignant mesothelioma cell lines. 28 of the interactors of MPM proteins are targets of 147 U.S. Food and Drug Administration (FDA)-approved drugs. By comparing disease-associated versus drug-induced differential expression profiles, we identified five potentially repurposable drugs, namely cabazitaxel, primaquine, pyrimethamine, trimethoprim and gliclazide. Preclinical studies may be con-ducted in vitro to validate these computational results. Interactome analysis of disease-associated genes is a powerful approach with high translational impact. It shows how MPM-associated genes identified by various high throughput studies are functionally linked, leading to clinically translatable results such as repurposed drugs. The PPIs are made available on a webserver with interactive user interface, visualization and advanced search capabilities.

Benchmark Evaluation of Protein-Protein Interaction Prediction Algorithms.

Protein-protein interactions (PPIs) perform various functions and regulate processes throughout cells. Knowledge of the full network of PPIs is vital to biomedical research, but most of the PPIs are still unknown. As it is infeasible to discover all of them experimentally due to technical and resource limitations, computational prediction of PPIs is essential and accurately assessing the performance of algorithms is required before further application or translation. However, many published methods compose their evaluation datasets incorrectly, using a higher proportion of positive class data than occuring naturally, leading to exaggerated performance. We re-implemented various published algorithms and evaluated them on datasets with realistic data compositions and found that their performance is overstated in original publications; with several methods outperformed by our control models built on 'illogical' and random number features. We conclude that these methods are influenced by an over-characterization of some proteins in the literature and due to scale-free nature of PPI network and that they fail when tested on all possible protein pairs. Additionally, we found that sequence-only-based algorithms performed worse than those that employ functional and expression features. We present a benchmark evaluation of many published algorithms for PPI prediction. The source code of our implementations and the benchmark datasets created here are made available in open source.

A reference catalog of DNA palindromes in the human genome and their variations in 1000 Genomes.

A palindrome in DNA is like a palindrome in language, but when read backwards, it is a complement of the forward sequence; effectively, the two halves of a sequence complement each other from its midpoint like in a double strand of DNA. Palindromes are distributed throughout the human genome and play significant roles in gene expression and regulation. Palindromic mutations are linked to many human diseases, such as neuronal disorders, mental retardation, and various cancers. In this work, we computed and analyzed the palindromic sequences in the human genome and studied their conservation in personal genomes using 1000 Genomes data. We found that ~30% of the palindromes exhibit variation, some of which are caused by rare variants. The analysis of disease/trait-associated single-nucleotide polymorphisms in palindromic regions showed that disease-associated risk variants are 14 times more likely to be present in palindromic regions than in other regions. The catalog of palindromes in the reference genome and 1000 Genomes is being made available here with details on their variations in each individual genome to serve as a resource for future and retrospective whole-genome studies identifying statistically significant palindrome variations associated with diseases or traits and their roles in disease mechanisms.

Cilia interactome with predicted protein-protein interactions reveals connections to Alzheimer's disease, aging and other neuropsychiatric processes.

Cilia are dynamic microtubule-based organelles present on the surface of many eukaryotic cell types and can be motile or non-motile primary cilia. Cilia defects underlie a growing list of human disorders, collectively called ciliopathies, with overlapping phenotypes such as developmental delays and cognitive and memory deficits. Consistent with this, cilia play an important role in brain development, particularly in neurogenesis and neuronal migration. These findings suggest that a deeper systems-level understanding of how ciliary proteins function together may provide new mechanistic insights into the molecular etiologies of nervous system defects. Towards this end, we performed a protein-protein interaction (PPI) network analysis of known intraflagellar transport, BBSome, transition zone, ciliary membrane and motile cilia proteins. Known PPIs of ciliary proteins were assembled from online databases. Novel PPIs were predicted for each ciliary protein using a computational method we developed, called High-precision PPI Prediction (HiPPIP) model. The resulting cilia "interactome" consists of 165 ciliary proteins, 1,011 known PPIs, and 765 novel PPIs. The cilia interactome revealed interconnections between ciliary proteins, and their relation to several pathways related to neuropsychiatric processes, and to drug targets. Approximately 184 genes in the cilia interactome are targeted by 548 currently approved drugs, of which 103 are used to treat various diseases of nervous system origin. Taken together, the cilia interactome presented here provides novel insights into the relationship between ciliary protein dysfunction and neuropsychiatric disorders, for e.g. interconnections of Alzheimer's disease, aging and cilia genes. These results provide the framework for the rational design of new therapeutic agents for treatment of ciliopathies and neuropsychiatric disorders.

Interactome of SARS-CoV-2 / nCoV19 modulated host proteins with computationally predicted PPIs.

World over, people are looking for solutions to tackle the pandemic coronavirus disease (COVID-19) caused by the virus SARS-CoV-2/nCoV-19. Notable contributions in biomedical field have been characterizing viral genomes, host transcriptomes and proteomes, repurposable drugs and vaccines. In one such study, 332 human proteins targeted by nCoV19 were identified. We expanded this set of host proteins by constructing their protein interactome, including in it not only the known protein-protein interactions (PPIs) but also novel, hitherto unknown PPIs predicted with our High-precision Protein-Protein Interaction Prediction (HiPPIP) model that was shown to be highly accurate. In fact, one of the earliest discoveries made possible by HiPPIP is related to activation of immunity upon viral infection. We found that several interactors of the host proteins are differentially expressed upon viral infection, are related to highly relevant pathways, and that the novel interaction of NUP98 with CHMP5 may activate an antiviral mechanism leading to disruption of viral budding. We are making the interactions available as downloadable files to facilitate future systems biology studies and also on a web-server at http://hagrid.dbmi.pitt.edu/corona that allows not only keyword search but also queries such as "PPIs where one protein is associated with 'virus' and the interactors with 'pulmonary'".

Potentially repurposable drugs for schizophrenia identified from its interactome.

We previously presented the protein-protein interaction network of schizophrenia associated genes, and from it, the drug-protein interactome which showed the drugs that target any of the proteins in the interactome. Here, we studied these drugs further to identify whether any of them may potentially be repurposable for schizophrenia. In schizophrenia, gene expression has been described as a measurable aspect of the disease reflecting the action of risk genes. We studied each of the drugs from the interactome using the BaseSpace Correlation Engine, and shortlisted those that had a negative correlation with differential gene expression of schizophrenia. This analysis resulted in 12 drugs whose differential gene expression (drug versus normal) had an anti-correlation with differential expression for schizophrenia (disorder versus normal). Some of these drugs were already being tested for their clinical activity in schizophrenia and other neuropsychiatric disorders. Several proteins in the protein interactome of the targets of several of these drugs were associated with various neuropsychiatric disorders. The network of genes with opposite drug-induced versus schizophrenia-associated expression profiles were significantly enriched in pathways relevant to schizophrenia etiology and GWAS genes associated with traits or diseases that had a pathophysiological overlap with schizophrenia. Drugs that targeted the same genes as the shortlisted drugs, have also demonstrated clinical activity in schizophrenia and other related disorders. This integrated computational analysis will help translate insights from the schizophrenia drug-protein interactome to clinical research - an important step, especially in the field of psychiatric drug development which faces a high failure rate.

Towards Extracting Supporting Information About Predicted Protein-Protein Interactions.

One of the goals of relation extraction is to identify protein-protein interactions (PPIs) in biomedical literature. Current systems are capturing binary relations and also the direction and type of an interaction. Besides assisting in the curation PPIs into databases, there has been little real-world application of these algorithms. We describe UPSITE, a text mining tool for extracting evidence in support of a hypothesized interaction. Given a predicted PPI, UPSITE uses a binary relation detector to check whether a PPI is found in abstracts in PubMed. If it is not found, UPSITE retrieves documents relevant to each of the two proteins separately, and extracts contextual information about biological events surrounding each protein, and calculates semantic similarity of the two proteins to provide evidential support for the predicted PPI. In evaluations, relation extraction achieved an Fscore of 0.88 on the HPRD50 corpus, and semantic similarity measured with angular distance was found to be statistically significant. With the development of PPI prediction algorithms, the burden of interpreting the validity and relevance of novel PPIs is on biologists. We suggest that presenting annotations of the two proteins in a PPI side-by-side and a score that quantifies their similarity lessens this burden to some extent.

Generating testable hypotheses for schizophrenia and rheumatoid arthritis pathogenesis by integrating epidemiological, genomic, and protein interaction data.

Patients with schizophrenia and their relatives have reduced prevalence of rheumatoid arthritis. Schizophrenia and rheumatoid arthritis genome-wide association studies also indicate negative genetic correlations, suggesting that there may be shared pathogenesis at the DNA level or downstream. A portion of the inverse prevalence could be attributed to pleiotropy, i.e., variants of a single nucleotide polymorphism that could confer differential risk for these disorders. To study the basis for such an interrelationship, we initially compared lists of single nucleotide polymorphisms with significant genetic associations (p < 1e-8) for schizophrenia or rheumatoid arthritis, evaluating patterns of linkage disequilibrium and apparent pleiotropic risk profiles. Single nucleotide polymorphisms that conferred risk for both schizophrenia and rheumatoid arthritis were localized solely to the extended HLA region. Among single nucleotide polymorphisms that conferred differential risk for schizophrenia and rheumatoid arthritis, the majority were localized to HLA-B, TNXB, NOTCH4, HLA-C, HCP5, MICB, PSORS1C1, and C6orf10; published functional data indicate that HLA-B and HLA-C have the most plausible pathogenic roles in both disorders. Interactomes of these eight genes were constructed from protein-protein interaction information using publicly available databases and novel computational predictions. The genes harboring apparently pleiotropic single nucleotide polymorphisms are closely connected to rheumatoid arthritis and schizophrenia associated genes through common interacting partners. A separate and independent analysis of the interactomes of rheumatoid arthritis and schizophrenia genes showed a significant overlap between the two interactomes and that they share several common pathways, motivating functional studies suggesting a relationship in the pathogenesis of schizophrenia/rheumatoid arthritis.

The complex genetics of hypoplastic left heart syndrome.

Congenital heart disease (CHD) affects up to 1% of live births. Although a genetic etiology is indicated by an increased recurrence risk, sporadic occurrence suggests that CHD genetics is complex. Here, we show that hypoplastic left heart syndrome (HLHS), a severe CHD, is multigenic and genetically heterogeneous. Using mouse forward genetics, we report what is, to our knowledge, the first isolation of HLHS mutant mice and identification of genes causing HLHS. Mutations from seven HLHS mouse lines showed multigenic enrichment in ten human chromosome regions linked to HLHS. Mutations in Sap130 and Pcdha9, genes not previously associated with CHD, were validated by CRISPR-Cas9 genome editing in mice as being digenic causes of HLHS. We also identified one subject with HLHS with SAP130 and PCDHA13 mutations. Mouse and zebrafish modeling showed that Sap130 mediates left ventricular hypoplasia, whereas Pcdha9 increases penetrance of aortic valve abnormalities, both signature HLHS defects. These findings show that HLHS can arise genetically in a combinatorial fashion, thus providing a new paradigm for the complex genetics of CHD.

Antibody Exchange: Information extraction of biological antibody donation and a web-portal to find donors and seekers.

Bio-molecular reagents like antibodies required in experimental biology are expensive and their effectiveness, among other things, is critical to the success of the experiment. Although such resources are sometimes donated by one investigator to another through personal communication between the two, there is no previous study to our knowledge on the extent of such donations, nor a central platform that directs resource seekers to donors. In this paper, we describe, to our knowledge, a first attempt at building a web-portal titled Antibody Exchange (or more general 'Bio-Resource Exchange') that attempts to bridge this gap between resource seekers and donors in the domain of experimental biology. Users on this portal can request for or donate antibodies, cell-lines and DNA Constructs. This resource could also serve as a crowd-sourced database of resources for experimental biology. Further, we also studied the extent of antibody donations by mining the acknowledgement sections of scientific articles. Specifically, we extracted the name of the donor, his/her affiliation and the name of the antibody for every donation by parsing the acknowledgements sections of articles. To extract annotations at this level, we adopted two approaches - a rule based algorithm and a bootstrapped pattern learning algorithm. The algorithms extracted donor names, affiliations and antibody names with average accuracies of 57% and 62% respectively. We also created a dataset of 50 expert-annotated acknowledgements sections that will serve as a gold standard dataset to evaluate extraction algorithms in the future.

Schizophrenia interactome with 504 novel protein-protein interactions.

Genome-wide association studies of schizophrenia (GWAS) have revealed the role of rare and common genetic variants, but the functional effects of the risk variants remain to be understood. Protein interactome-based studies can facilitate the study of molecular mechanisms by which the risk genes relate to schizophrenia (SZ) genesis, but protein-protein interactions (PPIs) are unknown for many of the liability genes. We developed a computational model to discover PPIs, which is found to be highly accurate according to computational evaluations and experimental validations of selected PPIs. We present here, 365 novel PPIs of liability genes identified by the SZ Working Group of the Psychiatric Genomics Consortium (PGC). Seventeen genes that had no previously known interactions have 57 novel interactions by our method. Among the new interactors are 19 drug targets that are targeted by 130 drugs. In addition, we computed 147 novel PPIs of 25 candidate genes investigated in the pre-GWAS era. While there is little overlap between the GWAS genes and the pre-GWAS genes, the interactomes reveal that they largely belong to the same pathways, thus reconciling the apparent disparities between the GWAS and prior gene association studies. The interactome including 504 novel PPIs overall, could motivate other systems biology studies and trials with repurposed drugs. The PPIs are made available on a webserver, called Schizo-Pi at http://severus.dbmi.pitt.edu/schizo-pi with advanced search capabilities.

A pilot study on the prevalence of DNA palindromes in breast cancer genomes.

BACKGROUND: DNA palindromes are a unique pattern of repeat sequences that are present in the human genome. It consists of a sequence of nucleotides in which the second half is the complement of the first half but appearing in reverse order. These palindromic sequences may have a significant role in DNA replication, transcription and gene regulation processes. They occur frequently in human cancers by clustering at specific locations of the genome that undergo gene amplification and tumorigenesis. Moreover, some studies showed that palindromes are clustered in amplified regions of breast cancer genomes especially in chromosomes (chr) 8 and 11. With the large number of personal genomes and cancer genomes becoming available, it is now possible to study their association to diseases using computational methods. Here, we conducted a pilot study on chromosomes 8 and 11 of cancer genomes to identify computationally the differentially occurring palindromes. METHODS: We processed 69 breast cancer genomes from The Cancer Genome Atlas including serum-normal and tumor genomes, and 1000 Genomes to serve as control group. The Biological Language Modelling Toolkit (BLMT) computes palindromes in whole genomes. We developed a computational pipeline integrating BLMT to compute and compare prevalence of palindromes in personal genomes. RESULTS: We carried out a pilot study on chr 8 and chr 11 taking into account single nucleotide polymorphisms, insertions and deletions. Of all the palindromes that showed any variation in cancer genomes, 38% of what were near breast cancer genes happened to be the most differentiated palindromes in tumor (i.e. they ranked among the top 25% by our heuristic measure). CONCLUSIONS: These results will shed light on the prevalence of palindromes in oncogenes and the mutations that are present in the palindromic regions that could contribute to genomic rearrangements, and breast cancer progression.