AFS: identification and quantification of species composition by metagenomic sequencing.

Summary: DNA-based methods to detect and quantify taxon composition in biological materials are often based on species-specific polymerase chain reaction, limited to detecting species targeted by the assay. Next-generation sequencing overcomes this drawback by untargeted shotgun sequencing of whole metagenomes at affordable cost. Here we present AFS, a software pipeline for quantification of species composition in food. AFS uses metagenomic shotgun sequencing and sequence read counting to infer species proportions. Using Illumina data from a reference sausage comprising four species, we reveal that AFS is independent of the sequencing assay and library preparation protocol. Cost-saving short (50-bp) single-end reads and Nextera ® library preparation yield reliable results. Availability and Implementation: Datasets, binaries and usage instructions are available under http://all-food-seq.sourceforge.net. Raw data is available at NCBI's SRA with accession number PRJNA271645. Contact: hankeln@uni-mainz.de. Supplementary information: Supplementary data are available at Bioinformatics online.

Chironomus riparius (Diptera) genome sequencing reveals the impact of minisatellite transposable elements on population divergence.

Active transposable elements (TEs) may result in divergent genomic insertion and abundance patterns among conspecific populations. Upon secondary contact, such divergent genetic backgrounds can theoretically give rise to classical Dobzhansky-Muller incompatibilities (DMI), thus contributing to the evolution of endogenous genetic barriers and eventually causing population divergence. We investigated differential TE abundance among conspecific populations of the nonbiting midge Chironomus riparius and evaluated their potential role in causing endogenous genetic incompatibilities between these populations. We focussed on a Chironomus-specific TE, the minisatellite-like Cla-element, whose activity is associated with speciation in the genus. Using a newly generated and annotated draft genome for a genomic study with five natural C. riparius populations, we found highly population-specific TE insertion patterns with many private insertions. A significant correlation of the pairwise FST estimated from genomewide single-nucleotide polymorphisms (SNPs) and the FST estimated from TEs is consistent with drift as the major force driving TE population differentiation. However, the significantly higher Cla-element FST level due to a high proportion of differentially fixed Cla-element insertions also indicates selection against segregating (i.e. heterozygous) insertions. With reciprocal crossing experiments and fluorescent in situ hybridization of Cla-elements to polytene chromosomes, we documented phenotypic effects on female fertility and chromosomal mispairings. We propose that the inferred negative selection on heterozygous Cla-element insertions may cause endogenous genetic barriers and therefore acts as DMI among C. riparius populations. The intrinsic genomic turnover exerted by TEs may thus have a direct impact on population divergence that is operationally different from drift and local adaptation.

All-Food-Seq (AFS): a quantifiable screen for species in biological samples by deep DNA sequencing.

BACKGROUND: DNA-based methods like PCR efficiently identify and quantify the taxon composition of complex biological materials, but are limited to detecting species targeted by the choice of the primer assay. We show here how untargeted deep sequencing of foodstuff total genomic DNA, followed by bioinformatic analysis of sequence reads, facilitates highly accurate identification of species from all kingdoms of life, at the same time enabling quantitative measurement of the main ingredients and detection of unanticipated food components. RESULTS: Sequence data simulation and real-case Illumina sequencing of DNA from reference sausages composed of mammalian (pig, cow, horse, sheep) and avian (chicken, turkey) species are able to quantify material correctly at the 1% discrimination level via a read counting approach. An additional metagenomic step facilitates identification of traces from animal, plant and microbial DNA including unexpected species, which is prospectively important for the detection of allergens and pathogens. CONCLUSIONS: Our data suggest that deep sequencing of total genomic DNA from samples of heterogeneous taxon composition promises to be a valuable screening tool for reference species identification and quantification in biosurveillance applications like food testing, potentially alleviating some of the problems in taxon representation and quantification associated with targeted PCR-based approaches.

Divergent roles of the Drosophila melanogaster globins.

In contrast to long-held assumptions, the gene repertoire of most insects includes hemoglobins. Analyses of the genome of the fruitfly Drosophila melanogaster identified three distinct hemoglobin genes (glob1, glob2, and glob3). While glob1 is predominantly associated with the tracheal system and fat body, glob2 and glob3 are almost exclusively expressed in the testis. The physiological role of globins in Drosophila is uncertain. Here, we studied the functions of the three globins in a cell culture system. Drosophila Schneider 2 (S2) cells were stably transfected with each of the three globins and the empty vector as control. Under hypoxia (1% atmospheric O2), only glob1 overexpression enhanced the activity of mitochondrial oxidases and the ATP content. However, the positive effect of glob1 expression disappeared after 24h hypoxia, suggesting metabolic adaptations of the S2 cells. glob2 and glob3 had no positive effect on hypoxia-survival. After application of oxidative stress by H2O2, glob2 dramatically enhanced the viability of S2 cells. Evaluation of the intracellular localization of the globins using specific antibodies and green fluorescent protein-fusion constructs suggested that glob1 and glob2 most likely reside in the cytoplasm, while glob3 is associated with structures that may represent parts of the intracellular transport machinery. In silico analyses of public RNA-Seq data from different developmental stages provided that glob1 is co-expressed with genes of the aerobic energy metabolism, while glob2 and glob3 expression can be related to spermatogenesis and reproduction. Together, the results indicate divergent functions of the Drosophila globins: glob1 may play a role in the O2-dependent metabolism while glob2 may protect spermatogenesis from reactive oxygen species.

Knockdown of Drosophila hemoglobin suggests a role in O2 homeostasis.

Almost all insects are equipped with a tracheal system, which appears to be sufficient for O2 supply even in phases of high metabolic activity. Therefore, with the exception of a few species dwelling in hypoxic habitats, specialized respiratory proteins had been considered unnecessary in insects. The recent discovery and apparently universal presence of intracellular hemoglobins in insects has remained functionally unexplained. The fruitfly Drosophila melanogaster harbors three different globin genes (referred to as glob1-3). Glob1 is the most highly expressed globin and essentially occurs in the tracheal system and the fat body. To better understand the functions of insect globins, the levels of glob1 were modulated in Drosophila larvae and adults by RNAi-mediated knockdown and transgenic over-expression. No effects on the development were observed in flies with manipulated glob1 levels. However, the knockdown of glob1 led to a significantly reduced survival rate of adult flies under hypoxia (5% and 1.5% O2). Surprisingly, the glob1 knockdown flies also displayed increased resistance towards the reactive oxygen species-forming agent paraquat, which may be explained by a restricted availability of O2 resulting in decreased formation of harmful O2(-). In summary, our results suggest an important functional role of glob1 in O2 homeostasis, possibly by enhancing O2 supply.