OUP accepted manuscript

Chromosome-level assembly of Lindenbergia philippensis and comparative genomic analyses shed light on genome evolution in Lamiales

Lamiales, comprising over 23,755 species across 24 families, stands as a highly diverse and prolific plant group, playing a significant role in the cultivation of horticultural, ornamental, and medicinal plant varieties. Whole-genome duplication (WGD) and its subsequent post-polyploid diploidization (PPD) process represent the most drastic type of karyotype evolution, injecting significant potential for promoting the diversity of this lineage. However, polyploidization histories, as well as genome and subgenome fractionation following WGD events in Lamiales species, are still not well investigated. In this study, we constructed a chromosome-level genome assembly of Lindenbergia philippensis (Orobanchaceae) and conducted comparative genomic analyses with 14 other Lamiales species. L. philippensis is positioned closest to the parasitic lineage within Orobanchaceae and has a conserved karyotype. Through a combination of Ks analysis and syntenic depth analysis, we reconstructed and validated polyploidization histories of Lamiales species. Our results indicated that Primulina huaijiensis underwent three rounds of diploidization events following the γ-WGT event, rather than two rounds as reported. Besides, we reconfirmed that most Lamiales species shared a common diploidization event (L-WGD). Subsequently, we constructed the Lamiales Ancestral Karyotype (LAK), comprising 11 proto-chromosomes, and elucidated its evolutionary trajectory, highlighting the highly flexible reshuffling of the Lamiales paleogenome. We identified biased fractionation of subgenomes following the L-WGD event across eight species, and highlighted the positive impacts of non-WGD genes on gene family expansion. This study provides novel genomic resources and insights into polyploidy and karyotype remodeling of Lamiales species, essential for advancing our understanding of species diversification and genome evolution.

Evolutionary Dynamics of Satellite DNA Repeats across the Tettigoniidae Family: Insights from Genomic Analysis

Satellite DNA repeats are repetitive DNA sequences found in eukaryotic genomes, typically consisting of short DNA motifs repeated in tandem arrays. Despite the vast body of literature on satellite DNA repeats in other taxa, investigations specifically targeting Tettigoniidae remain conspicuously absent. Our study aims to fill a critical gap in our understanding of satellitome evolutionary processes shaping Tettigoniidae genomes. Repeatome analysis revealed that the Meconema thalassinum genome comprises 92%, and Phryganogryllacris superangulata had the lowest value of 34%, with an average of 67% in other Tettigoniidae species. The analysis reveals significant variation in the number of satellite DNA repeats across species of the Tettigoniidae family, with M. thalassinum exhibiting the highest count, 246, reported in insects to date and the lowest count, 10, in Pholidoptera griseoptera. Ruspolia dubia and Ruspolia yunnana, which are congeneric species, showcase distinct counts of 104 and 84 families, respectively. Satellite DNA repeats in R. dubia exhibit the highest abundance, constituting 17.2% of the total genome, while the lowest abundance was reported in P. griseoptera, at 5.65%. The genome size correlates weakly with the satellite DNA family count (rs = 0.42, p = 0.29), but a strong correlation exists between satellite abundance and family number (rs = 0.73, p = 0.03). Moreover, the analysis of satellite DNA gain and loss patterns provides insights into the amplification and homogenization of satellite DNA families within the genome, with species-specific repeats exhibiting a positive trend toward amplification. The chromosomal distribution in M. thalassinum displayed that the highest accumulation was observed on Chr12, Chr01, and Chr04, constituting 17.79%, 17.4%, and 17.22% of the total chromosome size, respectively. The chromosome-specific propagation of satellite DNA families was evident, with MthSat01 solely on chromosome 1 and MthSat170 on chromosome 2, sharing 1.64% and 2.33%. The observed conservation and variations in satellite DNA number and abundances, along with distinct patterns of gain and loss, indicate the influence of potentially diverse evolutionary processes shaping the genomic landscape of these insects, which requires further investigation. Furthermore, the differential accumulation of satellite DNA on specific chromosomes implies that potential chromosome-specific functions or structural features influence the retention and proliferation of satellite sequences.

