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.

Taxonomic analysis of metagenomic data with kASA.

The taxonomic analysis of sequencing data has become important in many areas of life sciences. However, currently available tools for that purpose either consume large amounts of RAM or yield insufficient quality and robustness. Here, we present kASA, a k-mer based tool capable of identifying and profiling metagenomic DNA or protein sequences with high computational efficiency and a user-definable memory footprint. We ensure both high sensitivity and precision by using an amino acid-like encoding of k-mers together with a range of multiple k's. Custom algorithms and data structures optimized for external memory storage enable a full-scale taxonomic analysis without compromise on laptop, desktop, and HPCC.

NoPeak: k-mer-based motif discovery in ChIP-Seq data without peak calling.

MOTIVATION: The discovery of sequence motifs mediating DNA-protein binding usually implies the determination of binding sites using high-throughput sequencing and peak calling. The determination of peaks, however, depends strongly on data quality and is susceptible to noise. RESULTS: Here, we present a novel approach to reliably identify transcription factor-binding motifs from ChIP-Seq data without peak detection. By evaluating the distributions of sequencing reads around the different k-mers in the genome, we are able to identify binding motifs in ChIP-Seq data that yield no results in traditional pipelines. AVAILABILITY AND IMPLEMENTATION: NoPeak is published under the GNU General Public License and available as a standalone console-based Java application at https://github.com/menzel/nopeak. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

RNA-guided retargeting of Sleeping Beauty transposition in human cells.

An ideal tool for gene therapy would enable efficient gene integration at predetermined sites in the human genome. Here we demonstrate biased genome-wide integration of the Sleeping Beauty (SB) transposon by combining it with components of the CRISPR/Cas9 system. We provide proof-of-concept that it is possible to influence the target site selection of SB by fusing it to a catalytically inactive Cas9 (dCas9) and by providing a single guide RNA (sgRNA) against the human Alu retrotransposon. Enrichment of transposon integrations was dependent on the sgRNA, and occurred in an asymmetric pattern with a bias towards sites in a relatively narrow, 300 bp window downstream of the sgRNA targets. Our data indicate that the targeting mechanism specified by CRISPR/Cas9 forces integration into genomic regions that are otherwise poor targets for SB transposition. Future modifications of this technology may allow the development of methods for specific gene insertion for precision genetic engineering.

The Two Prevalent Genotypes of an Emerging Infectious Disease, Deformed Wing Virus, Cause Equally Low Pupal Mortality and Equally High Wing Deformities in Host Honey Bees.

Deformed wing virus (DWV) is an emerging infectious disease of the honey bee (Apis mellifera) that is considered a major cause of elevated losses of honey bee colonies. DWV comprises two widespread genotypes: the originally described genotype A, and genotype B. In adult honey bees, DWV-B has been shown to be more virulent than DWV-A. However, their comparative effects on earlier host developmental stages are unknown. Here, we experimentally inoculated honey bee pupae and tested for the relative impact of DWV-A versus DWV-B on mortality and wing deformities in eclosing adults. DWV-A and DWV-B caused similar, and only slightly elevated, pupal mortality (mean 18% greater mortality than control). Both genotypes caused similarly high wing deformities in eclosing adults (mean 60% greater wing deformities than control). Viral titer was high in all of the experimentally inoculated eclosing adults, and was independent of wing deformities, suggesting that the phenotype 'deformed wings' is not directly related to viral titer or viral genotype. These viral traits favor the emergence of both genotypes of DWV by not limiting the reproduction of its vector, the ectoparasitic Varroa destructor mite, in infected pupae, and thereby facilitating the spread of DWV in honey bees infested by the mite.

Enhort: a platform for deep analysis of genomic positions.

The rise of high-throughput methods in genomic research greatly expanded our knowledge about the functionality of the genome. At the same time, the amount of available genomic position data increased massively, e.g., through genome-wide profiling of protein binding, virus integration or DNA methylation. However, there is no specialized software to investigate integration site profiles of virus integration or transcription factor binding sites by correlating the sites with the diversity of available genomic annotations. Here we present Enhort, a user-friendly software tool for relating large sets of genomic positions to a variety of annotations. It functions as a statistics based genome browser, not focused on a single locus but analyzing many genomic positions simultaneously. Enhort provides comprehensive yet easy-to-use methods for statistical analysis, visualization, and the adjustment of background models according to experimental conditions and scientific questions. Enhort is publicly available online at enhort.mni.thm.de and published under GNU General Public License.

Interactive Responses of Solanum Dulcamara to Drought and Insect Feeding are Herbivore Species-Specific.

