RNA-Seq time-course analysis of neural precursor cell transcriptome in response to herpes simplex Virus-1 infection.

The neurogenic niches within the central nervous system serve as essential reservoirs for neural precursor cells (NPCs), playing a crucial role in neurogenesis. However, these NPCs are particularly vulnerable to infection by the herpes simplex virus 1 (HSV-1). In the present study, we investigated the changes in the transcriptome of NPCs in response to HSV-1 infection using bulk RNA-Seq, compared to those of uninfected samples, at different time points post infection and in the presence or absence of antivirals. The results showed that NPCs upon HSV-1 infection undergo a significant dysregulation of genes playing a crucial role in aspects of neurogenesis, including genes affecting NPC proliferation, migration, and differentiation. Our analysis revealed that the CREB signaling, which plays a crucial role in the regulation of neurogenesis and memory consolidation, was the most consistantly downregulated pathway, even in the presence of antivirals. Additionally, cholesterol biosynthesis was significantly downregulated in HSV-1-infected NPCs. The findings from this study, for the first time, offer insights into the intricate molecular mechanisms that underlie the neurogenesis impairment associated with HSV-1 infection.

Peripheral Ischemia Imprints Epigenetic Changes in Hematopoietic Stem Cells to Propagate Inflammation and Atherosclerosis.

BACKGROUND: Peripheral ischemia caused by peripheral artery disease is associated with systemic inflammation, which may aggravate underlying comorbidities such as atherosclerosis and heart failure. However, the mechanisms of increased inflammation and inflammatory cell production in patients with peripheral artery disease remain poorly understood. METHODS: We used peripheral blood collected from patients with peripheral artery disease and performed hind limb ischemia (HI) in Apoe-/- mice fed a Western diet and C57BL/6J mice with a standard laboratory diet. Bulk and single-cell RNA sequencing analysis, whole-mount microscopy, and flow cytometry were performed to analyze hematopoietic stem and progenitor cell (HSPC) proliferation, differentiation, and relocation. RESULTS: We observed augmented numbers of leukocytes in the blood of patients with peripheral artery disease and Apoe-/- mice with HI. RNA sequencing and whole-mount imaging of the bone marrow revealed HSPC migration into the vascular niche from the osteoblastic niche and their exaggerated proliferation and differentiation. Single-cell RNA sequencing demonstrated alterations in the genes responsible for inflammation, myeloid cell mobilization, and HSPC differentiation after HI. Heightened inflammation in Apoe-/- mice after HI aggravated atherosclerosis. Surprisingly, bone marrow HSPCs expressed higher amounts of the receptors for IL (interleukin)-1 and IL-3 after HI. Concomitantly, the promoters of Il1r1 and Il3rb had augmented H3K4me3 and H3K27ac marks after HI. Genetic and pharmacological inhibition of these receptors resulted in suppressed HSPC proliferation, reduced leukocyte production, and ameliorated atherosclerosis. CONCLUSIONS: Our findings demonstrate increased inflammation, HSPC abundance in the vascular niches of the bone marrow, and elevated IL-3Rb and IL-1R1 (IL-1 receptor 1) expression in HSPC following HI. Furthermore, the IL-3Rb and IL-1R1 signaling plays a pivotal role in HSPC proliferation, leukocyte abundance, and atherosclerosis aggravation after HI.

Identifying Molecular Signatures of Distinct Modes of Collective Migration in Response to the Microenvironment Using Three-Dimensional Breast Cancer Models.