Prevalent fast evolution of genes involved in heterochromatin functions

Heterochromatin is a gene-poor and repeat-rich genomic compartment universally found in eukaryotes. Despite its low transcriptional activity, heterochromatin plays important roles in maintaining genome stability, organizing chromosomes, and suppressing transposable elements (TEs). Given the importance of these functions, it is expected that the genes involved in heterochromatin regulation would be highly conserved. Yet, a handful of these genes were found to evolve rapidly. To investigate whether these previous findings are anecdotal or general to genes modulating heterochromatin, we compile an exhaustive list of 106 candidate genes involved in heterochromatin functions and investigate their evolution over short and long evolutionary time scales in Drosophila. Our analyses find that these genes exhibit significantly more frequent evolutionary changes, both in the forms of amino acid substitutions and gene copy number change, when compared to genes involved in Polycomb-based repressive chromatin. While positive selection drives amino acid changes within both structured domains with diverse functions and intrinsically disordered regions (IDRs), purifying selection may have maintained the proportions of IDRs of these proteins. Together with the observed negative associations between evolutionary rates of these genes and genomic TE abundance, we propose an evolutionary model where the fast evolution of genes involved in heterochromatin functions is an inevitable outcome of the unique functional roles of heterochromatin, while the rapid evolution of TEs may be an effect rather than cause. Our study provides an important global view of the evolution of genes involved in this critical cellular domain and provides insights into the factors driving the distinctive evolution of heterochromatin.

Comprehensive analysis of the Xya riparia genome uncovers the dominance of DNA transposons, LTR/Gypsy elements, and their evolutionary dynamics

Transposable elements (TEs) are DNA sequences that can move or replicate within a genome, and their study has become increasingly important in understanding genome evolution and function. The Tridactylidae family, including Xya riparia (pygmy mole cricket), harbors a variety of transposable elements (TEs) that have been insufficiently investigated. Further research is required to fully understand their diversity and evolutionary characteristics. Hence, we conducted a comprehensive repeatome analysis of X. riparia species using the chromosome-level assembled genome. The study aimed to comprehensively analyze the abundance, distribution, and age of transposable elements (TEs) in the genome. The results indicated that the genome was 1.67 Gb, with 731.63 Mb of repetitive sequences, comprising 27% of Class II (443.25 Mb), 16% of Class I (268.45 Mb), and 1% of unknown TEs (19.92 Mb). The study found that DNA transposons dominate the genome, accounting for approximately 60% of the total repeat size, with retrotransposons and unknown elements accounting for 37% and 3% of the genome, respectively. The members of the Gypsy superfamily were the most abundant amongst retrotransposons, accounting for 63% of them. The transposable superfamilies (LTR/Gypsy, DNA/nMITE, DNA/hAT, and DNA/Helitron) collectively constituted almost 70% of the total repeat size of all six chromosomes. The study further unveiled a significant linear correlation (Pearson correlation: r = 0.99, p-value = 0.00003) between the size of the chromosomes and the repetitive sequences. The average age of DNA transposon and retrotransposon insertions ranges from 25 My (million years) to 5 My. The satellitome analysis discovered 13 satellite DNA families that comprise about 0.15% of the entire genome. In addition, the transcriptional analysis of TEs found that DNA transposons were more transcriptionally active than retrotransposons. Overall, the study suggests that the genome of X. riparia is complex, characterized by a substantial portion of repetitive elements. These findings not only enhance our understanding of TE evolution within the Tridactylidae family but also provide a foundation for future investigations into the genomic intricacies of related species.

Recurrent Duplication and Diversification of a Vital DNA Repair Gene Family Across Drosophila

