High-throughput sequencing reveals principles of adeno-associated virus serotype 2 integration.

Viral integrations are important in human biology, yet genome-wide integration profiles have not been determined for many viruses. Adeno-associated virus (AAV) infects most of the human population and is a prevalent gene therapy vector. AAV integrates into the human genome with preference for a single locus, termed AAVS1. However, the genome-wide integration of AAV has not been defined, and the principles underlying this recombination remain unclear. Using a novel high-throughput approach, integrant capture sequencing, nearly 12 million AAV junctions were recovered from a human cell line, providing five orders of magnitude more data than were previously available. Forty-five percent of integrations occurred near AAVS1, and several thousand novel integration hotspots were identified computationally. Most of these occurred in genes, with dozens of hotspots targeting known oncogenes. Viral replication protein binding sites (RBS) and transcriptional activity were major factors favoring integration. In a first for eukaryotic viruses, the data reveal a unique asymmetric integration profile with distinctive directional orientation of viral genomes. These studies provide a new understanding of AAV integration biology through the use of unbiased high-throughput data acquisition and bioinformatics.

Comprehensive immunophenotyping of solid tumor-infiltrating immune cells reveals the expression characteristics of LAG-3 and its ligands.

Background: Immune cell expression profiling from patient samples is critical for the successful development of immuno-oncology agents and is useful to understand mechanism-of-action, to identify exploratory biomarkers predictive of response, and to guide treatment selection and combination therapy strategies. LAG-3 is an inhibitory immune checkpoint that can suppress antitumor T-cell responses and targeting LAG-3, in combination with PD-1, is a rational approach to enhance antitumor immunity that has recently demonstrated clinical success. Here, we sought to identify human immune cell subsets that express LAG-3 and its ligands, to characterize the marker expression profile of these subsets, and to investigate the potential relationship between LAG-3 expressing subsets and clinical outcomes to immuno-oncology therapies. Methods: Comprehensive high-parameter immunophenotyping was performed using mass and flow cytometry of tumor-infiltrating lymphocytes (TILs) and peripheral blood mononuclear cells (PBMCs) from two independent cohorts of samples from patients with various solid tumor types. Profiling of circulating immune cells by single cell RNA-seq was conducted on samples from a clinical trial cohort of melanoma patients treated with immunotherapy. Results: LAG-3 was most highly expressed by subsets of tumor-infiltrating CD8 T central memory (TCM) and effector memory (TEM) cells and was frequently co-expressed with PD-1. We determined that these PD-1+ LAG-3+ CD8 memory T cells exhibited a unique marker profile, with greater expression of activation (CD69, HLA-DR), inhibitory (TIM-3, TIGIT, CTLA-4) and stimulatory (4-1BB, ICOS) markers compared to cells that expressed only PD-1 or LAG-3, or that were negative for both checkpoints. In contrast to tumors, LAG-3 expression was more limited in circulating immune cells from healthy donors and solid tumor patients. Additionally, we found abundant expression of the LAG-3 ligands MHC-II and galectin-3 in diverse immune cell types, whereas FGL1 and LSECtin were minimally expressed by immune cells in the tumor microenvironment (TME). Lastly, we found an inverse relationship between baseline and on-treatment levels of circulating LAG3 transcript-expressing CD8 memory T cells and response to combination PD-1 and CTLA-4 blockade in a clinical trial cohort of melanoma patients profiled by scRNAseq. Conclusions: These results provide insights into the nature of LAG-3- and ligand-expressing immune cells within the TME, and suggest a biological basis for informing mechanistic hypotheses, treatment selection strategies, and combination immunotherapy approaches to support continued development of dual PD-1 and LAG-3 blockade.

Changes in DNA methylation accompany changes in gene expression during chondrocyte hypertrophic differentiation in vitro.

