Protective role of the vulture facial skin and gut microbiomes aid adaptation to scavenging.

BACKGROUND: Vultures have adapted the remarkable ability to feed on carcasses that may contain microorganisms that would be pathogenic to most other animals. The holobiont concept suggests that the genetic basis of such adaptation may not only lie within their genomes, but additionally in their associated microbes. To explore this, we generated shotgun DNA sequencing datasets of the facial skin and large intestine microbiomes of the black vulture (Coragyps atratus) and the turkey vulture (Cathartes aura). We characterized the functional potential and taxonomic diversity of their microbiomes, the potential pathogenic challenges confronted by vultures, and the microbial taxa and genes that could play a protective role on the facial skin and in the gut. RESULTS: We found microbial taxa and genes involved in diseases, such as dermatitis and pneumonia (more abundant on the facial skin), and gas gangrene and food poisoning (more abundant in the gut). Interestingly, we found taxa and functions with potential for playing beneficial roles, such as antilisterial bacteria in the gut, and genes for the production of antiparasitics and insecticides on the facial skin. Based on the identified phages, we suggest that phages aid in the control and possibly elimination, as in phage therapy, of microbes reported as pathogenic to a variety of species. Interestingly, we identified Adineta vaga in the gut, an invertebrate that feeds on dead bacteria and protozoans, suggesting a defensive predatory mechanism. Finally, we suggest a colonization resistance role through biofilm formation played by Fusobacteria and Clostridia in the gut. CONCLUSIONS: Our results highlight the importance of complementing genomic analyses with metagenomics in order to obtain a clearer understanding of the host-microbial alliance and show the importance of microbiome-mediated health protection for adaptation to extreme diets, such as scavenging.

Do Vertebrate Gut Metagenomes Confer Rapid Ecological Adaptation?

During times of rapid environmental change, survival of most vertebrate populations depends on their phenomic plasticity. Although differential gene-expression and post-transcriptional processes of the host genome receive focus as the main molecular mechanisms, growing evidence points to the gut microbiota as a key driver defining hosts' phenotypes. We propose that the plasticity of the gut microbiota might be an essential factor determining phenomic plasticity of vertebrates, and that it might play a pivotal role when vertebrates acclimate and adapt to fast environmental variation. We contemplate some key questions and suggest methodological approaches and experimental designs that can be used to evaluate whether gut microorganisms provide a boost of plasticity to vertebrates' phenomes, thereby increasing their acclimation and adaptation capacity.

DAMe: a toolkit for the initial processing of datasets with PCR replicates of double-tagged amplicons for DNA metabarcoding analyses.

BACKGROUND: DNA metabarcoding is an approach for identifying multiple taxa in an environmental sample using specific genetic loci and taxa-specific primers. When combined with high-throughput sequencing it enables the taxonomic characterization of large numbers of samples in a relatively time- and cost-efficient manner. One recent laboratory development is the addition of 5'-nucleotide tags to both primers producing double-tagged amplicons and the use of multiple PCR replicates to filter erroneous sequences. However, there is currently no available toolkit for the straightforward analysis of datasets produced in this way. RESULTS: We present DAMe, a toolkit for the processing of datasets generated by double-tagged amplicons from multiple PCR replicates derived from an unlimited number of samples. Specifically, DAMe can be used to (i) sort amplicons by tag combination, (ii) evaluate PCR replicates dissimilarity, and (iii) filter sequences derived from sequencing/PCR errors, chimeras, and contamination. This is attained by calculating the following parameters: (i) sequence content similarity between the PCR replicates from each sample, (ii) reproducibility of each unique sequence across the PCR replicates, and (iii) copy number of the unique sequences in each PCR replicate. We showcase the insights that can be obtained using DAMe prior to taxonomic assignment, by applying it to two real datasets that vary in their complexity regarding number of samples, sequencing libraries, PCR replicates, and used tag combinations. Finally, we use a third mock dataset to demonstrate the impact and importance of filtering the sequences with DAMe. CONCLUSIONS: DAMe allows the user-friendly manipulation of amplicons derived from multiple samples with PCR replicates built in a single or multiple sequencing libraries. It allows the user to: (i) collapse amplicons into unique sequences and sort them by tag combination while retaining the sample identifier and copy number information, (ii) identify sequences carrying unused tag combinations, (iii) evaluate the comparability of PCR replicates of the same sample, and (iv) filter tagged amplicons from a number of PCR replicates using parameters of minimum length, copy number, and reproducibility across the PCR replicates. This enables an efficient analysis of complex datasets, and ultimately increases the ease of handling datasets from large-scale studies.

