Resolving the Taxonomic Status of the Marbled Toad (Bufonidae: Incilius marmoreus): 2RAD-based Phylogeography Including an Isolated Population in Veracruz, Mexico

This study utilizes multidisciplinary morphological and genome-wide SNP analyses to confirm that the geographically disjunct Marbled Toad (*Incilius marmoreus*) populations in Veracruz and the Pacific Coast of Mexico represent a single species with a north-south split occurring approximately 0.86 million years ago, followed by a more recent east-west divergence around 0.33 million years ago likely driven by Late Pleistocene climatic changes.

The Gist

Imagine a family reunion where everyone looks a little different, and some relatives live in a completely different town, separated by a massive mountain range. You start to wonder: Are these people actually part of the same family, or are they two different families that just happen to look similar?

This is exactly the mystery scientists tackled in this paper regarding the Marbled Toad (Incilius marmoreus).

The Mystery of the "Lost" Toads

Most of these toads live along the sunny, dry Pacific Coast of Mexico, stretching for hundreds of miles. But there's a weird twist: a small group of them also lives far away on the Gulf Coast in the state of Veracruz.

Between these two groups lies the Isthmus of Tehuantepec, a narrow strip of land that acts like a giant wall. For a long time, scientists were confused. Because the Veracruz toads live in a different place and look slightly smaller, some thought they might be a different species entirely. It was like finding a group of cousins in a distant village who were shorter than the rest of the family and wondering if they were actually adopted.

The Detective Work: Three Tools for the Job

To solve this family feud, the researchers didn't just guess; they used three high-tech "detective tools":

1. The Tape Measure (Morphology): They grabbed a digital caliper and measured 171 toads from museums. They checked everything from snout length to eye size.

   • The Analogy: Imagine measuring the height and shoe size of everyone at a family reunion. They found that the Veracruz toads were indeed a bit smaller (like a "dwarf" branch of the family), but they weren't that different. It's like comparing a tall basketball player to a slightly shorter one; they're still clearly the same type of athlete.

2. The DNA Fingerprint (Genetics): This was the big gun. They used a method called 2RAD to read thousands of tiny genetic markers (SNPs) from the toads' DNA.

   • The Analogy: Think of DNA as a family recipe book. If two groups have different recipes, they are different families. If they share the same recipe but just added a pinch of extra salt, they are the same family with a local variation. The genetic "recipe" showed that the Veracruz toads and the Pacific toads share the exact same core family history.

3. The Family Tree (Phylogeny): They built a digital family tree to see who is related to whom.

   • The Result: The tree showed that the Veracruz toads are actually nested inside the family of the toads from Oaxaca (a state on the Pacific side). They aren't a separate branch; they are a twig growing off the same branch.

The Verdict: One Family, Two Neighborhoods

The study concluded that the Marbled Toad is just one single species, despite living in two separated areas.

• The "Dwarf" Effect: The Veracruz toads are smaller, but this is likely just a local adaptation to their specific environment, not a sign of a new species.
• The Timeline: The genetic clock ticked back about 330,000 years. That's when the Veracruz group split off from the southern Pacific group.
• The "How": How did they get separated? The scientists think that during the last Ice Age, the climate changed. The land bridge connecting them might have turned into a wet, swampy jungle that the dry-loving toads couldn't cross. It's like a river suddenly flooding a road, trapping a group of hikers on one side while the rest of the group stayed on the other.

Why Does This Matter?

In the world of biology, naming things correctly is crucial. If we think these are two different species, we might protect them differently or manage their habitats in separate ways.

This paper is like a peace treaty. It tells us: "Stop worrying! These are all the same family. The small ones in Veracruz are just the 'city cousins' who live in a different neighborhood and eat a slightly different diet, but they are still your kin."

The Bottom Line: The Marbled Toad is a single, widespread species that has managed to survive in two very different corners of Mexico, separated by a natural barrier that formed a few hundred thousand years ago. They are different in size and location, but genetically, they are one and the same.

Technical Summary

1. Problem Statement

The Marbled Toad (Incilius marmoreus) exhibits an unusual disjunct distribution in Mexico: a broad range along the Pacific Coast (from Sonora to Chiapas) and a small, isolated (allopatric) population in the state of Veracruz on the Gulf Coast. Historically, Günther (1901) synonymized Pacific coastal taxa (Bufo argillaceus and B. lateralis) with I. marmoreus, noting the type locality is in Veracruz. However, the validity of treating these geographically separated populations as a single species remained untested. Previous studies on Incilius suggest that allopatric montane populations often represent distinct species. Furthermore, the morphological similarity between I. marmoreus and the sympatric I. perplexus (particularly in males) complicates identification. The study aims to resolve the taxonomic status of the Veracruz population and the Pacific coastal populations using a multidisciplinary approach to determine if they represent independent evolutionary lineages or a single species with geographic variation.

2. Methodology

The authors employed a comprehensive, multidisciplinary framework combining morphometrics, phylogenomics, and population genetics:

• Morphological Analysis:

  • Data: 171 ethanol-preserved specimens were measured from museum collections (AMNH, CAS, KU, etc.), including 36 from Veracruz and 135 from the Pacific Coast.
  • Metrics: 10 linear measurements (e.g., Snout-Vent Length [SVL], Tibial Length, Head Width, Parotoid dimensions) were taken.
  • Statistical Approach: Welch's t-tests compared SVL between regions. Linear models (log-transformed) tested for allometric differences (interaction between SVL and region) and main effects. Principal Component Analysis (PCA) was used to visualize morphological variation, separating males and females.

• Genomic Data Generation (2RAD):

  • Sequencing: 2RAD (a reduced-representation library method similar to ddRAD) was used to generate genome-wide Single Nucleotide Polymorphisms (SNPs).
  • Samples: 16 I. marmoreus (including 2 from Veracruz) and outgroups (I. perplexus, I. canaliferus, I. mazatlanensis, I. valliceps, and Rhinella horribilis).
  • Assembly: Data were assembled using ipyrad (v0.9.94) with parameters including a 0.85 clustering threshold and minimum depth of 6 reads.

• Phylogenomic and Population Genetic Analyses:

  • Phylogeny: Maximum-likelihood (RAxML) and species-tree inference (SVDQuartets via tetrad) were performed. A Neighbor-Net network was constructed to visualize reticulation.
  • Population Structure: PCA and Bayesian clustering (Structure v2.3.4) were used to identify genetic clusters ($K$).
  • Coalescent Dating: Bayesian Phylogenetics & Phylogeography (BPP) was used under a multi-species coalescent model to estimate divergence times ($\tau$) and effective population sizes ($N_e$) for two scenarios: (1) Veracruz vs. all others, and (2) Veracruz vs. Southern Pacific vs. Northern Pacific.

3. Key Contributions

• Integration of Data: The study successfully integrates morphological data with high-resolution genome-wide SNP data to address a taxonomic question where morphology alone is ambiguous.
• Refinement of Phylogeny: It provides a robust, well-supported phylogeny for the Incilius clade, clarifying relationships between I. marmoreus, I. perplexus, and I. canaliferus.
• Biogeographic Insight: It offers a detailed demographic history of the species, specifically estimating the timing of the divergence between the Pacific and Gulf populations.
• Taxonomic Resolution: It formally evaluates the "General Lineage Concept" of species for this specific case, providing evidence to support or reject the synonymy of historical taxa.

4. Key Results

• Morphology:

  • Veracruz toads are significantly smaller in SVL than Pacific Coast toads (Males: ~49mm vs. ~57mm; Females: ~53mm vs. ~61mm).
  • Some allometric differences were found (e.g., head width and eye size scaling), but these were not consistent after Bonferroni correction.
  • PCA showed separation along PC1 (driven by size), but 95% confidence ellipses overlapped considerably, indicating no discrete morphological break.

• Phylogenetics:

  • Interspecific: I. marmoreus is sister to I. perplexus (Bootstrap support = 100), and this clade is sister to I. canaliferus. This topology differs from a previous study (Mendelson et al., 2011) but aligns better with morphological similarities between marmoreus and perplexus.
  • Intraspecific: The species splits into a Northern Pacific clade (Sinaloa, Nayarit, Jalisco) and a Southern Pacific clade (Guerrero, Oaxaca, Veracruz).
  • Veracruz Position: Crucially, the Veracruz samples are nested within the Oaxaca (Southern Pacific) samples in phylogenetic trees, rather than forming a distinct basal lineage.

• Population Genetics:

  • PCA & Structure: Genetic clustering ($K=2$) separates Northern and Southern populations. At $K=4$, Veracruz forms a distinct cluster, but assignment probabilities show a gradient, suggesting isolation-by-distance rather than a hard barrier.
  • Divergence Times (BPP):
    • The split between Northern and Southern Pacific lineages occurred ~0.86 Mya (Late Pleistocene).
    • The divergence between the Veracruz population and the nearest Southern Pacific (Oaxaca) population occurred much more recently, ~0.33 Mya.

5. Significance and Conclusions

• Taxonomic Conclusion: The authors conclude that Incilius marmoreus constitutes a single species. The Veracruz population, while geographically isolated and morphologically smaller (dwarfed), does not represent an independent evolutionary lineage. The genetic data show it is nested within the southern Pacific lineage, and there is no evidence of reproductive isolation or distinct ecological divergence sufficient to warrant species status.
• Synonymy: The study confirms Günther's (1901) recommendation that Bufo argillaceus and B. lateralis are junior synonyms of I. marmoreus.
• Biogeographic Hypothesis: The ~0.33 Mya divergence suggests the Veracruz population was isolated relatively recently, likely due to Late Pleistocene climatic and vegetation changes that created an ecological barrier (wet Gulf forest vs. dry Pacific forest) across the Isthmus of Tehuantepec, rather than a geological event like a Pliocene seaway (which would imply an older divergence).
• Conservation: While currently listed as "Least Concern," the study highlights the species' vulnerability to climate change, which may further fragment its already disjunct distribution.

In summary, this paper demonstrates that despite significant geographic isolation and morphological variation (size), the Marbled Toad populations across Mexico represent a single, continuously evolving metapopulation lineage, resolving a century-old taxonomic uncertainty through modern genomic techniques.