Critical considerations for communicating environmental DNA science.

The economic and methodological efficiencies of environmental DNA (eDNA) based survey approaches provide an unprecedented opportunity to assess and monitor aquatic environments. However, instances of inadequate communication from the scientific community about confidence levels, knowledge gaps, reliability, and appropriate parameters of eDNA-based methods have hindered their uptake in environmental monitoring programs and, in some cases, has created misperceptions or doubts in the management community. To help remedy this situation, scientists convened a session at the Second National Marine eDNA Workshop to discuss strategies for improving communications with managers. These include articulating the readiness of different eDNA applications, highlighting the strengths and limitations of eDNA tools for various applications or use cases, communicating uncertainties associated with specified uses transparently, and avoiding the exaggeration of exploratory and preliminary findings. Several key messages regarding implementation, limitations, and relationship to existing methods were prioritized. To be inclusive of the diverse managers, practitioners, and researchers, we and the other workshop participants propose the development of communication workflow plans, using RACI (Responsible, Accountable, Consulted, Informed) charts to clarify the roles of all pertinent individuals and parties and to minimize the chance for miscommunications. We also propose developing decision support tools such as Structured Decision-Making (SDM) to help balance the benefits of eDNA sampling with the inherent uncertainty, and developing an eDNA readiness scale to articulate the technological readiness of eDNA approaches for specific applications. These strategies will increase clarity and consistency regarding our understanding of the utility of eDNA-based methods, improve transparency, foster a common vision for confidently applying eDNA approaches, and enhance their benefit to the monitoring and assessment community.

Contrasting outcomes of genome reduction in mikrocytids and microsporidians.

BACKGROUND: Intracellular symbionts often undergo genome reduction, losing both coding and non-coding DNA in a process that ultimately produces small, gene-dense genomes with few genes. Among eukaryotes, an extreme example is found in microsporidians, which are anaerobic, obligate intracellular parasites related to fungi that have the smallest nuclear genomes known (except for the relic nucleomorphs of some secondary plastids). Mikrocytids are superficially similar to microsporidians: they are also small, reduced, obligate parasites; however, as they belong to a very different branch of the tree of eukaryotes, the rhizarians, such similarities must have evolved in parallel. Since little genomic data are available from mikrocytids, we assembled a draft genome of the type species, Mikrocytos mackini, and compared the genomic architecture and content of microsporidians and mikrocytids to identify common characteristics of reduction and possible convergent evolution. RESULTS: At the coarsest level, the genome of M. mackini does not exhibit signs of extreme genome reduction; at 49.7 Mbp with 14,372 genes, the assembly is much larger and gene-rich than those of microsporidians. However, much of the genomic sequence and most (8075) of the protein-coding genes code for transposons, and may not contribute much of functional relevance to the parasite. Indeed, the energy and carbon metabolism of M. mackini share several similarities with those of microsporidians. Overall, the predicted proteome involved in cellular functions is quite reduced and gene sequences are extremely divergent. Microsporidians and mikrocytids also share highly reduced spliceosomes that have retained a strikingly similar subset of proteins despite having reduced independently. In contrast, the spliceosomal introns in mikrocytids are very different from those of microsporidians in that they are numerous, conserved in sequence, and constrained to an exceptionally narrow size range (all 16 or 17 nucleotides long) at the shortest extreme of known intron lengths. CONCLUSIONS: Nuclear genome reduction has taken place many times and has proceeded along different routes in different lineages. Mikrocytids show a mix of similarities and differences with other extreme cases, including uncoupling the actual size of a genome with its functional reduction.

Trade-offs between reducing complex terminology and producing accurate interpretations from environmental DNA: Comment on "Environmental DNA: What's behind the term?" by Pawlowski et al., (2020).

