Genome structure of Brachionus asplanchnoidis, a Eukaryote with intrapopulation variation in genome size

Abstract

Eukaryotic genomes vary greatly in size due to variation in the proportion of non-coding DNA, a pattern that emerges both in comparisons at a larger taxonomic scale and at the level of individuals within a species. The rotifer Brachionus asplanchnoidis represents one of the most extreme cases of intraspecific genome size variation among Eukaryotes, displaying almost 2-fold variation within a geographic population. Here we used a whole-genome sequencing approach to identify the underlying DNA sequence differences by assembling a high-quality reference genome draft for one individual of the population and aligning short-reads of 15 individuals from the same geographic population. We identified large, contiguous copy number variable regions (CNVs), which exhibited significant coverage differences among individuals, and whose coverage overall scaled with genome size. CNVs were mainly composed of tandemly repeated satellite DNA, with only few interspersed genes or other sequences, and were characterized by an elevated GC-content. Judging from their distributions across contigs, some CNVs are fragments of accessory (B-)chromosomes while others resemble large insertions to normal chromosomes. CNV patterns in offspring of two parents with divergent genome size, and CNV patterns in several individuals from an inbred line differing in genome size demonstrated inheritance and accumulation of CNVs across generations. Our study provides unprecedented insights into genome size evolution at microevolutionary time scales and thus paves the way for studying genome size evolution in contemporary populations rather than inferring patterns and processes a posteriori from species comparisons.