Hybridization ddRAD‐sequencing for population genomics of nonmodel plants using highly degraded historical specimen DNA

Species’ responses at the genetic level are key to understanding the long‐term consequences of anthropogenic global change. Herbaria document such responses, and, with contemporary sampling, provide high‐resolution time‐series of plant evolutionary change. Characterizing genetic diversity is straightforward for model species with small genomes and a reference sequence. For nonmodel species—with small or large genomes—diversity is traditionally assessed using restriction‐enzyme‐based sequencing. However, age‐related DNA damage and fragmentation preclude the use of this approach for ancient herbarium DNA. Here, we combine reduced‐representation sequencing and hybridization‐capture to overcome this challenge and efficiently compare contemporary and historical specimens. Specifically, we describe how homemade DNA baits can be produced from reduced‐representation libraries of fresh samples, and used to efficiently enrich historical libraries for the same fraction of the genome to produce compatible sets of sequence data from both types of material. Applying this approach to both Arabidopsis thaliana and the nonmodel plant Cardamine bulbifera, we discovered polymorphisms de novo in an unbiased, reference‐free manner. We show that the recovered genetic variation recapitulates known genetic diversity in A. thaliana, and recovers geographical origin in both species and over time, independent of bait diversity. Hence, our method enables fast, cost‐efficient, large‐scale integration of contemporary and historical specimens for assessment of genome‐wide genetic trends over time, independent of genome size and presence of a reference genome.     

Evaluating multiplexed next‐generation sequencing as a method in palynology for mixed pollen samples

The identification of pollen plays an important role in ecology, palaeo‐climatology, honey quality control and other areas. Currently, expert knowledge and reference collections are essential to identify pollen origin through light microscopy. Pollen identification through molecular sequencing and DNA barcoding has been proposed as an alternative approach, but the assessment of mixed pollen samples originating from multiple plant species is still a tedious and error‐prone task. Next‐generation sequencing has been proposed to avoid this hindrance. In this study we assessed mixed pollen probes through next‐generation sequencing of amplicons from the highly variable, species‐specific internal transcribed spacer 2 region of nuclear ribosomal DNA. Further, we developed a bioinformatic workflow to analyse these high‐throughput data with a newly created reference database. To evaluate the feasibility, we compared results from classical identification based on light microscopy from the same samples with our sequencing results. We assessed in total 16 mixed pollen samples, 14 originated from honeybee colonies and two from solitary bee nests. The sequencing technique resulted in higher taxon richness (deeper assignments and more identified taxa) compared to light microscopy. Abundance estimations from sequencing data were significantly correlated with counted abundances through light microscopy. Simulation analyses of taxon specificity and sensitivity indicate that 96% of taxa present in the database are correctly identifiable at the genus level and 70% at the species level. Next‐generation sequencing thus presents a useful and efficient workflow to identify pollen at the genus and species level without requiring specialised palynological expert knowledge.     

The phylogenetic systematics of blue‐tailed skinks (Plestiodon) and the family Scincidae

Blue‐tailed skinks (genus Plestiodon) are a common component of the terrestrial herpetofauna throughout their range in eastern Eurasia and North and Middle America. Plestiodon species are also frequent subjects of ecological and evolutionary research, yet a comprehensive, well‐supported phylogenetic framework does not yet exist for this genus. We construct a comprehensive molecular phylogeny of Plestiodon using Bayesian phylogenetic analyses of a nine‐locus data set comprising 8308 base pairs of DNA, sampled from 38 of the 43 species in the genus. We evaluate potential gene tree/species tree discordance by conducting phylogenetic analyses of the concatenated and individual locus data sets, as well as employing coalescent‐based methods. Specifically, we address the placement of Plestiodon within the evolutionary tree of Scincidae, as well as the phylogenetic relationships between Plestiodon species, and their taxonomy. Given our sampling of major Scincidae lineages, we also re‐evaluate ‘deep’ relationships within the family, with the goal of resolving relationships that have been ambiguous in recent molecular phylogenetic analyses. We infer strong support for several scincid relationships, including a major clade of ‘scincines’ and the inter‐relationships of major Mediterranean and southern African genera. Although we could not estimate the precise phylogenetic affinities of Plestiodon with statistically significant support, we nonetheless infer significant support for its inclusion in a large ‘scincine’ clade exclusive of Acontinae, Lygosominae, Brachymeles, and Ophiomorus. Plestiodon comprises three major geographically cohesive clades. One of these clades is composed of mostly large‐bodied species inhabiting northern Indochina, south‐eastern China (including Taiwan), and the southern Ryukyu Islands of Japan. The second clade comprises species inhabiting central China (including Taiwan) and the entire Japanese archipelago. The third clade exclusively inhabits North and Middle America and the island of Bermuda. A vast majority of interspecific relationships are strongly supported in the concatenated data analysis, but there is nonetheless significant conflict amongst the individual gene trees. Coalescent‐based gene tree/species tree analyses indicate that incongruence amongst the nuclear loci may severely obscure the phylogenetic inter‐relationships of the primarily small‐bodied Plestiodon species that inhabit the central Mexican highlands. These same analyses do support the sister relationship between Plestiodon marginatus Hallowell, 1861 and Plestiodon stimpsonii (Thompson, 1912), and differ with the mitochondrial DNA analysis that supports Plestiodon elegans (Boulenger, 1887) + P. stimpsonii. Finally, because the existing Plestiodon taxonomy is a poor representation of evolutionary relationships, we replace the existing supraspecific taxonomy with one congruent with our phylogenetic results. © 2012 The Linnean Society of London, Zoological Journal of the Linnean Society, 2012, 165 , 163–189.     

