Ancient metapopulations and extreme sex-biased demography revealed by ABC inference and X-Chromosome diversity in Indo-Pacific reef sharks

By integrating Approximate Bayesian Computation with X-chromosome diversity analyses, this study reveals that grey reef sharks in the Indo-Pacific comprise two ancient, deeply separated metapopulations shaped by extreme sex-biased dispersal and recent secondary contact, supporting their designation as distinct Evolutionary Significant Units.

The Gist

Imagine the grey reef shark as a traveler who loves coral reefs but hates crossing the open ocean. For a long time, scientists have been trying to figure out how these sharks are related to each other across the vast Indo-Pacific Ocean. Are they all one big, connected family, or are they separated into distinct clans?

This paper is like a genetic detective story that uses three main clues to solve the mystery: ancient history, family trees, and a special "sex-based" genetic signature.

The Mystery: One Big Family or Two Separate Clans?

Scientists have debated whether grey reef sharks across the Indian and Pacific Oceans are one continuous population (like a long line of people holding hands across a room) or two separate groups that got cut off from each other long ago.

The Detective Work (The "Time Machine"):
The researchers used a powerful computer simulation called ABC (Approximate Bayesian Computation). Think of this as a "What If?" machine. They ran thousands of simulations of shark history:

• Scenario A: Sharks slowly drifted apart like a line of people walking in a circle.
• Scenario B: Two separate groups formed long ago, lived apart, and only recently met again.

The Verdict: The computer said Scenario B is the winner. The sharks aren't one big family. Instead, there are two ancient, deeply separated clans:

1. The Western/Central Indian Ocean Clan (around places like the Maldives and Chagos).
2. The Eastern Indian Ocean/Pacific Clan (around Australia, Indonesia, and the Pacific islands).

These two groups split apart roughly 1.3 million years ago (a long time in shark years!) and only recently started mixing a little bit again.

The "Sex-Selective" Clue: The X-Chromosome

Here is where it gets really interesting. Sharks, like humans, have sex chromosomes. Females have two X chromosomes, and males have one X and one Y.

In a perfectly balanced world where males and females move around equally, the genetic diversity on the X chromosome should be about 75% of the diversity found on the other chromosomes (autosomes).

The Twist:
The researchers found that for these sharks, the X chromosome diversity was drastically lower—sometimes only 18% to 30% of the normal amount!

The Analogy: The "Bachelors' Club" vs. The "Stay-at-Home Moms"
Imagine a town where:

• The Men (Males) are like adventurous travelers. They leave their hometowns, travel far and wide, and have kids with women in many different towns.
• The Women (Females) are like stay-at-home moms. They are very loyal to their birthplace (natal philopatry). They rarely leave their home reef to have babies.

Because the women stay put, the "female line" of the family tree gets stuck in one spot. If a local reef gets damaged or loses sharks, that specific female lineage disappears forever. The men, however, are constantly bringing in fresh genetic material from far away.

This creates a genetic bottleneck for the X chromosome (which comes from the mother). It's like if you only had one source of water for a whole city, but the pipes were constantly being replaced by travelers bringing water from elsewhere. The "local water" (female genes) gets very scarce and diverse-less, while the "traveler water" (male genes) keeps things mixed.

The study found this "stay-at-home mom" behavior is extreme in the Indian Ocean sharks, leading to a massive loss of genetic variety on the X chromosome.

The "Family Tree" Clue: Inbreeding

The team also looked for "Runs of Homozygosity" (ROH). Think of this as looking for long stretches of identical DNA, which happens when parents are closely related (inbreeding).

• The Finding: The sharks in the Maldives (the Western Indian Ocean) had more of these identical stretches than the Australian sharks.
• The Meaning: The Maldivian sharks are more isolated and have had a smaller population size for a long time. They aren't necessarily suffering from recent inbreeding (like siblings mating), but they have been living in a "small pond" for thousands of years, which makes their genetic pool shallower.

Why Does This Matter? (The "So What?")

This isn't just academic trivia; it changes how we need to protect these sharks.