Collective cell migration is a key feature of transition of ductal carcinoma in situ (DCIS) to invasive ductal carcinoma (IDC) among many other cancers, yet the microenvironmental factors and underlying mechanisms that trigger collective migration remain poorly understood. Here, we investigated two microenvironmental factors, tumor-intrinsic hypoxia and tumor-secreted factors (secretome), as triggers of collective migration using three-dimensional (3D) discrete-sized microtumor models that recapitulate hallmarks of DCIS-IDC transition. Interestingly, the two factors induced two distinct modes of collective migration: directional and radial migration in the 3D microtumors generated from the same breast cancer cell line model, T47D. Without external stimulus, large (600 µm) T47D microtumors exhibited tumor-intrinsic hypoxia and directional migration, while small (150 µm), non-hypoxic microtumors exhibited radial migration only when exposed to the secretome of large microtumors. To investigate the mechanisms underlying hypoxia- and secretome-induced directional vs. radial migration modes, we performed differential gene expression analysis of hypoxia- and secretome-induced migratory microtumors compared with non-hypoxic, non-migratory small microtumors as controls. We propose unique gene signature sets related to tumor-intrinsic hypoxia, hypoxia-induced epithelial-mesenchymal transition (EMT), as well as hypoxia-induced directional migration and secretome-induced radial migration. Gene Set Enrichment Analysis (GSEA) and protein-protein interaction (PPI) network analysis revealed enrichment and potential interaction between hypoxia, EMT, and migration gene signatures for the hypoxia-induced directional migration. In contrast, hypoxia and EMT were not enriched in the secretome-induced radial migration, suggesting that complete EMT may not be required for radial migration. Survival analysis identified unique genes associated with low survival rate and poor prognosis in TCGA-breast invasive carcinoma dataset from our tumor-intrinsic hypoxia gene signature (CXCR4, FOXO3, LDH, NDRG1), hypoxia-induced EMT gene signature (EFEMP2, MGP), and directional migration gene signature (MAP3K3, PI3K3R3). NOS3 was common between hypoxia and migration gene signature. Survival analysis from secretome-induced radial migration identified ATM, KCNMA1 (hypoxia gene signature), and KLF4, IFITM1, EFNA1, TGFBR1 (migration gene signature) to be associated with poor survival rate. In conclusion, our unique 3D cultures with controlled microenvironments respond to different microenvironmental factors, tumor-intrinsic hypoxia, and secretome by adopting distinct collective migration modes and their gene expression analysis highlights the phenotypic heterogeneity and plasticity of epithelial cancer cells.

Activation of endothelial Wnt/ß-catenin signaling by protective astrocytes repairs BBB damage in ischemic stroke.

The role of astrocytes in dysregulation of blood-brain barrier (BBB) function following ischemic stroke is not well understood. Here, we investigate the effects of restoring the repair properties of astrocytes on the BBB after ischemic stroke. Mice deficient for NHE1, a pH-sensitive Na+/H+ exchanger 1, in astrocytes have reduced BBB permeability after ischemic stroke, increased angiogenesis and cerebral blood flow perfusion, in contrast to wild-type mice. Bulk RNA-sequencing transcriptome analysis of purified astrocytes revealed that ∼177 genes were differentially upregulated in mutant astrocytes, with Wnt7a mRNA among the top genes. Using a Wnt reporter line, we confirmed that the pathway was upregulated in cerebral vessels of mutant mice after ischemic stroke. However, administration of the Wnt/ß-catenin inhibitor, XAV-939, blocked the reparative effects of Nhe1-deficient astrocytes. Thus, astrocytes lacking pH-sensitive NHE1 protein are transformed from injurious to "protective" by inducing Wnt production to promote BBB repair after ischemic stroke.

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.

ELL2 Is Required for the Growth and Survival of AR-Negative Prostate Cancer Cells.