Universal target‐enrichment baits for anthozoan (Cnidaria) phylogenomics: New approaches to long‐standing problems

Anthozoans (e.g., corals, anemones) are an ecologically important and diverse group of marine metazoans that occur from shallow to deep waters worldwide. However, our understanding of the evolutionary relationships among the ~7,500 species within this class is hindered by the lack of phylogenetically informative markers that can be reliably sequenced across a diversity of taxa. We designed and tested 16,306 RNA baits to capture 720 ultraconserved element loci and 1,071 exon loci. Library preparation and target enrichment were performed on 33 taxa from all orders within the class Anthozoa. Following Illumina sequencing and Trinity assembly, we recovered 1,774 of 1,791 targeted loci. The mean number of loci recovered from each species was 638 ± 222, with more loci recovered from octocorals (783 ± 138 loci) than hexacorals (475 ± 187 loci). Parsimony informative sites ranged from 26 to 49% for alignments at differing hierarchical taxonomic levels (e.g., Anthozoa, Octocorallia, Hexacorallia). The per cent of variable sites within each of three genera (Acropora, Alcyonium, and Sinularia) for which multiple species were sequenced ranged from 4.7% to 30%. Maximum‐likelihood analyses recovered highly resolved trees with topologies matching those supported by other studies, including the monophyly of the order Scleractinia. Our results demonstrate the utility of this target‐enrichment approach to resolve phylogenetic relationships from relatively old to recent divergences. Redesigning the baits with improved affinities to capture loci within each subclass will provide a valuable toolset to address systematic questions, further our understanding of the timing of diversifications and help resolve long‐standing controversial relationships in the class Anthozoa.     

Sequence capture phylogenomics of historical ethanol‐preserved museum specimens: Unlocking the rest of the vault

Natural history collections play a crucial role in biodiversity research, and museum specimens are increasingly being incorporated into modern genetics‐based studies. Sequence capture methods have proven incredibly useful for phylogenomics, providing the additional ability to sequence historical museum specimens with highly degraded DNA, which until recently have been deemed less valuable for genetic work. The successful sequencing of ultraconserved elements (UCEs) from historical museum specimens has been demonstrated on multiple tissue types including dried bird skins, formalin‐fixed squamates and pinned insects. However, no study has thoroughly demonstrated this approach for historical ethanol‐preserved museum specimens. Alongside sequencing of “fresh” specimens preserved in >95% ethanol and stored at ?80°C, we used extraction techniques specifically designed for degraded DNA coupled with sequence capture protocols to sequence UCEs from historical museum specimens preserved in 70%–80% ethanol and stored at room temperature, the standard for such ethanol‐preserved museum collections. Across 35 fresh and 15 historical museum samples of the arachnid order Opiliones, an average of 345 UCE loci were included in phylogenomic matrices, with museum samples ranging from six to 495 loci. We successfully demonstrate the inclusion of historical ethanol‐preserved museum specimens in modern sequence capture phylogenomic studies, show a high frequency of variant bases at the species and population levels, and from off‐target reads successfully recover multiple loci traditionally sequenced in multilocus studies including mitochondrial loci and nuclear rRNA loci. The methods detailed in this study will allow researchers to potentially acquire genetic data from millions of ethanol‐preserved museum specimens held in collections worldwide.     

