What's in a genome? The C-value enigma and the evolution of eukaryotic genome content

Repeatome diversity in sea anemone genomics (Cnidaria: Actiniaria) based on the Actiniaria-REPlib library

BACKGROUND

Genomic repetitive DNA sequences (Repeatomes, REPs) are widespread in eukaryotes, influencing biological form and function. In Cnidaria, an early-diverging animal lineage, these sequences remain largely uncharacterized. This study investigates sea anemone REPs (Cnidaria: Actiniaria) in a phylogenetic context. We sequenced and assembled de novo the genome of Actinostella flosculifera and analyzed a total of 38 nuclear genomes to create the first ActiniariaREP library (Actiniaria-REPlib). We compared Actiniaria-REPlib with Repbase and RepeatModeler2 libraries, and used dnaPipeTE to annotate REPs from genomic short-read datasets of 36 species for divergence landscapes.

RESULTS

Our study assembled and annotated the mitochondrial genomes, including 27 newly assembled ones. We re-annotated ~92% of the unknown sequences from the initial nuclear genome library, finding that 6.4-30.6% were DNA transposons, 2.1-11.6% retrotransposons, 1-28.4% tandem repeat sequences, and 1.2-7% unclassifiable sequences. Actiniaria-REPlib recovered 9.4x more REP sequences from actiniarian genomes than Dfam and 10.4x more than Repbase. It yielded 79,903 annotated TE consensus sequences (74,643 known, 5,260 unknown), compared to Dfam with 7,697 (3,742 known, 3,944 unknown) and Repbae (763 known).

CONCLUSIONS

Our study significantly enhances the characterization of sea anemone repetitive DNA, assembling mitochondrial genomes, re-annotating nuclear sequences, and identifying diverse repeat elements. Actiniaria-REPlib vastly outperforms existing databases, recovering significantly more REP sequences and providing a comprehensive resource for future genomic and evolutionary studies in Actiniaria.

Repeat competition and ecological shifts drive the evolution of the mobilome in Rhynchospora Vahl (Cyperaceae), the holocentric beaksedges

BACKGROUND AND AIMS

Genomic changes triggered by polyploidy, chromosomal rearrangements and/ or environmental stress are among factors that affect the activity of mobile elements, particularly long terminal repeat retrotransposons (LTR-RTs) and DNA transposons. Because these elements can proliferate and move throughout host genomes, altering the genetic, epigenetic and nucleotypic landscape, they have been recognized as a relevant evolutionary force. Beaksedges (Rhynchospora) stand out for their wide cosmopolitan distribution, high diversity (~400 species) and holocentric chromosomes related to high karyotypic diversity and a centromere-specific satDNA, Tyba. This makes the genus an interesting model to investigate the interactions between repetitive elements, phylogenetic relationships and ecological variables.

METHODS

Here we used comparative phylogenetic methods to investigate the forces driving the evolution of the entire set of mobile elements (mobilome) in the holocentric genus Rhynchospora. We statistically tested the impact of phylogenetic relationships, abundance of holocentromeric satDNA Tyba, diversity of repeatome composition, ecological variables and chromosome number in mobile element diversification.

KEY RESULTS

Tyba abundance was found to be inversely correlated with LTR-RT content. Decrease of LTR abundance and diversity was also related to increase in chromosome number (likely due to fission events) and colonization of dry environments in the northern hemisphere. In contrast, we found constant LTR insertions throughout time in species with lower chromosome numbers in rainier environments in South America. A multivariate model showed that different traits drive LTR abundance, especially repeat diversity and Tyba abundance. Other mobile elements, such as non-LTR RTs and DNA transposons, had insufficient abundance to be included in our models.

CONCLUSIONS

Our findings suggest that LTR evolution is strongly impacted by the holocentric characteristics of Rhynchospora chromosomes, correlating with species diversification and biome shifts, and supporting a holokinetic drive model of evolution and a competitive scenario with Tyba. Altogether, our results present evidence of multi-trait influence on LTR-RT dynamics and provide a broader understanding of transposable element evolution in a macroevolutionary context.

Comparative Genomics Points to Ecological Drivers of Genomic Divergence Among Intertidal Limpets

Comparative genomic studies of closely related taxa are important for our understanding of the causes of divergence on a changing Earth. This being said, the genomic resources available for marine intertidal molluscs are limited and currently, there are few publicly available high-quality annotated genomes for intertidal species and for molluscs in general. Here we report transcriptome assemblies for six species of Patellogastropoda and genome assemblies and annotations for three of these species (Scurria scurra, Scurria viridula and Scurria zebrina). Comparative analysis using these genomic resources suggest that and recently diverging lineages (10-20 Mya) have experienced similar amounts of contractions and expansions but across different gene families. Furthermore, differences among recently diverged species are reflected in variation in the amount of coding and noncoding material in genomes, such as amount of repetitive elements and lengths of transcripts and introns and exons. Additionally, functional ontologies of species-specific and duplicated genes together with demographic inference support the finding that recent divergence among members of the genus Scurria aligns with their unique ecological characteristics. Overall, the resources presented here will be valuable for future studies of adaptation in molluscs and in intertidal habitats as a whole.

Exploring mitonuclear interactions in the regulation of cell physiology: Insights from interspecies cybrids

Brown hares (Lepus europaeus) and mountain hares (Lepus timidus) frequently hybridize in regions where their range overlaps, producing fertile offspring and enabling gene flow between the species. Despite this, no hybrid species has emerged, suggesting that hybrid backcrosses may incur fitness costs. One potential mechanism for such costs involves the interactions between mitochondrial and nuclear gene products, where incompatibilities between species-specific alleles may reinforce species barriers and lead to hybrid breakdown. However, direct experimental evidence for this hypothesis remains limited. In this study, we used fibroblasts derived from skin biopsies of wild-caught hares to generate cytoplasmic hybrid (cybrid) cell lines, wherein mitochondria and mtDNA from one species were transferred to mitochondria-depleted cells of the other species, creating novel mitonuclear gene combinations while preserving the original diploid nuclear background. Employing a range of techniques - including transcriptomics, metabolomics, microscopy, and respirometry - we explored the consequences of mitochondrial transfer between these hare species. Our results reveal that in the studied species mitonuclear incompatibilities exhibit strong effects on cellular fitness but are limited to specific genotypes. We propose mechanisms of cellular-level incompatibility and their potential consequences for interspecific hybrids, offering new insights into the complexity of mitonuclear interactions.

Single specimen genome assembly of Culicoides stellifer shows evidence of a non-retroviral endogenous viral element

BACKGROUND

Advancing our knowledge of vector species genomes is a key step in our battle against the spread of diseases. Biting midges of the genus Culicoides are vectors of arboviruses that significantly affect livestock worldwide. Culicoides stellifer is a suspected vector with a wide range distribution in North America, for which cryptic diversity has been described.

RESULTS

With just one specimen of C. stellifer, we assembled and annotated the nuclear and mitochondrial genome using the ultra-low input DNA PacBio protocol. The genome assembly is 119 Mb in length with a contig N50 value of 479.3 kb, contains 11% repeat sequences and 18,895 annotated protein-coding genes. To further elucidate the role of this species as a vector, we provide genomic evidence of a non-retroviral endogenous viral element integrated into the genome that corresponds to rhabdovirus nucleocapsid proteins, the same family as the vesicular stomatitis virus.

CONCLUSIONS

This genomic information will pave the way for future investigations into this species's putative vector role. We also demonstrate the practicability of completing genomic studies in small dipterans using single specimens preserved in ethanol as well as introduce a workflow for data analysis that considers the challenges of insect genome assembly.

Comparing Small and Large Genomes Within Monogonont Rotifers