1. They are different species of "culture": Because the Indian Ocean sharks and the Pacific sharks have been separated for over a million years and have different genetic histories, they should be treated as separate conservation units. You can't manage them as one big group.
2. The "Stay-at-Home" risk: Since the females don't move much, if you destroy a specific reef in the Indian Ocean, you might wipe out a unique female lineage forever. The traveling males can't easily "rescue" that lost lineage because the females won't move to new areas to repopulate.
3. Urgent Protection Needed: The study suggests that the Western and Central Indian Ocean populations are Evolutionarily Significant Units (ESUs). This is a fancy way of saying: "These are unique, irreplaceable groups that need their own specific protection plans." Currently, very little of their habitat is fully protected.

The Bottom Line

The grey reef shark isn't just one big, happy, swimming family. It's actually two ancient clans that split up a million years ago. The females are very loyal to their home reefs, while the males are the wanderers. This behavior has left the Indian Ocean sharks with a very fragile genetic makeup, making them highly vulnerable to habitat loss. To save them, we need to stop treating them as a single group and start protecting their specific, isolated homes.

Technical Summary

1. Problem Statement

Reconstructing the demographic history of marine species like the grey reef shark (Carcharhinus amblyrhynchos) is challenging due to their vast geographic ranges, high site fidelity, and sex-biased dispersal behaviors.

• The Conflict: Previous studies on C. amblyrhynchos have yielded competing hypotheses: one suggesting a continuous metapopulation with isolation-by-distance across the Indo-Pacific, and another proposing a hierarchical structure with deep divergence between the Central/Western Indian Ocean and the Eastern Indian Ocean/Pacific.
• The Gap: Traditional genetic markers (e.g., mtDNA vs. nuclear DNA) often show "mito-nuclear discordance," which has been used to infer sex-biased dispersal (female philopatry). However, recent critiques suggest these inferences may be unreliable due to methodological limitations (e.g., low resolution, reliance on single non-recombining loci).
• The Goal: To resolve the true population structure and quantify the impact of sex-biased dispersal on genomic diversity using a novel integration of Approximate Bayesian Computation (ABC) and sex-chromosome specific diversity metrics.

2. Methodology

The authors employed a multi-faceted approach combining simulation-based inference with comparative genomic analysis.

A. Data Sources

• Samples: Integrated data from DArTseq (reduced representation) and Whole Genome Sequencing (WGS) across the Indo-Pacific.
• Locations: Included populations from the Western Indian Ocean (WIO), Central Indian Ocean (CIO), Eastern Indian Ocean (EIO), and Western/Central Pacific.
• Reference: Utilized a newly published high-quality grey reef shark reference genome (C. amblyrhynchos, GCA_965287045).

B. Analytical Frameworks

1. Approximate Bayesian Computation (ABC):

   • Developed a custom pipeline (MetapopABC) using msprime to simulate three competing demographic models:
     • SST (Stepping Stone): A single continuous metapopulation.
     • HIER (Hierarchical): Two distinct, ancient metapopulations with no current gene flow.
     • HIERSC (Hierarchical + Secondary Contact): Two distinct metapopulations with recent, limited secondary contact.
   • Used ABC Random Forest (ABC-RF) for model selection and parameter estimation, comparing observed summary statistics (Watterson's $\Theta$, $\pi$, Tajima's D, $F_{ST}$, $d_{xy}$) against 10,000–50,000 simulated datasets.

2. Sex-Biased Dispersal Metrics (X-Chromosome Analysis):

   • Calculated the ratio of effective population sizes ($N_e$) between the X chromosome and autosomes ($Q$).
   • $Q_\pi$ (Nucleotide Diversity): $\pi_X / \pi_A$. Under neutrality and equal sex ratios, $Q \approx 0.75$. Deviations indicate sex-biased processes in ancestral populations.
   • $Q_{FST}$ (Differentiation): Derived from $F_{ST}$ ratios. Deviations indicate sex-biased processes occurring after population splits.
   • Comparison: Compared C. amblyrhynchos against other carcharhinids with known dispersal patterns (C. limbatus, Galeocerdo cuvier vs. C. acronotus).

3. Demographic History & Inbreeding:

   • PSMC (Pairwise Sequentially Markovian Coalescent): Reconstructed historical $N_e$ trajectories for Maldives, Ningaloo, and Rowley Shoals populations.
   • Runs of Homozygosity (ROH): Analyzed ROH patterns to assess recent inbreeding and historical bottlenecks.

