The Role of Satellite DNAs in Genome Architecture and Sex Chromosome Evolution in Crambidae Moths

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.

Dominance of transposable element-related satDNAs results in great complexity of "satDNA library" and invokes the extension towards "repetitive DNA library"

Research on bivalves is fast-growing, including genome-wide analyses and genome sequencing. Several characteristics qualify oysters as a valuable model to explore repetitive DNA sequences and their genome organization. Here we characterize the satellitomes of five species in the family Ostreidae (Crassostrea angulata, C. virginica, C. hongkongensis, C. ariakensis, Ostrea edulis), revealing a substantial number of satellite DNAs (satDNAs) per genome (ranging between 33 and 61) and peculiarities in the composition of their satellitomes. Numerous satDNAs were either associated to or derived from transposable elements, displaying a scarcity of transposable element-unrelated satDNAs in these genomes. Due to the non-conventional satellitome constitution and dominance of Helitron-associated satDNAs, comparative satellitomics demanded more in-depth analyses than standardly employed. Comparative analyses (including C. gigas, the first bivalve species with a defined satellitome) revealed that 13 satDNAs occur in all six oyster genomes, with Cg170/HindIII satDNA being the most abundant in all of them. Evaluating the "satDNA library model" highlighted the necessity to adjust this term when studying tandem repeat evolution in organisms with such satellitomes. When repetitive sequences with potential variation in the organizational form and repeat-type affiliation are examined across related species, the introduction of the terms "TE library" and "repetitive DNA library" becomes essential.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-024-00218-0.

Satellitome Analysis in the Southern Lapwing (Vanellus chilensis) Genome: Implications for SatDNA Evolution in Charadriiform Birds

Vanellus (Charadriidae; Charadriiformes) comprises around 20 species commonly referred to as lapwings. In this study, by integrating cytogenetic and genomic approaches, we assessed the satellite DNA (satDNA) composition of one typical species, Vanellus chilensis, with a highly conserved karyotype. We additionally underlined its role in the evolution, structure, and differentiation process of the present ZW sex chromosome system. Seven distinct satellite DNA families were identified within its genome, accumulating on the centromeres, microchromosomes, and the W chromosome. However, these identified satellite DNA families were not found in two other Charadriiformes members, namely Jacana jacana and Calidris canutus. The hybridization of microsatellite sequences revealed the presence of a few repetitive sequences in V. chilensis, with only two out of sixteen displaying positive hybridization signals. Overall, our results contribute to understanding the genomic organization and satDNA evolution in Charadriiform birds.

The spread of satellite DNAs in euchromatin and insights into the multiple sex chromosome evolution in Hemiptera revealed by repeatome analysis of the bug Oxycarenus hyalinipennis

Satellite DNAs (satDNAs) are highly repeated tandem sequences primarily located in heterochromatin, although their occurrence in euchromatin has been reported. Here, our aim was to advance the understanding of satDNA and multiple sex chromosome evolution in heteropterans. We combined cytogenetic and genomic approaches to study, for the first time, the satDNA composition of the genome in an Oxycarenidae bug, Oxycarenus hyalinipennis. The species exhibits a male karyotype of 2n = 19 (14A + 2 m + X1 X2 Y), with a highly differentiated Y chromosome, as demonstrated by C-banding and comparative genomic hybridization, revealing an enrichment of repeats from the male genome. Additionally, comparative analysis between males and females revealed that the 26 identified satDNA families are significantly biased towards male genome, accumulating in discrete regions in the Y chromosome. Exceptionally, the OhyaSat04-125 family was found to be distributed virtually throughout the entire extension of the Y chromosome. This suggests an important role of satDNA in Y chromosome differentiation, in comparison of other repeats, which collectively shows similar abundance between sexes, about 50%. Furthermore, chromosomal mapping of all satDNA families revealed an unexpected high spread in euchromatic regions, covering the entire extension, irrespective of their abundance. Only discrete regions of heterochromatin on the Y chromosome and of the m-chromosomes (peculiar chromosomes commonly observed in heteropterans) were enriched with satDNAs. The putative causes of the intense enrichment of satDNAs in euchromatin are discussed, including the possible existence of burst cycles similar to transposable elements and as a result of holocentricity. These data challenge the classical notion that euchromatin is not enriched with satDNAs.

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.

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.

Holocentric repeat landscapes: From micro-evolutionary patterns to macro-evolutionary associations with karyotype evolution

