Mitochondrial genomic data are effective at placing mycoheterotrophic lineages in plant phylogeny

Fully mycoheterotrophic plants can be difficult to place in plant phylogeny due to elevated substitution rates associated with photosynthesis loss. This potentially limits the effectiveness of downstream analyses of mycoheterotrophy that depend on accurate phylogenetic inference. Although mitochondrial genomic data sets are rarely used in plant phylogenetics, theory predicts that they should be resilient to long-branch artefacts, thanks to their generally slow evolution, coupled with limited rate elevation in heterotrophs. We examined the utility of mitochondrial genomes for resolving contentious higher-order placements of mycoheterotrophic lineages in two test cases: monocots (focusing on Dioscoreales) and Ericaceae. We find Thismiaceae to be distantly related to Burmanniaceae in the monocot order Dioscoreales, conflicting with current classification schemes based on few gene data sets. We confirm that the unusual Afrothismia is related to Taccaceae-Thismiaceae, with a corresponding independent loss of photosynthesis. In Ericaceae we recovered the first well supported relationships among its five major lineages: mycoheterotrophic Ericaceae are not monophyletic, as pyroloids are inferred to be sister to core Ericaceae, and monotropoids to arbutoids. Genes recovered from mitochondrial genomes collectively resolved previously ambiguous mycoheterotroph higher-order relationships. We propose that mitochondrial genomic data should be considered in standardised gene panels for inferring overall plant phylogeny.

Message in a Bottle-Metabarcoding enables biodiversity comparisons across ecoregions

Background

Traditional biomonitoring approaches have delivered a basic understanding of biodiversity, but they cannot support the large-scale assessments required to manage and protect entire ecosystems. This study used DNA metabarcoding to assess spatial and temporal variation in species richness and diversity in arthropod communities from 52 protected areas spanning 3 Canadian ecoregions.

Results

This study revealed the presence of 26,263 arthropod species in the 3 ecoregions and indicated that at least another 3,000–5,000 await detection. Results further demonstrate that communities are more similar within than between ecoregions, even after controlling for geographical distance. Overall α-diversity declined from east to west, reflecting a gradient in habitat disturbance. Shifts in species composition were high at every site, with turnover greater than nestedness, suggesting the presence of many transient species.

Conclusions

Differences in species composition among their arthropod communities confirm that ecoregions are a useful synoptic for biogeographic patterns and for structuring conservation efforts. The present results also demonstrate that metabarcoding enables large-scale monitoring of shifts in species composition, making it possible to move beyond the biomass measurements that have been the key metric used in prior efforts to track change in arthropod communities.

Assessment of stream macroinvertebrate communities with eDNA is not congruent with tissue-based metabarcoding

Freshwater biomonitoring programmes routinely sample aquatic macroinvertebrates. These samples are time-consuming to collect, as well as challenging and costly to identify reliably genus or species. Environmental DNA (eDNA) metabarcoding has emerged as a surrogate to traditional collection techniques and has been used in whole-community approaches across several taxa and ecosystems. However, the usefulness of eDNA-based detection of freshwater macroinvertebrates has not been extensively explored. Few studies have directly compared bulk sample and eDNA metabarcoding at a local scale to assess how effective each method is at characterizing aquatic macroinvertebrate communities. Here, we collected both eDNA and kicknet samples at the same sample transect locations across nine different streams in southern Ontario, Canada. We observed minimal overlap in community composition between these paired samples. Bulk tissue metabarcoding resulted in a greater proportion of sequences belonging to metazoan taxa (over 99%) than eDNA (12%) and had higher OTU richness for macroinvertebrate taxa. We suggest that degenerate primers are not effective for eDNA metabarcoding due to the high degree of nontarget amplification and subsequently low yield of target DNA. While both bulk sample and eDNA metabarcoding had the power to detect differences between stream communities, eDNA did not represent local communities. Bulk tissue metabarcoding thus provides a more accurate representation of local stream macroinvertebrate communities and is the preferred method if smaller-scale spatial resolution is an important factor in data analyses.

