Jan Grašič

Genome structural variations within species are rare. How selective constraints preserve gene order and chromosome structure is a central question in evolutionary biology that remains unsolved. Our sequencing of several genomes of the appendicularian tunicateOikopleura dioicaaround the globe reveals extreme genome scrambling caused by thousands of chromosomal rearrangements, although showing no obvious morphological differences between these animals. The breakpoint accumulation rate is an order of magnitude higher than in ascidian tunicates, nematodes,Drosophila,or mammals. Chromosome arms and sex-specific regions appear to be the primary unit of macrosynteny conservation. At the microsyntenic level, scrambling did not preserve operon structures, suggesting an absence of selective pressure to maintain them. The uncoupling of the genome scrambling with morphological conservation inO. dioicasuggests the presence of previously unnoticed cryptic species and provides a new biological system that challenges our previous vision of speciation in which similar animals always share similar genome structures.

SUMMARY: Genes are not randomly distributed throughout chromosomes. How gene order evolves and how selective constraints act to preserve or vary gene order, both at the macrosyntenic level of whole chromosomes or microsyntenic level of gene blocks, are central questions of evolutionary biology and genomics that remain largely unsolved. Here, after sequencing several genomes of the appendicularian tunicate Oikopleura dioica from different locations around the globe, we show an unprecedented amount of genome scrambling in animals with no obvious morphological differences, consistent with cryptic speciation. Our assemblies suggest that all members of this clade possess a common 3-chromosome karyotype, and that different species largely preserve gene content, despite the presence of thousands of rearrangements in gene order. The movements of genes are largely restricted to chromosome arms and sex-specific regions, which appear to be the primary unit of macrosynteny conservation, and examples of these within-arm movements can be seen in the Hox and Fgf gene families. Our approach employing whole-genome alignments demonstrates that segments containing protein-coding elements tend to be preserved at the microsyntenic scale, consistent with strong purifying selection, with appreciably less preservation of non-coding elements. Unexpectedly, scrambling did not preserve operon structure across species, suggesting an absence of selective pressure to maintain operon structure. As well, genome scrambling does not occur uniformly across all chromosomes, as short chromosome arms possess shorter genes, smaller operons, more breakpoints, and elevated dN/dS values compared to long chromosome arms. Estimation of divergence times among the cryptic O. dioica lineages yielded an estimated breakpoint accumulation rate of 6 to 25 breakpoints per megabase per million years, which is an order of magnitude higher than the rates for other ascidian tunicates or Drosophila species. Therefore, O. dioica appears to be an attractive animal system to unravel the mechanisms that underlie gene order and synteny conservation, as well as exploring the limits of genome scrambling without an apparent impact on phenotypic evolution.