In a recent paper, "Environmental DNA: What's behind the term? Clarifying the terminology and recommendations for its future use in biomonitoring," Pawlowski et al. argue that the term eDNA should be used to refer to the pool of DNA isolated from environmental samples, as opposed to only extra-organismal DNA from macro-organisms. We agree with this view. However, we are concerned that their proposed two-level terminology specifying sampling environment and targeted taxa is overly simplistic and might hinder rather than improve clear communication about environmental DNA and its use in biomonitoring. This terminology is based on categories that are often difficult to assign and uninformative, and it overlooks a fundamental distinction within eDNA: the type of DNA (organismal or extra-organismal) from which ecological interpretations are derived.

Evaluation of histopathology, PCR, and qPCR to detect Mikrocytos mackini in oysters Crassostrea gigas using Bayesian latent class analysis.

Latent class analysis (LCA) is a common method to evaluate the diagnostic sensitivity (DSe) and specificity (DSp) for pathogen detection assays in the absence of a perfect reference standard. Here we used LCA to evaluate the diagnostic accuracy of 3 tests for the detection of Mikrocytos mackini in Pacific oysters Crassostrea gigas: conventional polymerase chain reaction (PCR), real-time quantitative PCR (qPCR), and histopathology. A total of 802 Pacific oysters collected over 12 sampling events from 9 locations were assessed. Preliminary investigations indicated that standard LCA assumptions of test independence and constant detection accuracy across locations were likely unrealistic. This was mitigated by restructuring the LCA in a Bayesian framework to include test-derived knowledge about pathogen prevalence and load for categorizing populations into 2 classes of infection severity (low or high) and assessing separate DSe and DSp estimates for each class. Median DSp estimates were high (>96%) for all 3 tests in both population classes. DSe estimates varied between tests and population classes but were consistently highest for qPCR (87-99%) and lowest for histopathology (21-51%). Acknowledging that detection of M. mackini may be fitted to multiple diagnostic and management purposes, qPCR had the highest DSe while maintaining similar DSp to both conventional PCR and histopathology and thus is generally well-suited to most applications.

Biotic signals associated with benthic impacts of salmon farms from eDNA metabarcoding of sediments.

Environmental DNA (eDNA) metabarcoding can rapidly characterize the composition and diversity of benthic communities, thus it has high potential utility for routine assessments of benthic impacts of marine finfish farming. In this study, 126 sediment grab samples from 42 stations were collected at six salmon farms in British Columbia, Canada. Benthic community changes were assessed by both eDNA metabarcoding of metazoans and macrofaunal polychaete surveys. The latter was done by analysing 11,466 individuals using a combination of morphology-based taxonomy and DNA barcoding. Study objectives were to: (i) compare biotic signals associated with benthic impacts of salmon farming in the two data sources, and (ii) identify potential eDNA indicators to facilitate monitoring in Canada. Alpha diversity parameters were consistently reduced near fish cage edge and negatively correlated with pore-water sulphide concentration, with coefficients ranging from -0.62 to -0.48. Although Polychaeta are a common indicator group, the negative correlation with pore-water sulphide concentration was much stronger for Nematoda OTU richness (correlation coefficient: -0.86) than for Polychaeta (correlation coefficient: -0.38). Presence/absence of Capitella generally agreed well between the two methods despite that they differed in the volume of sediments sampled and the molecular marker used. Multiple approaches were used to identify OTUs related to organic enrichment statuses. We demonstrate that eDNA metabarcoding generates biotic signals that could be leveraged for environmental assessment of benthic impacts of fish farms in multiple ways: both alpha diversity and Nematoda OTU richness could be used to assess the spatial extent of impact, and OTUs related to organic enrichment could be used to develop local biotic indices.

The risks of using molecular biodiversity data for incidental detection of species of concern.

