Variant profiling of evolving prokaryotic populations.

Genomic heterogeneity of bacterial species is observed and studied in experimental evolution experiments and clinical diagnostics, and occurs as micro-diversity of natural habitats. The challenge for genome research is to accurately capture this heterogeneity with the currently used short sequencing reads. Recent advances in NGS technologies improved the speed and coverage and thus allowed for deep sequencing of bacterial populations. This facilitates the quantitative assessment of genomic heterogeneity, including low frequency alleles or haplotypes. However, false positive variant predictions due to sequencing errors and mapping artifacts of short reads need to be prevented. We therefore created VarCap, a workflow for the reliable prediction of different types of variants even at low frequencies. In order to predict SNPs, InDels and structural variations, we evaluated the sensitivity and accuracy of different software tools using synthetic read data. The results suggested that the best sensitivity could be reached by a union of different tools, however at the price of increased false positives. We identified possible reasons for false predictions and used this knowledge to improve the accuracy by post-filtering the predicted variants according to properties such as frequency, coverage, genomic environment/localization and co-localization with other variants. We observed that best precision was achieved by using an intersection of at least two tools per variant. This resulted in the reliable prediction of variants above a minimum relative abundance of 2%. VarCap is designed for being routinely used within experimental evolution experiments or for clinical diagnostics. The detected variants are reported as frequencies within a VCF file and as a graphical overview of the distribution of the different variant/allele/haplotype frequencies. The source code of VarCap is available at https://github.com/ma2o/VarCap. In order to provide this workflow to a broad community, we implemeted VarCap on a Galaxy webserver, which is accessible at http://galaxy.csb.univie.ac.at.

Expressional and functional variation of horizontally acquired cellulases in the nematode Pristionchus pacificus.

Various whole genome-sequencing projects in the nematode phylum have revealed the widespread occurrence of horizontal gene transfer from different sources. Pristionchus pacificus was the first non-plant parasitic nematode that was found to contain cellulase genes in its genome and to have cellulolytic activity when grown in a monoxenic culture on Escherichia coli. The P. pacificus reference strain PS312 has seven cellulase genes, all of which were acquired by horizontal gene transfer. Previous phylogenomic studies indicated that the acquisition occurred at the base of the genus Pristionchus and was followed by rapid gene duplications and gene turnover. However, little was known about the protein domain architecture, gene expression and the functionality of individual proteins. Here, we analyzed the protein domain architecture, studied the expression at various developmentally stages and tried to induce cellulase gene expression by feeding nematodes with different polysaccharides. Only two of the encoded proteins, Ppa-CEL-2 and Ppa-CEL-3, have a carbohydrate-binding module. Interestingly, these were also the ones with developmental gene expression regulation. Ppa-cel-2 shows high expression in larval and adult stages but is only low expressed in eggs, while Ppa-cel-3 is highly upregulated in adult worms but hardly detectable in any other stage. Ppa-CEL-1, has a catalytic domain similar to Ppa-CEL-2 and Ppa-CEL-3, but lacks a carbohydrate-binding module. The other four cellulases have a very low transcriptional expression correlating with putative incompleteness of their catalytic domain. While, the expression of none of the genes is inducible by polysaccharides, zymographic studies and mass spectrometry indicate that Ppa-CEL-2 and Ppa-CEL-3 are the only two cellulases contributing to cellulase activity in carboxymethylcellulose. Thus, the cellulases of P. pacificus differ in their protein domain architecture, gene expression and functionality. These results indicate that horizontal gene transfer-acquired genes undergo rapid evolutionary changes that affect all aspects of their molecular biology.

Horizontal gene transfer of microbial cellulases into nematode genomes is associated with functional assimilation and gene turnover.

BACKGROUND: Natural acquisition of novel genes from other organisms by horizontal or lateral gene transfer is well established for microorganisms. There is now growing evidence that horizontal gene transfer also plays important roles in the evolution of eukaryotes. Genome-sequencing and EST projects of plant and animal associated nematodes such as Brugia, Meloidogyne, Bursaphelenchus and Pristionchus indicate horizontal gene transfer as a key adaptation towards parasitism and pathogenicity. However, little is known about the functional activity and evolutionary longevity of genes acquired by horizontal gene transfer and the mechanisms favoring such processes. RESULTS: We examine the transfer of cellulase genes to the free-living and beetle-associated nematode Pristionchus pacificus, for which detailed phylogenetic knowledge is available, to address predictions by evolutionary theory for successful gene transfer. We used transcriptomics in seven Pristionchus species and three other related diplogastrid nematodes with a well-defined phylogenetic framework to study the evolution of ancestral cellulase genes acquired by horizontal gene transfer. We performed intra-species, inter-species and inter-genic analysis by comparing the transcriptomes of these ten species and tested for cellulase activity in each species. Species with cellulase genes in their transcriptome always exhibited cellulase activity indicating functional integration into the host's genome and biology. The phylogenetic profile of cellulase genes was congruent with the species phylogeny demonstrating gene longevity. Cellulase genes show notable turnover with elevated birth and death rates. Comparison by sequencing of three selected cellulase genes in 24 natural isolates of Pristionchus pacificus suggests these high evolutionary dynamics to be associated with copy number variations and positive selection. CONCLUSION: We could demonstrate functional integration of acquired cellulase genes into the nematode's biology as predicted by theory. Thus, functional assimilation, remarkable gene turnover and selection might represent key features of horizontal gene transfer events in nematodes.

The Pristionchus pacificus genome provides a unique perspective on nematode lifestyle and parasitism.

Here we present a draft genome sequence of the nematode Pristionchus pacificus, a species that is associated with beetles and is used as a model system in evolutionary biology. With 169 Mb and 23,500 predicted protein-coding genes, the P. pacificus genome is larger than those of Caenorhabditis elegans and the human parasite Brugia malayi. Compared to C. elegans, the P. pacificus genome has more genes encoding cytochrome P450 enzymes, glucosyltransferases, sulfotransferases and ABC transporters, many of which were experimentally validated. The P. pacificus genome contains genes encoding cellulase and diapausin, and cellulase activity is found in P. pacificus secretions, indicating that cellulases can be found in nematodes beyond plant parasites. The relatively higher number of detoxification and degradation enzymes in P. pacificus is consistent with its necromenic lifestyle and might represent a preadaptation for parasitism. Thus, comparative genomics analysis of three ecologically distinct nematodes offers a unique opportunity to investigate the association between genome structure and lifestyle.