The research team, led by Marina Brasó-Vives and Diego A. Hartasánchez and including colleagues working in France and Switzerland, set out to resolve a puzzle that had been building in the literature for years. Earlier genomic surveys of amphioxus species had hinted at extraordinarily high diversity, but such estimates are notoriously sensitive to methodological choices. Assembling and reading the genome of an organism with so many heterozygous positions is technically challenging, because the two haplotypes inherited from an individual’s parents can differ so much that standard assembly pipelines may confuse them or collapse variants. The new work therefore aimed not only to measure diversity but to do so with an accuracy that could withstand scrutiny, and to explain the biological processes that could generate and maintain such extreme variation in a natural population.

To achieve this, the researchers leveraged whole-genome sequencing data from multiple individuals sampled at two geographically distant locations, one in the Atlantic and one in the Mediterranean Sea. This sampling design was deliberate. Adult amphioxus are sedentary burrowers that live partially buried in sandy substrates, filter-feeding in the water column only briefly, and they form localized settlements with very limited mobility as adults. If genetic diversity were shaped primarily by local processes, populations separated by hundreds of kilometers of coastline and an entire sea basin should show measurable differences. Instead, the team found that population differentiation between the Atlantic and Mediterranean samples was minimal, a result that immediately pointed to extensive gene flow connecting even the most distant adult settlements.

The measurement of heterozygosity itself deserves technical attention, because it is the central quantity of the study. Heterozygosity describes the proportion of positions in the genome where the two copies carried by an individual differ. An average of 2.73 percent means that tens of millions of single nucleotide differences exist between the maternal and paternal haplotypes of a single amphioxus. Such values rival or exceed those reported for many celebrated hyper-diverse organisms, including certain marine fishes and invertebrates with enormous census sizes. The authors confirmed that their estimates were consistent with previous measurements in other amphioxus species, suggesting that extreme diversity is a shared property of the lineage rather than a peculiarity of one population or one sequencing effort.

Explaining this diversity required distinguishing between two broad classes of evolutionary explanation. High heterozygosity can arise when two deeply diverged lineages hybridize, a scenario in which an individual inherits chromosomes from ancestral populations that separated long ago, producing a mosaic genome with ancient differences locked into every cell. Alternatively, high diversity can simply reflect a very large effective population size, the number of breeding individuals that governs how quickly genetic variation is lost through random drift. In small populations, drift purges variation; in vast ones, mutations accumulate over millions of generations and persist, inflating heterozygosity without any need for admixture between diverged lineages.

To separate these hypotheses, the team combined phylogenetic analyses with population genetic simulations. The phylogenetic approach examined the relationships among genome sequences to determine whether amphioxus genomes show the signatures expected from past population structure or admixture, while forward and backward simulations modeled how heterozygosity, differentiation, and genealogical patterns would evolve under different demographic scenarios. The verdict was clear: the elevated genomic diversity of the European amphioxus is primarily driven by a large effective population size. In other words, this is not a hybrid swarm or a relic of ancient population splits but the natural consequence of an enormous number of reproducing individuals maintaining variation across the species’ broad ecological range.

This conclusion raises an obvious question. If adults barely move, how does genetic variation stay evenly distributed across a range spanning the Atlantic and the Mediterranean? The answer lies in the animal’s early life history. Amphioxus reproduce by releasing gametes into the water, and the resulting larvae are planktonic, drifting with ocean currents for extended periods before settling into the sediment as juveniles. The study’s findings indicate that this long-distance larval dispersal is sufficient to generate a near-panmictic population structure, a term describing a population in which mating occurs essentially at random across the entire range, with no meaningful genetic barriers between distant localities. The ocean itself, in effect, functions as the amphioxus’s long-distance transport system, mixing genes on a scale that the sedentary adults never could.

The phrase in the study’s title, at the edge of panmixia, captures the nuance of the result. The population is not perfectly panmictic, and subtle structure may exist, but differentiation is so low that the species across its sampled range behaves, for most population genetic purposes, as a single interbreeding unit. This has practical implications for how the species should be studied and managed. Conservation assessments that treat Atlantic and Mediterranean settlements as separate units would misread the biology, since the larvae continuously homogenize the gene pool. It also matters for laboratory science: amphioxus has become an important model organism for evolutionary developmental biology, and knowing that wild individuals carry extraordinarily diverse genomes informs how researchers should interpret crosses, assemblies, and comparisons between laboratory strains and wild populations.

Beyond amphioxus itself, the study speaks to a broader theme in modern evolutionary genetics: the growing recognition that effective population size, sustained over deep time, is one of the most powerful determinants of genomic diversity. Species with vast, stable marine populations can accumulate and retain levels of variation that land-based animals with smaller or fluctuating populations rarely approach. The amphioxus result also demonstrates the value of pairing empirical sequencing with simulation-based inference, because raw measurements of heterozygosity alone cannot distinguish a large stable population from an admixed one. By combining both approaches, the authors converted a striking number, 2.73 percent, into a coherent biological story about larvae, currents, and drift.

