In Vitro Insertional Mutagenesis Screen Identifies Novel Genes Driving Breast Cancer Metastasis.

Metastasis, a complex, multistep process, is responsible for the overwhelming majority of cancer-related deaths. Despite its devastating consequences, it is not possible to effectively treat cancer that has spread to vital organs, the mechanisms leading to metastasis are still poorly understood, and the catalog of metastasis promoting genes is still incomprehensive. To identify new driver genes of metastasis development, we performed an in vitro Sleeping Beauty transposon-based forward genetic screen in nonmetastatic SKBR3 human breast cancer cells. Boyden chamber-based matrix invasion assays were used to harvest cells that acquired a de novo invasive phenotype. Using targeted RNA sequencing data from 18 pools of invasive cells, we carried out a gene-centric candidate gene prediction and identified established and novel metastasis driver genes. Analysis of these genes revealed their association with metastasis related processes and we further established their clinical relevance in metastatic breast cancer. Two novel candidate genes, G protein-coupled receptor kinase interacting ArfGAP 2 (GIT2) and muscle-associated receptor tyrosine kinase (MUSK), were functionally validated as metastasis driver genes in a series of in vitro and in vivo experimental metastasis models. We propose that our robust and scalable approach will be a useful addition to the toolkit of methodologic resources used to identify genes driving cancer metastasis. IMPLICATIONS: Novel metastasis drivers were identified in a human breast cancer cell line by performing an in vitro, Sleeping Beauty transposon-based forward genetic screen and an RNA fusion-based candidate gene prediction.

Ecological plant epigenetics: Evidence from model and non-model species, and the way forward.

Growing evidence shows that epigenetic mechanisms contribute to complex traits, with implications across many fields of biology. In plant ecology, recent studies have attempted to merge ecological experiments with epigenetic analyses to elucidate the contribution of epigenetics to plant phenotypes, stress responses, adaptation to habitat, and range distributions. While there has been some progress in revealing the role of epigenetics in ecological processes, studies with non-model species have so far been limited to describing broad patterns based on anonymous markers of DNA methylation. In contrast, studies with model species have benefited from powerful genomic resources, which contribute to a more mechanistic understanding but have limited ecological realism. Understanding the significance of epigenetics for plant ecology requires increased transfer of knowledge and methods from model species research to genomes of evolutionarily divergent species, and examination of responses to complex natural environments at a more mechanistic level. This requires transforming genomics tools specifically for studying non-model species, which is challenging given the large and often polyploid genomes of plants. Collaboration among molecular geneticists, ecologists and bioinformaticians promises to enhance our understanding of the mutual links between genome function and ecological processes.

Comprehensive Small RNA-Seq of Adeno-Associated Virus (AAV)-Infected Human Cells Detects Patterns of Novel, Non-Coding AAV RNAs in the Absence of Cellular miRNA Regulation.

Most DNA viruses express small regulatory RNAs, which interfere with viral or cellular gene expression. For adeno-associated virus (AAV), a small ssDNA virus with a complex biphasic life cycle miRNAs or other small regulatory RNAs have not yet been described. This is the first comprehensive Illumina-based RNA-Seq analysis of small RNAs expressed by AAV alone or upon co-infection with helper adenovirus or HSV. Several hotspots of AAV-specific small RNAs were detected mostly close to or within the AAV-ITR and apparently transcribed from the newly identified anti-p5 promoter. An additional small RNA hotspot was located downstream of the p40 promoter, from where transcription of non-coding RNAs associated with the inhibition of adenovirus replication were recently described. Parallel detection of known Ad and HSV miRNAs indirectly validated the newly identified small AAV RNA species. The predominant small RNAs were analyzed on Northern blots and by human argonaute protein-mediated co-immunoprecipitation. None of the small AAV RNAs showed characteristics of bona fide miRNAs, but characteristics of alternative RNA processing indicative of differentially regulated AAV promoter-associated small RNAs. Furthermore, the AAV-induced regulation of cellular miRNA levels was analyzed at different time points post infection. In contrast to other virus groups AAV infection had virtually no effect on the expression of cellular miRNA, which underscores the long-established concept that wild-type AAV infection is apathogenic.

Genome-wide Profiling Reveals Remarkable Parallels Between Insertion Site Selection Properties of the MLV Retrovirus and the piggyBac Transposon in Primary Human CD4(+) T Cells.

The inherent risks associated with vector insertion in gene therapy need to be carefully assessed. We analyzed the genome-wide distributions of Sleeping Beauty (SB) and piggyBac (PB) transposon insertions as well as MLV retrovirus and HIV lentivirus insertions in human CD4(+) T cells with respect to a panel of 40 chromatin states. The distribution of SB transposon insertions displayed the least deviation from random, while the PB transposon and the MLV retrovirus showed unexpected parallels across all chromatin states. Both MLV and PB insertions are enriched at transcriptional start sites (TSSs) and co-localize with BRD4-associated sites. We demonstrate physical interaction between the PB transposase and bromodomain and extraterminal domain proteins (including BRD4), suggesting convergent evolution of a tethering mechanism that directs integrating genetic elements into TSSs. We detect unequal biases across the four systems with respect to targeting genes whose deregulation has been previously linked to serious adverse events in gene therapy clinical trials. The SB transposon has the highest theoretical chance of targeting a safe harbor locus in the human genome. The data underscore the significance of vector choice to reduce the mutagenic load on cells in clinical applications.