Metagenomic Analysis from the Interior of a Speleothem in Tjuv-Ante's Cave, Northern Sweden.

Speleothems are secondary mineral deposits normally formed by water supersaturated with calcium carbonate percolating into underground caves, and are often associated with low-nutrient and mostly non-phototrophic conditions. Tjuv-Ante's cave is a shallow-depth cave formed by the action of waves, with granite and dolerite as major components, and opal-A and calcite as part of the speleothems, making it a rare kind of cave. We generated two DNA shotgun sequencing metagenomic datasets from the interior of a speleothem from Tjuv-Ante's cave representing areas of old and relatively recent speleothem formation. We used these datasets to perform i) an evaluation of the use of these speleothems as past biodiversity archives, ii) functional and taxonomic profiling of the speleothem's different formation periods, and iii) taxonomic comparison of the metagenomic results to previous microscopic analyses from a nearby speleothem of the same cave. Our analyses confirm the abundance of Actinobacteria and fungi as previously reported by microscopic analyses on this cave, however we also discovered a larger biodiversity. Interestingly, we identified photosynthetic genes, as well as genes related to iron and sulphur metabolism, suggesting the presence of chemoautotrophs. Furthermore, we identified taxa and functions related to biomineralization. However, we could not confidently establish the use of this type of speleothems as biological paleoarchives due to the potential leaching from the outside of the cave and the DNA damage that we propose has been caused by the fungal chemical etching.

Odintifier--A computational method for identifying insertions of organellar origin from modern and ancient high-throughput sequencing data based on haplotype phasing.

BACKGROUND: Cellular organelles with genomes of their own (e.g. plastids and mitochondria) can pass genetic sequences to other organellar genomes within the cell in many species across the eukaryote phylogeny. The extent of the occurrence of these organellar-derived inserted sequences (odins) is still unknown, but if not accounted for in genomic and phylogenetic studies, they can be a source of error. However, if correctly identified, these inserted sequences can be used for evolutionary and comparative genomic studies. Although such insertions can be detected using various laboratory and bioinformatic strategies, there is currently no straightforward way to apply them as a standard organellar genome assembly on next-generation sequencing data. Furthermore, most current methods for identification of such insertions are unsuitable for use on non-model organisms or ancient DNA datasets. RESULTS: We present a bioinformatic method that uses phasing algorithms to reconstruct both source and inserted organelle sequences. The method was tested in different shotgun and organellar-enriched DNA high-throughput sequencing (HTS) datasets from ancient and modern samples. Specifically, we used datasets from lions (Panthera leo ssp. and Panthera leo leo) to characterize insertions from mitochondrial origin, and from common grapevine (Vitis vinifera) and bugle (Ajuga reptans) to characterize insertions derived from plastid genomes. Comparison of the results against other available organelle genome assembly methods demonstrated that our new method provides an improvement in the sequence assembly. CONCLUSION: Using datasets from a wide range of species and different levels of complexity we showed that our novel bioinformatic method based on phasing algorithms can be used to achieve the next two goals: i) reference-guided assembly of chloroplast/mitochondrial genomes from HTS data and ii) identification and simultaneous assembly of odins. This method represents the first application of haplotype phasing for automatic detection of odins and reference-based organellar genome assembly.

Environmental genes and genomes: understanding the differences and challenges in the approaches and software for their analyses.