In nature, plants are frequently subjected to multiple biotic and abiotic stresses, resulting in a convergence of adaptive responses. We hypothesised that hormonal signalling regulating defences to different herbivores may interact with drought responses, causing distinct resistance phenotypes. To test this, we studied the hormonal and transcriptomic responses of Solanum dulcamara subjected to drought and herbivory by the generalist Spodoptera exigua (beet armyworm; BAW) or the specialist Leptinotarsa decemlineata (Colorado potato beetle; CPB). Bioassays showed that the performance of BAW, but not CPB, decreased on plants under drought compared to controls. While drought did not alter BAW-induced hormonal responses, it enhanced the CPB-induced accumulation of jasmonic acid and salicylic acid (SA), and suppressed ethylene (ET) emission. Microarray analyses showed that under drought, BAW herbivory enhanced several herbivore-induced responses, including cell-wall remodelling and the metabolism of carbohydrates, lipids, and secondary metabolites. In contrast, CPB herbivory enhanced several photosynthesis-related and pathogen responses in drought-stressed plants. This may divert resources away from defence production and increase leaf nutritive value. In conclusion, while BAW suffers from the drought-enhanced defences, CPB may benefit from the effects of enhanced SA and reduced ET signalling. This suggests that the fine-tuned interaction between the plant and its specialist herbivore is sustained under drought.

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.

Ambient temperature and genotype differentially affect developmental and phenotypic plasticity in Arabidopsis thaliana.

BACKGROUND: Global increase in ambient temperatures constitute a significant challenge to wild and cultivated plant species. Forward genetic analyses of individual temperature-responsive traits have resulted in the identification of several signaling and response components. However, a comprehensive knowledge about temperature sensitivity of different developmental stages and the contribution of natural variation is still scarce and fragmented at best. RESULTS: Here, we systematically analyze thermomorphogenesis throughout a complete life cycle in ten natural Arabidopsis thaliana accessions grown under long day conditions in four different temperatures ranging from 16 to 28 °C. We used Q10, GxE, phenotypic divergence and correlation analyses to assess temperature sensitivity and genotype effects of more than 30 morphometric and developmental traits representing five phenotype classes. We found that genotype and temperature differentially affected plant growth and development with variing strengths. Furthermore, overall correlations among phenotypic temperature responses was relatively low which seems to be caused by differential capacities for temperature adaptations of individual accessions. CONCLUSION: Genotype-specific temperature responses may be attractive targets for future forward genetic approaches and accession-specific thermomorphogenesis maps may aid the assessment of functional relevance of known and novel regulatory components.

Genome-wide profiling of S/MAR-based replicon contact sites.

Autonomously replicating vectors represent a simple and versatile model system for genetic modifications, but their localization in the nucleus and effect on endogenous gene expression is largely unknown. Using circular chromosome conformation capture we mapped genomic contact sites of S/MAR-based replicons in HeLa cells. The influence of cis-active sequences on genomic localization was assessed using replicons containing either an insulator sequence or an intron. While the original and the insulator-containing replicons displayed distinct contact sites, the intron-containing replicon showed a rather broad genomic contact pattern. Our results indicate a preference for certain chromatin structures and a rather non-dynamic behaviour during mitosis. Independent of inserted cis-active elements established vector molecules reside preferentially within actively transcribed regions, especially within promoter sequences and transcription start sites. However, transcriptome analyses revealed that established S/MAR-based replicons do not alter gene expression profiles of host genome. Knowledge of preferred contact sites of exogenous DNA, e.g. viral or non-viral episomes, contribute to our understanding of episome behaviour in the nucleus and can be used for vector improvement and guiding of DNA sequences to specific subnuclear sites.

Unity in defence: honeybee workers exhibit conserved molecular responses to diverse pathogens.

BACKGROUND: Organisms typically face infection by diverse pathogens, and hosts are thought to have developed specific responses to each type of pathogen they encounter. The advent of transcriptomics now makes it possible to test this hypothesis and compare host gene expression responses to multiple pathogens at a genome-wide scale. Here, we performed a meta-analysis of multiple published and new transcriptomes using a newly developed bioinformatics approach that filters genes based on their expression profile across datasets. Thereby, we identified common and unique molecular responses of a model host species, the honey bee (Apis mellifera), to its major pathogens and parasites: the Microsporidia Nosema apis and Nosema ceranae, RNA viruses, and the ectoparasitic mite Varroa destructor, which transmits viruses. RESULTS: We identified a common suite of genes and conserved molecular pathways that respond to all investigated pathogens, a result that suggests a commonality in response mechanisms to diverse pathogens. We found that genes differentially expressed after infection exhibit a higher evolutionary rate than non-differentially expressed genes. Using our new bioinformatics approach, we unveiled additional pathogen-specific responses of honey bees; we found that apoptosis appeared to be an important response following microsporidian infection, while genes from the immune signalling pathways, Toll and Imd, were differentially expressed after Varroa/virus infection. Finally, we applied our bioinformatics approach and generated a gene co-expression network to identify highly connected (hub) genes that may represent important mediators and regulators of anti-pathogen responses. CONCLUSIONS: Our meta-analysis generated a comprehensive overview of the host metabolic and other biological processes that mediate interactions between insects and their pathogens. We identified key host genes and pathways that respond to phylogenetically diverse pathogens, representing an important source for future functional studies as well as offering new routes to identify or generate pathogen resilient honey bee stocks. The statistical and bioinformatics approaches that were developed for this study are broadly applicable to synthesize information across transcriptomic datasets. These approaches will likely have utility in addressing a variety of biological questions.