Genome size is an important correlate of many biological features including body size, metabolic rate, and developmental rate and can vary due to a variety of mechanisms, including incorporation of repetitive elements, duplication events, or reduction due to selective constraints. Our ability to understand the causes of genome size variation is hampered by limited sampling of many nonmodel taxa, including monogonont rotifers. Here, we used high-throughput Nanopore sequencing and flow cytometry to estimate genome sizes of nine species of monogonont rotifers representing seven families, including three representatives of Superorder Gnesiotrocha. We annotated the genomes and classified the repetitive elements. We also compared genome size with two biological features: body size and metabolic rate. Body sizes were obtained from the literature and our estimates. Oxygen consumption was used as a proxy for metabolic rate and was determined using a respirometer. We obtained similar genome size estimates from genome assemblies and flow cytometry, which were positively correlated with body size and size-specific respiration rate. Importantly, we determined that genome size variation is not due to increased numbers of repetitive elements or large regions of duplication. Instead, we observed higher numbers of predicted proteins as genome size increased, but currently many have no known function. Our results substantially expand the taxonomic scope of available genomes for Rotifera and provide opportunities for addressing genetic mechanisms underlying evolutionary and ecological processes in the phylum.

Comparative analysis using a chromosome-scale genome assembly for Funaria hygrometrica suggests greater collinearity in mosses than in seed plants

Mosses, the largest lineage of seed-free plants, have smaller and less variable genome sizes than flowering plants. Nevertheless, whether this difference results from divergent genome dynamics is poorly known. Here, we use newly generated chromosome-scale genome assemblies for Funaria hygrometrica and comparative analysis with other moss and seed plant genomes to investigate moss genome dynamics. Although some aspects of moss genome dynamics are seed plant-like, such as the mechanism of genome size change and de novo gain/loss of genes, moss genomes retain higher synteny, and collinearity over evolutionary time than seed plant genomes. Furthermore, transposable elements and genes are more evenly distributed along chromosomes in mosses than in seed plants, a feature shared with other sequenced seed-free plant genomes. Overall, our findings support the hypothesis that large-scale genome structure and dynamics of mosses and seed plants differ. In particular, our data suggest a lower rate of gene order reshuffling along chromosomes in mosses compared to seed plants. We speculate that such lower rate of structural genomic variation and unique chromosome structure in mosses may contribute to their relatively smaller and less variable genome sizes.

Short INDELs and SNPs as markers of evolutionary processes in hybrid zones

Polymorphic short insertions and deletions (INDELs ≤ 50 bp) are abundant, although less common than single nucleotide polymorphisms (SNPs). Evidence from model organisms shows INDELs to be more strongly influenced by purifying selection than SNPs. Partly for this reason, INDELs are rarely used as markers for demographic processes or to detect divergent selection. Here, we compared INDELs and SNPs in the intertidal snail Littorina saxatilis, focussing on hybrid zones between ecotypes, in order to test the utility of INDELs in the detection of divergent selection. We computed INDEL and SNP site frequency spectra using capture sequencing data. We assessed the impact of divergent selection by analyzing allele frequency clines across habitat boundaries. We also examined the influence of GC-biased gene conversion because it may be confounded with signatures of selection. We show evidence that short INDELs are affected more by purifying selection than SNPs, but part of the observed site frequency spectra difference can be attributed to GC-biased gene conversion. We did not find a difference in the impact of divergent selection between short INDELs and SNPs. Short INDELs and SNPs were similarly distributed across the genome and so are likely to respond to indirect selection in the same way. A few regions likely affected by divergent selection were revealed by INDELs and not by SNPs. Short INDELs can be useful (additional) genetic markers helping to identify genomic regions important for adaptation and population divergence.

Organismal complexity strongly correlates with the number of protein families and domains

In the pregenomic era, scientists were puzzled by the observation that haploid genome size (the C-value) did not correlate well with organismal complexity. This phenomenon, called the "C-value paradox," is mostly explained by the fact that protein-coding genes occupy only a small fraction of eukaryotic genomes. When the first genome sequences became available, scientists were even more surprised by the fact that the number of genes (G-value) was also a poor predictor of complexity, which gave rise to the "G-value paradox." The proposed explanations usually invoke mechanisms that increase the information content of each individual gene (protein-protein interactions, intrinsic disorder, posttranslational modifications, alternative splicing, etc.). Less attention has been paid to mechanisms that increase the amount of genetic material but do not increase (or not to the same extent) the amount of information encoded in the genome, such as gene duplication and domain shuffling. Proteins belonging to the same family and/or sharing the same domains often carry out similar or even redundant functions. We thus hypothesized that an organism's number of different protein families and domains should be suitable predictors of organismal complexity. In agreement with our hypothesis, we observed that the number of protein families, clans, domains, and motifs increases from simple to progressively more complex organisms. In addition, these metrics correlate with the number of cell types better than and independently of the number of protein-coding genes and several previously proposed predictors of organismal complexity. Our observations have the potential to represent a resolution to the G-value paradox.

Expansion of three types of transposon superfamilies within 25 Mya lead to large genome size of a rice insect pest

The brown planthoppers (BPH, Nilaparvata lugens), white backed planthopper (WBPH, Sogatella furcifera) and small brown planthopper (SBPH, Laodelphax striatellus) are widely distributed rice insect pests, causing huge annual yield loss of rice production. Though these three planthoppers belong to the same family, Delphacidae of Hemiptera, their genome sizes (GS) are very different, ranging from 541 to 1088 Mb. To uncover the main factors driving GS changes of three planthoppers, we first estimated the GS of their ancestor Fulgoroidea, to be 794.33 Mb, indicating GS expansion in BPH but contraction in SBPH and WBPH. Next, we identified repetitive sequences and compared the TE landscapes, showed that three types of transposon superfamilies, hAT, Tc1-Mariner and Gypsy, expanded within 25 Mya in BPH. In addition, BPH kept ancient TEs of Fulgoroidea dated back to 175 Mya, while SBPH and WBPH have lost most of these ancient TEs. Here, we present evidence that the gain of recently expanded TEs driving the GS expansion and loss of ancient TEs leading to the GS contraction, providing new insights into the mechanism of GS variation.

Multiple losses of aKRAB from PRDM9 coincide with a teleost-specific intron size distribution

BACKGROUND

Primary transcripts are largely comprised of intronic sequences that are excised and discarded shortly after synthesis. In vertebrates, the shape of the intron size distribution is largely constant; however, most teleost fish have a diverged log-bimodal 'teleost distribution' (TD) that is seen only in teleosts. How the TD evolved and to what extent this was affected by adaptative or non-adaptive mechanisms is unknown.

RESULTS

Here, we show that the TD has evolved independently at least six times and that its appearance is linked to the loss of the aKRAB domain from PRDM9. We determined intron size distributions and identified PRDM9 orthologues from annotated genomes in addition to scanning 1193 teleost assemblies for the aKRAB domain. We show that a diverged form of PRDM9 ( β ) is predominant in teleosts whereas the α version is absent from most species. Only a subset of PRDM9- α proteins contain aKRAB, and hence, it is present only in a small number of teleost lineages. Almost all lineages lacking aKRAB (but no species with) had TDs.

CONCLUSIONS

In mammals, PRDM9 defines the sites of meiotic recombination through a mechanism that increases structural variance and depends on aKRAB. The loss of aKRAB is likely to have shifted the locations of both recombination and structural variance hotspots. Our observations suggest that the TD evolved as a side-effect of these changes and link recombination to the evolution of intron size illustrating how genome architectures can evolve in the absence of selection.

The genome of the simultaneously hermaphroditic snail Lymnaea stagnalis reveals an evolutionary expansion of FMRFamide-like receptors