During osteoarthritis (OA), articular chondrocytes undergo phenotypic changes that resemble developmental patterns characteristic of growth plate chondrocytes. These phenotypic alterations lead to a hypertrophy-like phenotype characterized by altered production of extracellular matrix constituents and increased collagenase activity, which, in turn, results in cartilage destruction in OA disease. Recent studies have shown that the phenotypic instability and dysregulated gene expression in OA are associated with changes in DNA methylation patterns. Subsequent efforts have aimed to identify changes in DNA methylation with functional impact in OA disease, to potentially uncover therapeutic targets. Here, we paired an in vitro 3D/pellet culture system that mimics chondrocyte hypertrophy with RNA sequencing (RNA-Seq) and enhanced reduced representation of bisulfite sequencing (ERRBS) to identify transcriptomic and epigenomic changes in murine primary articular chondrocytes undergoing hypertrophy-like differentiation. We identified hypertrophy-associated changes in DNA methylation patterns in vitro. Integration of RNA-Seq and ERRBS datasets identified associations between changes in methylation and gene expression. Our integrative analyses showed that hypertrophic differentiation of articular chondrocytes is accompanied by transcriptomic and epigenomic changes in vitro. We believe that our integrative approaches have the potential to uncover new targets for therapeutic intervention.

Transcriptomic and epigenomic analyses uncovered Lrrc15 as a contributing factor to cartilage damage in osteoarthritis.

In osteoarthritis (OA), articular chondrocytes display phenotypic and functional changes associated with epigenomic alterations. These changes contribute to the disease progression, which is characterized by dysregulated reparative processes and abnormal extracellular matrix remodeling leading to cartilage degradation. Recent studies using a murine model of posttraumatic OA highlighted the contribution of changes in DNA hydroxymethylation (5hmC) to OA progression. Here, we integrated transcriptomic and epigenomic analyses in cartilage after induction of OA to show that the structural progression of OA is accompanied by early transcriptomic and pronounced DNA methylation (5mC) changes in chondrocytes. These changes accumulate over time and are associated with recapitulation of developmental processes, including cartilage development, chondrocyte hypertrophy, and ossification. Our integrative analyses also uncovered that Lrrc15 is differentially methylated and expressed in OA cartilage, and that it may contribute to the functional and phenotypic alterations of chondrocytes, likely coordinating stress responses and dysregulated extracellular matrix remodeling.

The long noncoding RNA H19 regulates tumor plasticity in neuroendocrine prostate cancer.

Neuroendocrine (NE) prostate cancer (NEPC) is a lethal subtype of castration-resistant prostate cancer (PCa) arising either de novo or from transdifferentiated prostate adenocarcinoma following androgen deprivation therapy (ADT). Extensive computational analysis has identified a high degree of association between the long noncoding RNA (lncRNA) H19 and NEPC, with the longest isoform highly expressed in NEPC. H19 regulates PCa lineage plasticity by driving a bidirectional cell identity of NE phenotype (H19 overexpression) or luminal phenotype (H19 knockdown). It contributes to treatment resistance, with the knockdown of H19 re-sensitizing PCa to ADT. It is also essential for the proliferation and invasion of NEPC. H19 levels are negatively regulated by androgen signaling via androgen receptor (AR). When androgen is absent SOX2 levels increase, driving H19 transcription and facilitating transdifferentiation. H19 facilitates the PRC2 complex in regulating methylation changes at H3K27me3/H3K4me3 histone sites of AR-driven and NEPC-related genes. Additionally, this lncRNA induces alterations in genome-wide DNA methylation on CpG sites, further regulating genes associated with the NEPC phenotype. Our clinical data identify H19 as a candidate diagnostic marker and predictive marker of NEPC with elevated H19 levels associated with an increased probability of biochemical recurrence and metastatic disease in patients receiving ADT. Here we report H19 as an early upstream regulator of cell fate, plasticity, and treatment resistance in NEPC that can reverse/transform cells to a treatable form of PCa once therapeutically deactivated.

Complement anaphylatoxin C5a neuroprotects through regulation of glutamate receptor subunit 2 in vitro and in vivo.