Maintaining genome integrity is vital for organismal survival and reproduction. Essential, broadly conserved DNA repair pathways actively preserve genome integrity. However, many DNA repair proteins evolve adaptively. Ecological forces like UV exposure are classically cited drivers of DNA repair evolution. Intrinsic forces like repetitive DNA, which also imperil genome integrity, have received less attention. We recently reported that a Drosophila melanogaster-specific DNA satellite array triggered species-specific, adaptive evolution of a DNA repair protein called Spartan/MH. The Spartan family of proteases cleave hazardous, covalent crosslinks that form between DNA and proteins ("DNA-protein crosslink repair"). Appreciating that DNA satellites are both ubiquitous and universally fast-evolving, we hypothesized that satellite DNA turnover spurs adaptive evolution of DNA-protein crosslink repair beyond a single gene and beyond the D. melanogaster lineage. This hypothesis predicts pervasive Spartan gene family diversification across Drosophila species. To study the evolutionary history of the Drosophila Spartan gene family, we conducted population genetic, molecular evolution, phylogenomic, and tissue-specific expression analyses. We uncovered widespread signals of positive selection across multiple Spartan family genes and across multiple evolutionary timescales. We also detected recurrent Spartan family gene duplication, divergence, and gene loss. Finally, we found that ovary-enriched parent genes consistently birthed functionally diverged, testis-enriched daughter genes. To account for Spartan family diversification, we introduce a novel mechanistic model of antagonistic coevolution that links DNA satellite evolution and adaptive regulation of Spartan protease activity. This framework promises to accelerate our understanding of how DNA repeats drive recurrent evolutionary innovation to preserve genome integrity.

The implications of satellite DNA instability on cellular function and evolution

Abundant tandemly repeated satellite DNA is present in most eukaryotic genomes. Previous limitations including a pervasive view that it was uninteresting junk DNA, combined with challenges in studying it, are starting to dissolve - and recent studies have found important functions for satellite DNAs. The observed rapid evolution and implied instability of satellite DNA now has important significance for their functions and maintenance within the genome. In this review, we discuss the processes that lead to satellite DNA copy number instability, and the importance of mechanisms to manage the potential negative effects of instability. Satellite DNA is vulnerable to challenges during replication and repair, since it forms difficult-to-process secondary structures and its homology within tandem arrays can result in various types of recombination. Satellite DNA instability may be managed by DNA or chromatin-binding proteins ensuring proper nuclear localization and repair, or by proteins that process aberrant structures that satellite DNAs tend to form. We also discuss the pattern of satellite DNA mutations from recent mutation accumulation (MA) studies that have tracked changes in satellite DNA for up to 1000 generations with minimal selection. Finally, we highlight examples of satellite evolution from studies that have characterized satellites across millions of years of Drosophila fruit fly evolution, and discuss possible ways that selection might act on the satellite DNA composition.

Contribution of the satellitome to the exceptionally large genome of dinoflagellates: The case of the harmful alga Alexandrium minutum

Dinoflagellates are known to possess an exceptionally large genome organized in permanently condensed chromosomes. Focusing on the contribution of satellite DNA (satDNA) to the whole DNA content of genomes and its potential role in the architecture of the chromosomes, we present the characterization of the satellitome of Alexandriun minutum strain VGO577. To achieve this, we analyzed Illumina reads using graph-based clustering and performed complementary bioinformatic analyses. In this way, we discovered 180 satDNAs occupying 17.38 % of the genome. The 12 most abundant satDNAs represent the half of the satellitome but no satDNA is overrepresented, with the most abundant contributing ∼1.56 % of the genome. The largest repeat unit is 517 bp long but more than the half of the satDNAs (101) have repeat units shorter than 20 bp. We used FISH to map a selected set of 26 satDNAs. Although some satDNAs generate discrete hybridization signals at specific chromosomal locations (hybridization sites, HS), our cytological analysis showed that most satDNAs are dispersed throughout the genome, probably forming short arrays. Two satDNAs co-localize with the 45S rDNA. With the exception of telomeric DNA, no other satDNA yields HS on all chromosomes. In addition, we analyzed nine satDNAs yielding HS in VGO577 in four other A. minutum strains. Polymorphism at the intraspecific level was found for the presence/absence and/or abundance of some satDNAs, suggesting the amplification/deletion of these satDNAs following geographic separation or during culture maintenance of the strains. We also discuss how these results contribute to the understanding of chromosome architecture and evolution of dinoflagellate genomes.