The great pond snail Lymnaea stagnalis has served as a model organism for over a century in diverse disciplines such as neurophysiology, evolution, ecotoxicology and developmental biology. To support both established uses and newly emerging research interests we have performed whole genome sequencing (avg.176 × depth), assembly and annotation of a single individual derived from an inbred line. These efforts resulted in a final assembly of 943 Mb (L50 = 257; N50 = 957,215) with a total of 22,499 predicted gene models. The mitogenome was found to be 13,834 bp long and similarly organized as in other lymnaeid species, with minor differences in location of tRNA genes. As a first step towards understanding the hermaphroditic reproductive biology of L. stagnalis, we identified molecular receptors, specifically nuclear receptors (including newly discovered 2xDNA binding domain-NRs), G protein-coupled receptors, and receptor tyrosine kinases, that may be involved in the cellular specification and maintenance of simultaneously active male and female reproductive systems. A phylogenetic analysis of one particular family of GPCRs (Rhodopsin neuropeptide FMRFamide-receptor-like genes) shows a remarkable expansion that coincides with the occurrence of simultaneous hermaphroditism in the Euthyneura gastropods. As some GPCRs and NRs also showed qualitative differences in expression in female (albumen gland) and male (prostate gland) organs, it is possible that separate regulation of male and female reproductive processes may in part have been enabled by an increased abundance of receptors in the transition from a separate-sexed state to a hermaphroditic condition. These findings will support efforts to pair receptors with their activating ligands, and more generally stimulate deeper insight into the mechanisms that underlie the modes of action of compounds involved in neuroendocrine regulation of reproduction, induced toxicity, and development in L. stagnalis, and molluscs in general.

Chromosome-level genome assembly from a single planthopper Nilaparvata muiri (Hemiptera: Delphacidae)

The planthopper Nilaparvata muiri is a sister species to N. lugens (Hemiptera: Delphacidae), a notorious insect pest in Asian rice fields. N. muiri and N. lugens have a different host preference despite the similarities in many biological features. To better understand the adaptive evolution of planthoppers, comprehensive genomic information on N. muiri and N. lugens are urgently needed. In this study, we used ultra-low input PacBio HiFi libraries and Hi-C sequencing technologies to assemble a reference genome of a single N. muiri at the chromosomal level. The genome size was determined to be 531.62 Mb with a contig N50 size of 2.47 Mb and scaffold N50 size of 38.37 Mb. Totally, 96.61% assembled sequences were anchored to the 15 pseudo-chromosomes. BUSCO analysis yielded an Insecta completeness score of 98.6%. A total of 22,057 protein-coding genes were annotated, and 168.16 Mb repetitive sequences occupying 31.63% of genome were pinpointed. The assembled genome is valuable for evolutionary and genetic studies of planthoppers, and may provide sights to pest control.

Genome size and chromosome number are critical metrics for accurate genome assembly assessment in Eukaryota

The number of genome assemblies has rapidly increased in recent history, with NCBI databases reaching over 41,000 eukaryotic genome assemblies across about 2,300 species. Increases in read length and improvements in assembly algorithms have led to increased contiguity and larger genome assemblies. While this number of assemblies is impressive, only about a third of these assemblies have corresponding genome size estimations for their respective species on publicly available databases. In this paper, genome assemblies are assessed regarding their total size compared to their respective publicly available genome size estimations. These deviations in size are assessed related to genome size, kingdom, sequencing platform, and standard assembly metrics, such as N50 and BUSCO values. A large proportion of assemblies deviate from their estimated genome size by more than 10%, with increasing deviations in size with increased genome size, suggesting nonprotein coding and structural DNA may be to blame. Modest differences in performance of sequencing platforms are noted as well. While standard metrics of genome assessment are more likely to indicate an assembly approaching the estimated genome size, much of the variation in this deviation in size is not explained with these raw metrics. A new, proportional N50 metric is proposed, in which N50 values are made relative to the average chromosome size of each species. This new metric has a stronger relationship with complete genome assemblies and, due to its proportional nature, allows for a more direct comparison across assemblies for genomes with variation in sizes and architectures.

A 160 Gbp fork fern genome shatters size record for eukaryotes

Vascular plants are exceptional among eukaryotes due to their outstanding genome size diversity which ranges ∼2,400-fold, including the largest genome so far recorded in the angiosperm Paris japonica (148.89 Gbp/1C). Despite available data showing that giant genomes are restricted across the Tree of Life, the biological limits to genome size expansion remain to be established. Here, we report the discovery of an even larger eukaryotic genome in Tmesipteris oblanceolata, a New Caledonian fork fern. At 160.45 Gbp/1C, this record-breaking genome challenges current understanding and opens new avenues to explore the evolutionary dynamics of genomic gigantism.

Diversity of genome size and chromosome number in homothallic and heterothallic strains of the Closterium peracerosum-strigosum-littorale complex (Desmidiales, Zygnematophyceae, Streptophyta)

The evolutionary transitions of mating systems between outcrossing and self-fertilization are often suggested to associate with the cytological and genomic changes, but the empirical reports are limited in multicellular organisms. Here we used the unicellular zygnematophycean algae, the Closterium peracerosum-strigosum-littorale (C. psl.) complex, to address whether genomic properties such as genome sizes and chromosome numbers are associated with mating system transitions between homothallism (self-fertility) and heterothallism (self-sterility). Phylogenetic analysis revealed the polyphyly of homothallic strains, suggesting multiple transitions between homothallism and heterothallism in the C. psl. complex. Flow cytometry analysis identified a more than 2-fold genome size variation, ranging from 0.53 to 1.42 Gbp, which was positively correlated with chromosome number variation between strains. Although we did not find consistent trends in genome size change and mating system transitions, the mean chromosome sizes tend to be smaller in homothallic strains than in their relative heterothallic strains. This result suggests that homothallic strains possibly have more fragmented chromosomes, which is consistent with the argument that self-fertilizing populations may tolerate more chromosomal rearrangements.

Chromosome-level genome assembly and characterization of the Calophaca sinica genome

Calophaca sinica is a rare plant endemic to northern China which belongs to the Fabaceae family and possesses rich nutritional value. To support the preservation of the genetic resources of this plant, we have successfully generated a high-quality genome of C. sinica (1.06 Gb). Notably, transposable elements (TEs) constituted ~73% of the genome, with long terminal repeat retrotransposons (LTR-RTs) dominating this group of elements (~54% of the genome). The average intron length of the C. sinica genome was noticeably longer than what has been observed for closely related species. The expansion of LTR-RTs and elongated introns emerged had the largest influence on the enlarged genome size of C. sinica in comparison to other Fabaceae species. The proliferation of TEs could be explained by certain modes of gene duplication, namely, whole genome duplication (WGD) and dispersed duplication (DSD). Gene family expansion, which was found to enhance genes associated with metabolism, genetic maintenance, and environmental stress resistance, was a result of transposed duplicated genes (TRD) and WGD. The presented genomic analysis sheds light on the genetic architecture of C. sinica, as well as provides a starting point for future evolutionary biology, ecology, and functional genomics studies centred around C. sinica and closely related species.

Parallel gene size and isoform expansion of ancient neuronal genes

How nervous systems evolved is a central question in biology. A diversity of synaptic proteins is thought to play a central role in the formation of specific synapses leading to nervous system complexity. The largest animal genes, often spanning hundreds of thousands of base pairs, are known to be enriched for expression in neurons at synapses and are frequently mutated or misregulated in neurological disorders and diseases. Although many of these genes have been studied independently in the context of nervous system evolution and disease, general principles underlying their parallel evolution remain unknown. To investigate this, we directly compared orthologous gene sizes across eukaryotes. By comparing relative gene sizes within organisms, we identified a distinct class of large genes with origins predating the diversification of animals and, in many cases, the emergence of neurons as dedicated cell types. We traced this class of ancient large genes through evolution and found orthologs of the large synaptic genes potentially driving the immense complexity of metazoan nervous systems, including in humans and cephalopods. Moreover, we found that while these genes are evolving under strong purifying selection, as demonstrated by low dN/dS ratios, they have simultaneously grown larger and gained the most isoforms in animals. This work provides a new lens through which to view this distinctive class of large and multi-isoform genes and demonstrates how intrinsic genomic properties, such as gene length, can provide flexibility in molecular evolution and allow groups of genes and their host organisms to evolve toward complexity.