Shotgun mitogenomics across body size classes in a local assemblage of tropical Diptera: Phylogeny,species diversity and mitochondrial abundance spectrum

Mitochondrial genomes can be assembled readily from shotgun‐sequenced DNA mixtures of mass‐trapped arthropods (“mitochondrial metagenomics”), speeding up the taxonomic characterization. Bulk sequencing was conducted on some 800 individuals of Diptera obtained by canopy fogging of a single tree in Borneo dominated by small (<1.5 mm) individuals. Specimens were split into five body size classes for DNA extraction, to equalize read numbers across specimens and to study how body size, a key ecological trait, interacts with species and phylogenetic diversity. Genome assembly produced 304 orthologous mitochondrial contigs presumed to each represent a different species. The small‐bodied fraction was the by far most species‐rich (187 contigs). Identification of contigs was through phylogenetic analysis together with 56 reference mitogenomes, which placed most of the Bornean community into seven clades of small‐bodied species, indicating phylogenetic conservation of body size. Mapping of shotgun reads against the mitogenomes showed wide ranges of read abundances within each size class. Ranked read abundance plots were largely log‐linear, indicating a uniformly filled abundance spectrum, especially for small‐bodied species. Small‐bodied species differed greatly from other size classes in neutral metacommunity parameters, exhibiting greater levels of immigration, besides greater total community size. We suggest that the established uses of mitochondrial metagenomics for analysis of species and phylogenetic diversity can be extended to parameterize recent theories of community ecology and biodiversity, and by focusing on the number mitochondria, rather than individuals, a new theoretical framework for analysis of mitochondrial abundance spectra can be developed that incorporates metabolic activity approximated by the count of mitochondria.     

Next‐generation sequencing workflow for assembly of nonmodel mitogenomes exemplified with North Pacific albatrosses (Phoebastria spp.)

Use of complete mitochondrial genomes (mitogenomes) can greatly increase the resolution achievable in phylogeographic and historical demographic studies. Using next‐generation sequencing methods, it is now feasible to efficiently sequence mitogenomes of large numbers of individuals once a reference mitogenome is available. However, assembling the initial mitogenomes of nonmodel organisms can present challenges, for example, in birds, where mtDNA is often subject to gene rearrangements and duplications. We developed a workflow based on Illumina paired‐end, whole‐genome shotgun sequencing, which we used to generate complete 19‐kilobase mitogenomes for each of three species of North Pacific albatross, a group of birds known to carry a tandem duplication. Although this duplication had been described previously, our procedure did not depend on this prior knowledge, nor did it require a closely related reference mitogenome (e.g. a mammalian mitogenome was sufficient). We employed an iterative process including de novo assembly, reference‐guided assembly and gap closing, which enabled us to detect duplications, determine gene order and identify sequence for primer positioning to resolve any mitogenome ambiguity (via minimal targeted Sanger sequencing). We present full mtDNA annotations, including 22 tRNAs, 2 rRNAs, 13 protein‐coding genes, a control region and a duplicated feature for all three species. Pairwise comparisons supported previous hypotheses regarding the phylogenetic relationships within this group and occurrence of a shared tandem duplication. The resulting mitogenome sequences will enable rapid, high‐throughput NGS mitogenome sequencing of North Pacific albatrosses via direct reference‐guided assembly. Moreover, our approach to assembling mitogenomes should be applicable to any taxon.     

A from‐benchtop‐to‐desktop workflow for validating HTS data and for taxonomic identification in diet metabarcoding studies

The main objective of this work was to develop and validate a robust and reliable “from‐benchtop‐to‐desktop” metabarcoding workflow to investigate the diet of invertebrate‐eaters. We applied our workflow to faecal DNA samples of an invertebrate‐eating fish species. A fragment of the cytochrome c oxidase I (COI) gene was amplified by combining two minibarcoding primer sets to maximize the taxonomic coverage. Amplicons were sequenced by an Illumina MiSeq platform. We developed a filtering approach based on a series of nonarbitrary thresholds established from control samples and from molecular replicates to address the elimination of cross‐contamination, PCR/sequencing errors and mistagging artefacts. This resulted in a conservative and informative metabarcoding data set. We developed a taxonomic assignment procedure that combines different approaches and that allowed the identification of ~75% of invertebrate COI variants to the species level. Moreover, based on the diversity of the variants, we introduced a semiquantitative statistic in our diet study, the minimum number of individuals, which is based on the number of distinct variants in each sample. The metabarcoding approach described in this article may guide future diet studies that aim to produce robust data sets associated with a fine and accurate identification of prey items.     