Auxin-induced expression divergence between Arabidopsis species may originate within the TIR1/AFB-AUX/IAA-ARF module.

Auxin is an essential regulator of plant growth and development, and auxin signaling components are conserved among land plants. Yet, a remarkable degree of natural variation in physiological and transcriptional auxin responses has been described among Arabidopsis thaliana accessions. As intraspecies comparisons offer only limited genetic variation, we here inspect the variation of auxin responses between A. thaliana and A. lyrata. This approach allowed the identification of conserved auxin response genes including novel genes with potential relevance for auxin biology. Furthermore, promoter divergences were analyzed for putative sources of variation. De novo motif discovery identified novel and variants of known elements with potential relevance for auxin responses, emphasizing the complex, and yet elusive, code of element combinations accounting for the diversity in transcriptional auxin responses. Furthermore, network analysis revealed correlations of interspecies differences in the expression of AUX/IAA gene clusters and classic auxin-related genes. We conclude that variation in general transcriptional and physiological auxin responses may originate substantially from functional or transcriptional variations in the TIR1/AFB, AUX/IAA, and ARF signaling network. In that respect, AUX/IAA gene expression divergence potentially reflects differences in the manner in which different species transduce identical auxin signals into gene expression responses.

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.

Elevated virulence of an emerging viral genotype as a driver of honeybee loss.

Emerging infectious diseases (EIDs) have contributed significantly to the current biodiversity crisis, leading to widespread epidemics and population loss. Owing to genetic variation in pathogen virulence, a complete understanding of species decline requires the accurate identification and characterization of EIDs. We explore this issue in the Western honeybee, where increasing mortality of populations in the Northern Hemisphere has caused major concern. Specifically, we investigate the importance of genetic identity of the main suspect in mortality, deformed wing virus (DWV), in driving honeybee loss. Using laboratory experiments and a systematic field survey, we demonstrate that an emerging DWV genotype (DWV-B) is more virulent than the established DWV genotype (DWV-A) and is widespread in the landscape. Furthermore, we show in a simple model that colonies infected with DWV-B collapse sooner than colonies infected with DWV-A. We also identify potential for rapid DWV evolution by revealing extensive genome-wide recombination in vivo The emergence of DWV-B in naive honeybee populations, including via recombination with DWV-A, could be of significant ecological and economic importance. Our findings emphasize that knowledge of pathogen genetic identity and diversity is critical to understanding drivers of species decline.

A Helitron transposon reconstructed from bats reveals a novel mechanism of genome shuffling in eukaryotes.

Helitron transposons capture and mobilize gene fragments in eukaryotes, but experimental evidence for their transposition is lacking in the absence of an isolated active element. Here we reconstruct Helraiser, an ancient element from the bat genome, and use this transposon as an experimental tool to unravel the mechanism of Helitron transposition. A hairpin close to the 3'-end of the transposon functions as a transposition terminator. However, the 3'-end can be bypassed by the transposase, resulting in transduction of flanking sequences to new genomic locations. Helraiser transposition generates covalently closed circular intermediates, suggestive of a replicative transposition mechanism, which provides a powerful means to disseminate captured transcriptional regulatory signals across the genome. Indeed, we document the generation of novel transcripts by Helitron promoter capture both experimentally and by transcriptome analysis in bats. Our results provide mechanistic insight into Helitron transposition, and its impact on diversification of gene function by genome shuffling.

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.

Genomic analysis of Sleeping Beauty transposon integration in human somatic cells.

The Sleeping Beauty (SB) transposon is a non-viral integrating vector system with proven efficacy for gene transfer and functional genomics. However, integration efficiency is negatively affected by the length of the transposon. To optimize the SB transposon machinery, the inverted repeats and the transposase gene underwent several modifications, resulting in the generation of the hyperactive SB100X transposase and of the high-capacity "sandwich" (SA) transposon. In this study, we report a side-by-side comparison of the SA and the widely used T2 arrangement of transposon vectors carrying increasing DNA cargoes, up to 18 kb. Clonal analysis of SA integrants in human epithelial cells and in immortalized keratinocytes demonstrates stability and integrity of the transposon independently from the cargo size and copy number-dependent expression of the cargo cassette. A genome-wide analysis of unambiguously mapped SA integrations in keratinocytes showed an almost random distribution, with an overrepresentation in repetitive elements (satellite, LINE and small RNAs) compared to a library representing insertions of the first-generation transposon vector and to gammaretroviral and lentiviral libraries. The SA transposon/SB100X integrating system therefore shows important features as a system for delivering large gene constructs for gene therapy applications.

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.