Integrated epigenomic profiling reveals endogenous retrovirus reactivation in renal cell carcinoma.

BACKGROUND: Transcriptional dysregulation drives cancer formation but the underlying mechanisms are still poorly understood. Renal cell carcinoma (RCC) is the most common malignant kidney tumor which canonically activates the hypoxia-inducible transcription factor (HIF) pathway. Despite intensive study, novel therapeutic strategies to target RCC have been difficult to develop. Since the RCC epigenome is relatively understudied, we sought to elucidate key mechanisms underpinning the tumor phenotype and its clinical behavior. METHODS: We performed genome-wide chromatin accessibility (DNase-seq) and transcriptome profiling (RNA-seq) on paired tumor/normal samples from 3 patients undergoing nephrectomy for removal of RCC. We incorporated publicly available data on HIF binding (ChIP-seq) in a RCC cell line. We performed integrated analyses of these high-resolution, genome-scale datasets together with larger transcriptomic data available through The Cancer Genome Atlas (TCGA). FINDINGS: Though HIF transcription factors play a cardinal role in RCC oncogenesis, we found that numerous transcription factors with a RCC-selective expression pattern also demonstrated evidence of HIF binding near their gene body. Examination of chromatin accessibility profiles revealed that some of these transcription factors influenced the tumor's regulatory landscape, notably the stem cell transcription factor POU5F1 (OCT4). Elevated POU5F1 transcript levels were correlated with advanced tumor stage and poorer overall survival in RCC patients. Unexpectedly, we discovered a HIF-pathway-responsive promoter embedded within a endogenous retroviral long terminal repeat (LTR) element at the transcriptional start site of the PSOR1C3 long non-coding RNA gene upstream of POU5F1. RNA transcripts are induced from this promoter and read through PSOR1C3 into POU5F1 producing a novel POU5F1 transcript isoform. Rather than being unique to the POU5F1 locus, we found that HIF binds to several other transcriptionally active LTR elements genome-wide correlating with broad gene expression changes in RCC. INTERPRETATION: Integrated transcriptomic and epigenomic analysis of matched tumor and normal tissues from even a small number of primary patient samples revealed remarkably convergent shared regulatory landscapes. Several transcription factors appear to act downstream of HIF including the potent stem cell transcription factor POU5F1. Dysregulated expression of POU5F1 is part of a larger pattern of gene expression changes in RCC that may be induced by HIF-dependent reactivation of dormant promoters embedded within endogenous retroviral LTRs.

Integrative detection and analysis of structural variation in cancer genomes.

Structural variants (SVs) can contribute to oncogenesis through a variety of mechanisms. Despite their importance, the identification of SVs in cancer genomes remains challenging. Here, we present a framework that integrates optical mapping, high-throughput chromosome conformation capture (Hi-C), and whole-genome sequencing to systematically detect SVs in a variety of normal or cancer samples and cell lines. We identify the unique strengths of each method and demonstrate that only integrative approaches can comprehensively identify SVs in the genome. By combining Hi-C and optical mapping, we resolve complex SVs and phase multiple SV events to a single haplotype. Furthermore, we observe widespread structural variation events affecting the functions of noncoding sequences, including the deletion of distal regulatory sequences, alteration of DNA replication timing, and the creation of novel three-dimensional chromatin structural domains. Our results indicate that noncoding SVs may be underappreciated mutational drivers in cancer genomes.

Polymyxin resistance of Pseudomonas aeruginosa phoQ mutants is dependent on additional two-component regulatory systems.

Pseudomonas aeruginosa can develop resistance to polymyxin as a consequence of mutations in the PhoPQ regulatory system, mediated by covalent lipid A modification. Transposon mutagenesis of a polymyxin-resistant phoQ mutant defined 41 novel loci required for resistance, including two regulatory systems, ColRS and CprRS. Deletion of the colRS genes, individually or in tandem, abrogated the polymyxin resistance of a ΔphoQ mutant, as did individual or tandem deletion of cprRS. Individual deletion of colR or colS in a ΔphoQ mutant also suppressed 4-amino-L-arabinose addition to lipid A, consistent with the known role of this modification in polymyxin resistance. Surprisingly, tandem deletion of colRS or cprRS in the ΔphoQ mutant or individual deletion of cprR or cprS failed to suppress 4-amino-L-arabinose addition to lipid A, indicating that this modification alone is not sufficient for PhoPQ-mediated polymyxin resistance in P. aeruginosa. Episomal expression of colRS or cprRS in tandem or of cprR individually complemented the Pm resistance phenotype in the ΔphoQ mutant, while episomal expression of colR, colS, or cprS individually did not. Highly polymyxin-resistant phoQ mutants of P. aeruginosa isolated from polymyxin-treated cystic fibrosis patients harbored mutant alleles of colRS and cprS; when expressed in a ΔphoQ background, these mutant alleles enhanced polymyxin resistance. These results define ColRS and CprRS as two-component systems regulating polymyxin resistance in P. aeruginosa, indicate that addition of 4-amino-L-arabinose to lipid A is not the only PhoPQ-regulated biochemical mechanism required for resistance, and demonstrate that colRS and cprS mutations can contribute to high-level clinical resistance.