Shark genome size evolution and its relationship with cellular, life-history, ecological, and diversity traits

Among vertebrates, sharks exhibit both large and heterogeneous genome sizes ranging from 2.86 to 17.05 pg. Aiming for a better understanding of the patterns and causalities of shark genome size evolution, we applied phylogenetic comparative methods to published genome-size estimates for 71 species representing the main phylogenetic lineages, life-histories and ecological traits. The sixfold range of genome size variation was strongly traceable throughout the phylogeny, with a major expansion preceding shark diversification during the late Paleozoic and an ancestral state (6.33 pg) close to the present-day average (6.72 pg). Subsequent deviations from this average occurred at higher rates in squalomorph than in galeomorph sharks and were unconnected to evolutionary changes in the karyotype architecture, which were dominated by descending disploidy events. Genome size was positively correlated with cell and nucleus sizes and negatively with metabolic rate. The metabolic constraints on increasing genome size also manifested at higher phenotypic scales, with large genomes associated with slow lifestyles and purely marine waters. Moreover, large genome sizes were also linked to non-placental reproductive modes, which may entail metabolically less demanding embryological developments. Contrary to ray-finned fishes, large genome size was associated neither with the taxonomic diversity of affected clades nor with low genetic diversity.

Remarkably High Repeat Content in the Genomes of Sparrows: The Importance of Genome Assembly Completeness for Transposable Element Discovery

Transposable elements (TE) play critical roles in shaping genome evolution. Highly repetitive TE sequences are also a major source of assembly gaps making it difficult to fully understand the impact of these elements on host genomes. The increased capacity of long-read sequencing technologies to span highly repetitive regions promises to provide new insights into patterns of TE activity across diverse taxa. Here we report the generation of highly contiguous reference genomes using PacBio long-read and Omni-C technologies for three species of Passerellidae sparrow. We compared these assemblies to three chromosome-level sparrow assemblies and nine other sparrow assemblies generated using a variety of short- and long-read technologies. All long-read based assemblies were longer (range: 1.12 to 1.41 Gb) than short-read assemblies (0.91 to 1.08 Gb) and assembly length was strongly correlated with the amount of repeat content. Repeat content for Bell's sparrow (31.2% of genome) was the highest level ever reported within the order Passeriformes, which comprises over half of avian diversity. The highest levels of repeat content (79.2% to 93.7%) were found on the W chromosome relative to other regions of the genome. Finally, we show that proliferation of different TE classes varied even among species with similar levels of repeat content. These patterns support a dynamic model of TE expansion and contraction even in a clade where TEs were once thought to be fairly depauperate and static. Our work highlights how the resolution of difficult-to-assemble regions of the genome with new sequencing technologies promises to transform our understanding of avian genome evolution.

Satellite DNAs, heterochromatin, and sex chromosomes of the wattled jacana (Charadriiformes; Jacanidae): a species with highly rearranged karyotype

Charadriiformes, which comprises shorebirds and their relatives, is one of the most diverse avian orders, with over 390 species showing a wide range of karyotypes. Here, we isolated and characterized the whole collection of satellite DNAs (satDNAs) at both molecular and cytogenetic levels of one of its representative species, named the wattled jacana (Jacana jacana), a species that contains a typical ZZ/ZW sex chromosome system and a highly rearranged karyotype. In addition, we also investigate the in situ location of telomeric and microsatellite repeats. A small catalog of 11 satDNAs was identified that typically accumulated on microchromosomes and on the W chromosome. The latter also showed a significant accumulation of telomeric signals, being (GA)10 the only microsatellite with positive hybridization signals among all the 16 tested ones. These current findings contribute to our understanding of the genomic organization of repetitive DNAs in a bird species with high degree of chromosomal reorganization contrary to the majority of bird species that have stable karyotypes.

A Small Genome amidst the Giants: Evidence of Genome Reduction in a Small Tubulinid Free-Living Amoeba

This study investigates the genomic characteristics of Echinamoeba silvestris, a small-sized amoeba within the Tubulinea clade of the Amoebozoa supergroup. Despite Tubulinea's significance in various fields, genomic data for this clade have been scarce. E. silvestris presents the smallest free-living amoeba genome within Tubulinea and Amoebozoa to date. Comparative analysis reveals intriguing parallels with parasitic lineages in terms of genome size and predicted gene numbers, emphasizing the need to understand the consequences of reduced genomes in free-living amoebae. Functional categorization of predicted genes in E. silvestris shows similar percentages of ortholog groups to other amoebae in various categories, but a distinctive feature is the extensive gene contraction in orphan (ORFan) genes and those involved in biological processes. Notably, among the few genes that underwent expansion, none are related to cellular components, suggesting adaptive processes that streamline biological processes and cellular components for efficiency and energy conservation. Additionally, our investigation into noncoding and repetitive elements sheds light on the evolution of genome size in amoebae, with E. silvestris distinguished by low percentage of repetitive elements. Furthermore, the analysis reveals that E. silvestris has the lowest mean number of introns per gene among the species studied, providing further support for its observed compact genome. Overall, this research underscores the diversity within Tubulinea, highlights knowledge gaps in Amoebozoa genomics, and positions E. silvestris as a valuable addition to genomic data sets, prompting further exploration of complexities in Amoebozoa diversity and genome evolution.

Anciently duplicated genes continuously recruited to heart expression in vertebrate evolution are associated with heart chamber increase

Although gene/genome duplications in the early stage of vertebrates have been thought to provide major resources of raw genetic materials for evolutionary innovations, it is unclear whether they continuously contribute to the evolution of morphological complexity during the course of vertebrate evolution, such as the evolution from two heart chambers (fishes) to four heart chambers (mammals and birds). We addressed this issue by our heart RNA-Seq experiments combined with published data, using 13 vertebrates and one invertebrate (sea squirt, as an outgroup). Our evolutionary transcriptome analysis showed that number of ancient paralogous genes expressed in heart tends to increase with the increase of heart chamber number along the vertebrate phylogeny, in spite that most of them were duplicated at the time near to the origin of vertebrates or even more ancient. Moreover, those paralogs expressed in heart exert considerably different functions from heart-expressed singletons: the former are functionally enriched in cardiac muscle and muscle contraction-related categories, whereas the latter play more basic functions of energy generation like aerobic respiration. These findings together support the notion that recruiting anciently paralogous genes that are expressed in heart is associated with the increase of chamber number in vertebrate evolution.

Half a Century of Controversy: The Neutralist/Selectionist Debate in Molecular Evolution

The neutral and nearly neutral theories, introduced more than 50 yr ago, have raised and still raise passionate discussion regarding the forces governing molecular evolution and their relative importance. The debate, initially focused on the amount of within-species polymorphism and constancy of the substitution rate, has spread, matured, and now underlies a wide range of topics and questions. The neutralist/selectionist controversy has structured the field and influences the way molecular evolutionary scientists conceive their research.

Moderating the neutralist-selectionist debate: exactly which propositions are we debating, and which arguments are valid?

Half a century after its foundation, the neutral theory of molecular evolution continues to attract controversy. The debate has been hampered by the coexistence of different interpretations of the core proposition of the neutral theory, the 'neutral mutation-random drift' hypothesis. In this review, we trace the origins of these ambiguities and suggest potential solutions. We highlight the difference between the original, the revised and the nearly neutral hypothesis, and re-emphasise that none of them equates to the null hypothesis of strict neutrality. We distinguish the neutral hypothesis of protein evolution, the main focus of the ongoing debate, from the neutral hypotheses of genomic and functional DNA evolution, which for many species are generally accepted. We advocate a further distinction between a narrow and an extended neutral hypothesis (of which the latter posits that random non-conservative amino acid substitutions can cause non-ecological phenotypic divergence), and we discuss the implications for evolutionary biology beyond the domain of molecular evolution. We furthermore point out that the debate has widened from its initial focus on point mutations, and also concerns the fitness effects of large-scale mutations, which can alter the dosage of genes and regulatory sequences. We evaluate the validity of neutralist and selectionist arguments and find that the tested predictions, apart from being sensitive to violation of underlying assumptions, are often derived from the null hypothesis of strict neutrality, or equally consistent with the opposing selectionist hypothesis, except when assuming molecular panselectionism. Our review aims to facilitate a constructive neutralist-selectionist debate, and thereby to contribute to answering a key question of evolutionary biology: what proportions of amino acid and nucleotide substitutions and polymorphisms are adaptive?