3. Key Results

A. Population Structure (ABC Inference)

• Model Selection: The HIERSC model (two distinct metapopulations with recent secondary contact) received overwhelming support (Posterior Probability > 0.97).
• Divergence: The Central/Western Indian Ocean metapopulation (including Chagos) split from the Eastern Indian Ocean/Pacific metapopulation approximately 1.28 million years ago (~78,000 generations).
• Connectivity: While deep divergence exists, there is evidence of recent, asymmetric, and limited gene flow (secondary contact) between the two metapopulations.
• Demography: The Central Indian Ocean metapopulation has significantly smaller effective deme sizes ($N_e \approx 1,629$) compared to the Pacific/Eastern Indian Ocean ($N_e \approx 5,047$), implying lower connectivity and higher isolation in the Indian Ocean.

B. Sex-Biased Dispersal ($Q$ Metrics)

• Extreme Reduction in X Diversity: $Q_\pi$ values were drastically lower than the neutral expectation of 0.75 across the species range.
  • WIO/CIO: $Q_\pi \approx 0.18–0.30$ (Extreme reduction).
  • EIO/Pacific: $Q_\pi \approx 0.37–0.46$.
• Spatial Gradient: In the Indian Ocean, $Q_\pi$ decreased with distance from the inferred expansion origin (West Papua), suggesting male-biased dispersal during sequential colonization waves.
• $F_{ST}$ Patterns: $Q_{FST}$ was also below 0.75, indicating that male-biased dispersal has shaped genetic differentiation both historically and in contemporary populations.

C. Historical Demography and Inbreeding

• PSMC Trajectories: Showed a sustained decline in effective population size for the Maldives population, contrasting with stable or increasing trends in Australian populations.
• ROH Analysis: The Maldivian population exhibited higher numbers of ROH (>100kb) and higher $F_{ROH}$ (15.2%) compared to Australian populations (5–7.4%).
• Inbreeding Status: While ROH counts were elevated, no ROH > 1 Mb were detected. This suggests the low diversity is due to long-term historical $N_e$ reduction rather than recent severe bottlenecks or acute inbreeding depression.

4. Key Contributions

1. Resolution of Metapopulation Structure: Provided robust statistical evidence (via ABC-RF) that C. amblyrhynchos consists of two ancient, deeply separated metapopulations rather than a single continuous population, resolving previous conflicting hypotheses.
2. Novel Application of Sex-Chromosome Metrics: Successfully applied $Q_\pi$ and $Q_{FST}$ metrics to a non-model marine organism to quantify sex-biased dispersal, moving beyond the limitations of mito-nuclear discordance.
3. Temporal Scaling of Sex Bias: Demonstrated that sex-biased dispersal operates at different temporal scales: ancestral sex bias (reflected in $Q_\pi$) drove the extreme depletion of X-chromosome diversity, while contemporary processes (reflected in $Q_{FST}$) continue to shape differentiation.
4. Methodological Advancement: Developed and open-sourced a scalable ABC framework (MetapopABC) capable of simulating complex, spatially structured metapopulations with secondary contact.

5. Significance and Conservation Implications

• Evolutionary Significant Units (ESUs): The study provides convergent evidence (deep divergence, reduced diversity, extreme X-chromosome depletion) to designate the Western and Central Indian Ocean populations as distinct Evolutionary Significant Units (ESUs).
• Management Strategy: Current range-wide management is insufficient. The extreme isolation and low genetic diversity of the Indian Ocean populations make them highly vulnerable.
• Sex-Biased Vulnerability: The findings highlight that female philopatry (site fidelity) creates a "genetic bottleneck" for females, making local populations disproportionately vulnerable to habitat loss. Conservation efforts must prioritize protecting specific breeding and nursery grounds in the Indian Ocean, as these populations cannot be easily rescued by immigration from the Pacific.
• Genomic Health: Despite low diversity, the lack of long ROH suggests these populations are not currently suffering from severe inbreeding depression, though their long-term evolutionary potential is compromised by small effective population sizes.

In summary, this study leverages advanced genomic inference to reveal a complex, ancient history of isolation and recent contact in grey reef sharks, driven by strong male-biased dispersal and female philopatry, necessitating a shift toward region-specific conservation strategies.