Methods in Ecology and Evolution · 2019 · 53 citations · 37 references
Life SciencesEcologyPlant SciencesGenomic Data AnalysisGenomicsPlant Molecular BiologyLong-read MetagenomicsMetagenomic AnalysisPollen AnalysisShotgun Metagenomic SequencingMetabarcoding AnalysisPlant Species AbundanceShort-read Metagenomics
Abstract The ability to identify and quantify the constituent plant species that make up a mixed‐species sample of pollen has important applications in ecology, conservation, and agriculture. Recently, metabarcoding protocols have been developed for pollen that can identify constituent plant species, but there are strong reasons to doubt that metabarcoding can accurately quantify their relative abundances. A PCR‐free, shotgun metagenomics approach has greater potential for accurately quantifying species relative abundances, but applying metagenomics to eukaryotes is challenging due to low numbers of reference genomes. We have developed a pipeline, RevMet (Reverse Metagenomics) that allows reliable and semi‐quantitative characterization of the species composition of mixed‐species eukaryote samples, such as bee‐collected pollen, without requiring reference genomes. Instead, reference species are represented only by ‘genome skims’: low‐cost, low‐coverage, short‐read sequence datasets. The skims are mapped to individual long reads sequenced from mixed‐species samples using the MinION, a portable nanopore sequencing device, and each long read is uniquely assigned to a plant species. We genome‐skimmed 49 wild UK plant species, validated our pipeline with mock DNA mixtures of known composition, and then applied RevMet to pollen loads collected from wild bees. We demonstrate that RevMet can identify plant species present in mixed‐species samples at proportions of DNA ≥ 1%, with few false positives and false negatives, and reliably differentiate species represented by high versus low amounts of DNA in a sample. RevMet could readily be adapted to generate semi‐quantitative datasets for a wide range of mixed eukaryote samples. Our per‐sample costs were £90 per genome skim and £60 per pollen sample, and new versions of sequencers available now will further reduce these costs.
37
R: A Language and Environment for Statistical Computing
R Core Team · 2000 · 353K citations · Full text
Natural Sciences, Computer Science, Statistical Computing +10
The Sequence Alignment/Map format and SAMtools
Heng Li, Alec Wysoker, Tim Fennell et al. · Bioinformatics · 2009 · 67.5K citations · Full text
Cutadapt removes adapter sequences from high-throughput sequencing reads
Marcel Martin · EMBnet journal · 2011 · 35.8K citations · Full text
Minimap2: pairwise alignment for nucleotide sequences
Heng Li · Bioinformatics · 2018 · 16.8K citations · Full text