A Comprehensive RNA Sequencing Analysis of the Adeno-Associated Virus (AAV) Type 2 Transcriptome Reveals Novel AAV Transcripts, Splice Variants, and Derived Proteins.

UNLABELLED: Adeno-associated virus (AAV) is recognized for its bipartite life cycle with productive replication dependent on coinfection with adenovirus (Ad) and AAV latency being established in the absence of a helper virus. The shift from latent to Ad-dependent AAV replication is mostly regulated at the transcriptional level. The current AAV transcription map displays highly expressed transcripts as found upon coinfection with Ad. So far, AAV transcripts have only been characterized on the plus strand of the AAV single-stranded DNA genome. The AAV minus strand is assumed not to be transcribed. Here, we apply Illumina-based RNA sequencing (RNA-Seq) to characterize the entire AAV2 transcriptome in the absence or presence of Ad. We find known and identify novel AAV transcripts, including additional splice variants, the most abundant of which leads to expression of a novel 18-kDa Rep/VP fusion protein. Furthermore, we identify for the first time transcription on the AAV minus strand with clustered reads upstream of the p5 promoter, confirmed by 5' rapid amplification of cDNA ends and RNase protection assays. The p5 promoter displays considerable activity in both directions, a finding indicative of divergent transcription. Upon infection with AAV alone, low-level transcription of both AAV strands is detectable and is strongly stimulated upon coinfection with Ad. IMPORTANCE: Next-generation sequencing (NGS) allows unbiased genome-wide analyses of transcription profiles, used here for an in depth analysis of the AAV2 transcriptome during latency and productive infection. RNA-Seq analysis led to the discovery of novel AAV transcripts and splice variants, including a derived, novel 18-kDa Rep/VP fusion protein. Unexpectedly, transcription from the AAV minus strand was discovered, indicative of divergent transcription from the p5 promoter. This finding opens the door for novel concepts of the switch between AAV latency and productive replication. In the absence of a suitable animal model to study AAV in vivo, combined in cellulae and in silico studies will help to forward the understanding of the unique, bipartite AAV life cycle.

Adeno-associated virus type 2 wild-type and vector-mediated genomic integration profiles of human diploid fibroblasts analyzed by third-generation PacBio DNA sequencing.

UNLABELLED: Genome-wide analysis of adeno-associated virus (AAV) type 2 integration in HeLa cells has shown that wild-type AAV integrates at numerous genomic sites, including AAVS1 on chromosome 19q13.42. Multiple GAGY/C repeats, resembling consensus AAV Rep-binding sites are preferred, whereas rep-deficient AAV vectors (rAAV) regularly show a random integration profile. This study is the first study to analyze wild-type AAV integration in diploid human fibroblasts. Applying high-throughput third-generation PacBio-based DNA sequencing, integration profiles of wild-type AAV and rAAV are compared side by side. Bioinformatic analysis reveals that both wild-type AAV and rAAV prefer open chromatin regions. Although genomic features of AAV integration largely reproduce previous findings, the pattern of integration hot spots differs from that described in HeLa cells before. DNase-Seq data for human fibroblasts and for HeLa cells reveal variant chromatin accessibility at preferred AAV integration hot spots that correlates with variant hot spot preferences. DNase-Seq patterns of these sites in human tissues, including liver, muscle, heart, brain, skin, and embryonic stem cells further underline variant chromatin accessibility. In summary, AAV integration is dependent on cell-type-specific, variant chromatin accessibility leading to random integration profiles for rAAV, whereas wild-type AAV integration sites cluster near GAGY/C repeats. IMPORTANCE: Adeno-associated virus type 2 (AAV) is assumed to establish latency by chromosomal integration of its DNA. This is the first genome-wide analysis of wild-type AAV2 integration in diploid human cells and the first to compare wild-type to recombinant AAV vector integration side by side under identical experimental conditions. Major determinants of wild-type AAV integration represent open chromatin regions with accessible consensus AAV Rep-binding sites. The variant chromatin accessibility of different human tissues or cell types will have impact on vector targeting to be considered during gene therapy.

Retargeting transposon insertions by the adeno-associated virus Rep protein.

The Sleeping Beauty (SB), piggyBac (PB) and Tol2 transposons are promising instruments for genome engineering. Integration site profiling of SB, PB and Tol2 in human cells showed that PB and Tol2 insertions were enriched in genes, whereas SB insertions were randomly distributed. We aimed to introduce a bias into the target site selection properties of the transposon systems by taking advantage of the locus-specific integration system of adeno-associated virus (AAV). The AAV Rep protein binds to Rep recognition sequences (RRSs) in the human genome, and mediates viral integration into nearby sites. A series of fusion constructs consisting of the N-terminal DNA-binding domain of Rep and the transposases or the N57 domain of SB were generated. A plasmid-based transposition assay showed that Rep/SB yielded a 15-fold enrichment of transposition at a particular site near a targeted RRS. Genome-wide insertion site analysis indicated that an approach based on interactions between the SB transposase and Rep/N57 enriched transgene insertions at RRSs. We also provide evidence of biased insertion of the PB and Tol2 transposons. This study provides a comparative insight into target site selection properties of transposons, as well as proof-of-principle for targeted chromosomal transposition by composite protein-protein and protein-DNA interactions.