A Small Genome Amidst the Giants: Evidence of Genome Reduction in a Small Tubulinid Free-Living Amoeba

This study investigates the genomic characteristics of Echinamoeba silvestris , a small-sized amoeba within the Tubulinea clade of the Amoebozoa supergroup. Despite Tubulinea's significance in various fields, genomic data for this clade have been scarce . E. silvestris presents the smallest free-living amoeba genome within Tubulinea and Amoebozoa to date. Comparative analysis reveals intriguing parallels with parasitic lineages in terms of genome size and predicted gene numbers, emphasizing the need to understand the consequences of reduced genomes in free-living amoebae. Functional categorization of predicted genes in E. silvestris shows similar percentages of ortholog groups to other amoebae in various categories, but a distinctive feature is the extensive gene contraction in orphan (ORFan) genes and those involved in biological processes. Notably, among the few genes that underwent expansion, none are related to cellular components, suggesting adaptive processes that streamline biological processes and cellular components for efficiency and energy conservation. The investigation delves into genomic structural evidence, including gene content and repetitive elements, illuminating the distinctive genomic traits of E. silvestris and providing reinforcement for its compact genome size. Overall, this research underscores the diversity within Tubulinea, highlights knowledge gaps in Amoebozoa genomics, and positions E. silvestris as a valuable addition to genomic datasets, prompting further exploration of complexities in Amoebozoa diversity and genome evolution.

The MetaInvert soil invertebrate genome resource provides insights into below-ground biodiversity and evolution

Soil invertebrates are among the least understood metazoans on Earth. Thus far, the lack of taxonomically broad and dense genomic resources has made it hard to thoroughly investigate their evolution and ecology. With MetaInvert we provide draft genome assemblies for 232 soil invertebrate species, representing 14 common groups and 94 families. We show that this data substantially extends the taxonomic scope of DNA- or RNA-based taxonomic identification. Moreover, we confirm that theories of genome evolution cannot be generalised across evolutionarily distinct invertebrate groups. The soil invertebrate genomes presented here will support the management of soil biodiversity through molecular monitoring of community composition and function, and the discovery of evolutionary adaptations to the challenges of soil conditions.

Ecological factors and parity mode correlate with genome size variation in squamate reptiles

BACKGROUND

Evidence of correlation between genome size, the nuclear haploid DNA content of a cell, environmental factors and life-history traits have been reported in many animal species. Genome size, however, spans over three orders of magnitude across taxa and such a correlation does not seem to follow a universal pattern. In squamate reptiles, the second most species-rich order of vertebrates, there are currently no studies investigating drivers of genome size variability. We run a series of phylogenetic generalized least-squares models on 227 species of squamates to test for possible relationships between genome size and ecological factors including latitudinal distribution, bioclimatic variables and microhabitat use. We also tested whether genome size variation can be associated with parity mode, a highly variable life history trait in this order of reptiles.

RESULTS

The best-fitting model showed that the interaction between microhabitat use and parity mode mainly accounted for genome size variation. Larger genome sizes were found in live-bearing species that live in rock/sand ecosystems and in egg-laying arboreal taxa. On the other hand, smaller genomes were found in fossorial live-bearing species.

CONCLUSIONS

Environmental factors and species parity mode appear to be among the main parameters explaining genome size variation in squamates. Our results suggest that genome size may favour adaptation of some species to certain environments or could otherwise result from the interaction between environmental factors and parity mode. Integration of genome size and genome sequencing data could help understand the role of differential genome content in the evolutionary process of genome size variation in squamates.

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.

Long-read genome assemblies for the study of chromosome expansion: Drosophila kikkawai, Drosophila takahashii, Drosophila bipectinata, and Drosophila ananassae

Flow cytometry estimates of genome sizes among species of Drosophila show a 3-fold variation, ranging from ∼127 Mb in Drosophila mercatorum to ∼400 Mb in Drosophila cyrtoloma. However, the assembled portion of the Muller F element (orthologous to the fourth chromosome in Drosophila melanogaster) shows a nearly 14-fold variation in size, ranging from ∼1.3 Mb to >18 Mb. Here, we present chromosome-level long-read genome assemblies for 4 Drosophila species with expanded F elements ranging in size from 2.3 to 20.5 Mb. Each Muller element is present as a single scaffold in each assembly. These assemblies will enable new insights into the evolutionary causes and consequences of chromosome size expansion.

The draft genome of the microscopic Nemertoderma westbladi sheds light on the evolution of Acoelomorpha genomes

Background: Xenacoelomorpha is a marine clade of microscopic worms that is an important model system for understanding the evolution of key bilaterian novelties, such as the excretory system. Nevertheless, Xenacoelomorpha genomics has been restricted to a few species that either can be cultured in the lab or are centimetres long. Thus far, no genomes are available for Nemertodermatida, one of the group's main clades and whose origin has been dated more than 400 million years ago. Methods: DNA was extracted from a single specimen and sequenced with HiFi following the PacBio Ultra-Low DNA Input protocol. After genome assembly, decontamination, and annotation, the genome quality was benchmarked using two acoel genomes and one Illumina genome as reference. The gene content of three cnidarians, three acoelomorphs, four deuterostomes, and eight protostomes was clustered in orthogroups to make inferences of gene content evolution. Finally, we focused on the genes related to the ultrafiltration excretory system to compare patterns of presence/absence and gene architecture among these clades. Results: We present the first nemertodermatid genome sequenced from a single specimen of Nemertoderma westbladi. Although genome contiguity remains challenging (N50: 60 kb), it is very complete (BUSCO: 80.2%, Metazoa; 88.6%, Eukaryota) and the quality of the annotation allows fine-detail analyses of genome evolution. Acoelomorph genomes seem to be relatively conserved in terms of the percentage of repeats, number of genes, number of exons per gene and intron size. In addition, a high fraction of genes present in both protostomes and deuterostomes are absent in Acoelomorpha. Interestingly, we show that all genes related to the excretory system are present in Xenacoelomorpha except Osr, a key element in the development of these organs and whose acquisition seems to be interconnected with the origin of the specialised excretory system. Conclusion: Overall, these analyses highlight the potential of the Ultra-Low Input DNA protocol and HiFi to generate high-quality genomes from single animals, even for relatively large genomes, making it a feasible option for sequencing challenging taxa, which will be an exciting resource for comparative genomics analyses.

A highly contiguous reference genome for the Steller's jay (Cyanocitta stelleri)

The Steller's jay is a familiar bird of western forests from Alaska south to Nicaragua. Here, we report a draft reference assembly for the species generated from PacBio HiFi long-read and Omni-C chromatin-proximity sequencing data as part of the California Conservation Genomics Project (CCGP). Sequenced reads were assembled into 352 scaffolds totaling 1.16 Gb in length. Assembly metrics indicate a highly contiguous and complete assembly with a contig N50 of 7.8 Mb, scaffold N50 of 25.8 Mb, and BUSCO completeness score of 97.2%. Repetitive elements span 16.6% of the genome including nearly 90% of the W chromosome. Compared with high-quality assemblies from other members of the family Corvidae, the Steller's jay genome contains a larger proportion of repetitive elements than 4 crow species (Corvus), but a lower proportion of repetitive elements than the California scrub-jay (Aphelocoma californica). This reference genome will serve as an essential resource for future studies on speciation, local adaptation, phylogeography, and conservation genetics in this species of significant biological interest.

