Comparative analysis reveals within-population genome size variation in a rotifer is driven by large genomic elements with highly abundant satellite DNA repeat elements

Genome sizes and repeatome evolution in zoantharians (Cnidaria: Hexacorallia: Zoantharia)

Across eukaryotes, large variations of genome sizes have been observed even between closely related species. Transposable elements as part of the repeated DNA have been proposed and confirmed as one of the most important contributors to genome size variation. However, the evolutionary implications of genome size variation and transposable element dynamics are not well understood. Together with phenotypic traits, they are commonly referred to as the "C-value enigma". The order Zoantharia are benthic cnidarians found from intertidal zones to the deep sea, and some species are particularly abundant in coral reefs. Despite their high ecological relevance, zoantharians have yet to be largely studied from the genomic point of view. This study aims at investigating the role of the repeatome (total content of repeated elements) in genome size variations across the order Zoantharia. To this end, whole-genomes of 32 zoantharian species representing five families were sequenced. Genome sizes were estimated and the abundances of different repeat classes were assessed. In addition, the repeat overlap between species was assessed by a sequence clustering method. The genome sizes in the dataset varied up to 2.4 fold magnitude. Significant correlations between genome size, repeated DNA content and transposable elements, respectively (Pearson's correlation test R2 = 0.47, p = 0.0016; R2 = 0.22, p = 0.05) were found, suggesting their involvement in the dynamics of genome expansion and reduction. In all species, long interspersed nuclear elements and DNA transposons were the most abundant identified elements. These transposable elements also appeared to have had a recent expansion event. This was in contrast to the comparative clustering analysis which revealed species-specific patterns of satellite elements' amplification. In summary, the genome sizes of zoantharians likely result from the complex dynamics of repeated elements. Finally, the majority of repeated elements (up to 70%) could not be annotated to a known repeat class, highlighting the need to further investigate non-model cnidarian genomes. More research is needed to understand how repeated DNA dynamics relate to zoantharian evolution and their biology.

Meiotic transmission patterns of additional genomic elements in Brachionus asplanchnoidis, a rotifer with intraspecific genome size variation

Intraspecific genome size (GS) variation in Eukaryotes is often mediated by additional, nonessential genomic elements. Physically, such additional elements may be represented by supernumerary (B-)chromosomes or by large heterozygous insertions into the regular chromosome set. Here we analyze meiotic transmission patterns of Megabase-sized, independently segregating genomic elements (ISEs) in Brachionus asplanchnoidis, a planktonic rotifer that displays an up to two-fold intraspecific GS variation due to variation in size and number of these elements. To gain insights into the meiotic transmission patterns of ISEs, we measured GS distributions of haploid males produced by individual mother clones using flow cytometry and compared these distributions to theoretical distributions expected under a range of scenarios. These scenarios considered transmission biases resembling (meiotic) drive, or cosegregation biases, e.g., if pairs of ISEs preferentially migrated towards the same pole during meiosis. We found that the inferred transmission patterns were diverse and ranged from positive biases (suggesting drive) to negative biases (suggesting drag), depending on rotifer clone and its ISE composition. Additionally, we obtained evidence for a negative cosegregation bias in some of the rotifer clones, i.e., pairs of ISEs exhibited an increased probability of migrating towards opposite poles during meiosis. Strikingly, these transmission and segregation patterns were more similar among members of a genetically homogeneous inbred line than among outbred members of the population. Comparisons between early and late stages of haploid male embryonic development (e.g., young synchronized male eggs vs. hatched males) showed very similar GS distributions, suggesting that transmission biases occur very early in male development, or even during meiosis. Very large genome size was associated with reduced male embryonic survival, suggesting that excessive amounts of ISEs might be detrimental to male fitness. Altogether, our results indicate considerable functional diversity of ISEs in B. asplanchnoidis, with consequences on meiotic transmission and embryonic survival.

Transposons and Non-coding Regions Drive the Intrafamily Differences of Genome Size in Insects

Genome size (GS) can vary considerably between phylogenetically close species, but the landscape of GS changes in insects remain largely unclear. To better understand the specific evolutionary factors that determine GS in insects, we examined flow cytometry-based published GS data from 1,326 insect species, spanning 700 genera, 155 families, and 21 orders. Model fitting showed that GS generally followed an Ornstein–Uhlenbeck adaptive evolutionary model in Insecta overall. Ancestral reconstruction indicated a likely GS of 1,069 Mb, suggesting that most insect clades appeared to undergo massive genome expansions or contractions. Quantification of genomic components in 56 species from nine families in four insect orders revealed that the proliferation of transposable elements contributed to high variation in GS between close species, such as within Coleoptera. This study sheds lights on the pattern of GS variation in insects and provides a better understanding of insect GS evolution.

•

The most comprehensive variation pattern of insects genome size (GS) to date

•

GS evolution of the Insecta was reflected by the adaptive Ornstein–Uhlenbeck model

•

Ancestral state of insect GS was estimated to be ∼1 Gb

•

Intrafamily GS variations were driven by recent transpositions and non-coding regions