BACKGROUND: The complement system is thought to be involved in the pathogenesis of numerous neurological diseases. We previously reported that pre-treatment of murine cortico-hippocampal neuronal cultures with the complement derived anaphylatoxin C5a, protects against glutamate mediated apoptosis. Our present study with C5a receptor knock out (C5aRKO) mice corroborates that the deficiency of C5a renders C5aRKO mouse more susceptible to apoptotic injury in vivo. In this study we explored potential upstream mechanisms involved in C5a mediated neuroprotection in vivo and in vitro. METHODS: Based on evidence suggesting that reduced expression of glutamate receptor subunit 2 (GluR2) may influence apoptosis in neurons, we studied the effect of human recombinant C5a on GluR2 expression in response to glutamate neurotoxicity. Glutamate analogs were injected into C5aRKO mice or used to treat in vitro neuronal culture and GluR2 expression were assessed in respect with cell death. RESULTS: In C5aRKO mice we found that the neurons are more susceptible to excitotoxicity resulting in apoptotic injury in the absence of the C5a receptor compared to WT control mice. Our results suggest that C5a protects against apoptotic pathways in neurons in vitro and in vivo through regulation of GluR2 receptor expression. CONCLUSION: Complement C5a neuroprotects through regulation of GluR2 receptor subunit.

Gene expression profiling separates chromophobe renal cell carcinoma from oncocytoma and identifies vesicular transport and cell junction proteins as differentially expressed genes.

PURPOSE: To compare gene expression profiles of chromophobe renal cell carcinoma (RCC) and benign oncocytoma, aiming at identifying differentially expressed genes. EXPERIMENTAL DESIGN: Nine cases each of chromophobe RCC and oncocytoma were analyzed by oligonucleotide microarray. Candidate genes that showed consistent differential expression were validated by reverse transcription-PCR using 25 fresh-frozen and 15 formalin-fixed, paraffin-embedded tumor samples. Immunohistochemical analysis was also done for two selected gene products, claudin 8 and MAL2. RESULTS: Unsupervised hierarchical clustering separated the chromophobe RCC and oncocytoma into two distinct groups. By a combination of data analysis approaches, we identified 11 candidate genes showing consistent differential expression between chromophobe RCC and oncocytoma. Five of these genes, AP1M2, MAL2, PROM2, PRSS8, and FLJ20171, were shown to effectively separate these two tumor groups by quantitative reverse transcription-PCR using fresh tissue samples, with similar trends seen on formalin-fixed tissues. Immunohistochemical analysis revealed selective expression of MAL2 and claudin 8 in distal renal tubules, with MAL2 antibody showing differential expression between chromophobe RCC and oncocytoma. Functional analyses suggest that genes encoding tight junction proteins and vesicular membrane trafficking proteins, normally expressed in distal nephrons, are retained in chromophobe RCC and lost or consistently down-regulated in oncocytoma, indicating that these two tumor types, believed to be both derived from distal tubules, are likely distinctive in their histogenesis. CONCLUSIONS: We showed that chromophobe RCC and oncocytoma are distinguishable by mRNA expression profiles and a panel of gene products potentially useful as diagnostic markers were identified.

Lowered H3K27me3 and DNA hypomethylation define poorly prognostic pediatric posterior fossa ependymomas.

Childhood posterior fossa (PF) ependymomas cause substantial morbidity and mortality. These tumors lack recurrent genetic mutations, but a subset of these ependymomas exhibits CpG island (CpGi) hypermethylation [PF group A (PFA)], implicating epigenetic alterations in their pathogenesis. Further, histological grade does not reliably predict prognosis, highlighting the importance of developing more robust prognostic markers. We discovered global H3K27me3 reduction in a subset of these tumors (PF-ve ependymomas) analogous to H3K27M mutant gliomas. PF-ve tumors exhibited many clinical and biological similarities with PFA ependymomas. Genomic H3K27me3 distribution showed an inverse relationship with CpGi methylation, suggesting that CpGi hypermethylation drives low H3K27me3 in PF-ve ependymomas. Despite CpGi hypermethylation and global H3K27me3 reduction, these tumors showed DNA hypomethylation in the rest of the genome and exhibited increased H3K27me3 genomic enrichment at limited genomic loci similar to H3K27M mutant gliomas. Combined integrative analysis of PF-ve ependymomas with H3K27M gliomas uncovered common epigenetic deregulation of select factors that control radial glial biology, and PF radial glia in early human development exhibited reduced H3K27me3. Finally, H3K27me3 immunostaining served as a biomarker of poor prognosis and delineated radiologically invasive tumors, suggesting that reduced H3K27me3 may be a prognostic indicator in PF ependymomas.