Repetitive elements can cause large-scale chromosomal rearrangements, for example through ectopic recombination, potentially promoting reproductive isolation and speciation. Species with holocentric chromosomes, that lack a localized centromere, might be more likely to retain chromosomal rearrangements that lead to karyotype changes such as fusions and fissions. This is because chromosome segregation during cell division should be less affected than in organisms with a localized centromere. The relationships between repetitive elements and chromosomal rearrangements and how they may translate to patterns of speciation in holocentric organisms are though poorly understood. Here, we use a reference-free approach based on low-coverage short-read sequencing data to characterize the repeat landscape of two independently evolved holocentric groups: Erebia butterflies and Carex sedges. We consider both micro- and macro-evolutionary scales to investigate the repeat landscape differentiation between Erebia populations and the association between repeats and karyotype changes in a phylogenetic framework for both Erebia and Carex. At a micro-evolutionary scale, we found population differentiation in repeat landscape that increases with overall intraspecific genetic differentiation among four Erebia species. At a macro-evolutionary scale, we found indications for an association between repetitive elements and karyotype changes along both Erebia and Carex phylogenies. Altogether, our results suggest that repetitive elements are associated with the level of population differentiation and chromosomal rearrangements in holocentric clades and therefore likely play a role in adaptation and potentially species diversification.

The Role of Repetitive Sequences in Repatterning of Major Ribosomal DNA Clusters in Lepidoptera

Genes for major ribosomal RNAs (rDNA) are present in multiple copies mainly organized in tandem arrays. The number and position of rDNA loci can change dynamically and their repatterning is presumably driven by other repetitive sequences. We explored a peculiar rDNA organization in several representatives of Lepidoptera with either extremely large or numerous rDNA clusters. We combined molecular cytogenetics with analyses of second- and third-generation sequencing data to show that rDNA spreads as a transcription unit and reveal association between rDNA and various repeats. Furthermore, we performed comparative long read analyses among the species with derived rDNA distribution and moths with a single rDNA locus, which is considered ancestral. Our results suggest that satellite arrays, rather than mobile elements, facilitate homology-mediated spread of rDNA via either integration of extrachromosomal rDNA circles or ectopic recombination. The latter arguably better explains preferential spread of rDNA into terminal regions of lepidopteran chromosomes as efficiency of ectopic recombination depends on the proximity of homologous sequences to telomeres.

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.

Accumulation of retrotransposons contributes to W chromosome differentiation in the willow beauty Peribatodes rhomboidaria (Lepidoptera: Geometridae)

The W chromosome of Lepidoptera is typically gene-poor, repeat-rich and composed of heterochromatin. Pioneering studies investigating this chromosome reported an abundance of mobile elements. However, the actual composition of the W chromosome varies greatly between species, as repeatedly demonstrated by comparative genomic hybridization (CGH) or fluorescence in situ hybridization (FISH). Here we present an analysis of repeats on the W chromosome in the willow beauty, Peribatodes rhomboidaria (Geometridae), a species in which CGH predicted an abundance of W-enriched or W-specific sequences. Indeed, comparative analysis of male and female genomes using RepeatExplorer identified ten putative W chromosome-enriched repeats, most of which are LTR or LINE mobile elements. We analysed the two most abundant: PRW LINE-like and PRW Bel-Pao. The results of FISH mapping and bioinformatic analysis confirmed their enrichment on the W chromosome, supporting the hypothesis that mobile elements are the driving force of W chromosome differentiation in Lepidoptera. As the W chromosome is highly underrepresented in chromosome-level genome assemblies of Lepidoptera, this recently introduced approach, combining bioinformatic comparative genome analysis with molecular cytogenetics, provides an elegant tool for studying this elusive and rapidly evolving part of the genome.

A step forward in the genome characterization of the sugarcane borer, Diatraea saccharalis: karyotype analysis, sex chromosome system and repetitive DNAs through a cytogenomic approach

Moths of the family Crambidae include a number of pests that cause economic losses to agricultural crops. Despite their economic importance, little is known about their genome architecture and chromosome evolution. Here, we characterized the chromosomes and repetitive DNA of the sugarcane borer Diatraea saccharalis using a combination of low-pass genome sequencing, bioinformatics, and cytogenetic methods, focusing on the sex chromosomes. Diploid chromosome numbers differed between the sexes, i.e., 2n = 33 in females and 2n = 34 in males. This difference was caused by the occurrence of a WZ₁Z₂ trivalent in female meiosis, indicating a multiple sex-chromosome system WZ₁Z₂/Z₁Z₁Z₂Z₂. A strong interstitial telomeric signal was observed on the W chromosome, indicating a fusion of the ancestral W chromosome with an autosome. Among repetitive DNAs, transposable elements (TEs) accounted for 39.18% (males) to 41.35% (females), while satDNAs accounted for only 0.214% (males) and 0.215% (females) of the genome. FISH mapping revealed different chromosomal organization of satDNAs, such as single localized clusters, spread repeats, and non-clustered repeats. Two TEs mapped by FISH were scattered. Although we found a slight enrichment of some satDNAs in the female genome, they were not differentially enriched on the W chromosome. However, we found enriched FISH signals for TEs on the W chromosome, suggesting their involvement in W chromosome degeneration and differentiation. These data shed light on karyotype and repetitive DNA dynamics due to multiple chromosome fusions in D. saccharalis, contribute to the understanding of genome structure in Lepidoptera and are important for future genomic studies.