Sequencing historical specimens: successful preparation of small specimens with low amounts of degraded DNA

Despite advances that allow DNA sequencing of old museum specimens, sequencing small‐bodied, historical specimens can be challenging and unreliable as many contain only small amounts of fragmented DNA. Dependable methods to sequence such specimens are especially critical if the specimens are unique. We attempt to sequence small‐bodied (3–6 mm) historical specimens (including nomenclatural types) of beetles that have been housed, dried, in museums for 58–159 years, and for which few or no suitable replacement specimens exist. To better understand ideal approaches of sample preparation and produce preparation guidelines, we compared different library preparation protocols using low amounts of input DNA (1–10 ng). We also explored low‐cost optimizations designed to improve library preparation efficiency and sequencing success of historical specimens with minimal DNA, such as enzymatic repair of DNA. We report successful sample preparation and sequencing for all historical specimens despite our low‐input DNA approach. We provide a list of guidelines related to DNA repair, bead handling, reducing adapter dimers and library amplification. We present these guidelines to facilitate more economical use of valuable DNA and enable more consistent results in projects that aim to sequence challenging, irreplaceable historical specimens.     

High‐throughput sequencing of fecal DNA to identify insects consumed by wild Weddell's saddleback tamarins (Saguinus weddelli,Cebidae, Primates) in Bolivia

The genus Saguinus represents a successful radiation of over 20 species of small‐bodied New World monkeys. Studies of the tamarin diet indicate that insects and small vertebrates account for ～16–45% of total feeding and foraging time, and represent an important source of lipids, protein, and metabolizable energy. Although tamarins are reported to commonly consume large‐bodied insects such as grasshoppers and walking sticks (Orthoptera), little is known concerning the degree to which smaller or less easily identifiable arthropod prey comprises an important component of their diet. To better understand tamarin arthropod feeding behavior, fecal samples from 20 wild Bolivian saddleback tamarins (members of five groups) were collected over a 3 week period in June 2012, and analyzed for the presence of arthropod DNA. DNA was extracted using a Qiagen stool extraction kit, and universal insect primers were created and used to amplify a ～280 bp section of the COI mitochondrial gene. Amplicons were sequenced on the Roche 454 sequencing platform using high‐throughput sequencing techniques. An analysis of these samples indicated the presence of 43 taxa of arthropods including 10 orders, 15 families, and 12 identified genera. Many of these taxa had not been previously identified in the tamarin diet. These results highlight molecular analysis of fecal DNA as an important research tool for identifying anthropod feeding patterns in primates, and reveal broad diversity in the taxa, foraging microhabitats, and size of arthropods consumed by tamarin monkeys. Am J Phys Anthropol 156:474–481, 2015. © 2014 Wiley Periodicals, Inc.     

Phylogenetic analysis of cultivation‐resistant terrestrial cyanobacteria with massive sheaths (Stigonema spp. and Petalonema alatum,Nostocales, Cyanobacteria) using single‐cell and filament sequencing of environmental samples

Molecular assessment of a large portion of traditional cyanobacterial taxa has been hindered by the failure to isolate and grow them in culture. In this study, we developed an optimized protocol for single cell/filament isolation and 16S rRNA gene sequencing of terrestrial cyanobacteria with large mucilaginous sheaths, and applied it to determine the phylogenetic position of typical members of the genera Petalonema and Stigonema. A methodology based on a glass‐capillary isolation technique and a semi‐nested PCR protocol enabled reliable sequencing of the 16S rRNA gene from all samples analyzed. Ten samples covering seven species of Stigonema from Europe, North and Central America, and Hawaii, and the type species of Petalonema from Slovakia were sequenced. Contrary to some previous studies, which proposed a relationship with heteropolar nostocalean cyanobacteria, Petalonema appeared to belong to the family Scytonemataceae. Analysis of Stigonema specimens recovered a unique coherent phylogenetic cluster, substantially broadening our knowledge of the molecular diversity within this genus. Neither the uni‐ to biseriate species nor the multiseriate species formed monophyletic subclusters within the genus. Typical multiseriate species of Stigonema clustered in a phylogenetic branch derived from uni‐ to biseriate S. ocellatum Thuret ex Bornet & Flahault in our analysis, suggesting that species with more complex thalli may have evolved from the more simple ones. We propose the technique tested in this study as a promising tool for a future revision of the molecular taxonomy in cyanobacteria.     