Ancient origins of complex neuronal genes

How nervous systems evolved is a central question in biology. An increasing diversity of synaptic proteins is thought to play a central role in the formation of specific synapses leading to nervous system complexity. The largest animal genes, often spanning millions of base pairs, are known to be enriched for expression in neurons at synapses and are frequently mutated or misregulated in neurological disorders and diseases. While many of these genes have been studied independently in the context of nervous system evolution and disease, general principles underlying their parallel evolution remain unknown. To investigate this, we directly compared orthologous gene sizes across eukaryotes. By comparing relative gene sizes within organisms, we identified a distinct class of large genes with origins predating the diversification of animals and in many cases the emergence of dedicated neuronal cell types. We traced this class of ancient large genes through evolution and found orthologs of the large synaptic genes driving the immense complexity of metazoan nervous systems, including in humans and cephalopods. Moreover, we found that while these genes are evolving under strong purifying selection as demonstrated by low dN/dS scores, they have simultaneously grown larger and gained the most isoforms in animals. This work provides a new lens through which to view this distinctive class of large and multi-isoform genes and demonstrates how intrinsic genomic properties, such as gene length, can provide flexibility in molecular evolution and allow groups of genes and their host organisms to evolve toward complexity.

The State of Squamate Genomics: Past, Present, and Future of Genome Research in the Most Speciose Terrestrial Vertebrate Order

Squamates include more than 11,000 extant species of lizards, snakes, and amphisbaenians, and display a dazzling diversity of phenotypes across their over 200-million-year evolutionary history on Earth. Here, we introduce and define squamates (Order Squamata) and review the history and promise of genomic investigations into the patterns and processes governing squamate evolution, given recent technological advances in DNA sequencing, genome assembly, and evolutionary analysis. We survey the most recently available whole genome assemblies for squamates, including the taxonomic distribution of available squamate genomes, and assess their quality metrics and usefulness for research. We then focus on disagreements in squamate phylogenetic inference, how methods of high-throughput phylogenomics affect these inferences, and demonstrate the promise of whole genomes to settle or sustain persistent phylogenetic arguments for squamates. We review the role transposable elements play in vertebrate evolution, methods of transposable element annotation and analysis, and further demonstrate that through the understanding of the diversity, abundance, and activity of transposable elements in squamate genomes, squamates can be an ideal model for the evolution of genome size and structure in vertebrates. We discuss how squamate genomes can contribute to other areas of biological research such as venom systems, studies of phenotypic evolution, and sex determination. Because they represent more than 30% of the living species of amniote, squamates deserve a genome consortium on par with recent efforts for other amniotes (i.e., mammals and birds) that aim to sequence most of the extant families in a clade.

Human Pangenomics: Promises and Challenges of a Distributed Genomic Reference

A pangenome is a collection of the common and unique genomes that are present in a given species. It combines the genetic information of all the genomes sampled, resulting in a large and diverse range of genetic material. Pangenomic analysis offers several advantages compared to traditional genomic research. For example, a pangenome is not bound by the physical constraints of a single genome, so it can capture more genetic variability. Thanks to the introduction of the concept of pangenome, it is possible to use exceedingly detailed sequence data to study the evolutionary history of two different species, or how populations within a species differ genetically. In the wake of the Human Pangenome Project, this review aims at discussing the advantages of the pangenome around human genetic variation, which are then framed around how pangenomic data can inform population genetics, phylogenetics, and public health policy by providing insights into the genetic basis of diseases or determining personalized treatments, targeting the specific genetic profile of an individual. Moreover, technical limitations, ethical concerns, and legal considerations are discussed.

Protein length distribution is remarkably uniform across the tree of life

BACKGROUND

In every living species, the function of a protein depends on its organization of structural domains, and the length of a protein is a direct reflection of this. Because every species evolved under different evolutionary pressures, the protein length distribution, much like other genomic features, is expected to vary across species but has so far been scarcely studied.

RESULTS

Here we evaluate this diversity by comparing protein length distribution across 2326 species (1688 bacteria, 153 archaea, and 485 eukaryotes). We find that proteins tend to be on average slightly longer in eukaryotes than in bacteria or archaea, but that the variation of length distribution across species is low, especially compared to the variation of other genomic features (genome size, number of proteins, gene length, GC content, isoelectric points of proteins). Moreover, most cases of atypical protein length distribution appear to be due to artifactual gene annotation, suggesting the actual variation of protein length distribution across species is even smaller.

CONCLUSIONS

These results open the way for developing a genome annotation quality metric based on protein length distribution to complement conventional quality measures. Overall, our findings show that protein length distribution between living species is more uniform than previously thought. Furthermore, we also provide evidence for a universal selection on protein length, yet its mechanism and fitness effect remain intriguing open questions.

The evolution of early-life telomere length, pace-of-life and telomere-chromosome length dynamics in birds

Telomeres, the short DNA sequences that protect chromosome ends, are an ancient molecular structure, which is highly conserved across most eukaryotes. Species differ in their telomere lengths, but the causes of this variation are not well understood. Here, we demonstrate that mean early-life telomere length is an evolutionary labile trait across 57 bird species (representing 35 families in 12 orders) with the greatest trait diversity found among passerines. Among these species, telomeres are significantly shorter in fast-lived than in slow-lived species, suggesting that telomere length may have evolved to mediate trade-offs between physiological requirements underlying the diversity of pace-of-life strategies in birds. This association was attenuated when excluding studies that may include interstitial telomeres in the estimation of mean telomere length. Curiously, within some species, larger individual chromosome size predicts longer telomere lengths on that chromosome, leading to the hypothesis that telomere length also covaries with chromosome length across species. We show that longer mean chromosome length or genome size tends to be associated with longer mean early-life telomere length (measured across all chromosomes) within a phylogenetic framework constituting up to 31 bird species. These associations were strengthened when excluding highly influential outliers. However, sensitivity analyses suggested that they were susceptible to sample size effects and not robust to the exclusion of studies that may include interstitial telomeres. Combined, our analyses generalize patterns previously found within a few species and provide potential adaptive explanations for the 10-fold variation in telomere lengths observed among birds.

Long-read genome assemblies for the study of chromosome expansion: Drosophila kikkawai , Drosophila takahashii , Drosophila bipectinata , and Drosophila ananassae

Flow cytometry estimates of genome sizes among species of Drosophila show a 3-fold variation, ranging from ∼127 Mb in Drosophila mercatorum to ∼400 Mb in Drosophila cyrtoloma . However, the assembled portion of the Muller F Element (orthologous to the fourth chromosome in Drosophila melanogaster ) shows a nearly 14-fold variation in size, ranging from ∼1.3 Mb to > 18 Mb. Here, we present chromosome-level long read genome assemblies for four Drosophila species with expanded F Elements ranging in size from 2.3 Mb to 20.5 Mb. Each Muller Element is present as a single scaffold in each assembly. These assemblies will enable new insights into the evolutionary causes and consequences of chromosome size expansion.

Information-energy equivalence and the emergence of self-replicating biological systems

Biological processes are characterized by a decrease in entropy in apparent violation of the second law of thermodynamics. Information stored in genomes help to solve this paradox when interpreted under the relationship between information and energy stated by Brillouin in the 1950's. However, the origins of living forms from inanimate matter which have no information storage device remains an open question. In this paper, a theoretical approach is developed on this issue. The replication of a simple entity with a binary genome is assumed to require an information-equivalent energy in addition to the standard activation energy. It is found that, in some conditions, a decrease in entropy can be accomplished together with a decrease in Gibbs free energy. An equation of the total energy for the replication of this entity is derived. Three factors are predicted to lower this energy: a small number of states of the coding sequence, a lower temperature, and a high ratio of the reaction on diffusion coefficients. These factors may have favoured the emergence of evolutionary demons-information storage devices that are able to decrease entropy. It is evaluated that some short, single-stranded RNA sequences made only of G and of C may conform to this model. The consequences of this model and its predictions on the origins of life on Earth and on other planets are discussed.