A chromosome-level genome assembly of Paracymoriza distinctalis (Lepidoptera: Crambidae: Acentropinae)

Paracymoriza distinctalis is a semiaquatic lepidopteran insect, which is of great value for studying the differentiation of the Pyraloidea superfamily. However, the understanding of heredity, evolution, and functional genomics of P. distinctalis are limited by few genome-wide resources. Here, we applied PacBio sequencing and the chromosome capture technique to assemble the first P. distinctalis genome from a single female individual. The genome size is 1.2 Gb with 32 chromosomes and the N50 is 38.91 Mb. Approximately 576.37 Mb, accounting for 48.93% of the genome, was identified as repeats. The genome comprises 39,003 protein-coding genes, 66.56% of which were functionally annotated. Comparative genomics analysis suggested that the common ancestor of P. distinctalis and Chilo suppressalis lived ~83.5 million years ago. This chromosome-level genome assembly work is not only conducive to the understanding of P. distinctalis, but also may promote the study of the genomes of other lepidopteran species.

Satellitome Analysis and Transposable Elements Comparison in Geographically Distant Populations of Spodoptera frugiperda

Spodoptera frugiperda (fall armyworm) is a member of the superfamily Noctuoidea that accounts for more than a third of all Lepidoptera and includes a considerable number of agricultural and forest pest species. Spodoptera frugiperda is a polyphagous species that is a significant agricultural pest worldwide, emphasizing its economic importance. Spodoptera frugiperda's genome size, assembly, phylogenetic classification, and transcriptome analysis have all been previously described. However, the different studies reported different compositions of repeated DNA sequences that occupied the whole assembled genome, and the Spodoptera frugiperda genome also lacks the comprehensive study of dynamic satellite DNA. We conducted a comparative analysis of repetitive DNA across geographically distant populations of Spodoptera frugiperda, particularly satellite DNA, using publicly accessible raw genome data from eight different geographical regions. Our results showed that most transposable elements (TEs) were commonly shared across all geographically distant samples, except for the Maverick and PIF/Harbinger elements, which have divergent repeat copies. The TEs age analysis revealed that most TEs families consist of young copies 1-15 million years old; however, PIF/Harbinger has some older/degenerated copies of 30-35 million years old. A total of seven satellite DNA families were discovered, accounting for approximately 0.65% of the entire genome of the Spodoptera frugiperda fall armyworm. The repeat profiling analysis of satellite DNA families revealed differential read depth coverage or copy numbers. The satellite DNA families range in size from the lowest 108 bp SfrSat06-108 families to the largest (1824 bp) SfrSat07-1824 family. We did not observe a statistically significant correlation between monomer length and K2P divergence, copy number, or abundance of each satellite family. Our findings suggest that the satellite DNA families identified in Spodoptera frugiperda account for a considerable proportion of the genome's repetitive fraction. The satellite DNA families' repeat profiling revealed a point mutation along the reference sequences. Limited TEs differentiation exists among geographically distant populations of Spodoptera frugiperda.

Satellitome analysis illuminates the evolution of ZW sex chromosomes of Triportheidae fishes (Teleostei: Characiformes)