Incidental detection of species of concern (e.g., invasive species, pathogens, threatened and endangered species) during biodiversity assessments based on high-throughput DNA sequencing holds significant risks in the absence of rigorous, fit-for-purpose data quality and reporting standards. Molecular biodiversity data are predominantly collected for ecological studies and thus are generated to common quality assurance standards. However, the detection of certain species of concern in these data would likely elicit interest from end users working in biosecurity or other surveillance contexts (e.g., pathogen detection in health-related fields), for which more stringent quality control standards are essential to ensure that data are suitable for informing decision-making and can withstand legal or political challenges. We suggest here that data quality and reporting criteria are urgently needed to enable clear identification of those studies that may be appropriately applied to surveillance contexts. In the interim, more pointed disclaimers on uncertainties associated with the detection and identification of species of concern may be warranted in published studies. This is not only to ensure the utility of molecular biodiversity data for consumers, but also to protect data generators from uncritical and potentially ill-advised application of their science in decision-making.

A new approach to molecular biosurveillance of invasive species using DNA metabarcoding.

Non-indigenous species (NIS) reach every corner of the world, at times wreaking havoc on ecosystems and costing the global economy billions of dollars. A rapid and accurate biosurveillance tool tailored to a particular biogeographic region is needed to detect NIS when they are first introduced into an area as traditional detection methods are expensive and require specialized expertise. Metabarcoding of environmental and community DNA meets those biosurveillance requirements; a novel tool tailored to the Northwest Pacific Ocean is presented here using an approach that could revolutionize early detection of NIS. Eight newly designed genetic markers for multiple gene regions were implemented to meet the stringent taxonomic requirements for the detection of NIS across four major marine phyla. The tool was considered highly successful because it identified 12 known NIS in the study area and a further seven species representing potential new records. Overall community composition detected here was statistically different between substrate types; zooplankton sampling accounted for significantly higher species richness than filtered sea water in most cases, but this was dominated by mollusk and arthropod species. Both substrate types sampled were required to identify the wide taxonomic breadth of known NIS in the study area. Intensive sampling is known to be paramount for the detection of rare species, including new incursions of NIS, thus it is recommended to include diverse DNA sampling protocols based on species' life-history characteristics for broad detection capacity. Application of a metabarcoding-based molecular biosurveillance tool optimized for biogeographic regions enables rapid and accurate early detection across a wide taxonomic range to allow quick implementation of eradication or control efforts and potentially mitigate some of the devastating effects of NIS worldwide.

Responses of foraminifera communities to aquaculture-derived organic enrichment as revealed by environmental DNA metabarcoding.

Current monitoring methods to assess benthic impacts of marine finfish aquaculture are based on complex biological indices and/or geochemistry data. The former requires benthic macrofauna morpho-taxonomic characterization that is time- and cost-intensive, while the latter provides rapid assessment of the organic enrichment status of sediments but does not directly measure biotic impacts. In this study, sediment samples were collected from seven stations at six salmon farms in British Columbia, Canada, and analyzed for geochemical parameters and by eDNA metabarcoding to investigate linkages between geochemistry and foraminifera. Sediment texture across farm sites ranged from sand to silty loam, while the maximum sediment pore-water sulphide concentration at each site ranged from 1,000 to 13,000 µM. Foraminifera alpha diversity generally increased with distance from cage edge. Adonis analyses revealed that farm site explained the most variation in foraminifera community, followed by sediment type, enrichment status, and distance from cage edge. Farm-specific responses were observed in diversity analyses, taxonomic difference analyses, and correlation analyses. Results demonstrated that species diversity and composition of foraminifera characterized by eDNA metabarcoding generated signals consistent with benthic biodiversity being impacted by finfish farming activities. This substantiates the validity of eDNA metabarcoding for augmenting current approaches to benthic impact assessments by providing more cost-effective and practicable biotic measures than traditional morpho-taxonomy. To capitalize on this potential, further work is needed to design a new nomogram that combines eDNA metabarcoding data and geochemistry data to enable accurate monitoring of benthic impacts of fish farming in a time- and cost-efficient way.

Metabarcoding using multiplexed markers increases species detection in complex zooplankton communities.