Metabarcoding a diverse arthropod mock community

Although DNA metabarcoding is an attractive approach for monitoring biodiversity, it is often difficult to detect all the species present in a bulk sample. In particular, sequence recovery for a given species depends on its biomass and mitome copy number as well as the primer set employed for PCR. To examine these variables, we constructed a mock community of terrestrial arthropods comprised of 374 species. We used this community to examine how species recovery was impacted when amplicon pools were constructed in four ways. The first two protocols involved the construction of bulk DNA extracts from different body segments (Bulk Abdomen, Bulk Leg). The other protocols involved the production of DNA extracts from single legs which were then merged prior to PCR (Composite Leg) or PCR‐amplified separately (Single Leg) and then pooled. The amplicons generated by these four treatments were then sequenced on three platforms (Illumina MiSeq, Ion Torrent PGM and Ion Torrent S5). The choice of sequencing platform did not substantially influence species recovery, although the Miseq delivered the highest sequence quality. As expected, species recovery was most efficient from the Single Leg treatment because amplicon abundance varied little among taxa. Among the three treatments where PCR occurred after pooling, the Bulk Abdomen treatment produced a more uniform read abundance than the Bulk Leg or Composite Leg treatment. Primer choice also influenced species recovery and evenness. Our results reveal how variation in protocols can have substantial impacts on perceived diversity unless sequencing coverage is sufficient to reach an asymptote.

Finding the pond through the weeds: eDNA reveals underestimated diversity of pondweeds

PREMISE OF THE STUDY

The detection of environmental DNA (eDNA) using high-throughput sequencing has rapidly emerged as a method to detect organisms from environmental samples. However, eDNA studies of aquatic biomes have focused on surveillance of animal species with less emphasis on plants. Pondweeds are important bioindicators of freshwater ecosystems, although their diversity is underestimated due to difficulties in morphological identification and monitoring.

METHODS

A protocol was developed to detect pondweeds in water samples using atpB-rbcL and ITS2 markers. The water samples were collected from the Grand River within the rare Charitable Research Reserve, Ontario (RARE). Short fragments were amplified using primers targeting pondweeds, sequenced on an Ion Torrent Personal Genome Machine, and assigned to the taxonomy using a local DNA reference library and GenBank.

RESULTS

We detected two species earlier documented at the experimental site during ecological surveys (Potamogeton crispus and Stuckenia pectinata) and three species new to the RARE checklist (P. foliosus, S. filiformis, and Zannichellia palustris).

DISCUSSION

Our targeted approach to track the species composition of pondweeds in freshwater ecosystems revealed underestimation of their diversity. This result suggests that eDNA is an effective tool for monitoring plant diversity in aquatic habitats.

A Sequel to Sanger: amplicon sequencing that scales

BACKGROUND

Although high-throughput sequencers (HTS) have largely displaced their Sanger counterparts, the short read lengths and high error rates of most platforms constrain their utility for amplicon sequencing. The present study tests the capacity of single molecule, real-time (SMRT) sequencing implemented on the SEQUEL platform to overcome these limitations, employing 658 bp amplicons of the mitochondrial cytochrome c oxidase I gene as a model system.

RESULTS

By examining templates from more than 5000 species and 20,000 specimens, the performance of SMRT sequencing was tested with amplicons showing wide variation in GC composition and varied sequence attributes. SMRT and Sanger sequences were very similar, but SMRT sequencing provided more complete coverage, especially for amplicons with homopolymer tracts. Because it can characterize amplicon pools from 10,000 DNA extracts in a single run, the SEQUEL can reduce greatly reduce sequencing costs in comparison to first (Sanger) and second generation platforms (Illumina, Ion).

CONCLUSIONS