Long non-coding RNAs reveal new regulatory mechanisms controlling gene expression

A plethora of non-coding RNAs have been found in eukaryotes, notably with the advent of modern sequencing technologies to analyze the transcriptome. Apart from the well-known housekeeping RNA genes (such as the ribosomal RNA or the transfer RNA), many thousands of transcripts detected are not evidently linked to a protein-coding gene. These, so called non-coding RNAs, may code for crucial regulators of gene expression, the small si/miRNAs, for small peptides (translated under specific conditions) or may act as long RNA molecules (antisense, intronic or intergenic long non-coding RNAs or lncRNAs). The lncRNAs interact with members of multiple machineries involved in gene regulation. In this review, we discussed about how plant lncRNAs permitted to discover new regulatory mechanisms acting in epigenetic control, chromatin 3D structure and alternative splicing. These novel regulations diversified the expression patterns and protein variants of target protein-coding genes and are an important element of the response of plants to environmental stresses and their adaptation to changing conditions.

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.

A Reference Genome Assembly of the Bobcat, Lynx rufus

The bobcat (Lynx rufus) is a medium-sized carnivore well adapted to various environments and an indicator species for landscape connectivity. It is one of the 4 species within the extant Lynx genus in the family Felidae. Because of its broad geographic distribution and central role in food webs, the bobcat is important for conservation. Here we present a high-quality de novo genome assembly of a male bobcat located in Mendocino County, CA, as part of the California Conservation Genomics Project (CCGP). The assembly was generated using the standard CCGP pipeline from a combination of Omni-C and HiFi technologies. The primary assembly comprises 76 scaffolds spanning 2.4 Gb, represented by a scaffold N50 of 142 Mb, a contig N50 of 66.2 Mb, and a BUSCO completeness score of 95.90%. The bobcat genome will be an important resource for the effective management and conservation of this species and comparative genomics exploration.

Comparative genomics reveals insight into the evolutionary origin of massively scrambled genomes

Ciliates are microbial eukaryotes that undergo extensive programmed genome rearrangement, a natural genome editing process that converts long germline chromosomes into smaller gene-rich somatic chromosomes. Three well-studied ciliates include Oxytricha trifallax, Tetrahymena thermophila, and Paramecium tetraurelia, but only the Oxytricha lineage has a massively scrambled genome, whose assembly during development requires hundreds of thousands of precisely programmed DNA joining events, representing the most complex genome dynamics of any known organism. Here we study the emergence of such complex genomes by examining the origin and evolution of discontinuous and scrambled genes in the Oxytricha lineage. This study compares six genomes from three species, the germline and somatic genomes for Euplotes woodruffi, Tetmemena sp., and the model ciliate O. trifallax. We sequenced, assembled, and annotated the germline and somatic genomes of E. woodruffi, which provides an outgroup, and the germline genome of Tetmemena sp. We find that the germline genome of Tetmemena is as massively scrambled and interrupted as Oxytricha's: 13.6% of its gene loci require programmed translocations and/or inversions, with some genes requiring hundreds of precise gene editing events during development. This study revealed that the earlier diverged spirotrich, E. woodruffi, also has a scrambled genome, but only roughly half as many loci (7.3%) are scrambled. Furthermore, its scrambled genes are less complex, together supporting the position of Euplotes as a possible evolutionary intermediate in this lineage, in the process of accumulating complex evolutionary genome rearrangements, all of which require extensive repair to assemble functional coding regions. Comparative analysis also reveals that scrambled loci are often associated with local duplications, supporting a gradual model for the origin of complex, scrambled genomes via many small events of DNA duplication and decay.

Transposable element expansion and low-level piRNA silencing in grasshoppers may cause genome gigantism

Background

Transposable elements (TEs) have been likened to parasites in the genome that reproduce and move ceaselessly in the host, continuously enlarging the host genome. However, the Piwi-interacting RNA (piRNA) pathway defends animal genomes against the harmful consequences of TE invasion by imposing small-RNA-mediated silencing. Here we compare the TE activity of two grasshopper species with different genome sizes in Acrididae (Locusta migratoria manilensis♀1C = 6.60 pg, Angaracris rhodopa♀1C = 16.36 pg) to ascertain the influence of piRNAs.

Results

We discovered that repetitive sequences accounted for 74.56% of the genome in A. rhodopa, more than 56.83% in L. migratoria, and the large-genome grasshopper contained a higher TEs proportions. The comparative analysis revealed that 41 TEs (copy number > 500) were shared in both species. The two species exhibited distinct “landscapes” of TE divergence. The TEs outbreaks in the small-genome grasshopper occurred at more ancient times, while the large-genome grasshopper maintains active transposition events in the recent past. Evolutionary history studies on TEs suggest that TEs may be subject to different dynamics and resistances in these two species. We found that TE transcript abundance was higher in the large-genome grasshopper and the TE-derived piRNAs abundance was lower than in the small-genome grasshopper. In addition, we found that the piRNA methylase HENMT, which is underexpressed in the large-genome grasshopper, impedes the piRNA silencing to a lower level.

Conclusions

Our study revealed that the abundance of piRNAs is lower in the gigantic genome grasshopper than in the small genome grasshopper. In addition, the key gene HENMT in the piRNA biogenesis pathway (Ping-Pong cycle) in the gigantic genome grasshopper is underexpressed. We hypothesize that low-level piRNA silencing unbalances the original positive correlation between TEs and piRNAs, and triggers TEs to proliferate out of control, which may be one of the reasons for the gigantism of grasshopper genomes.

Species-specific chromatin landscape determines how transposable elements shape genome evolution

Transposable elements (TEs) are selfish genetic parasites that increase their copy number at the expense of host fitness. The ‘success’, or genome-wide abundance, of TEs differs widely between species. Deciphering the causes for this large variety in TE abundance has remained a central question in evolutionary genomics. We previously proposed that species-specific TE abundance could be driven by the inadvertent consequences of host-direct epigenetic silencing of TEs—the spreading of repressive epigenetic marks from silenced TEs into adjacent sequences. Here, we compared this TE-mediated local enrichment of repressive marks, or ‘the epigenetic effect of TEs’, in six species in the Drosophila melanogaster subgroup to dissect step-by-step the role of such effect in determining genomic TE abundance. We found that TE-mediated local enrichment of repressive marks is prevalent and substantially varies across and even within species. While this TE-mediated effect alters the epigenetic states of adjacent genes, we surprisingly discovered that the transcription of neighboring genes could reciprocally impact this spreading. Importantly, our multi-species analysis provides the power and appropriate phylogenetic resolution to connect species-specific host chromatin regulation, TE-mediated epigenetic effects, the strength of natural selection against TEs, and genomic TE abundance unique to individual species. Our findings point toward the importance of host chromatin landscapes in shaping genome evolution through the epigenetic effects of a selfish genetic parasite.

All the instructions required for life are encoded in the set of DNA present in a cell. It therefore seems natural to think that every bit of this genetic information should serve the organism. And yet most species carry parasitic ‘transposable’ sequences, or transposons, whose only purpose is to multiply and insert themselves at other positions in the genome.

It is possible for cells to suppress these selfish elements. Chemical marks can be deposited onto the DNA to temporarily ‘silence’ transposons and prevent them from being able to move and replicate. However, this sometimes comes at a cost: the repressive chemical modifications can spread to nearby genes that are essential for the organism and perturb their function.

Strangely, the prevalence of transposons varies widely across the tree of life. These sequences form the majority of the genome of certain species – in fact, they represent about half of the human genetic information. But their abundance is much lower in other organisms, forming a measly 6% of the genome of puffer fish for instance. Even amongst fruit fly species, the prevalence of transposable elements can range between 2% and 25%. What explains such differences?