Systematic localization of common disease-associated variation in regulatory DNA.

Genome-wide association studies have identified many noncoding variants associated with common diseases and traits. We show that these variants are concentrated in regulatory DNA marked by deoxyribonuclease I (DNase I) hypersensitive sites (DHSs). Eighty-eight percent of such DHSs are active during fetal development and are enriched in variants associated with gestational exposure-related phenotypes. We identified distant gene targets for hundreds of variant-containing DHSs that may explain phenotype associations. Disease-associated variants systematically perturb transcription factor recognition sequences, frequently alter allelic chromatin states, and form regulatory networks. We also demonstrated tissue-selective enrichment of more weakly disease-associated variants within DHSs and the de novo identification of pathogenic cell types for Crohn's disease, multiple sclerosis, and an electrocardiogram trait, without prior knowledge of physiological mechanisms. Our results suggest pervasive involvement of regulatory DNA variation in common human disease and provide pathogenic insights into diverse disorders.

Large-insert genome analysis technology detects structural variation in Pseudomonas aeruginosa clinical strains from cystic fibrosis patients.

Large-insert genome analysis (LIGAN) is a broadly applicable, high-throughput technology designed to characterize genome-scale structural variation. Fosmid paired-end sequences and DNA fingerprints from a query genome are compared to a reference sequence using the Genomic Variation Analysis (GenVal) suite of software tools to pinpoint locations of insertions, deletions, and rearrangements. Fosmids spanning regions that contain new structural variants can then be sequenced. Clonal pairs of Pseudomonas aeruginosa isolates from four cystic fibrosis patients were used to validate the LIGAN technology. Approximately 1.5 Mb of inserted sequences were identified, including 743 kb containing 615 ORFs that are absent from published P. aeruginosa genomes. Six rearrangement breakpoints and 220 kb of deleted sequences were also identified. Our study expands the "genome universe" of P. aeruginosa and validates a technology that complements emerging, short-read sequencing methods that are better suited to characterizing single-nucleotide polymorphisms than structural variation.

Multiarm high-throughput integration site detection: limitations of LAM-PCR technology and optimization for clonal analysis.

Retroviral integration provides a unique and heritable genomic tag for a target cell and its progeny, enabling studies of clonal composition and repopulation kinetics after gene transfer into hematopoietic stem cells. The clonal tracking method, linear amplification-mediated polymerase chain reaction (LAM-PCR) is widely employed to follow the hematopoietic output of retrovirally marked stem cells. Here we examine the capabilities and limitations of conventional LAM-PCR to track individual clones in a complex multiclonal mix. Using artificial mixtures of retrovirally marked, single-cell-derived clones, we demonstrate that LAM-PCR fails to detect 30-40% of the clones, even after exhaustive analysis. Furthermore, the relative abundance of specific clones within a mix is not accurately represented, deviating by as much as 60-fold from their true abundance. We describe an optimized, multiarm, high-throughput modification of LAM-PCR that improves the global detection capacity to greater than 90% with exhaustive sampling, facilitates accurate estimates of the total pool size from smaller samplings, and provides a rapid, cost-effective approach to the generation of large insertion-site data bases required for evaluation of vector integration preferences. The inability to estimate the abundance of individual clones within mixtures remains a serious limitation. Thus, although LAM-PCR is a powerful tool for identification of integration sites and for estimations of clonal complexity, it fails to provide the semiquantitative information necessary for direct, reliable tracking of individual clones in a chimeric background.

Unique integration profiles in a canine model of long-term repopulating cells transduced with gammaretrovirus, lentivirus, or foamy virus.