Integrated epigenome profiling of repressive histone modifications, DNA methylation and gene expression in normal and malignant urothelial cells.

Epigenetic regulation of gene expression is commonly altered in human cancer. We have observed alterations of DNA methylation and microRNA expression that reflect the biology of bladder cancer. This common disease arises by distinct pathways with low and high-grade differentiation. We hypothesized that epigenetic gene regulation reflects an interaction between histone and DNA modifications, and differences between normal and malignant urothelial cells represent carcinogenic events within bladder cancer. To test this we profiled two repressive histone modifications (H3K9m3 and H3K27m3) using ChIP-Seq, cytosine methylation using MeDIP and mRNA expression in normal and malignant urothelial cell lines. In genes with low expression we identified H3K27m3 and DNA methylation each in 20-30% of genes and both marks in 5% of genes. H3K9m3 was detected in 5-10% of genes but was not associated with overall expression. DNA methylation was more closely related to gene expression in malignant than normal cells. H3K27m3 was the epigenetic mark most specifically correlated to gene silencing. Our data suggest that urothelial carcinogenesis is accompanied by a loss of control of both DNA methylation and H3k27 methylation. From our observations we identified a panel of genes with cancer specific-epigenetic mediated aberrant expression including those with reported carcinogenic functions and members potentially mediating a positive epigenetic feedback loop. Pathway enrichment analysis revealed genes marked by H3K9m3 were involved with cell homeostasis, those marked by H3K27m3 mediated pro-carcinogenic processes and those marked with cytosine methylation were mixed in function. In 150 normal and malignant urothelial samples, our gene panel correctly estimated expression in 65% of its members. Hierarchical clustering revealed that this gene panel stratified samples according to the presence and phenotype of bladder cancer.

Integration preferences of wildtype AAV-2 for consensus rep-binding sites at numerous loci in the human genome.

Adeno-associated virus type 2 (AAV) is known to establish latency by preferential integration in human chromosome 19q13.42. The AAV non-structural protein Rep appears to target a site called AAVS1 by simultaneously binding to Rep-binding sites (RBS) present on the AAV genome and within AAVS1. In the absence of Rep, as is the case with AAV vectors, chromosomal integration is rare and random. For a genome-wide survey of wildtype AAV integration a linker-selection-mediated (LSM)-PCR strategy was designed to retrieve AAV-chromosomal junctions. DNA sequence determination revealed wildtype AAV integration sites scattered over the entire human genome. The bioinformatic analysis of these integration sites compared to those of rep-deficient AAV vectors revealed a highly significant overrepresentation of integration events near to consensus RBS. Integration hotspots included AAVS1 with 10% of total events. Novel hotspots near consensus RBS were identified on chromosome 5p13.3 denoted AAVS2 and on chromsome 3p24.3 denoted AAVS3. AAVS2 displayed seven independent junctions clustered within only 14 bp of a consensus RBS which proved to bind Rep in vitro similar to the RBS in AAVS3. Expression of Rep in the presence of rep-deficient AAV vectors shifted targeting preferences from random integration back to the neighbourhood of consensus RBS at hotspots and numerous additional sites in the human genome. In summary, targeted AAV integration is not as specific for AAVS1 as previously assumed. Rather, Rep targets AAV to integrate into open chromatin regions in the reach of various, consensus RBS homologues in the human genome.

Comparative analysis of transposable element vector systems in human cells.

Transposon-based gene vectors have become indispensable tools in vertebrate genetics for applications ranging from insertional mutagenesis and transgenesis in model species to gene therapy in humans. The transposon toolkit is expanding, but a careful, side-by-side characterization of the diverse transposon systems has been lacking. Here we compared the Sleeping Beauty (SB), piggyBac (PB), and Tol2 transposons with respect to overall activity, overproduction inhibition (OPI), target site selection, transgene copy number as well as long-term expression in human cells. SB was the most efficient system under conditions where the availability of the transposon DNA is limiting the transposition reaction including hard-to-transfect hematopoietic stem/progenitor cells (HSCs), and the most sensitive to OPI, underpinning the need for careful optimization of the transposon components. SB and PB were about equally active, and both more efficient than Tol2, under nonrestrictive conditions. All three systems provided long-term transgene expression in human cells with minimal signs of silencing. Indeed, mapping of Tol2 insertion sites revealed significant underrepresentation within chromosomal regions with H3K27me3 histone marks typically associated with transcriptionally repressed heterochromatin. SB, Tol2, and PB constitute complementary research tools for gene transfer in mammalian cells with important implications for fundamental and translational research.