Fast satellite DNA evolution in Nothobranchius annual killifishes

Satellite DNA (satDNA) is a rapidly evolving class of tandem repeats, with some monomers being involved in centromere organization and function. To identify repeats associated with (peri)centromeric regions, we investigated satDNA across Southern and Coastal clades of African annual killifishes of the genus Nothobranchius. Molecular cytogenetic and bioinformatic analyses revealed that two previously identified satellites, designated here as NkadSat01-77 and NfurSat01-348, are associated with (peri)centromeres only in one lineage of the Southern clade. NfurSat01-348 was, however, additionally detected outside centromeres in three members of the Coastal clade. We also identified a novel satDNA, NrubSat01-48, associated with (peri)centromeres in N. foerschi, N. guentheri, and N. rubripinnis. Our findings revealed fast turnover of satDNA associated with (peri)centromeres and different trends in their evolution in two clades of the genus Nothobranchius.

Chromosomal Rearrangements and Satellite DNAs: Extensive Chromosome Reshuffling and the Evolution of Neo-Sex Chromosomes in the Genus Pyrrhulina (Teleostei; Characiformes)

Chromosomal rearrangements play a significant role in the evolution of fish genomes, being important forces in the rise of multiple sex chromosomes and in speciation events. Repetitive DNAs constitute a major component of the genome and are frequently found in heterochromatic regions, where satellite DNA sequences (satDNAs) usually represent their main components. In this work, we investigated the association of satDNAs with chromosome-shuffling events, as well as their potential relevance in both sex and karyotype evolution, using the well-known Pyrrhulina fish model. Pyrrhulina species have a conserved karyotype dominated by acrocentric chromosomes present in all examined species up to date. However, two species, namely P. marilynae and P. semifasciata, stand out for exhibiting unique traits that distinguish them from others in this group. The first shows a reduced diploid number (with 2n = 32), while the latter has a well-differentiated multiple X1X2Y sex chromosome system. In addition to isolating and characterizing the full collection of satDNAs (satellitomes) of both species, we also in situ mapped these sequences in the chromosomes of both species. Moreover, the satDNAs that displayed signals on the sex chromosomes of P. semifasciata were also mapped in some phylogenetically related species to estimate their potential accumulation on proto-sex chromosomes. Thus, a large collection of satDNAs for both species, with several classes being shared between them, was characterized for the first time. In addition, the possible involvement of these satellites in the karyotype evolution of P. marilynae and P. semifasciata, especially sex-chromosome formation and karyotype reduction in P. marilynae, could be shown.

Following the Pathway of W Chromosome Differentiation in Triportheus (Teleostei: Characiformes)

In this work, we trace the dynamics of satellite DNAs (SatDNAs) accumulation and elimination along the pathway of W chromosome differentiation using the well-known Triportheus fish model. Triportheus stands out due to a conserved ZZ/ZW sex chromosome system present in all examined species. While the Z chromosome is conserved in all species, the W chromosome is invariably smaller and exhibits differences in size and morphology. The presumed ancestral W chromosome is comparable to that of T. auritus, and contains 19 different SatDNA families. Here, by examining five additional Triportheus species, we showed that the majority of these repetitive sequences were eliminated as speciation was taking place. The W chromosomes continued degeneration, while the Z chromosomes of some species began to accumulate some TauSatDNAs. Additional species-specific SatDNAs that made up the heterochromatic region of both Z and W chromosomes were most likely amplified in each species. Therefore, the W chromosomes of the various Triportheus species have undergone significant evolutionary changes in a short period of time (15-25 Myr) after their divergence.

Three recent sex chromosome-to-autosome fusions in a Drosophila virilis strain with high satellite DNA content