Metabarcoding combines DNA barcoding with high-throughput sequencing, often using one genetic marker to understand complex and taxonomically diverse samples. However, species-level identification depends heavily on the choice of marker and the selected primer pair, often with a trade-off between successful species amplification and taxonomic resolution. We present a versatile metabarcoding protocol for biomonitoring that involves the use of two barcode markers (COI and 18S) and four primer pairs in a single high-throughput sequencing run, via sample multiplexing. We validate the protocol using a series of 24 mock zooplanktonic communities incorporating various levels of genetic variation. With the use of a single marker and single primer pair, the highest species recovery was 77%. With all three COI fragments, we detected 62%-83% of species across the mock communities, while the use of the 18S fragment alone resulted in the detection of 73%-75% of species. The species detection level was significantly improved to 89%-93% when both markers were used. Furthermore, multiplexing did not have a negative impact on the proportion of reads assigned to each species and the total number of species detected was similar to when markers were sequenced alone. Overall, our metabarcoding approach utilizing two barcode markers and multiple primer pairs per barcode improved species detection rates over a single marker/primer pair by 14% to 35%, making it an attractive and relatively cost-effective method for biomonitoring natural zooplankton communities. We strongly recommend combining evolutionary independent markers and, when necessary, multiple primer pairs per marker to increase species detection (i.e., reduce false negatives) in metabarcoding studies.

Seawater detection and biological assessments regarding transmission of the oyster parasite Mikrocytos mackini using qPCR.

Mikrocytos mackini is an intracellular parasite of oysters and causative agent of Denman Island disease in Pacific oysters Crassostrea gigas. Although M. mackini has been investigated for decades, its natural mode of transmission, mechanism for host entry, and environmental stability are largely unknown. We explored these biological characteristics of M. mackini using a recently described quantitative PCR (qPCR) assay. We detected M. mackini in the flow-through tank water of experimentally infected oysters and during disease remission in host tissues following 6 wk of elevated water temperature. Waterborne exposure of oysters to M. mackini further confirmed the potential for extracellular seawater transmission of this parasite and also identified host gill to have the highest early and continued prevalence for M. mackini DNA compared to stomach, mantle, labial palps, or adductor muscle samples. However, infections following waterborne challenge were slow to develop despite a substantial exposure (>106 M. mackini l-1 for 24 h), and further investigation demonstrated that M. mackini occurrence and infectivity severely declined following extracellular seawater incubation of more than 24 h. This study demonstrates a potential for using qPCR to monitor M. mackini in wild or farmed oyster populations during periods of disease remission or from environmental seawater samples. This work also suggests that gill tissues may provide a primary site for waterborne entry and possibly shedding of M. mackini in oysters. Further, although extracellular seawater transmission of M. mackini was possible, poor environmental stability and infection efficiency likely restricts the geographic transmission of M. mackini between oysters in natural environs and may help to explain localized areas of infection.

Barcode-based species delimitation in the marine realm: a test using Hexanauplia (Multicrustacea: Thecostraca and Copepoda).

DNA barcoding has been used successfully for identifying specimens belonging to marine planktonic groups. However, the ability to delineate species within taxonomically diverse and widely distributed marine groups, such as the Copepoda and Thecostraca, remains largely untested. We investigate whether a cytochrome c oxidase subunit I (COI-5P) global pairwise sequence divergence threshold exists between intraspecific and interspecific divergences in the copepods plus the thecostracans (barnacles and allies). Using publicly accessible sequence data, we applied a graphical method to determine an optimal threshold value. With these thresholds, and using a newly generated planktonic marine data set, we quantify the degree of concordance using a bidirectional analysis and discuss different analytical methods for sequence-based species delimitation (e.g., BIN, ABGD, jMOTU, UPARSE, Mothur, PTP, and GMYC). Our results support a COI-5P threshold between 2.1% and 2.6% p-distance across methods for these crustacean taxa, yielding molecular groupings largely concordant with traditional, morphologically defined species. The adoption of internal methods for clustering verification enables rapid biodiversity studies and the exploration of unknown faunas using DNA barcoding. The approaches taken here for concordance assessment also provide a more quantitative comparison of clustering results (as contrasted with "success/failure" of barcoding), and we recommend their further consideration for barcoding studies.