Satellites are an abundant source of repetitive DNAs that play an essential role in the chromosomal organization and are tightly linked with the evolution of sex chromosomes. Among fishes, Triportheidae stands out as the only family where almost all species have a homeologous ZZ/ZW sex chromosomes system. While the Z chromosome is typically conserved, the W is always smaller, with variations in size and morphology between species. Here, we report an analysis of the satellitome of Triportheus auritus (TauSat) by integrating genomic and chromosomal data, with a special focus on the highly abundant and female-biased satDNAs. In addition, we investigated the evolutionary trajectories of the ZW sex chromosomes in the Triportheidae family by mapping satDNAs in selected representative species of this family. The satellitome of T. auritus comprised 53 satDNA families of which 24 were also hybridized by FISH. Most satDNAs differed significantly between sexes, with 19 out of 24 being enriched on the W chromosome of T. auritus. The number of satDNAs hybridized into the W chromosomes of T. signatus and T. albus decreased to six and four, respectively, in accordance with the size of their W chromosomes. No TauSat probes produced FISH signals on the chromosomes of Agoniates halecinus. Despite its apparent conservation, our results indicate that each species differs in the satDNA accumulation on the Z chromosome. Minimum spanning trees (MSTs), generated for three satDNA families with different patterns of FISH mapping data, revealed different homogenization rates between the Z and W chromosomes. These results were linked to different levels of recombination between them. The most abundant satDNA family (TauSat01) was exclusively hybridized in the centromeres of all 52 chromosomes of T. auritus, and its putative role in the centromere evolution was also highlighted. Our results identified a high differentiation of both ZW chromosomes regarding satellites composition, highlighting their dynamic role in the sex chromosomes evolution.

Genomic Differences Between the Sexes in a Fish Species Seen Through Satellite DNAs

Neotropical fishes have highly diversified karyotypic and genomic characteristics and present many diverse sex chromosome systems, with various degrees of sex chromosome differentiation. Knowledge on their sex-specific composition and evolution, however, is still limited. Satellite DNAs (satDNAs) are tandemly repeated sequences with pervasive genomic distribution and distinctive evolutionary pathways, and investigating satDNA content might shed light into how genome architecture is organized in fishes and in their sex chromosomes. The present study investigated the satellitome of Megaleporinus elongatus, a freshwater fish with a proposed Z₁Z₁Z₂Z₂/Z₁W₁Z₂W₂ multiple sex chromosome system that encompasses a highly heterochromatic and differentiated W₁ chromosome. The species satellitome comprises of 140 different satDNA families, including previously isolated sequences and new families found in this study. This diversity is remarkable considering the relatively low proportion that satDNAs generally account for the M. elongatus genome (around only 5%). Differences between the sexes in regards of satDNA content were also evidenced, as these sequences are 14% more abundant in the female genome. The occurrence of sex-biased signatures of satDNA evolution in the species is tightly linked to satellite enrichment associated with W₁ in females. Although both sexes share practically all satDNAs, the overall massive amplification of only a few of them accompanied the W₁ differentiation. We also investigated the expansion and diversification of the two most abundant satDNAs of M. elongatus, MelSat01-36 and MelSat02-26, both highly amplified sequences in W₁ and, in MelSat02-26’s case, also harbored by Z₂ and W₂ chromosomes. We compared their occurrences in M. elongatus and the sister species M. macrocephalus (with a standard ZW sex chromosome system) and concluded that both satDNAs have led to the formation of highly amplified arrays in both species; however, they formed species-specific organization on female-restricted sex chromosomes. Our results show how satDNA composition is highly diversified in M. elongatus, in which their accumulation is significantly contributing to W₁ differentiation and not satDNA diversity per se. Also, the evolutionary behavior of these repeats may be associated with genome plasticity and satDNA variability between the sexes and between closely related species, influencing how seemingly homeologous heteromorphic sex chromosomes undergo independent satDNA evolution.

Satellitome Analysis of Rhodnius prolixus, One of the Main Chagas Disease Vector Species

The triatomine Rhodnius prolixus is the main vector of Chagas disease in countries such as Colombia and Venezuela, and the first kissing bug whose genome has been sequenced and assembled. In the repetitive genome fraction (repeatome) of this species, the transposable elements represented 19% of R. prolixus genome, being mostly DNA transposon (Class II elements). However, scarce information has been published regarding another important repeated DNA fraction, the satellite DNA (satDNA), or satellitome. Here, we offer, for the first time, extended data about satellite DNA families in the R. prolixus genome using bioinformatics pipeline based on low-coverage sequencing data. The satellitome of R. prolixus represents 8% of the total genome and it is composed by 39 satDNA families, including four satDNA families that are shared with Triatoma infestans, as well as telomeric (TTAGG)n and (GATA)n repeats, also present in the T. infestans genome. Only three of them exceed 1% of the genome. Chromosomal hybridization with these satDNA probes showed dispersed signals over the euchromatin of all chromosomes, both in autosomes and sex chromosomes. Moreover, clustering analysis revealed that most abundant satDNA families configured several superclusters, indicating that R. prolixus satellitome is complex and that the four most abundant satDNA families are composed by different subfamilies. Additionally, transcription of satDNA families was analyzed in different tissues, showing that 33 out of 39 satDNA families are transcribed in four different patterns of expression across samples.