The karyotype, or number and arrangement of chromosomes, has varying levels of stability across both evolution and disease. Karyotype changes often originate from DNA breaks near the centromeres of chromosomes, which generally contain long arrays of tandem repeats or satellite DNA. Drosophila virilis possesses among the highest relative satellite abundances of studied species, with almost half its genome composed of three related 7 bp satellites. We discovered a strain of D. virilis that we infer recently underwent three independent chromosome fusion events involving the X and Y chromosomes, in addition to one subsequent fission event. Here, we isolate and characterize the four different karyotypes we discovered in this strain which we believe demonstrates remarkable genome instability. We discovered that one of the substrains with an X-autosome fusion has an X-to-Y chromosome nondisjunction rate 20 × higher than the D. virilis reference strain (21% vs 1%). Finally, we found an overall higher rate of DNA breakage in the substrain with higher satellite DNA compared to a genetically similar substrain with less satellite DNA. This suggests that satellite DNA abundance may play a role in the risk of genome instability. Overall, we introduce a novel system consisting of a single strain with four different karyotypes, which we believe will be useful for future studies of genome instability, centromere function, and sex chromosome evolution.

Satellite DNAs-From Localized to Highly Dispersed Genome Components

According to the established classical view, satellite DNAs are defined as abundant non-coding DNA sequences repeated in tandem that build long arrays located in heterochromatin. Advances in sequencing methodologies and development of specialized bioinformatics tools enabled defining a collection of all repetitive DNAs and satellite DNAs in a genome, the repeatome and the satellitome, respectively, as well as their reliable annotation on sequenced genomes. Supported by various non-model species included in recent studies, the patterns of satellite DNAs and satellitomes as a whole showed much more diversity and complexity than initially thought. Differences are not only in number and abundance of satellite DNAs but also in their distribution across the genome, array length, interspersion patterns, association with transposable elements, localization in heterochromatin and/or in euchromatin. In this review, we compare characteristic organizational features of satellite DNAs and satellitomes across different animal and plant species in order to summarize organizational forms and evolutionary processes that may lead to satellitomes' diversity and revisit some basic notions regarding repetitive DNA landscapes in genomes.

The Satellite DNAs Populating the Genome of Trigona hyalinata and the Sharing of a Highly Abundant satDNA in Trigona Genus

Among Meliponini species, c-heterochromatin can occupy large portions of chromosomes. This characteristic could be useful for understanding evolutionary patterns of satellite DNAs (satDNAs), although few sequences have been characterized in these bees. In Trigona, phylogenetically represented by clades A and B, the c-heterochromatin is mostly located in one chromosome arm. Here we used different techniques, including restriction endonucleases and genome sequencing followed by chromosomal analysis, to identify satDNAs that may be contributing to the evolution of c-heterochromatin in Trigona. Our results revealed a highly abundant ThyaSat01-301 satDNA, corresponding to about 13.77% of the Trigona hyalinata genome. Another seven satDNAs were identified, one corresponding to 2.24%, and the other six corresponding to 0.545% of the genome. The satDNA ThyaSat01-301 was shown to be one of the main constituents of the c-heterochromatin of this species, as well as of other species belonging to clade B of Trigona. However, this satDNA was not observed on the chromosomes of species from clade A, demonstrating that the c-heterochromatin is evolving divergently between species of clade A and B, as a consequence of the evolution of repetitive DNA sequences. Finally, our data suggest the molecular diversification of the karyotypes, despite a conservated macrochromosomal structure on the genus.

Making the Genome Huge: The Case of Triatoma delpontei, a Triatominae Species with More than 50% of Its Genome Full of Satellite DNA

The genome of Triatoma delpontei Romaña & Abalos 1947 is the largest within Heteroptera, approximately two to three times greater than other evaluated Heteroptera genomes. Here, the repetitive fraction of the genome was determined and compared with its sister species Triatoma infestans Klug 1834, in order to shed light on the karyotypic and genomic evolution of these species. The T. delpontei repeatome analysis showed that the most abundant component in its genome is satellite DNA, which makes up more than half of the genome. The T. delpontei satellitome includes 160 satellite DNA families, most of them also present in T. infestans. In both species, only a few satellite DNA families are overrepresented on the genome. These families are the building blocks of the C-heterochromatic regions. Two of these satellite DNA families that form the heterochromatin are the same in both species. However, there are satellite DNA families highly amplified in the heterochromatin of one species that in the other species are in low abundance and located in the euchromatin. Therefore, the present results depicted the great impact of the satellite DNA sequences in the evolution of Triatominae genomes. Within this scenario, satellitome determination and analysis led to a hypothesis that explains how satDNA sequences have grown on T. delpontei to reach its huge genome size within true bugs.