Phylogeographic insights into the invasion history and secondary spread of the signal crayfish in Japan.

Successful invasion by nonindigenous species is often attributed to high propagule pressure, yet some foreign species become widespread despite showing reduced genetic variation due to founder effects. The signal crayfish (Pacifastacus leniusculus) is one such example, where rapid spread across Japan in recent decades is believed to be the result of only three founding populations. To infer the history and explore the success of this remarkable crayfish invasion, we combined detailed phylogeographical and morphological analyses conducted in both the introduced and native ranges. We sequenced 16S mitochondrial DNA of signal crayfish from across the introduced range in Japan (537 samples, 20 sites) and the native range in western North America (700 samples, 50 sites). Because chela size is often related to aggressive behavior in crayfish, and hence, their invasion success, we also measured chela size of a subset of specimens in both introduced and native ranges. Genetic diversity of introduced signal crayfish populations was as high as that of the dominant phylogeographic group in the native range, suggesting high propagule pressure during invasion. More recently established crayfish populations in Japan that originated through secondary spread from one of the founding populations exhibit reduced genetic diversity relative to older populations, probably as a result of founder effects. However, these newer populations also show larger chela size, consistent with expectations of rapid adaptations or phenotypic responses during the invasion process. Introduced signal crayfish populations in Japan originate from multiple source populations from a wide geographic range in the native range of western North America. A combination of high genetic diversity, especially for older populations in the invasive range, and rapid adaptation to colonization, manifested as larger chela in recent invasions, likely contribute to invasion success of signal crayfish in Japan.

Phylogenetic species delimitation for crayfishes of the genus Pacifastacus.

Molecular genetic approaches are playing an increasing role in conservation science by identifying biodiversity that may not be evident by morphology-based taxonomy and systematics. So-called cryptic species are particularly prevalent in freshwater environments, where isolation of dispersal-limited species, such as crayfishes, within dendritic river networks often gives rise to high intra- and inter-specific genetic divergence. We apply here a multi-gene molecular approach to investigate relationships among extant species of the crayfish genus Pacifastacus, representing the first comprehensive phylogenetic study of this taxonomic group. Importantly, Pacifastacus includes both the widely invasive signal crayfish Pacifastacus leniusculus, as well as several species of conservation concern like the Shasta crayfish Pacifastacus fortis. Our analysis used 83 individuals sampled across the four extant Pacifastacus species (omitting the extinct Pacifastacus nigrescens), representing the known taxonomic diversity and geographic distributions within this genus as comprehensively as possible. We reconstructed phylogenetic trees from mitochondrial (16S, COI) and nuclear genes (GAPDH), both separately and using a combined or concatenated dataset, and performed several species delimitation analyses (PTP, ABGD, GMYC) on the COI phylogeny to propose Primary Species Hypotheses (PSHs) within the genus. All phylogenies recovered the genus Pacifastacus as monophyletic, within which we identified a range of six to 21 PSHs; more abundant PSHs delimitations from GMYC and ABGD were always nested within PSHs delimited by the more conservative PTP method. Pacifastacus leniusculus included the majority of PSHs and was not monophyletic relative to the other Pacifastacus species considered. Several of these highly distinct P. leniusculus PSHs likely require urgent conservation attention. Our results identify research needs and conservation priorities for Pacifastacus crayfishes in western North America, and may inform better understanding and management of P. leniusculus in regions where it is invasive, such as Europe and Japan.

Molecular detection of Mikrocytos mackini in Pacific oysters using quantitative PCR.