Recent advances have allowed for improved retrovirus-mediated gene transfer, and therapeutic benefits have been described in patients. These successes have shown the potential of hematopoietic stem cell (HSC) gene therapy, but treatment-related leukemia and benign expansion of gene-modified clones have shifted the attention toward safety. The delayed onset of adverse events in gene therapy clinical trials emphasizes the importance of long-term integration site studies in large animal models. We have addressed safety by characterizing the genomic location of 555 integration sites of the three most commonly used integrating retroviral vectors, that is, gammaretrovirus, lentivirus, and foamy virus, in long-term repopulating cells from dogs. Gammaretroviral integrants showed the most significant frequency of occurrence very close (<2.5 kb) to transcription start sites, but a substantial portion of all three retroviral integrants were within 50 kb. Importantly, gammaretroviral integrants were found more frequently in and near proto-oncogenes, suggesting this retroviral system may be the most prone to adverse gene activation. These data suggest that gammaretroviral vectors may have the highest intrinsic risk, but also emphasize that no vector system can be defined as "safe" based solely on integration profile.

Comparison of HIV-derived lentiviral and MLV-based gammaretroviral vector integration sites in primate repopulating cells.

The potential for leukemia caused by retroviral vector integration has become a significant concern for hematopoietic stem cell gene therapy. We analyzed the distribution of vector integrants in pigtailed macaque and baboon repopulating cells for the two most commonly used retroviral vector systems, human immunodeficiency virus (HIV)-based lentiviral vectors and murine leukemia virus (MLV)-based gammaretroviral vectors, to help define their relative genotoxicity. All animals had polyclonal engraftment with no apparent adverse effects from transplantation with gene-modified cells. In all, 380 MLV and 235 HIV unique vector integration sites were analyzed and had distinct distribution patterns in relation to genes and CpG islands as observed in previous in vitro studies. Both vector types were found more frequently in and near proto-oncogenes in repopulating cells than in a random dataset. Analysis of functional classes of genes with integrants within 100 kilobases (kb) of their transcription start sites showed an over-representation of genes involved in growth or survival near both lentiviral and gammaretroviral integrants. Microarray analysis showed that both gammaretroviral and lentiviral vectors were found close to genes with high expression levels in primitive cells enriched for hematopoietic stem cells. These data help define the relative risk of insertional mutagenesis with MLV-, HIV-, and simian immunodeficiency virus (SIV)-based vectors in a highly relevant primate model.

Type IV pili-mediated secretion modulates Francisella virulence.

Francisella tularensis are the causative agent of the zoonotic disease, tularaemia. Among four F. tularensis subspecies, ssp. novicida (F. novicida) is pathogenic only for immunocompromised individuals, while all four subspecies are pathogenic for mice. This study utilized proteomic and bioinformatic approaches to identify seven F. novicida secreted proteins and the corresponding Type IV pilus (T4P) secretion system. The secreted proteins were predicted to encode two chitinases, a chitin binding protein, a protease (PepO), and a beta-glucosidase (BglX). The transcription of F. novicida pepO and bglX was regulated by the virulence regulator MglA. Intradermal infection of mice with F. novicida mutants defective in T4P secretion system or PepO resulted in enhanced F. novicida spread to systemic sites. Infection with F. novicida pepO mutants also resulted in increased neutrophil infiltration into the mouse airways. PepO is a zinc protease that is homologous to mammalian endothelin-converting enzyme ECE-1. Therefore, secretion of PepO likely results in increased production of endothelin and increased vasoconstriction at the infection site in skin that limits the F. novicida spread. Francisella human pathogenic strains contain a mutation in pepO predicted to abolish its secretion. Loss of PepO function may have contributed to evolution of highly virulent Francisellae.

Genetic adaptation by Pseudomonas aeruginosa to the airways of cystic fibrosis patients.

In many human infections, hosts and pathogens coexist for years or decades. Important examples include HIV, herpes viruses, tuberculosis, leprosy, and malaria. With the exception of intensively studied viral infections such as HIV/AIDs, little is known about the extent to which the clonal expansion that occurs during long-term infection by pathogens involves important genetic adaptations. We report here a detailed, whole-genome analysis of one such infection, that of a cystic fibrosis (CF) patient by the opportunistic bacterial pathogen Pseudomonas aeruginosa. The bacteria underwent numerous genetic adaptations during 8 years of infection, as evidenced by a positive-selection signal across the genome and an overwhelming signal in specific genes, several of which are mutated during the course of most CF infections. Of particular interest is our finding that virulence factors that are required for the initiation of acute infections are often selected against during chronic infections. It is apparent that the genotypes of the P. aeruginosa strains present in advanced CF infections differ systematically from those of "wild-type" P. aeruginosa and that these differences may offer new opportunities for treatment of this chronic disease.