Huang et al. set out to examine this question through the lens of transposon silencing, systematically comparing how this process impacts nearby regions in six species of fruit flies. This revealed variations in the strength of the side effects associated with transposon silencing, resulting in different levels of perturbation on neighbouring genes. A stronger impact was associated with the species having fewer transposons in its genome, suggesting that an evolutionary pressure is at work to keep the abundance of transposons at a low level in these species. Further analyses showed that the genes which determine how silencing marks are distributed may also be responsible for the variations in the impact of transposon silencing. They could therefore be the ones driving differences in the abundance of transposons between species.

Overall, this work sheds light on the complex mechanisms shaping the evolution of genomes, and it may help to better understand how transposons are linked to processes such as aging and cancer.

Repeat Dynamics across Timescales: A Perspective from Sibling Allotetraploid Marsh Orchids (Dactylorhiza majalis s.l.)

To provide insights into the fate of transposable elements (TEs) across timescales in a post-polyploidization context, we comparatively investigate five sibling Dactylorhiza allotetraploids (Orchidaceae) formed independently and sequentially between 500 and 100K generations ago by unidirectional hybridization between diploids D. fuchsii and D. incarnata. Our results first reveal that the paternal D. incarnata genome shows a marked increased content of LTR retrotransposons compared to the maternal species, reflected in its larger genome size and consistent with a previously hypothesized bottleneck. With regard to the allopolyploids, in the youngest D. purpurella both genome size and TE composition appear to be largely additive with respect to parents, whereas for polyploids of intermediate ages we uncover rampant genome expansion on a magnitude of multiple entire genomes of some plants such as Arabidopsis. The oldest allopolyploids in the series are not larger than the intermediate ones. A putative tandem repeat, potentially derived from a non-autonomous miniature inverted-repeat TE (MITE) drives much of the genome dynamics in the allopolyploids. The highly dynamic MITE-like element is found in higher proportions in the maternal diploid, D. fuchsii, but is observed to increase in copy number in both subgenomes of the allopolyploids. Altogether, the fate of repeats appears strongly regulated and therefore predictable across multiple independent allopolyploidization events in this system. Apart from the MITE-like element, we consistently document a mild genomic shock following the allopolyploidizations investigated here, which may be linked to their relatively large genome sizes, possibly associated with strong selection against further genome expansions.

Patterns of genome size variation in caridean shrimps: new estimates for non-gambarelloides Synalpheus species

Genome size (GS) or DNA nuclear content is considered a useful index for making inferences about evolutionary models and life history in animals, including taxonomic, biogeographical, and ecological scenarios. However, patterns of GS variation and their causes in crustaceans are still poorly understood. This study aimed to describe the GS of five Neotropical Synalpheus nongambarelloides shrimps (S. apioceros, S. minus, S. brevicarpus, S. fritzmueller, and S. scaphoceris) and compare the C-values of all Caridea Infraorder in terms of geography and phylogenetics. All animals were sampled in the coast of São Paulo State, Brazil and GS was assessed by flow cytometry analysis (FCA). The C-values ranged from 7.89 pg in S. apioceros to 12.24 pg in S. scaphoceris. Caridean shrimps had higher GS than other Decapoda crustaceans. The results reveal a tendency of obtaining larger genomes in species with direct development in Synalpheus shrimps. In addition, a tendency of positive biogeographical (latitudinal) correlation with Caridea Infraorder was also observed. This study provides remarkable and new protocol for FCA (using gating strategy for the analysis), which led to the discovery of new information regarding GS of caridean shrimps, especially for Neotropical Synalpheus, which represents the second-largest group in the Caridea Infraorder.

The origin of RNA interference: Adaptive or neutral evolution?

The origin of RNA interference (RNAi) is usually explained by a defense-based hypothesis, in which RNAi evolved as a defense against transposable elements (TEs) and RNA viruses and was already present in the last eukaryotic common ancestor (LECA). However, since RNA antisense regulation and double-stranded RNAs (dsRNAs) are ancient and widespread phenomena, the origin of defensive RNAi should have occurred in parallel with its regulative functions to avoid imbalances in gene regulation. Thus, we propose a neutral evolutionary hypothesis for the origin of RNAi in which qualitative system drift from a prokaryotic antisense RNA gene regulation mechanism leads to the formation of RNAi through constructive neutral evolution (CNE). We argue that RNAi was already present in the ancestor of LECA before the need for a new defense system arose and that its presence helped to shape eukaryotic genomic architecture and stability.

Where does RNA interference come from? This Essay describes a new step-by-step evolutionary model of how RNA interference might have originated in early eukaryotes through neutral events from the molecular machinery present in prokaryotes.

Transposable Element Diversity Remains High in Gigantic Genomes

Transposable elements (TEs) are repetitive sequences of DNA that replicate and proliferate throughout genomes. Taken together, all the TEs in a genome form a diverse community of sequences, which can be studied to draw conclusions about genome evolution. TE diversity can be measured using models for ecological community diversity that consider species richness and evenness. Several models predict TE diversity decreasing as genomes expand because of selection against ectopic recombination and/or competition among TEs to garner host replicative machinery and evade host silencing mechanisms. Salamanders have some of the largest vertebrate genomes and highest TE loads. Salamanders of the genus Plethodon, in particular, have genomes that range in size from 20 to 70 Gb. Here, we use Oxford Nanopore sequencing to generate low-coverage genomic sequences for four species of Plethodon that encompass two independent genome expansion events, one in the eastern clade (Plethodon cinereus, 29.3 Gb vs. Plethodon glutinosus, 38.9 Gb) and one in the western clade (Plethodon vehiculum, 46.4 Gb vs Plethodon idahoensis, 67.0 Gb). We classified the TEs in these genomes and found > 40 TE superfamilies, accounting for 22-27% of the genomes. We calculated Simpson's and Shannon's diversity indices to quantify overall TE diversity. In both pairwise comparisons, the diversity index values for the smaller and larger genome were almost identical. This result indicates that, when genomes reach extremely large sizes, they maintain high levels of TE diversity at the superfamily level, in contrast to predictions made by previous studies on smaller genomes.

Annotation-free delineation of prokaryotic homology groups

Phylogenomic studies of prokaryotic taxa often assume conserved marker genes are homologous across their length. However, processes such as horizontal gene transfer or gene duplication and loss may disrupt this homology by recombining only parts of genes, causing gene fission or fusion. We show using simulation that it is necessary to delineate homology groups in a set of bacterial genomes without relying on gene annotations to define the boundaries of homologous regions. To solve this problem, we have developed a graph-based algorithm to partition a set of bacterial genomes into Maximal Homologous Groups of sequences (MHGs) where each MHG is a maximal set of maximum-length sequences which are homologous across the entire sequence alignment. We applied our algorithm to a dataset of 19 Enterobacteriaceae species and found that MHGs cover much greater proportions of genomes than markers and, relatedly, are less biased in terms of the functions of the genes they cover. We zoomed in on the correlation between each individual marker and their overlapping MHGs, and show that few phylogenetic splits supported by the markers are supported by the MHGs while many marker-supported splits are contradicted by the MHGs. A comparison of the species tree inferred from marker genes with the species tree inferred from MHGs suggests that the increased bias and lack of genome coverage by markers causes incorrect inferences as to the overall relationship between bacterial taxa.

Assuming genes to be the basic evolutionary unit has been commonplace in bacterial genomics. For example, when quantifying the extent of horizontal gene transfer it is common to infer gene trees and reconcile them against a species tree to account for recombination-based processes. We have developed a new method which challenges this assumption by identifying contiguous regions of true homology without regards to gene boundaries and applied it to Enterobacteriaceae, a family of bacteria containing several important human pathogens. Our results show that genes are composed of distinct homologous regions with conflicting phylogenetic histories. We further demonstrate that failing to take account of this conflict, together with the functional biases we show exist among single-copy marker genes, significantly changes the consensus evolutionary tree of Enterobacteriaceae.