Mikrocytos mackini is an internationally regulated pathogen and causative agent of Denman Island disease in Pacific oysters Crassostrea gigas. Recent phylogenetic breakthroughs have placed this parasite within a highly divergent and globally distributed eukaryotic lineage that has been designated a new taxonomic order, Mikrocytida. The discovery of this new radiation of parasites is accompanied by a heightened awareness of the many knowledge gaps that exist with respect to the general biology, epizootiology, and potential impact of mikrocytid parasites on hosts, ecosystems, and commercial fisheries. It has also highlighted current shortcomings regarding our ability to detect these organisms. In this study, we developed a species-specific, sensitive, and quantitative method for detecting M. mackini DNA from host tissues using probe-based real-time qPCR technology. A limit of sensitivity between 2 and 5 genome copy equivalents was achieved in a reaction matrix containing ≥ 40 ng/µL host gDNA without inhibition. This detection proved superior to existing methods based on conventional PCR, histology or gross pathology and is the first species-specific diagnostic test for M. mackini. Quantitative assessment of parasite DNA using this assay remained accurate to between 10 and 50 copies identifying that during infection, M. mackini DNA was significantly more prevalent in hemolymph, labial palp, and mid-body cross-sections compared to mantle or adductor muscle. DNA extracted from a mid-body cross-section also provided the highest likelihood for detection during diagnostic screening of infected oysters. Taken together, these findings provide strong analytical evidence for the adoption of qPCR as the new reference standard for detecting M. mackini and give preliminary insight into the distribution of the parasite within host tissues. Standardised operating methodologies for sample collection and qPCR testing are provided to aid in the international regulatory diagnosis of M. mackini and serve as a useful platform for the future development of multiplexed or alternate mikrocytid species detection.

New primers for DNA barcoding of digeneans and cestodes (Platyhelminthes).

Digeneans and cestodes are species-rich taxa and can seriously impact human health, fisheries, aqua- and agriculture, and wildlife conservation and management. DNA barcoding using the COI Folmer region could be applied for species detection and identification, but both 'universal' and taxon-specific COI primers fail to amplify in many flatworm taxa. We found that high levels of nucleotide variation at priming sites made it unrealistic to design primers targeting all flatworms. We developed new degenerate primers that enabled acquisition of the COI barcode region from 100% of specimens tested (n = 46), representing 23 families of digeneans and 6 orders of cestodes. This high success rate represents an improvement over existing methods. Primers and methods provided here are critical pieces towards redressing the current paucity of COI barcodes for these taxa in public databases.

Comparison of whole mitochondrial genome sequences from two clades of the invasive ascidian, Didemnum vexillum.

The mitochondria are the main source of cellular energy production and have an important role in development, fertility, and thermal limitations. Adaptive mitochondrial DNA mutations have the potential to be of great importance in determining aspects of the life history of an organism. Phylogenetic analyses of the globally invasive marine ascidian Didemnum vexillum using the mitochondrial cytochrome c oxidase 1 (COX1) coding region, revealed two distinct clades. Representatives of one clade (denoted by 'B') are geographically restricted to D. vexillum's native region (north-west Pacific Ocean, including Japan), whereas members of the other clade (denoted by 'A') have been introduced and become invasive in temperate coastal areas around the world. Persistence of clade B's restricted distribution may reflect it being inherently less invasive than clade A. To investigate this we sought to determine if the two clades differ significantly in other mitochondrial genes of functional significance, specifically, alterations in amino acids encoded in mitochondrial enzyme subunits. Differences in functional mitochondrial genes could indicate an increased ability for clade A colonies to tolerate a wider range of environmental temperature. Full mitochondrial genomic sequences from D. vexillum clades A and B were obtained and they predict significant sequence differences in genes encoding for enzymes involved in oxidative phosphorylation. Diversity levels were relatively high and showed divergence across almost all genes, with p-distance values between the two clades indicating recent divergence. Both clades showed an excess of rare variants, which is consistent with balancing selection or a recent population expansion. Results presented here will inform future research focusing on examining the functional properties of the corresponding mitochondrial respiration enzymes, of A and B clade enzymes. By comparing closely related taxa that have differing distributions it is possible to identify genes and phenotypes suited to particular environments. The examination of mitochondrial genotypes, and associated enzyme functioning, across populations may aid in our understanding of thermal tolerance and environmental adaptation.