SMRT analysis generates high-fidelity sequences from amplicons with varying GC content and is resilient to homopolymer tracts. Analytical costs are low, substantially less than those for first or second generation sequencers. When implemented on the SEQUEL platform, SMRT analysis enables massive amplicon characterization because each instrument can recover sequences from more than 5 million DNA extracts a year.

Using herbarium-derived DNAs to assemble a large-scale DNA barcode library for the vascular plants of Canada

PREMISE OF THE STUDY

Constructing complete, accurate plant DNA barcode reference libraries can be logistically challenging for large-scale floras. Here we demonstrate the promise and challenges of using herbarium collections for building a DNA barcode reference library for the vascular plant flora of Canada.

METHODS

Our study examined 20,816 specimens representing 5076 of 5190 vascular plant species in Canada (98%). For 98% of the specimens, at least one of the DNA barcode regions was recovered from the plastid loci rbcL and matK and from the nuclear ITS2 region. We used beta regression to quantify the effects of age, type of preservation, and taxonomic affiliation (family) on DNA sequence recovery.

RESULTS

Specimen age and method of preservation had significant effects on sequence recovery for all markers, but influenced some families more (e.g., Boraginaceae) than others (e.g., Asteraceae).

DISCUSSION

Our DNA barcode library represents an unparalleled resource for metagenomic and ecological genetic research working on temperate and arctic biomes. An observed decline in sequence recovery with specimen age may be associated with poor primer matches, intragenomic variation (for ITS2), or inhibitory secondary compounds in some taxa.

On the brink: the highly reduced plastomes of nonphotosynthetic Ericaceae

Ericaceae (the heather family) is a large and diverse group of plants that forms elaborate symbiotic relationships with mycorrhizal fungi, and includes several nonphotosynthetic lineages. Using an extensive sample of fully mycoheterotrophic (MH) species, we explored inter- and intraspecific variation as well as selective constraints acting on the plastomes of these unusual plants. The plastomes of seven MH genera were analysed in a phylogenetic context with two geographically disparate individuals sequenced for Allotropa, Monotropa, and Pityopus. The plastomes of nonphotosynthetic Ericaceae are highly reduced in size (c. 33-41 kbp) and content, having lost all photosynthesis-related genes, and are reduced to encoding housekeeping genes as well as a protease subunit (clpP)-like and acetyl-CoA carboxylase subunit D (accD)-like open reading frames. Despite an increase in the rate of their nucleotide substitutions, the remaining protein-coding genes are typically under purifying selection in full MHs. We also identified ribosomal proteins under relaxed or neutral selection. These plastomes also exhibit striking structural rearrangements. Intraspecific variation within MH Ericaceae ranges from a few differences (Allotropa) to extensive population divergences (Monotropa, Hypopitys), which indicates that cryptic speciation may be occurring in several lineages. The pattern of gene loss within fully MH Ericaceae plastomes suggests an advanced state of degradation.

Escaping introns in COI through cDNA barcoding of mushrooms: Pleurotus as a test case

DNA barcoding involves the use of one or more short, standardized DNA fragments for the rapid identification of species. A 648-bp segment near the 5' terminus of the mitochondrial cytochrome c oxidase subunit I (COI) gene has been adopted as the universal DNA barcode for members of the animal kingdom, but its utility in mushrooms is complicated by the frequent occurrence of large introns. As a consequence, ITS has been adopted as the standard DNA barcode marker for mushrooms despite several shortcomings. This study employed newly designed primers coupled with cDNA analysis to examine COI sequence diversity in six species of Pleurotus and compared these results with those for ITS. The ability of the COI gene to discriminate six species of Pleurotus, the commonly cultivated oyster mushroom, was examined by analysis of cDNA. The amplification success, sequence variation within and among species, and the ability to design effective primers was tested. We compared ITS sequences to their COI cDNA counterparts for all isolates. ITS discriminated between all six species, but some sequence results were uninterpretable, because of length variation among ITS copies. By comparison, a complete COI sequences were recovered from all but three individuals of Pleurotus giganteus where only the 5' region was obtained. The COI sequences permitted the resolution of all species when partial data was excluded for P. giganteus. Our results suggest that COI can be a useful barcode marker for mushrooms when cDNA analysis is adopted, permitting identifications in cases where ITS cannot be recovered or where it offers higher resolution when fresh tissue is. The suitability of this approach remains to be confirmed for other mushrooms.