Effects of tobacco smoke on gene expression and cellular pathways in a cellular model of oral leukoplakia.

In addition to being causally linked to the formation of multiple tumor types, tobacco use has been associated with decreased efficacy of anticancer treatment and reduced survival time. A detailed understanding of the cellular mechanisms that are affected by tobacco smoke (TS) should facilitate the development of improved preventive and therapeutic strategies. We have investigated the effects of a TS extract on the transcriptome of MSK-Leuk1 cells, a cellular model of oral leukoplakia. Using Affymetrix HGU133 Plus 2 arrays, 411 differentially expressed probe sets were identified. The observed transcriptome changes were grouped according to functional information and translated into molecular interaction network maps and signaling pathways. Pathways related to cellular proliferation, inflammation, apoptosis, and tissue injury seemed to be perturbed. Analysis of networks connecting the affected genes identified specific modulated molecular interactions, hubs, and key transcription regulators. Thus, TS was found to induce several epidermal growth factor receptor (EGFR) ligands forming an EGFR-centered molecular interaction network, as well as several aryl hydrocarbon receptor-dependent genes, including the xenobiotic metabolizing enzymes CYP1A1 and CYP1B1. Notably, the latter findings in vitro are consistent with our parallel finding that CYP1A1 and CYP1B1 levels were increased in oral mucosa of smokers. Collectively, these results offer insights into the mechanisms underlying the procarcinogenic effects of TS and raise the possibility that inhibitors of EGFR or aryl hydrocarbon receptor signaling will prevent or delay the development of TS-related tumors. Moreover, the inductive effects of TS on xenobiotic metabolizing enzymes may help explain the reduced efficacy of chemotherapy, and suggest targets for chemopreventive agents in smokers.

Genomic resources for songbird research and their use in characterizing gene expression during brain development.

Vocal learning and neuronal replacement have been studied extensively in songbirds, but until recently, few molecular and genomic tools for songbird research existed. Here we describe new molecular/genomic resources developed in our laboratory. We made cDNA libraries from zebra finch (Taeniopygia guttata) brains at different developmental stages. A total of 11,000 cDNA clones from these libraries, representing 5,866 unique gene transcripts, were randomly picked and sequenced from the 3' ends. A web-based database was established for clone tracking, sequence analysis, and functional annotations. Our cDNA libraries were not normalized. Sequencing ESTs without normalization produced many developmental stage-specific sequences, yielding insights into patterns of gene expression at different stages of brain development. In particular, the cDNA library made from brains at posthatching day 30-50, corresponding to the period of rapid song system development and song learning, has the most diverse and richest set of genes expressed. We also identified five microRNAs whose sequences are highly conserved between zebra finch and other species. We printed cDNA microarrays and profiled gene expression in the high vocal center of both adult male zebra finches and canaries (Serinus canaria). Genes differentially expressed in the high vocal center were identified from the microarray hybridization results. Selected genes were validated by in situ hybridization. Networks among the regulated genes were also identified. These resources provide songbird biologists with tools for genome annotation, comparative genomics, and microarray gene expression analysis.

PDZBase: a protein-protein interaction database for PDZ-domains.

SUMMARY: PDZBase is a database that aims to contain all known PDZ-domain-mediated protein-protein interactions. Currently, PDZBase contains approximately 300 such interactions, which have been manually extracted from > 200 articles. The database can be queried through both sequence motif and keyword-based searches, and the sequences of interacting proteins can be visually inspected through alignments (for the comparison of several interactions), or as residue-based diagrams including schematic secondary structure information (for individual complexes).