Review of Mikrocytos microcell parasites at the dawn of a new age of scientific discovery.

The genus Mikrocytos is traditionally known for Mikrocytos mackini, the microcell parasite that typically infects Pacific oysters along the west coast of North America. Multiple factors have conspired to create difficulty for scientific research on Mikrocytos parasites. These include their tiny cell size, infections that are often of light intensity, lack of suitable cell lines and techniques for in vitro culture, and the seasonal nature of infections. The extreme rate of molecular evolution in Mikrocytos stymied new species discovery and confounded attempts to resolve its phylogenetic position for many years. Fortunately, 2 recent landmark studies have paved the way forward for future research by drastically changing our understanding of the evolution and diversity of these parasites. No longer an orphan eukaryotic lineage, the phylogenetic placement of Mikrocytos has been confidently resolved within Rhizaria and as sister taxon to Haplosporidia. The genus has also found a taxonomic home within the newly-discovered order, Mikrocytida - a globally distributed lineage of parasites infecting a wide range of invertebrate hosts. Here we review available scientific information on Mikrocytos parasites including their evolution and diversity, host and geographic ranges, epizootiology, and detection of the regulated pathogen, M. mackini. We also make recommendations towards a consistent taxonomic framework for this genus by minimally suggesting the use of 18S rDNA sequence, host species information, and histopathological presentation in new species descriptions. This is timely given that we are likely embarking on a new era of scientific advancements, including species discovery, in this genus and its relatives.

Molecular taxonomy of Mikrocytos boweri sp. nov. from Olympia oysters Ostrea lurida in British Columbia, Canada.

Mikrocytos mackini is a microcell parasite that usually infects Crassostrea gigas distributed along the Pacific Northwest coast of North America. For many years, M. mackini was the only known species in the genus, but there have been multiple recent findings of genetically divergent forms of Mikrocytos in different hosts and in distantly located geographic locations. This note describes M. boweri sp. nov. found in Olympia oysters Ostrea lurida collected from and native to British Columbia, Canada, primarily using a molecular taxonomic approach.

Evolution: hidden at the end of a very long branch.

DNA-based methods continue to unveil the diversity and evolutionary origins of life on Earth. 'Next generation' methods have just solved a long-standing puzzle by uncovering previously unseen yet globally distributed diversity within a lineage of amitochondriate parasites affecting commercially exploited aquatic hosts. This discovery will impact both pure and applied research fields.

Phylogenomics of the intracellular parasite Mikrocytos mackini reveals evidence for a mitosome in rhizaria.

Mikrocytos mackini is an intracellular protistan parasite of oysters whose position in the phylogenetic tree of eukaryotes has been a mystery for many years [1,2]. M. mackini is difficult to isolate, has not been cultured, and has no defining morphological feature. Furthermore, its only phylogenetic marker that has been successfully sequenced to date (the small subunit ribosomal RNA) is highly divergent and has failed to resolve its evolutionary position [2]. M. mackini is also one of the few eukaryotes that lacks mitochondria [1], making both its phylogenetic position and comparative analysis of mitochondrial function particularly important. Here, we have obtained transcriptomic data for M. mackini from enriched isolates and constructed a 119-gene phylogenomic data set. M. mackini proved to be among the fastest-evolving eukaryote lineages known to date, but, nevertheless, our analysis robustly placed it within Rhizaria. Searching the transcriptome for genetic evidence of a mitochondrion-related organelle (MRO) revealed only four mitochondrion-derived genes: IscS, IscU, mtHsp70, and FdxR. Interestingly, all four genes are involved in iron-sulfur cluster formation, a biochemical pathway common to other highly reduced "mitosomes" in unrelated MRO-containing lineages [7]. This is the first evidence of MRO in Rhizaria, and it suggests the parallel evolution of mitochondria to mitosomes in this supergroup.