Phylogenetic analysis of RAD‐seq data: examining the influence of gene genealogy conflict on analysis of concatenated data

One of the major issues in phylogenetic analysis is that gene genealogies from different gene regions may not reflect the true species tree or history of speciation. This has led to considerable debate about whether concatenation of loci is the best approach for phylogenetic analysis. The application of Next‐generation sequencing techniques such as RAD‐seq generates thousands of relatively short sequence reads from across the genomes of the sampled taxa. These data sets are typically concatenated for phylogenetic analysis leading to data sets that contain millions of base pairs per taxon. The influence of gene region conflict among so many loci in determining the phylogenetic relationships among taxa is unclear. We simulated RAD‐seq data by sampling 100 and 500 base pairs from alignments of over 6000 coding regions that each produce one of three highly supported alternative phylogenies of seven species of Drosophila. We conducted phylogenetic analyses on different sets of these regions to vary the sampling of loci with alternative gene trees to examine the effect on detecting the species tree. Irrespective of sequence length sampled per region and which subset of regions was used, phylogenetic analyses of the concatenated data always recovered the species tree. The results suggest that concatenated alignments of Next‐generation data that consist of many short sequences are robust to gene tree/species tree conflict when the goal is to determine the phylogenetic relationships among taxa.     

Phylogenomic species delimitation and host‐symbiont coevolution in the fungus‐farming ant genus Sericomyrmex Mayr (Hymenoptera: Formicidae): ultraconserved elements (UCEs) resolve a recent radiation

Ants in the Neotropical genus Sericomyrmex Mayr cultivate fungi for food. Both ants and fungi are obligate, coevolved symbionts. The taxonomy of Sericomyrmex is problematic because the morphology of the worker caste is generally homogeneous across all of the species within the genus, species limits are vague, and the relationships between them are unknown. We used ultraconserved elements (UCEs) as genome‐scale markers to reconstruct evolutionary history and to infer species boundaries in Sericomyrmex. We recovered an average of ～990 UCE loci for 88 Sericomyrmex samples from across the geographical range of the genus as well as for five outgroup taxa. Using maximum likelihood and species‐tree approaches, we recovered nearly identical topologies across datasets with 50–95% matrix completeness. We identify nine species‐level lineages in Sericomyrmex, including two new species. This is less than the previously described 19 species, even accounting for two species for which we had no UCE samples, which brings the total number of Sericomyrmex species to 11. Divergence‐dating analyses recovered 4.3 Ma as the crown‐group age estimates for Sericomyrmex, indicating a recent, rapid radiation. We also sequenced mitochondrial cytochrome oxidase subunit I (COI) for 125 specimens. Resolution and support for clades in our COI phylogeny are weak, indicating that COI is not an appropriate species‐delimitation tool. However, taxa within species consistently cluster together, suggesting that COI is useful as a species identification (‘DNA barcoding’) tool. We also sequenced internal transcribed spacer (ITS) and large subunit (LSU) for 32 Sericomyrmex fungal cultivars. The fungal phylogeny confirms that Sericomyrmex fungi are generalized higher‐attine cultivars, interspersed with Trachymyrmex‐associated fungal species, indicating cultivar sharing and horizontal transfer between these two genera. Our results indicate that UCEs offer immense potential for delimiting and resolving relationships of problematic, recently diverged species.     

Phylogenomics using formalin‐fixed and 100+ year‐old intractable natural history specimens

Museum specimens provide a wealth of information to biologists, but obtaining genetic data from formalin‐fixed and fluid‐preserved specimens remains challenging. While DNA sequences have been recovered from such specimens, most approaches are time‐consuming and produce low data quality and quantity. Here, we use a modified DNA extraction protocol combined with high‐throughput sequencing to recover DNA from formalin‐fixed and fluid‐preserved snakes that were collected over a century ago and for which little or no modern genetic materials exist in public collections. We successfully extracted DNA and sequenced ultraconserved elements ( = 2318 loci) from 10 fluid‐preserved snakes and included them in a phylogeny with modern samples. This phylogeny demonstrates the general use of such specimens in phylogenomic studies and provides evidence for the placement of enigmatic snakes, such as the rare and never‐before sequenced Indian Xylophis stenorhynchus. Our study emphasizes the relevance of museum collections in modern research and simultaneously provides a protocol that may prove useful for specimens that have been previously intractable for DNA sequencing.     

