Transcriptomes resolve phylogenetic relationships and reveal undescribed diversity in taildropper slugs (Genus Prophysaon)

By generating transcriptomes for six *Prophysaon* species, this study resolves previously unclear phylogenetic relationships, reveals undescribed diversity within the *P. andersonii-P. foliolatum* complex, suggests a role for introgression in their evolution, and proposes synonymizing the subgenus *Mimetarion* with *Prophysaon*.

The Gist

Imagine the Pacific Northwest as a grand, misty theater where the actors are taildropper slugs. These aren't just any slugs; they have a superpower. When a predator grabs them, they can snap off their own tails like a lizard losing a tail to escape, leaving the predator with a wriggling snack while the slug scurries away to safety.

For a long time, scientists were trying to figure out the family tree of these slugs (a group called Prophysaon). They knew there were nine different "characters" in the play, but they were confused about who was related to whom. The old script, written based on how the slugs looked (their shell-less bodies and reproductive organs), didn't quite match the story the DNA was trying to tell.

Here is the story of how this paper solved the mystery, explained simply:

1. The Old Map Was Wrong

Previously, scientists divided these slugs into two distinct "families" (subgenera) based on their anatomy. It was like sorting a deck of cards by color, assuming all red cards were one family and all black cards were another. But the new research, which used a high-tech "DNA microscope" called transcriptomics (reading the slug's entire instruction manual), showed that the old sorting method was a mess.

The new data revealed that the "families" were actually mixed up. One group, which scientists thought was a separate family, was actually just a branch of the main family that had lost a few specific features over time. To fix this, the authors decided to merge the two groups into one big, happy family, simplifying the name of the genus.

2. The "Lost" Slugs Were Found

While studying the slugs, the researchers found something exciting: two groups of slugs that didn't fit the existing descriptions.

Imagine you have a box of identical-looking blue marbles. Suddenly, you find two marbles that are slightly different shades of blue and have unique speckle patterns. You suspect they might be a different type of marble entirely, but you aren't 100% sure yet.

• The "Longview" Slugs: Found near Longview, Washington, these had solid dark brown bodies.
• The "Klahowya" Slugs: Found in the Olympic National Forest, these had a mottled, speckled pattern.

The DNA analysis confirmed these weren't just random variations; they were likely distinct populations, possibly even new species waiting to be officially named. It's like finding two new characters in a play who have been hiding in the wings this whole time.

3. The Family Reunion (and the Mix-Up)

The researchers built a massive family tree using thousands of genetic markers. Here's what they discovered:

• The Blue-Grey Ancestor: One species, the blue-grey slug (P. coeruleum), turned out to be the "great-aunt" of the whole group. It sits alone at the top of the tree, separate from everyone else.
• The Tangled Cousins: The other slugs were closely related, but their family tree was a bit tangled. Because they evolved so quickly and lived in the same areas, their DNA got mixed up a bit over time.
• The Great Mix-Up (Introgression): The study found evidence that these different slug groups had "borrowed" DNA from each other. Think of it like two neighboring villages swapping recipes. Even though they are distinct groups, they occasionally mated and shared genetic traits. This "genetic borrowing" helped shape who they are today.

4. The Tail-Dropping Mystery Solved

One of the big questions was: Why do some slugs have complex reproductive organs while others have simple ones?
The old theory said the "simple" ones were a separate family. The new DNA tree showed that the "simple" ones are actually descendants of the "complex" ones. It's like a family where the great-grandparents had a fancy, elaborate house, but the grandchildren moved into a tiny, simple cabin. The house got smaller, but they are still the same family. The "complex" house was the original, and the "simple" version was a later simplification.

The Bottom Line

This paper is like a genetic detective story. By using advanced technology to read the slugs' DNA, the authors:

1. Rewrote the family tree, proving the old way of grouping them was wrong.
2. Found new suspects (the Longview and Klahowya slugs) that might be new species.
3. Discovered that these slugs have been swapping DNA like neighbors sharing cookies, which helped them survive and evolve.

The authors are essentially saying, "We finally have the right map for these slugs, but there are still some hidden treasures (new species) waiting to be fully explored in the misty forests of the Pacific Northwest."

Technical Summary

1. Problem Statement

The genus Prophysaon (taildropper slugs), endemic to the temperate rainforests of the Pacific Northwest, presents significant taxonomic and phylogenetic challenges.

• Phylogenetic Uncertainty: Despite previous studies using mitochondrial DNA (COI) and reduced-representation genomic data (GBS), relationships among the nine recognized species remained unresolved. The COI gene tree failed to resolve deep nodes, and GBS data lacked sufficient shared loci across species for robust phylogenetic reconstruction.
• Taxonomic Instability: The genus is currently divided into two subgenera based on morphology: Prophysaon (nominotypical) and Mimetarion. This classification relies on genital morphology, specifically the presence of a massive, convoluted epiphallus in Prophysaon versus a simpler tract in Mimetarion. However, the monophyly of these groups and the validity of the subgeneric split have not been tested with genomic data.
• Undescribed Diversity: There are indications of cryptic diversity within the P. andersonii–P. foliolatum species complex, particularly in Washington, but molecular evidence for distinct species status has been lacking.

2. Methodology

The authors employed a phylotranscriptomic approach combined with population genetics and geometric morphometrics.

• Sampling and Sequencing:

  • Specimens of six valid species (P. andersonii, P. coeruleum, P. dubium, P. foliolatum, P. humile, P. vanattae) and two morphologically distinct populations (tentatively P. aff. andersonii "Longview" and "Klahowya") were collected.
  • Total RNA was extracted (focusing on albumin gland and epiphallus for quality) and sequenced using Illumina NextSeq500 (150 bp paired-end).
  • An outgroup (Arion vulgaris) transcriptome was downloaded and processed similarly.

• Transcriptome Assembly and Orthology:

  • Reads were error-corrected (RCorrector), trimmed (TrimGalore), and rRNA removed.
  • De novo assembly was performed using Trinity.
  • Ortholog clusters were identified via all-by-all BLASTN and the MCL algorithm.
  • Dataset Construction: Multiple datasets were created to test robustness:
    1. Single-copy orthologs only.
    2. Single-copy genes extracted from larger families using three branch-cutting methods (Maximum Inclusion, Monophyletic Outgroups, DISCO).
    3. All gene families (including paralogs).

• Phylogenetic Inference:

  • Species Tree: Inferred using ASTRAL-III (for single-copy) and ASTRAL-Pro (for all families) to account for Incomplete Lineage Sorting (ILS) and gene duplication/loss.
  • Concatenation: Gene trees were concatenated and analyzed using IQ-TREE with ModelFinder for substitution models and 1,000 ultrafast bootstrap replicates.
  • Support Metrics: Gene Concordance Factors (GCF) and Site Concordance Factors (SCF) were calculated to assess discordance.

• Population Genetics & Introgression:

  • Reads were mapped to the Arion vulgaris reference genome.
  • SNPs were called (bcftools) and filtered (vcftools).
  • Structure Analysis: Used to infer population structure within the P. andersonii–P. foliolatum complex (testing $K=1$ to $7$).
  • Introgression: D-statistics (ABBA-BABA) and $f$-branch statistics were calculated using Dsuite to detect historical gene flow.

• Morphological Analysis:

  • Radulae (teeth) were imaged via Scanning Electron Microscopy (SEM).
  • Geometric morphometrics (6 fixed landmarks + 4 curves) were applied to the central (rachidian) teeth.
  • Shape variation was analyzed using Generalized Procrustes Analysis (GPA), PCA, and Permutational Procrustes ANOVA.

3. Key Results

• Phylogenetic Resolution:

  • Transcriptomic data successfully resolved species relationships with high support (100% bootstrap, 1.0 posterior probability) across all inference methods.
  • Subgeneric Classification Overturned: The subgenus Mimetarion was recovered as a monophyletic clade nested within the nominotypical Prophysaon. Specifically, P. coeruleum (currently in Prophysaon) was resolved as the sister lineage to all other Prophysaon species. This renders the current subgenus Prophysaon paraphyletic.
  • Rapid Radiation: Low concordance factors in the backbone suggest rapid speciation and high levels of ILS.

• Taxonomic Revision:

  • The authors synonymize the subgenus Mimetarion with Prophysaon. The morphological distinction (epiphallus complexity) is interpreted as a derived trait lost in the Mimetarion lineage, rather than a synapomorphy defining a separate subgenus.

• Undescribed Diversity:

  • Two distinct populations from Washington ("Longview" and "Klahowya") were identified.
  • Phylogenetic Placement: Transcriptome data placed "Klahowya" as sister to P. andersonii and "Longview" as sister to P. foliolatum, rather than to each other.
  • Population Structure: Structure analysis ($K=4$) identified these two populations as distinct genetic clusters with minimal admixture, supporting their status as distinct evolutionary units.
  • Mitochondrial Discordance: The COI gene tree placed both populations as sisters nested within P. andersonii, highlighting the limitations of mitochondrial data for this group.
  • Morphology: Geometric morphometrics of radulae showed no significant shape differences between the populations after removing an outlier (likely a juvenile), suggesting the divergence is primarily genetic or involves non-radular traits.

• Introgression:

  • D-statistics indicated a moderate role of introgression in the genus history (26.8% of tests significant).
  • Evidence of gene flow was found between the P. andersonii–P. foliolatum complex and the Mimetarion clade (P. vanattae and P. humile), likely facilitated by range overlaps during Pleistocene glacial cycles.

4. Significance and Contributions

• First Robust Phylogeny: This study provides the first well-supported species-level phylogeny for Prophysaon, overcoming the limitations of previous mitochondrial and reduced-representation genomic studies.
• Nomenclatural Stability: By demonstrating that the morphological basis for the subgenus Mimetarion is not a reliable indicator of deep evolutionary splits, the authors propose a simplified taxonomy (synonymizing Mimetarion with Prophysaon), resolving a long-standing taxonomic conflict.
• Discovery of Cryptic Diversity: The study highlights the existence of potentially undescribed species within the P. andersonii complex, emphasizing the need for broader sampling and the integration of genomic data in malacology.
• Evolutionary Insights: The findings illustrate how major geological events (Cascades orogeny) and climatic fluctuations (Pleistocene glaciation) have shaped the diversification of PNW endemics through both isolation and secondary contact (introgression).
• Methodological Demonstration: The paper successfully demonstrates the utility of transcriptomics (even with moderate BUSCO completeness) for resolving difficult phylogenetic nodes in non-model organisms where single-copy genes are scarce or difficult to isolate.

5. Conclusion

The authors conclude that while transcriptomic data has resolved the deep phylogeny of Prophysaon and necessitated a taxonomic revision, the genus likely contains more species than currently recognized. Future work must focus on broader geographic sampling, particularly from type localities, and the examination of additional morphological traits (especially reproductive anatomy) to formally describe the newly detected lineages.