There is also an evolutionary irony worth savoring. Amphioxus are often described as living fossils, conservative in body plan and seemingly unchanged for hundreds of millions of years, yet their genomes reveal a dynamic, well-mixed, and enormously variable population architecture. Morphological stasis and genomic exuberance are not contradictory; a species can look the same while its gene pool churns with variation. As sequencing costs fall and accurate haplotype-aware analyses become routine, researchers expect more animals to join amphioxus in the hyper-diverse club, forcing a rethinking of what typical animal genomes look like. For now, the European amphioxus stands as a reminder that some of the richest genetic libraries on the planet are hidden in plain sand, stirred every generation by the patient mixing of the sea.

Subject of Research: Genomic diversity and population structure of the European amphioxus Branchiostoma lanceolatum

Article Title: The highly heterozygous European amphioxus ( Branchiostoma lanceolatum ) at the edge of panmixia

Article References: Brasó-Vives, M., Hartasánchez, D. A., Pulido, A., Bertrand, S., Escriva, H., & Robinson-Rechavi, M. (2026). The highly heterozygous European amphioxus (Branchiostoma lanceolatum) at the edge of panmixia. PLOS Genetics, 22(9), e1012293. https://doi.org/10.1371/journal.pgen.1012293

Image Credits: AI Generated

DOI: 10.1371/journal.pgen.1012293

Keywords: amphioxus, Branchiostoma lanceolatum, genetic diversity, heterozygosity, population genetics, panmixia, larval dispersal, effective population size, whole-genome sequencing, marine invertebrates, gene flow, PLOS Genetics

The study, published in Advanced Biotechnology, rested on an ambitious sequencing effort. The researchers generated more than fifty-fold short-read whole-genome sequencing coverage for each of 106 individuals from a two-generation family: two parents and 104 F1 offspring. Rather than aligning these reads to an existing reference genome, which would have failed because of the low sequence similarity between the laboratory animals and the published assembly, the team used Platanus-allee, a haplotype-aware assembler designed for highly heterozygous regions, to assemble each parent’s genome de novo. In these assemblies, heterozygous regions appear as paired ‘bubble’ contigs, with a longer primary sequence and a shorter secondary sequence representing the two alleles. Homozygous regions, by contrast, collapse into single non-bubble contigs. Merging the two parental assemblies produced a custom reference to which offspring reads could be mapped accurately, because offspring inherit their DNA directly from these parents.

The logic of the recombination detection is elegant. Each bubble contig pair acts as a biallelic marker: an offspring can inherit either the primary or the secondary contig from a given parent, but not both. Using the R package hapi, the researchers phased these markers and reconstructed chromosome-level haplotypes for each parent. Any switch in the parental haplotype along an offspring’s chromosome signals a crossover (CO) event that occurred during that parent’s meiosis. Because bubble contigs can themselves contain internal crossovers, the team also searched for breakpoints within contigs, where reads from an offspring align to complementary segments of both the primary and secondary contigs, generating four distinct haplotypes in the offspring population. Benchmarking with simulated meioses across a gradient of heterozygosity levels and family sizes showed that when parental heterozygosity exceeded 0.01 and more than 20 offspring were sequenced, the strategy recovered roughly 95 percent of simulated crossovers, with most misses confined to chromosome termini lacking marker coverage.

The results paint a vivid picture of amphioxus meiosis. Across all offspring, the researchers identified 2,329 paternal and 2,288 maternal crossovers, averaging 22.4 and 22.0 events per meiosis respectively. After adjusting for approximately 9 megabases of uncallable chromosome-end regions, this translates into recombination rates of 4.66 centimorgans per megabase in the father and 4.57 cM/Mb in the mother, substantially higher than the roughly 0.5 cM/Mb typical of mammals and the approximately 2 cM/Mb seen in birds. The authors note, however, that this elevated rate per unit of DNA partly reflects the compact size of the amphioxus genome: averaged per haploid chromosome, the count of about 1.2 crossovers per meiosis is close to that of vertebrates. Strikingly, unlike in humans and several other vertebrates where one sex shows markedly higher recombination, no significant sex bias emerged in either the quantity or the chromosomal distribution of crossovers, with paternal and maternal distributions correlating strongly across the genome.

One of the most intriguing findings concerns recombination coldspots. The paternal genome contained roughly 209 megabases of regions spanning at least five megabases with no crossovers in any offspring, while the maternal genome harbored about 142 megabases. These coldspots could not be explained by marker sparsity, since many showed dense marker coverage while crossovers were detected in marker-poor regions elsewhere. Instead, the evidence points to large-scale heterozygous inversions. Amphioxus genomes are known to carry structural variation rates approximately thirty times higher than those observed between humans and chimpanzees, and the team identified four large heterozygous inversions in the parental genomes, all situated within coldspots. In one maternal example on chromosome 3, one haplotype carried an inversion spanning roughly 1.4 to 13.5 megabases, likely preventing crossovers by disrupting pairing and recombination repair in that region. Because inversion polymorphisms appear common in wild amphioxus populations, individual genomes may differ substantially in where recombination is suppressed, which may explain why the pedigree-based crossover map correlated only weakly with historical recombination rates inferred from population genomic variation.