Evolution of satDNAs on holocentric chromosomes: insights from hemipteran insects of the genus Mahanarva

Satellite DNAs (satDNAs) constitute one of the main components of eukaryote genomes and are involved in chromosomal organization and diversification. Although largely studied, little information was gathered about their evolution on holocentric species, i.e., diffuse centromeres, which, due to differences in repeat organization, could result in different evolutionary patterns. Here, we combined bioinformatics and cytogenetic approaches to evaluate the evolution of the satellitomes in Mahanarva holocentric insects. In two species, de novo identification revealed a high number of satDNAs, 110 and 113, with an extreme monomer length range of 18-4228 bp. The overall abundance of satDNAs was observed to be 6.67% in M. quadripunctata and 1.98% in M. spectabilis, with different abundances for the shared satDNAs. Chromosomal mapping of the most abundant repeats of M. quadripunctata and M. spectabilis on other Mahanarva reinforced the dynamic nature of satDNAs. Variable patterns of chromosomal distribution for the satDNAs were noticed, with the occurrence of clusters on distinct numbers of chromosomes and at different positions and the occurrence of scattered signals or nonclustered satDNAs. Altogether, our data demonstrated the high dynamism of satDNAs in Mahanarva with the involvement of this genomic fraction in chromosome diversification of the genus. The general characteristics and patterns of evolution of satDNAs are similar to those observed on monocentric chromosomes, suggesting that the differential organization of genome compartments observed on holocentric chromosomes compared with monocentric chromosomes does not have a large impact on the evolution of satDNAs. Analysis of the satellitomes of other holocentric species in a comparative manner will shed light on this issue.

In Silico Identification and Characterization of Satellite DNAs in 23 Drosophila Species from the Montium Group

Satellite DNA (satDNA) is a class of tandemly repeated non-protein coding DNA sequences which can be found in abundance in eukaryotic genomes. They can be functional, impact the genomic architecture in many ways, and their rapid evolution has consequences for species diversification. We took advantage of the recent availability of sequenced genomes from 23 Drosophila species from the montium group to study their satDNA landscape. For this purpose, we used publicly available whole-genome sequencing Illumina reads and the TAREAN (tandem repeat analyzer) pipeline. We provide the characterization of 101 non-homologous satDNA families in this group, 93 of which are described here for the first time. Their repeat units vary in size from 4 bp to 1897 bp, but most satDNAs show repeat units < 100 bp long and, among them, repeats ≤ 10 bp are the most frequent ones. The genomic contribution of the satDNAs ranges from ~1.4% to 21.6%. There is no significant correlation between satDNA content and genome sizes in the 23 species. We also found that at least one satDNA originated from an expansion of the central tandem repeats (CTRs) present inside a Helitron transposon. Finally, some satDNAs may be useful as taxonomic markers for the identification of species or subgroups within the group.

The Satellite DNA Catalogues of Two Serrasalmidae (Teleostei, Characiformes): Conservation of General satDNA Features over 30 Million Years

Satellite DNAs (satDNAs) are tandemly repeated sequences that are usually located on the heterochromatin, and the entire collection of satDNAs within a genome is called satellitome. Primarily, these sequences are not under selective pressure and evolve by concerted evolution, resulting in elevated rates of divergence between the satDNA profiles of reproductive isolated species/populations. Here, we characterized two additional satellitomes of Characiformes fish (Colossoma macropomum and Piaractus mesopotamicus) that diverged approximately 30 million years ago, while still retaining conserved karyotype features. The results we obtained indicated that several satDNAs (50% of satellite sequences in P. mesopotamicus and 43% in C. macropomum) show levels of conservation between the analyzed species, in the nucleotide and chromosomal levels. We propose that long-life cycles and few genomic changes could slow down rates of satDNA differentiation.