DNA-based taxonomic and functional profiling is widely used for the characterization of organismal communities across a rapidly increasing array of research areas that include the role of microbiomes in health and disease, biomonitoring, and estimation of both microbial and metazoan species richness. Two principal approaches are currently used to assign taxonomy to DNA sequences: DNA metabarcoding and metagenomics. When initially developed, each of these approaches mandated their own particular methods for data analysis; however, with the development of high-throughput sequencing (HTS) techniques they have begun to share many aspects in data set generation and processing. In this review we aim to define the current characteristics, goals and boundaries of each field, and describe the different software used for their analysis. We argue that an appreciation of the potential and limitations of each method can help underscore the improvements required by each field so as to better exploit the richness of current HTS-based data sets.

DivA: detection of non-homologous and very divergent regions in protein sequence alignments.

BACKGROUND: Sequence alignments are used to find evidence of homology but sometimes contain regions that are difficult to align which can interfere with the quality of the subsequent analyses. Although it is possible to remove problematic regions manually, this is non-practical in large genome scale studies, and the results suffer from irreproducibility arising from subjectivity. Some automated alignment trimming methods have been developed to remove problematic regions in alignments but these mostly act by removing complete columns or complete sequences from the MSA, discarding a lot of informative sites. FINDINGS: Here we present a tool that identifies Divergent windows in protein sequence Alignments (DivA). DivA makes no assumptions on evolutionary models, and it is ideal for detecting incorrectly annotated segments within individual gene sequences. DivA works with a sliding-window approach to estimate four divergence-based parameters and their outlier values. It then classifies a window of a sequence of an alignment as very divergent (potentially non-homologous) if it presents a combination of outlier values for the four parameters it calculates. The windows classified as very divergent can optionally be masked in the alignment. CONCLUSIONS: DivA automatically identifies very divergent and incorrectly annotated genic regions in MSAs avoiding the subjective and time-consuming problem of manual annotation. The output is clear to interpret and allows the user to take more informed decisions for reducing the amount of sequence discarded but still finding the potentially erroneous and non-homologous regions.

Second generation sequencing and morphological faecal analysis reveal unexpected foraging behaviour by Myotis nattereri (Chiroptera, Vespertilionidae) in winter.

BACKGROUND: Temperate winters produce extreme energetic challenges for small insectivorous mammals. Some bat species inhabiting locations with mild temperate winters forage during brief inter-torpor normothermic periods of activity. However, the winter diet of bats in mild temperate locations is studied infrequently. Although microscopic analyses of faeces have traditionally been used to characterise bat diet, recently the coupling of PCR with second generation sequencing has offered the potential to further advance our understanding of animal dietary composition and foraging behaviour by allowing identification of a much greater proportion of prey items often with increased taxonomic resolution. We used morphological analysis and Illumina-based second generation sequencing to study the winter diet of Natterer's bat (Myotis nattereri) and compared the results obtained from these two approaches. For the first time, we demonstrate the applicability of the Illumina MiSeq platform as a data generation source for bat dietary analyses. RESULTS: Faecal pellets collected from a hibernation site in southern England during two winters (December-March 2009-10 and 2010-11), indicated that M. nattereri forages throughout winter at least in a location with a mild winter climate. Through morphological analysis, arthropod fragments from seven taxonomic orders were identified. A high proportion of these was non-volant (67.9% of faecal pellets) and unexpectedly included many lepidopteran larvae. Molecular analysis identified 43 prey species from six taxonomic orders and confirmed the frequent presence of lepidopteran species that overwinter as larvae. CONCLUSIONS: The winter diet of M. nattereri is substantially different from other times of the year confirming that this species has a wide and adaptable dietary niche. Comparison of DNA derived from the prey to an extensive reference dataset of potential prey barcode sequences permitted fine scale taxonomic resolution of prey species. The high occurrence of non-volant prey suggests that gleaning allows prey capture at low ambient temperatures when the abundance of flying insects may be substantially reduced. Interesting questions arise as to how M. nattereri might successfully locate and capture some of the non-volant prey species encountered in its faeces. The consumption of lepidopteran larvae such as cutworms suggests that M. nattereri eats agricultural pest species.