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.

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.

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.

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.

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.

Ultra-deep profiling of alternatively spliced Drosophila Dscam isoforms by circularization-assisted multi-segment sequencing.

The Drosophila melanogaster gene Dscam (Down syndrome cell adhesion molecule) can generate thousands of different ectodomains via mutual exclusive splicing of three large exon clusters. The isoform diversity plays a profound role in both neuronal wiring and pathogen recognition. However, the isoform expression pattern at the global level remained unexplored. Here, we developed a novel method that allows for direct quantification of the alternatively spliced exon combinations from over hundreds of millions of Dscam transcripts in one sequencing run. With unprecedented sequencing depth, we detected a total of 18,496 isoforms, out of 19,008 theoretically possible combinations. Importantly, we demonstrated that alternative splicing between different clusters is independent. Moreover, the isoforms were expressed across a broad dynamic range, with significant bias in cell/tissue and developmental stage-specific patterns. Hitherto underappreciated, such bias can dramatically reduce the ability of neurons to display unique surface receptor codes. Therefore, the seemingly excessive diversity encoded in the Dscam locus might nevertheless be essential for a robust self and non-self discrimination in neurons.

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.

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.

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.

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.

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.

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.

Global profiling of miRNAs and the hairpin precursors: insights into miRNA processing and novel miRNA discovery.

MicroRNAs (miRNAs) constitute an important class of small regulatory RNAs that are derived from distinct hairpin precursors (pre-miRNAs). In contrast to mature miRNAs, which have been characterized in numerous genome-wide studies of different organisms, research on global profiling of pre-miRNAs is limited. Here, using massive parallel sequencing, we have performed global characterization of both mouse mature and precursor miRNAs. In total, 87 369 704 and 252 003 sequencing reads derived from 887 mature and 281 precursor miRNAs were obtained, respectively. Our analysis revealed new aspects of miRNA/pre-miRNA processing and modification, including eight Ago2-cleaved pre-miRNAs, eight new instances of miRNA editing and exclusively 5' tailed mirtrons. Furthermore, based on the sequences of both mature and precursor miRNAs, we developed a miRNA discovery pipeline, miRGrep, which does not rely on the availability of genome reference sequences. In addition to 239 known mouse pre-miRNAs, miRGrep predicted 41 novel ones with high confidence. Similar as known ones, the mature miRNAs derived from most of these novel loci showed both reduced abundance following Dicer knockdown and the binding with Argonaute2. Evaluation on data sets obtained from Caenorhabditis elegans and Caenorhabditis sp.11 demonstrated that miRGrep could be widely used for miRNA discovery in metazoans, especially in those without genome reference sequences.

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.

Retargeting sleeping beauty transposon insertions by engineered zinc finger DNA-binding domains.

The Sleeping Beauty (SB) transposon is a nonviral, integrating vector system with proven efficacy in preclinical animal models, and thus holds promise for future clinical applications. However, SB has a close-to-random insertion profile that could lead to genotoxic effects, thereby presenting a potential safety issue. We evaluated zinc finger (ZF) DNA-binding domains (DBDs) for their abilities to introduce a bias into SB's insertion profile. E2C, that binds a unique site in the erbB-2 gene, mediated locus-specific transposon insertions at low frequencies. A novel ZF targeting LINE1 repeats, ZF-B, showed specific binding to an 18-bp site represented by ~12,000 copies in the human genome. We mapped SB insertions using linear-amplification (LAM)-PCR and Illumina sequencing. Targeted insertions with ZF-B peaked at approximately fourfold enrichment of transposition around ZF-B binding sites yielding ~45% overall frequency of insertion into LINE1. A decrease in the ZF-B dataset with respect to transposon insertions in genes was found, suggesting that LINE1 repeats act as a sponge that "soak up" a fraction of SB insertions and thereby redirect them away from genes. Improvements in ZF technology and a careful choice of targeted genomic regions may improve the safety profile of SB for future clinical applications.

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 genomic integration profiling of Sleeping Beauty transposons mobilized with high efficacy from integrase-defective lentiviral vectors in primary human cells.

It has been previously shown that integrase-defective HIV-1-based gene vectors can serve, with moderate efficiency, as substrate for DNA transposition by a transiently expressed Sleeping Beauty (SB) transposase. Here, we describe the enhanced gene transfer properties of a HIV-1/SB hybrid vector that allows efficient DNA transposition, facilitated by the hyperactive SB100X transposase, from vector DNA intermediates in primary human cells. Potent transposase-dependent integration of genetic cargo carried by the hybrid HIV-1/SB vector (up to 160-fold above background) is reported in human cell lines as well as in primary human fibroblasts and keratinocytes. The efficiency of transgene integration in context of the newly developed hybrid vector is comparable with that of conventional lentiviral vectors (LVs). Integration profiles of integrating HIV-1-derived vectors and SB transposons mobilized from LVs are investigated by deep sequencing of a large number of integration sites. A significant bias of lentiviral integrations in genes is reported, confirming that biological properties of the viral integration machinery facilitate preferred insertion into actively transcribed genomic regions. In sharp contrast, lentiviral insertions catalyzed by the SB100X transposase are far more random with respect to genes. Based on these properties, HIV-1/SB vectors may become valuable tools for genetic engineering and therapeutic gene transfer.

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.