Testing the Efficacy of DNA Barcodes for Identifying the Vascular Plants of Canada

Their relatively slow rates of molecular evolution, as well as frequent exposure to hybridization and introgression, often make it difficult to discriminate species of vascular plants with the standard barcode markers (rbcL, matK, ITS2). Previous studies have examined these constraints in narrow geographic or taxonomic contexts, but the present investigation expands analysis to consider the performance of these gene regions in discriminating the species in local floras at sites across Canada. To test identification success, we employed a DNA barcode reference library with sequence records for 96% of the 5108 vascular plant species known from Canada, but coverage varied from 94% for rbcL to 60% for ITS2 and 39% for matK. Using plant lists from 27 national parks and one scientific reserve, we tested the efficacy of DNA barcodes in identifying the plants in simulated species assemblages from six biogeographic regions of Canada using BLAST and mothur. Mean pairwise distance (MPD) and mean nearest taxon distance (MNTD) were strong predictors of barcode performance for different plant families and genera, and both metrics supported ITS2 as possessing the highest genetic diversity. All three genes performed strongly in assigning the taxa present in local floras to the correct genus with values ranging from 91% for rbcL to 97% for ITS2 and 98% for matK. However, matK delivered the highest species discrimination (~81%) followed by ITS2 (~72%) and rbcL (~44%). Despite the low number of plant taxa in the Canadian Arctic, DNA barcodes had the least success in discriminating species from this biogeographic region with resolution ranging from 36% with rbcL to 69% with matK. Species resolution was higher in the other settings, peaking in the Woodland region at 52% for rbcL and 87% for matK. Our results indicate that DNA barcoding is very effective in identifying Canadian plants to a genus, and that it performs well in discriminating species in regions where floristic diversity is highest.

Authentication of Herbal Supplements Using Next-Generation Sequencing

BACKGROUND

DNA-based testing has been gaining acceptance as a tool for authentication of a wide range of food products; however, its applicability for testing of herbal supplements remains contentious.

METHODS

We utilized Sanger and Next-Generation Sequencing (NGS) for taxonomic authentication of fifteen herbal supplements representing three different producers from five medicinal plants: Echinacea purpurea, Valeriana officinalis, Ginkgo biloba, Hypericum perforatum and Trigonella foenum-graecum. Experimental design included three modifications of DNA extraction, two lysate dilutions, Internal Amplification Control, and multiple negative controls to exclude background contamination. Ginkgo supplements were also analyzed using HPLC-MS for the presence of active medicinal components.

RESULTS

All supplements yielded DNA from multiple species, rendering Sanger sequencing results for rbcL and ITS2 regions either uninterpretable or non-reproducible between the experimental replicates. Overall, DNA from the manufacturer-listed medicinal plants was successfully detected in seven out of eight dry herb form supplements; however, low or poor DNA recovery due to degradation was observed in most plant extracts (none detected by Sanger; three out of seven-by NGS). NGS also revealed a diverse community of fungi, known to be associated with live plant material and/or the fermentation process used in the production of plant extracts. HPLC-MS testing demonstrated that Ginkgo supplements with degraded DNA contained ten key medicinal components.

CONCLUSION

Quality control of herbal supplements should utilize a synergetic approach targeting both DNA and bioactive components, especially for standardized extracts with degraded DNA. The NGS workflow developed in this study enables reliable detection of plant and fungal DNA and can be utilized by manufacturers for quality assurance of raw plant materials, contamination control during the production process, and the final product. Interpretation of results should involve an interdisciplinary approach taking into account the processes involved in production of herbal supplements, as well as biocomplexity of plant-plant and plant-fungal biological interactions.