All for one or one for all? Disentangling the Juncus bufonius complex through morphometrics, cytometry and genomics

By integrating genomic, cytometric, and morphological data across a broad latitudinal range, this study reveals that the *Juncus bufonius* complex consists of two distinct ploidy groups (diploids and polyploids) with continuous morphological variation and extensive geographic mixing driven by long-distance dispersal, leading to the recommendation to synonymize *J. ranarius* with *J. hybridus* and group all polyploids under *J. bufonius*.

The Gist

Imagine you are trying to sort a giant, messy pile of LEGO bricks. You have been told there are four distinct types of bricks: Red, Blue, Green, and Yellow. The instruction manual (the old scientific classification) says, "If it's small and has a smooth top, it's Red. If it's big and has a bumpy top, it's Blue."

But when you look at the pile, you realize the bricks are all jumbled together. Some "Red" bricks look exactly like "Blue" ones. Some "Green" bricks are the same size as "Yellow" ones. It's a chaotic mess, and nobody can agree on how to sort them.

This is exactly the problem scientists faced with the Toad Rush (Juncus bufonius), a common wetland plant. For decades, botanists argued over whether there were four different species of Toad Rush or just a few. They tried to sort them by looking at their size and shape (morphology), but the plants were too variable and confusing.

In this new study, the researchers decided to stop guessing and start using DNA detective work to solve the mystery. Here is what they found, explained simply:

1. The "DNA Fingerprint" vs. The "Look-Alike"

The researchers collected plants from all over Europe, from the sunny marshes of Spain to the wetlands of the UK. They used three tools to investigate:

• The Tape Measure (Morphometrics): Measuring the size of flowers and seeds.
• The DNA Counter (Cytometry): Counting how many "sets" of chromosomes (the plant's instruction manual) each plant had.
• The Genetic Scanner (Genomics): Reading the actual DNA code to see the family tree.

The Result: The tape measure was useless. A plant that looked like a "Red" brick could actually be a "Blue" brick genetically. The old way of sorting them by appearance was like trying to sort people by their hair color to figure out their family history—it just doesn't work because hair color varies too much.

2. The Real Family Tree: Two Groups, Not Four

When they looked at the DNA, the messy pile suddenly sorted itself into two clear groups:

• Group A (The Diploids): Plants with two sets of chromosomes.
• Group B (The Polyploids): Plants with extra sets of chromosomes (four or six sets).

The Big Surprise: The old names for the different species (like J. ranarius and J. hybridus) were completely wrong. The DNA showed that all the "two-set" plants were actually one big, mixed-up family. They couldn't tell them apart genetically.

• The Verdict: The scientists say, "Stop calling them different species. They are all the same." They recommend merging the names.

3. The "Genetic Smoothie" (How the Big Plants Were Made)

The researchers also figured out how the bigger plants (the ones with 4 or 6 sets of chromosomes) came to be.

• Imagine the Diploid (2 sets) and the Tetraploid (4 sets) plants as two different families.
• The Hexaploid (6 sets) plants are essentially a genetic smoothie made by mixing the DNA of the 2-set and 4-set families together.
• It's like if a parent from Family A and a parent from Family B had a baby, and that baby somehow got both parents' full instruction manuals plus a copy of one of them. This is called allopolyploidy.

4. The "Bird Taxi" Effect

One of the most fascinating parts of the story is how these plants traveled so far.

• The study found that a plant from Spain is genetically almost identical to a plant in the UK or the Netherlands.
• How did they get there? Birds.
• Migratory waterbirds (like geese and ducks) eat the seeds of these rushes. The seeds survive the trip through the bird's stomach and are dropped hundreds of miles away in a new wetland.
• The Analogy: Think of the birds as taxi drivers for the seeds. They pick up a "passenger" (a seed) in Spain and drop it off in the UK. Because the birds fly so far and so fast, the plants don't stay in one place. They get mixed up constantly, which is why you can't find a "pure" Spanish plant or a "pure" British plant—they are all part of one giant, traveling family.

The Final Conclusion

The paper concludes that we need to throw out the old instruction manual.

1. Forget the four species: There aren't four distinct types.
2. There are really just two main types: The "Small" ones (Diploids) and the "Big" ones (Polyploids).
3. The "Big" ones are hybrids: They are the result of the "Small" ones mixing with other versions of themselves.
4. Birds are the heroes: Migratory birds are the reason these plants are spread all over Europe and why they are so hard to tell apart.

In short: The Toad Rush isn't a group of four different species trying to hide from us. It's a single, chaotic, traveling family that has been mixing its DNA for a long time, carried across the continent by hungry birds. The scientists are saying, "Let's stop trying to sort the LEGO bricks by color and just admit they're all part of the same big box."

Technical Summary

1. Problem Statement

The Juncus bufonius L. species complex (toad rushes) is a taxonomically intricate group within the Juncaceae family, characterized by multiple ploidy levels (diploid, tetraploid, and hexaploid). Despite extensive morphological research, species delimitation remains unclear due to:

• Lack of Diagnostic Characters: Extreme morphological variability and continuous variation make it difficult to distinguish between proposed morphospecies (J. bufonius s.str., J. minutulus, J. ranarius, and J. hybridus).
• Taxonomic Controversy: Previous studies have yielded contradictory results regarding the validity of separating tetraploids (J. minutulus) from hexaploids (J. bufonius s.str.) or splitting diploids into J. ranarius and J. hybridus.
• Data Gaps: There is a paucity of molecular studies, particularly genomic ones, to resolve the reticulate evolutionary history caused by polyploidy and potential hybridization.

2. Methodology

The study employed an integrative approach combining morphometrics, flow cytometry, and phylogenomics across a broad latitudinal range (England to Spain).

• Sampling:

  • Morphometrics: 485 individuals from 17 locations and 60 herbarium specimens were measured. Quantitative variables included capsule length (CL), outer/inner tepal lengths (OTL/ITL), and their ratios. Qualitative traits included seed coat morphology and flower disposition.
  • Cytometry: Flow cytometry was performed on individuals from 31 populations to determine ploidy levels (2C DNA content) using Solanum lycopersicum as an internal standard.
  • Genomics (Hyb-Seq): Genomic data were generated for a subset of 36 individuals (representing all ploidy levels) using the Angiosperm353 kit. This involved target enrichment of nuclear and plastid genes.

• Data Analysis:

  • Morphological: Principal Component Analysis (PCA) and K-means clustering were used to test for morphological grouping. Non-parametric tests (Kruskal-Wallis, Dunn's test) compared morphological variables across ploidy levels.
  • Phylogenomics:
    • Nuclear: 272 nuclear supercontigs were aligned, trimmed, and used to infer a species tree under a multispecies coalescent model (ASTRAL-III) based on 272 gene trees.
    • Plastid: 57 concatenated plastid supercontigs were analyzed using Maximum-Likelihood (IQ-TREE 2).
  • Population Genomics: Single Nucleotide Polymorphisms (SNPs) were called (847 filtered SNPs from 28 individuals). Population structure was analyzed using ADMIXTURE (ancestry proportions), PCA, and Discriminant Analysis of Principal Components (DAPC).

3. Key Results

• Cytometric Distribution:

  • Three distinct ploidy levels were confirmed: Diploid (0.55–0.79 pg), Tetraploid (1.02–1.37 pg), and Hexaploid (1.49–2.08 pg).
  • Latitudinal Gradient: Ploidy levels correlated with latitude; diploids increased toward the south (Spain), while hexaploids increased toward the north (UK/Netherlands).
  • Sympatry: Most populations (14/23) contained mixed ploidy levels, primarily diploids and hexaploids.

• Morphological Findings:

  • No Discrete Groups: PCA and K-means clustering failed to separate individuals into distinct morphological groups corresponding to the proposed morphospecies. Individuals were distributed in a continuous cloud.
  • Ploidy vs. Morphology: While some significant differences existed between diploids and hexaploids (e.g., inner tepal length), tetraploids were statistically indistinguishable from both. Qualitative traits (seed ornamentation, flower clustering) showed extensive overlap and were unreliable for taxonomy.

• Phylogenomic and Population Structure:

  • Major Split: Both nuclear and plastid phylogenies strongly supported a dichotomy between diploids and polyploids (tetraploids + hexaploids).
  • Lack of Sub-structuring:
    • Diploids: J. ranarius and J. hybridus were completely intermixed in both nuclear and plastid trees, with no genetic support for separating them.
    • Polyploids: Tetraploids and hexaploids did not form monophyletic groups based on morphology or geography.
  • Hybrid Origin: ADMIXTURE analysis suggested that most hexaploids are allopolyploids resulting from the hybridization of diploid and tetraploid gene pools.
  • Geographic Mixing: Phylogenetic lineages were extensively mixed geographically (e.g., UK and Spanish samples clustering together), indicating frequent long-distance dispersal events.

4. Key Contributions

1. Taxonomic Revision: The study provides robust evidence to collapse the current complex taxonomy. It recommends:
   • Treating J. ranarius as a synonym of J. hybridus (the older name), as they represent a single diploid taxon.
   • Grouping tetraploids (J. minutulus) and hexaploids (J. bufonius s.str.) under a single taxon, J. bufonius s.l.
2. Mechanism of Polyploidization: Confirmed the allopolyploid origin of most hexaploids via hybridization between diploids and tetraploids, while noting the potential for autopolyploidy in specific lineages (e.g., hexaploid ES-DULC12 clustering with diploids).
3. Dispersal Dynamics: Demonstrated that migratory waterbirds facilitate extensive long-distance seed dispersal (>100 km), leading to high genetic connectivity and the breakdown of geographic barriers between cytotypes.
4. Ecological Differentiation: Identified potential local adaptation where diploids dominate stable ricefield environments in SW Spain, while polyploids are more prevalent in disturbed or variable natural marshes, suggesting niche differentiation based on ploidy.

5. Significance

This paper resolves a long-standing taxonomic debate in the Juncus bufonius complex by demonstrating that traditional morphological characters are insufficient for species delimitation in polyploid complexes. By integrating genomics, the authors show that the complex is best understood as two primary evolutionary units (diploids and polyploids) rather than multiple distinct species.

The findings highlight the role of endozoochory (seed dispersal by birds) in shaping the biogeography of polyploid plants, allowing them to overcome minority cytotype exclusion and colonize novel niches. This study serves as a model for disentangling reticulate evolution in other polyploid plant complexes, emphasizing the necessity of genomic data over morphology alone in such systems.