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