BACKGROUND: Elongation factor for RNA polymerase II 2 (ELL2) was reported as a putative tumor suppressor in the prostate. ELL2 is frequently down-regulated in prostatic adenocarcinoma specimens, and loss of ELL2 induced murine prostatic intraepithelial neoplasia and enhanced AR-positive prostate cancer cell proliferation. However, the ELL2 gene appears to be amplified in AR-negative neuroendocrine prostate tumors, suggesting a potential oncogenic role for ELL2 in AR-negative prostate cancer cells. In this study, we explored the potential function of ELL2 in PC-3 and DU145, two AR-negative prostate cancer cell lines. MATERIALS AND METHODS: The role of ELL2 in PC-3 and DU145 cells was studied using siRNA-mediated ELL2 knockdown. Genes regulated by ELL2 knockdown in PC-3 cells were identified and analyzed using RNA-Seq and bioinformatics. The expression of representative genes was confirmed by Western blot and/or quantitative PCR. Cell growth was determined by BrdU, MTT and colony formation assays. Cell death was analyzed by 7-AAD/Annexin V staining and trypan blue exclusion staining. Cell cycle was determined by PI staining and flow cytometry. RESULTS: ELL2 knockdown inhibited the proliferation of PC-3 and DU145 cells. RNA-Seq analysis showed an enrichment in genes associated with cell death and survival following ELL2 knockdown. The interferon-γ pathway was identified as the top canonical pathway comprising of 55.6% of the genes regulated by ELL2. ELL2 knockdown induced an increase in STAT1 and IRF1 mRNA and an induction of total STAT1 and phosphorylated STAT1 protein. Inhibition of cell proliferation by ELL2 knockdown was partly abrogated by STAT1 knockdown. ELL2 knockdown inhibited colony formation and induced apoptosis in both PC-3 and DU145 cells. Furthermore, knockdown of ELL2 caused S-phase cell cycle arrest, inhibition of CDK2 phosphorylation and cyclin D1 expression, and increased expression of cyclin E. CONCLUSION: ELL2 knockdown in PC-3 and DU145 cells induced S-phase cell cycle arrest and profound apoptosis, which was accompanied by the induction of genes associated with cell death and survival pathways. These observations suggest that ELL2 is a potential oncogenic protein required for survival and proliferation in AR-negative prostate cancer cells.

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.

Peptide vaccine immunotherapy biomarkers and response patterns in pediatric gliomas.

Low-grade gliomas (LGGs) are the most common brain tumor affecting children. We recently reported an early phase clinical trial of a peptide-based vaccine, which elicited consistent antigen-specific T cell responses in pediatric LGG patients. Additionally, we observed radiologic responses of stable disease (SD), partial response (PR), and near-complete/complete response (CR) following therapy. To identify biomarkers of clinical response in peripheral blood, we performed RNA sequencing on PBMC samples collected at multiple time points. Patients who showed CR demonstrated elevated levels of T cell activation markers, accompanied by a cytotoxic T cell response shortly after treatment initiation. At week 34, patients with CR demonstrated both IFN signaling and Poly-IC:LC adjuvant response patterns. Patients with PR demonstrated a unique, late monocyte response signature. Interestingly, HLA-V expression, before or during therapy, and an early monocytic hematopoietic response were strongly associated with SD. Finally, low IDO1 and PD-L1 expression before treatment and early elevated levels of T cell activation markers were associated with prolonged progression-free survival. Overall, our data support the presence of unique peripheral immune patterns in LGG patients associated with different radiographic responses to our peptide vaccine immunotherapy. Future clinical trials, including our ongoing phase II LGG vaccine immunotherapy, should monitor these response patterns.

EAF2 and p53 Co-Regulate STAT3 Activation in Prostate Cancer.

The tumor suppressor genes EAF2 and p53 are frequently dysregulated in prostate cancers. Recently, we reported that concurrent p53 nuclear staining and EAF2 downregulation were associated with high Gleason score. Combined loss of EAF2 and p53 in a murine model induced prostate tumors, and concurrent knockdown of EAF2 and p53 in prostate cancer cells enhanced proliferation and migration, further suggesting that EAF2 and p53 could functionally interact in the suppression of prostate tumorigenesis. Here, RNA-seq analyses identified differentially regulated genes in response to concurrent knockdown of p53 and EAF2. Several of these genes were associated with the STAT3 signaling pathway, and this was verified by significantly increased p-STAT3 immunostaining in the Eaf2-/-p53-/- mouse prostate. STAT3 knockdown abrogated the stimulation of C4-2 cell proliferation by concurrent knockdown of EAF2 and p53. Furthermore, immunostaining of p-STAT3 was increased in human prostate cancer specimens with EAF2 downregulation and/or p53 nuclear staining. Our findings suggest that simultaneous inactivation of EAF2 and p53 can act to activate STAT3 and drive prostate tumorigenesis.

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.

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.