The DNA sequence, annotation and analysis of human chromosome 3.

After the completion of a draft human genome sequence, the International Human Genome Sequencing Consortium has proceeded to finish and annotate each of the 24 chromosomes comprising the human genome. Here we describe the sequencing and analysis of human chromosome 3, one of the largest human chromosomes. Chromosome 3 comprises just four contigs, one of which currently represents the longest unbroken stretch of finished DNA sequence known so far. The chromosome is remarkable in having the lowest rate of segmental duplication in the genome. It also includes a chemokine receptor gene cluster as well as numerous loci involved in multiple human cancers such as the gene encoding FHIT, which contains the most common constitutive fragile site in the genome, FRA3B. Using genomic sequence from chimpanzee and rhesus macaque, we were able to characterize the breakpoints defining a large pericentric inversion that occurred some time after the split of Homininae from Ponginae, and propose an evolutionary history of the inversion.

Foamy virus vector integration sites in normal human cells.

Foamy viruses (FVs) or spumaviruses are retroviruses that have been developed as vectors, but their integration patterns have not been described. We have performed a large-scale analysis of FV integration sites in unselected human fibroblasts (n = 1,008) and human CD34(+) hematopoietic cells (n = 1,821) by using a bacterial shuttle vector and a comparable analysis of lentiviral vector integration sites in CD34(+) cells (n = 1,331). FV vectors had a distinct integration profile relative to other types of retroviruses. They did not integrate preferentially within genes, despite a modest preference for integration near transcription start sites and a significant preference for CpG islands. The genomewide distribution of FV vector proviruses was nonrandom, with both clusters and gaps. Transcriptional profiling showed that gene expression had little influence on integration site selection. Our findings suggest that FV vectors may have desirable integration properties for gene therapy applications.

Large-scale analysis of adeno-associated virus vector integration sites in normal human cells.

The integration sites of viral vectors used in human gene therapy can have important consequences for safety and efficacy. However, an extensive evaluation of adeno-associated virus (AAV) vector integration sites has not been completed, despite the ongoing use of AAV vectors in clinical trials. Here we have used a shuttle vector system to isolate and analyze 977 unique AAV vector-chromosome integration junctions from normal human fibroblasts and describe their genomic distribution. We found a significant preference for integrating within CpG islands and the first 1 kb of genes, but only a slight overall preference for transcribed sequences. Integration sites were clustered throughout the genome, including a major preference for integration in ribosomal DNA repeats, and 13 other hotspots that contained three or more proviruses within a 500-kb window. Both junctions were localized from 323 proviruses, allowing us to characterize the chromosomal deletions, insertions, and translocations associated with vector integration. These studies establish a profile of insertional mutagenesis for AAV vectors and provide unique insight into the chromosomal distribution of DNA strand breaks that may facilitate integration.

Fine-scale structural variation of the human genome.

Inversions, deletions and insertions are important mediators of disease and disease susceptibility. We systematically compared the human genome reference sequence with a second genome (represented by fosmid paired-end sequences) to detect intermediate-sized structural variants >8 kb in length. We identified 297 sites of structural variation: 139 insertions, 102 deletions and 56 inversion breakpoints. Using combined literature, sequence and experimental analyses, we validated 112 of the structural variants, including several that are of biomedical relevance. These data provide a fine-scale structural variation map of the human genome and the requisite sequence precision for subsequent genetic studies of human disease.

Whole-genome sequence variation among multiple isolates of Pseudomonas aeruginosa.

Whole-genome shotgun sequencing was used to study the sequence variation of three Pseudomonas aeruginosa isolates, two from clonal infections of cystic fibrosis patients and one from an aquatic environment, relative to the genomic sequence of reference strain PAO1. The majority of the PAO1 genome is represented in these strains; however, at least three prominent islands of PAO1-specific sequence are apparent. Conversely, approximately 10% of the sequencing reads derived from each isolate fail to align with the PAO1 backbone. While average sequence variation among all strains is roughly 0.5%, regions of pronounced differences were evident in whole-genome scans of nucleotide diversity. We analyzed two such divergent loci, the pyoverdine and O-antigen biosynthesis regions, by complete resequencing. A thorough analysis of isolates collected over time from one of the cystic fibrosis patients revealed independent mutations resulting in the loss of O-antigen synthesis alternating with a mucoid phenotype. Overall, we conclude that most of the PAO1 genome represents a core P. aeruginosa backbone sequence while the strains addressed in this study possess additional genetic material that accounts for at least 10% of their genomes. Approximately half of these additional sequences are novel.