The analysis also probed crossover interference, the phenomenon by which one crossover reduces the likelihood of another nearby. About a quarter of chromosomes in the offspring carried two or more crossovers in both parental genomes. After excluding coldspots and comparing observed inter-crossover distances with a shuffled null distribution, paternal crossovers showed no evidence of interference, while maternal crossovers exhibited significant negative interference, meaning one crossover actually increased the chance of another nearby (Wilcoxon signed-rank test, p = 7.07 x 10^-8). Negative interference remains rare among studied species and its molecular basis is largely unknown, making this maternal effect in amphioxus a curious exception that may shed light on how interference mechanisms evolve across chordates. The researchers also detected crossovers near the sex-determining region on the fused chromosome 16-20 in the paternal genome but not the maternal one, consistent with recombination suppression in the heterogametic ZW sex of this species.

Examining genomic features, the team found that crossover density correlated positively with transposable element density in both parents, with Pearson correlations of 0.43 in the father and 0.40 in the mother, and all major TE classes, including SINEs, LINEs, LTRs, MITEs and DNA transposons, individually showed significant associations. This positive relationship contrasts with the negative correlation often reported in other organisms, underscoring the complex and lineage-dependent interplay between transposable elements and recombination. Gene density also correlated with crossover rates, though more weakly, plausibly reflecting the more open chromatin configuration of gene-rich regions that facilitates binding of recombination proteins. GC content showed no significant correlation with paternal crossovers and only a weak negative correlation with maternal ones, suggesting that amphioxus recombination is not strongly shaped by base composition, unlike in many mammals where PRDM9-directed hotspots often coincide with GC-rich sequence.

Beyond crossovers, the high heterozygosity of amphioxus enabled an unusually sensitive search for non-crossover (NCO) gene conversion events, in which a short segment of one homolog is copied onto the other without a reciprocal exchange. Within analyzable bubble regions, the researchers identified 9,992 paternal and 5,850 maternal NCO events, averaging 96.1 and 56.3 per offspring. Most tracts were short, with 90.9 percent under 200 base pairs and an average length of 99 base pairs, and 83.6 percent covered only a single converted marker. Extrapolated genome-wide, the NCO-to-CO ratios of 28.8 in males and 20.5 in females far exceed those reported for humans and zebra finch, although the authors caution that the amphioxus figures may be inflated by the exceptional marker density and are likely underestimates overall since NCOs in marker-free regions cannot be detected. The conversion spectrum revealed GC-biased gene conversion, favoring C/G over A/T alleles, in most event classes, particularly in the paternal genome, indicating that this force, well documented in mammals, also operates in the basal chordate lineage.

The broader significance of this work is methodological as much as biological. Previous direct detections of recombination without gamete sequencing typically required pedigrees spanning three or more generations, and long-read sequencing approaches can cost more than ten times as much per unit of data. By exploiting bubble contigs as natural markers, the new pipeline achieves contig-level parental haplotyping from just two generations and affordable short-read data, offering a template for studying recombination in the many non-model organisms burdened with high heterozygosity. The findings illuminate both conserved features, such as short NCO tracts and GC-biased gene conversion, and divergent ones, including the absence of sex-biased crossover rates and the prevalence of inversion-driven coldspots, in a species that anchors our understanding of vertebrate origins. As the authors acknowledge, results from a single family cannot capture population-level diversity, and regions of lower heterozygosity remain under-sampled. Nevertheless, by transforming amphioxus’s genomic eccentricity into a high-resolution lens on meiosis, the study opens the door to recombination mapping across a wide swath of the living world that has until now remained largely invisible to geneticists.

Subject of Research: Direct genome-wide detection of meiotic crossover and non-crossover recombination events in the highly heterozygous amphioxus Branchiostoma floridae

Article Title: Direct detection of meiotic recombination events in the highly heterozygous amphioxus genome

Article References: Tao, L., Xue, J., Cao, J., Li, G., & Li, C. (2025). Direct detection of meiotic recombination events in the highly heterozygous amphioxus genome. Advanced Biotechnology, 3(4), Article 30. https://doi.org/10.1007/s44307-025-00083-7

Image Credits: AI Generated

DOI: 10.1007/s44307-025-00083-7

Keywords: amphioxus, meiotic recombination, crossover, non-crossover gene conversion, heterozygosity, Branchiostoma floridae, haplotype phasing, genome assembly, heterozygous inversions, GC-biased gene conversion, transposable elements, chordate evolution