Disparities in second‐generation DNA metabarcoding results exposed with accessible and repeatable workflows

Different second‐generation sequencing technologies may have taxon‐specific biases when DNA metabarcoding prey in predator faeces. Our major objective was to examine differences in prey recovery from bat guano across two different sequencing workflows using the same faecal DNA extracts. We compared results between the Ion Torrent PGM and the Illumina MiSeq with similar library preparations and the same analysis pipeline. We focus on repeatability and provide an R Notebook in an effort towards transparency for future methodological improvements. Full documentation of each step enhances the accessibility of our analysis pipeline. We tagged DNA from insectivorous bat faecal samples, targeted the arthropod cytochrome c oxidase I minibarcode region and sequenced the product on both second‐generation sequencing platforms. We developed an analysis pipeline with a high operational taxonomic unit (OTU) clustering threshold (i.e., ≥98.5%) followed by copy number filtering to avoid merging rare but genetically similar prey into the same OTUs. With this workflow, we detected 297 unique prey taxa, of which 74% were identified at the species level. Of these, 104 (35%) prey OTUs were detected by both platforms, 176 (59%) OTUs were detected by the Illumina MiSeq system only, and 17 (6%) OTUs were detected using the Ion Torrent system only. Costs were similar between platforms but the Illumina MiSeq recovered six times more reads and four additional insect orders than did Ion Torrent. The considerations we outline are particularly important for long‐term ecological monitoring; a more standardized approach will facilitate comparisons between studies and allow faster recognition of changes within ecological communities.     

Major radiation of cheilostome bryozoans: Triggered by the evolution of a new larval type?

A macroevolutionary model is developed to account for the “adaptive radiation”; of cheilostome bryozoans that commenced in the Cenomanian after a long phase of low diversity. Living cheilostome species possess one of two types of larvae; planktotrophic (cyphonautes) larvae of relatively long duration, and brooded non‐planktotrophic (coronate) larvae of short duration. Planktotrophic larvae characterize the paraphyletic “malacostegans”; from which “advanced”; cheilostomes with non‐planktotrophic larvae are thought to have evolved monophyletically. Research on other marine invertebrates suggests that gene flow within and between populations is likely to be poorer in species having non‐planktotrophic larvae, and hence the frequency of allopatric and quasi‐sympatric speciation may be greater. Skeletal evidence of larval brooding in the cheilostomes first appears in the late Albian, immediately before their adaptive radiation, and the evolution of non‐planktotrophy with associated increase in speciation rate is proposed to have triggered this radiation.     

Long‐read metabarcoding of the eukaryotic rDNA operon to phylogenetically and taxonomically resolve environmental diversity

High‐throughput DNA metabarcoding of amplicon sizes below 500 bp has revolutionized the analysis of environmental microbial diversity. However, these short regions contain limited phylogenetic signal, which makes it impractical to use environmental DNA in full phylogenetic inferences. This lesser phylogenetic resolution of short amplicons may be overcome by new long‐read sequencing technologies. To test this idea, we amplified soil DNA and used PacBio Circular Consensus Sequencing (CCS) to obtain an ~4500‐bp region spanning most of the eukaryotic small subunit (18S) and large subunit (28S) ribosomal DNA genes. We first treated the CCS reads with a novel curation workflow, generating 650 high‐quality operational taxonomic units (OTUs) containing the physically linked 18S and 28S regions. To assign taxonomy to these OTUs, we developed a phylogeny‐aware approach based on the 18S region that showed greater accuracy and sensitivity than similarity‐based methods. The taxonomically annotated OTUs were then combined with available 18S and 28S reference sequences to infer a well‐resolved phylogeny spanning all major groups of eukaryotes, allowing us to accurately derive the evolutionary origin of environmental diversity. A total of 1,019 sequences were included, of which a majority (58%) corresponded to the new long environmental OTUs. The long reads also allowed us to directly investigate the relationships among environmental sequences themselves, which represents a key advantage over the placement of short reads on a reference phylogeny. Together, our results show that long amplicons can be treated in a full phylogenetic framework to provide greater taxonomic resolution and a robust evolutionary perspective to environmental DNA.