RepeatProfiler: A pipeline for visualization and comparative analysis of repetitive DNA profiles

Satellitome Analysis of Adalia bipunctata (Coleoptera): Revealing Centromeric Turnover and Potential Chromosome Rearrangements in a Comparative Interspecific Study

Eukaryotic genomes exhibit a dynamic interplay between single-copy sequences and repetitive DNA elements, with satellite DNA (satDNA) representing a substantial portion, mainly situated at telomeric and centromeric chromosomal regions. We utilized Illumina next-generation sequencing data from Adalia bipunctata to investigate its satellitome. Cytogenetic mapping via fluorescence in situ hybridization was performed for the most abundant satDNA families. In silico localization of satDNAs was carried out using the CHRISMAPP (Chromosome In Silico Mapping) pipeline on the high-fidelity chromosome-level assembly already available for this species, enabling a meticulous characterization and localization of multiple satDNA families. Additionally, we analyzed the conservation of the satellitome at an interspecific scale. Specifically, we employed the CHRISMAPP pipeline to map the satDNAs of A. bipunctata onto the genome of Adalia decempunctata, which has also been sequenced and assembled at the chromosome level. This analysis, along with the creation of a synteny map between the two species, suggests a rapid turnover of centromeric satDNA between these species and the potential occurrence of chromosomal rearrangements, despite the considerable conservation of their satellitomes. Specific satDNA families in the sex chromosomes of both species suggest a role in sex chromosome differentiation. Our interspecific comparative study can provide a significant advance in the understanding of the repeat genome organization and evolution in beetles.

The genomic and cellular basis of biosynthetic innovation in rove beetles

How evolution at the cellular level potentiates macroevolutionary change is central to understanding biological diversification. The >66,000 rove beetle species (Staphylinidae) form the largest metazoan family. Combining genomic and cell type transcriptomic insights spanning the largest clade, Aleocharinae, we retrace evolution of two cell types comprising a defensive gland-a putative catalyst behind staphylinid megadiversity. We identify molecular evolutionary steps leading to benzoquinone production by one cell type via a mechanism convergent with plant toxin release systems, and synthesis by the second cell type of a solvent that weaponizes the total secretion. This cooperative system has been conserved since the Early Cretaceous as Aleocharinae radiated into tens of thousands of lineages. Reprogramming each cell type yielded biochemical novelties enabling ecological specialization-most dramatically in symbionts that infiltrate social insect colonies via host-manipulating secretions. Our findings uncover cell type evolutionary processes underlying the origin and evolvability of a beetle chemical innovation.

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.

Horizontal transfer and the widespread presence of Galileo transposons in Drosophilidae (Insecta: Diptera)

Galileo is a transposon notoriously involved with inversions in Drosophila buzzatii by ectopic recombination. Although widespread in Drosophila, little is known about this transposon in other lineages of Drosophilidae. Here, the abundance of the canonical Galileo and its evolutionary history in Drosophilidae genomes was estimated and reconstructed across genera within its two subfamilies. Sequences of this transposon were masked in these genomes and their transposase sequences were recovered using BLASTn. Phylogenetic analyses were employed to reconstruct their evolutionary history and compare it to that of host genomes. Galileo was found in nearly all 163 species, however, only 37 harbored nearly complete transposase sequences. In the remaining, Galileo was found highly fragmented. Copies from related species were clustered, however horizontal transfer events were detected between the melanogaster and montium groups of Drosophila, and between the latter and the Lordiphosa genus. The similarity of sequences found in the virilis and willistoni groups of Drosophila was found to be a consequence of lineage sorting. Therefore, the evolution of Galileo is primarily marked by vertical transmission and long-term inactivation, mainly through the deletion of open reading frames. The latter has the potential to lead copies of this transposon to become miniature inverted-repeat transposable elements.

Orthoptera-TElib: a library of Orthoptera transposable elements for TE annotation

Transposable elements (TEs) are a major component of eukaryotic genomes and are present in almost all eukaryotic organisms. TEs are highly dynamic between and within species, which significantly affects the general applicability of the TE databases. Orthoptera is the only known group in the class Insecta with a significantly enlarged genome (0.93-21.48 Gb). When analyzing the large genome using the existing TE public database, the efficiency of TE annotation is not satisfactory. To address this limitation, it becomes imperative to continually update the available TE resource library and the need for an Orthoptera-specific library as more insect genomes are publicly available. Here, we used the complete genome data of 12 Orthoptera species to de novo annotate TEs, then manually re-annotate the unclassified TEs to construct a non-redundant Orthoptera-specific TE library: Orthoptera-TElib. Orthoptera-TElib contains 24,021 TE entries including the re-annotated results of 13,964 unknown TEs. The naming of TE entries in Orthoptera-TElib adopts the same naming as RepeatMasker and Dfam and is encoded as the three-level form of "level1/level2-level3". Orthoptera-TElib can be directly used as an input reference database and is compatible with mainstream repetitive sequence analysis software such as RepeatMasker and dnaPipeTE. When analyzing TEs of Orthoptera species, Orthoptera-TElib performs better TE annotation as compared to Dfam and Repbase regardless of using low-coverage sequencing or genome assembly data. The most improved TE annotation result is Angaracris rhodopa, which has increased from 7.89% of the genome to 53.28%. Finally, Orthoptera-TElib is stored in Sqlite3 for the convenience of data updates and user access.

Analyses of 600+ insect genomes reveal repetitive element dynamics and highlight biodiversity-scale repeat annotation challenges

Repetitive elements (REs) are integral to the composition, structure, and function of eukaryotic genomes, yet remain understudied in most taxonomic groups. We investigated REs across 601 insect species and report wide variation in RE dynamics across groups. Analysis of associations between REs and protein-coding genes revealed dynamic evolution at the interface between REs and coding regions across insects, including notably elevated RE-gene associations in lineages with abundant long interspersed nuclear elements (LINEs). We leveraged this large, empirical data set to quantify impacts of long-read technology on RE detection and investigate fundamental challenges to RE annotation in diverse groups. In long-read assemblies, we detected ∼36% more REs than short-read assemblies, with long terminal repeats (LTRs) showing 162% increased detection, whereas DNA transposons and LINEs showed less respective technology-related bias. In most insect lineages, 25%-85% of repetitive sequences were "unclassified" following automated annotation, compared with only ∼13% in Drosophila species. Although the diversity of available insect genomes has rapidly expanded, we show the rate of community contributions to RE databases has not kept pace, preventing efficient annotation and high-resolution study of REs in most groups. We highlight the tremendous opportunity and need for the biodiversity genomics field to embrace REs and suggest collective steps for making progress toward this goal.

The genomic and cellular basis of biosynthetic innovation in rove beetles

How evolution at the cellular level potentiates change at the macroevolutionary level is a major question in evolutionary biology. With >66,000 described species, rove beetles (Staphylinidae) comprise the largest metazoan family. Their exceptional radiation has been coupled to pervasive biosynthetic innovation whereby numerous lineages bear defensive glands with diverse chemistries. Here, we combine comparative genomic and single-cell transcriptomic data from across the largest rove beetle clade, Aleocharinae. We retrace the functional evolution of two novel secretory cell types that together comprise the tergal gland-a putative catalyst behind Aleocharinae's megadiversity. We identify key genomic contingencies that were critical to the assembly of each cell type and their organ-level partnership in manufacturing the beetle's defensive secretion. This process hinged on evolving a mechanism for regulated production of noxious benzoquinones that appears convergent with plant toxin release systems, and synthesis of an effective benzoquinone solvent that weaponized the total secretion. We show that this cooperative biosynthetic system arose at the Jurassic-Cretaceous boundary, and that following its establishment, both cell types underwent ∼150 million years of stasis, their chemistry and core molecular architecture maintained almost clade-wide as Aleocharinae radiated globally into tens of thousands of lineages. Despite this deep conservation, we show that the two cell types have acted as substrates for the emergence of adaptive, biochemical novelties-most dramatically in symbiotic lineages that have infiltrated social insect colonies and produce host behavior-manipulating secretions. Our findings uncover genomic and cell type evolutionary processes underlying the origin, functional conservation and evolvability of a chemical innovation in beetles.

The nature and genomic landscape of repetitive DNA classes in Chrysanthemum nankingense shows recent genomic changes

BACKGROUND AND AIMS

Tandemly repeated DNA and transposable elements represent most of the DNA in higher plant genomes. High-throughput sequencing allows a survey of the DNA in a genome, but whole-genome assembly can miss a substantial fraction of highly repeated sequence motifs. Chrysanthemum nankingense (2n = 2x = 18; genome size = 3.07 Gb; Asteraceae), a diploid reference for the many auto- and allopolyploids in the genus, was considered as an ancestral species and serves as an ornamental plant and high-value food. We aimed to characterize the major repetitive DNA motifs, understand their structure and identify key features that are shaped by genome and sequence evolution.

METHODS

Graph-based clustering with RepeatExplorer was used to identify and classify repetitive motifs in 2.14 millions of 250-bp paired-end Illumina reads from total genomic DNA of C. nankingense. Independently, the frequency of all canonical motifs k-bases long was counted in the raw read data and abundant k-mers (16, 21, 32, 64 and 128) were extracted and assembled to generate longer contigs for repetitive motif identification. For comparison, long terminal repeat retrotransposons were checked in the published C. nankingense reference genome. Fluorescent in situ hybridization was performed to show the chromosomal distribution of the main types of repetitive motifs.

KEY RESULTS

Apart from rDNA (0.86 % of the total genome), a few microsatellites (0.16 %), and telomeric sequences, no highly abundant tandem repeats were identified. There were many transposable elements: 40 % of the genome had sequences with recognizable domains related to transposable elements. Long terminal repeat retrotransposons showed widespread distribution over chromosomes, although different sequence families had characteristic features such as abundance at or exclusion from centromeric or subtelomeric regions. Another group of very abundant repetitive motifs, including those most identified as low-complexity sequences (9.07 %) in the genome, showed no similarity to known sequence motifs or tandemly repeated elements.

CONCLUSIONS

The Chrysanthemum genome has an unusual structure with a very low proportion of tandemly repeated sequences (~1.02 %) in the genome, and a high proportion of low-complexity sequences, most likely degenerated remains of transposable elements. Identifying the presence, nature and genomic organization of major genome fractions enables inference of the evolutionary history of sequences, including degeneration and loss, critical to understanding biodiversity and diversification processes in the genomes of diploid and polyploid Chrysanthemum, Asteraceae and plants more widely.

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.

Epigenetic Stress and Long-Read cDNA Sequencing of Sunflower (Helianthus annuus L.) Revealed the Origin of the Plant Retrotranscriptome

Transposable elements (TEs) contribute not only to genome diversity but also to transcriptome diversity in plants. To unravel the sources of LTR retrotransposon (RTE) transcripts in sunflower, we exploited a recently developed transposon activation method ('TEgenesis') along with long-read cDNA Nanopore sequencing. This approach allows for the identification of 56 RTE transcripts from different genomic loci including full-length and non-autonomous RTEs. Using the mobilome analysis, we provided a new set of expressed and transpositional active sunflower RTEs for future studies. Among them, a Ty3/Gypsy RTE called SUNTY3 exhibited ongoing transposition activity, as detected by eccDNA analysis. We showed that the sunflower genome contains a diverse set of non-autonomous RTEs encoding a single RTE protein, including the previously described TR-GAG (terminal repeat with the GAG domain) as well as new categories, TR-RT-RH, TR-RH, and TR-INT-RT. Our results demonstrate that 40% of the loci for RTE-related transcripts (nonLTR-RTEs) lack their LTR sequences and resemble conventional eucaryotic genes encoding RTE-related proteins with unknown functions. It was evident based on phylogenetic analysis that three nonLTR-RTEs encode GAG (HadGAG1-3) fused to a host protein. These HadGAG proteins have homologs found in other plant species, potentially indicating GAG domestication. Ultimately, we found that the sunflower retrotranscriptome originated from the transcription of active RTEs, non-autonomous RTEs, and gene-like RTE transcripts, including those encoding domesticated proteins.

Transposon control as a checkpoint for tissue regeneration

Tissue regeneration requires precise temporal control of cellular processes such as inflammatory signaling, chromatin remodeling and proliferation. The combination of these processes forms a unique microenvironment permissive to the expression, and potential mobilization of, transposable elements (TEs). Here, we develop the hypothesis that TE activation creates a barrier to tissue repair that must be overcome to achieve successful regeneration. We discuss how uncontrolled TE activity may impede tissue restoration and review mechanisms by which TE activity may be controlled during regeneration. We posit that the diversification and co-evolution of TEs and host control mechanisms may contribute to the wide variation in regenerative competency across tissues and species.

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.

Genome diploidization associates with cladogenesis, trait disparity, and plastid gene evolution

Angiosperm genome evolution was marked by many clade-specific whole-genome duplication events. The Microlepidieae is one of the monophyletic clades in the mustard family (Brassicaceae) formed after an ancient allotetraploidization. Postpolyploid cladogenesis has resulted in the extant c. 17 genera and 60 species endemic to Australia and New Zealand (10 species). As postpolyploid genome diploidization is a trial-and-error process under natural selection, it may proceed with different intensity and be associated with speciation events. In Microlepidieae, different extents of homoeologous recombination between the two parental subgenomes generated clades marked by slow ("cold") versus fast ("hot") genome diploidization. To gain a deeper understanding of postpolyploid genome evolution in Microlepidieae, we analyzed phylogenetic relationships in this tribe using complete chloroplast sequences, entire 35S rDNA units, and abundant repetitive sequences. The four recovered intra-tribal clades mirror the varied diploidization of Microlepidieae genomes, suggesting that the intrinsic genomic features underlying the extent of diploidization are shared among genera and species within one clade. Nevertheless, even congeneric species may exert considerable morphological disparity (e.g. in fruit shape), whereas some species within different clades experience extensive morphological convergence despite the different pace of their genome diploidization. We showed that faster genome diploidization is positively associated with mean morphological disparity and evolution of chloroplast genes (plastid-nuclear genome coevolution). Higher speciation rates in perennials than in annual species were observed. Altogether, our results confirm the potential of Microlepidieae as a promising subject for the analysis of postpolyploid genome diploidization in Brassicaceae.

Idahoa and Subularia: Hidden polyploid origins of two enigmatic genera of crucifers

PREMISE

The monotypic Idahoa (I. scapigera) and the bispecific Subularia (S. aquatica and S. monticola) belong to Brassicaceae with unclear phylogenetic relationships and no tribal assignment. To fill this knowledge gap, we investigated these species and their closest relatives by combining cytogenomic and phylogenomic methods.

METHODS

We used whole plastome sequences in maximum likelihood and Bayesian inference analyses. We tested the phylogenetic informativeness of shared genomic repeats. We combined nuclear gene tree reconciliation and comparative chromosome painting (CCP) to examine the occurrence of past whole-genome duplications (WGDs).

RESULTS

The plastid data set corroborated the sister relationship between Idahoa and Subularia within the crucifer Lineage V but failed to resolve consistent topologies using both inference methods. The shared repetitive sequences provided conflicting pwhylogenetic signals. CCP analysis unexpectedly revealed that Idahoa (2n = 16) has a diploidized mesotetraploid genome, whereas two Subularia species (2n = 28 and 30) have diploidized mesoctoploid genomes. Several ancient allopolyploidy events have also been detected in closely related taxa (Chamira circaeoides, Cremolobeae, Eudemeae, and Notothlaspideae).

CONCLUSIONS

Our results suggest that the contentious phylogenetic placement of Idahoa and Subularia is best explained by two WGDs involving one or more shared parental genomes. The newly identified mesopolyploid genomes highlight the challenges of studying plant clades with complex polyploidy histories and provide a better framework for understanding genome evolution in the crucifer family.

Comparative repeatome analysis reveals new evidence on genome evolution in wild diploid Arachis (Fabaceae) species

Opposing changes in the abundance of satellite DNA and long terminal repeat (LTR) retroelements are the main contributors to the variation in genome size and heterochromatin amount in Arachis diploids. The South American genus Arachis (Fabaceae) comprises 83 species organized in nine taxonomic sections. Among them, section Arachis is characterized by species with a wide genome and karyotype diversity. Such diversity is determined mainly by the amount and composition of repetitive DNA. Here we performed computational analysis on low coverage genome sequencing to infer the dynamics of changes in major repeat families that led to the differentiation of genomes in diploid species (x = 10) of genus Arachis, focusing on section Arachis. Estimated repeat content ranged from 62.50 to 71.68% of the genomes. Species with different genome composition tended to have different landscapes of repeated sequences. Athila family retrotransposons were the most abundant and variable lineage among Arachis repeatomes, with peaks of transpositional activity inferred at different times in the evolution of the species. Satellite DNAs (satDNAs) were less abundant, but differentially represented among species. High rates of evolution of an AT-rich superfamily of satDNAs led to the differential accumulation of heterochromatin in Arachis genomes. The relationship between genome size variation and the repetitive content is complex. However, largest genomes presented a higher accumulation of LTR elements and lower contents of satDNAs. In contrast, species with lowest genome sizes tended to accumulate satDNAs in detriment of LTR elements. Phylogenetic analysis based on repetitive DNA supported the genome arrangement of section Arachis. Altogether, our results provide the most comprehensive picture on the repeatome dynamics that led to the genome differentiation of Arachis species.

Revealing the Satellite DNA History in Psalidodon and Astyanax Characid Fish by Comparative Satellitomics

Eukaryotic genomes are usually enriched in repetitive DNA sequences, which can be classified as dispersed or tandemly repeated elements. Satellite DNAs are noncoding monomeric sequences organized in a head-to-tail fashion that are generally located on the subtelomeric and/or pericentromeric heterochromatin. In general, a single species incorporates a diverse group of satellite DNA families, which collection is called satellitome. Here, we characterized three new satellitomes from distinct characid fish (Psalidodon fasciatus, P. bockmanni, and Astyanax lacustris) using a combination of genomic, cytogenetic, and bioinformatic protocols. We also compared our data with the available satellitome of P. paranae. We described 57 satellite DNA (satDNA) families of P. fasciatus (80 variants), 50 of P. bockmanni (77 variants), and 33 of A. lacustris (54 variants). Our analyses demonstrated that several sequences were shared among the analyzed species, while some were restricted to two or three species. In total, we isolated 104 distinctive satDNA families present in the four species, of which 10 were shared among all four. Chromosome mapping revealed that the clustered satDNA was mainly located in the subtelomeric and pericentromeric areas. Although all Psalidodon species demonstrated the same pattern of clusterization of satDNA, the number of clusters per genome was variable, indicating a high dynamism of these sequences. In addition, our results expand the knowledge of the As51 satellite DNA family, revealing that P. bockmanni and P. paranae exhibited an abundant variant of 39 bp, while P. fasciatus showed a variant of 43 bp. The majority of satDNAs in the satellitomes analyzed here presented a common library repetitive sequence in Psalidodon and Astyanax, with abundance variations in each species, as expected for closely related groups. In addition, we concluded that the most abundant satDNA in Psalidodon (As51) passed through a diversification process in this group, resulting in new variants exclusive of Psalidodon.

Genomes, repeatomes and interphase chromosome organization in the meadowfoam family (Limnanthaceae, Brassicales)

The meadowfoam family (Limnanthaceae) is one of the smallest and genomically underexplored families of the Brassicales. The Limnanthaceae harbor about seven species in the genus Limnanthes (meadowfoam) and Floerkea proserpinacoides (false mermaidweed), all native to North America. Because all Limnanthes and Floerkea species have only five chromosome pairs, i.e., a chromosome number rare in Brassicales and shared with Arabidopsis thaliana (Arabidopsis), we examined the Limnanthaceae genomes as a potential model system. Using low-coverage whole-genome sequencing data, we reexamined phylogenetic relationships and characterized the repeatomes of Limnanthaceae genomes. Phylogenies based on complete chloroplast and 35S rDNA sequences corroborated the sister relationship between Floerkea and Limnanthes and two major clades in the latter genus. The genome size of Limnanthaceae species ranges from 1.5 to 2.1 Gb, apparently due to the large increase in DNA repeats, which constitute 60-70% of their genomes. Repeatomes are dominated by long terminal repeat retrotransposons, while tandem repeats represent only less than 0.5% of the genomes. The average chromosome size in Limnanthaceae species (340-420 Mb) is more than 10 times larger than in Arabidopsis (32 Mb). A three-dimensional fluorescence in situ hybridization analysis demonstrated that the five chromosome pairs in interphase nuclei of Limnanthes species adopt the Rabl-like configuration.

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.

Genome size evolution in the diverse insect order Trichoptera

BACKGROUND

Genome size is implicated in the form, function, and ecological success of a species. Two principally different mechanisms are proposed as major drivers of eukaryotic genome evolution and diversity: polyploidy (i.e., whole-genome duplication) or smaller duplication events and bursts in the activity of repetitive elements. Here, we generated de novo genome assemblies of 17 caddisflies covering all major lineages of Trichoptera. Using these and previously sequenced genomes, we use caddisflies as a model for understanding genome size evolution in diverse insect lineages.

RESULTS

We detect a ∼14-fold variation in genome size across the order Trichoptera. We find strong evidence that repetitive element expansions, particularly those of transposable elements (TEs), are important drivers of large caddisfly genome sizes. Using an innovative method to examine TEs associated with universal single-copy orthologs (i.e., BUSCO genes), we find that TE expansions have a major impact on protein-coding gene regions, with TE-gene associations showing a linear relationship with increasing genome size. Intriguingly, we find that expanded genomes preferentially evolved in caddisfly clades with a higher ecological diversity (i.e., various feeding modes, diversification in variable, less stable environments).

CONCLUSION

Our findings provide a platform to test hypotheses about the potential evolutionary roles of TE activity and TE-gene associations, particularly in groups with high species, ecological, and functional diversities.

Power and Weakness of Repetition - Evaluating the Phylogenetic Signal From Repeatomes in the Family Rosaceae With Two Case Studies From Genera Prone to Polyploidy and Hybridization (Rosa and Fragaria)

Plant genomes consist, to a considerable extent, of non-coding repetitive DNA. Several studies showed that phylogenetic signals can be extracted from such repeatome data by using among-species dissimilarities from the RepeatExplorer2 pipeline as distance measures. Here, we advanced this approach by adjusting the read input for comparative clustering indirectly proportional to genome size and by summarizing all clusters into a main distance matrix subjected to Neighbor Joining algorithms and Principal Coordinate Analyses. Thus, our multivariate statistical method works as a "repeatomic fingerprint," and we proved its power and limitations by exemplarily applying it to the family Rosaceae at intrafamilial and, in the genera Fragaria and Rosa, at the intrageneric level. Since both taxa are prone to hybridization events, we wanted to show whether repeatome data are suitable to unravel the origin of natural and synthetic hybrids. In addition, we compared the results based on complete repeatomes with those from ribosomal DNA clusters only, because they represent one of the most widely used barcoding markers. Our results demonstrated that repeatome data contained a clear phylogenetic signal supporting the current subfamilial classification within Rosaceae. Accordingly, the well-accepted major evolutionary lineages within Fragaria were distinguished, and hybrids showed intermediate positions between parental species in data sets retrieved from both complete repeatomes and rDNA clusters. Within the taxonomically more complicated and particularly frequently hybridizing genus Rosa, we detected rather weak phylogenetic signals but surprisingly found a geographic pattern at a population scale. In sum, our method revealed promising results at larger taxonomic scales as well as within taxa with manageable levels of reticulation, but success remained rather taxon specific. Since repeatomes can be technically easy and comparably inexpensively retrieved even from samples of rather poor DNA quality, our phylogenomic method serves as a valuable alternative when high-quality genomes are unavailable, for example, in the case of old museum specimens.

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.

Comparative Analysis of Transposable Elements in Genus Calliptamus Grasshoppers Revealed That Satellite DNA Contributes to Genome Size Variation

Simple Summary

Calliptamus is a genus of grasshoppers belonging to the family Acrididae. The genus Calliptamus includes approximately 17 recognized species. Calliptamus abbreviatus, Calliptamus italicus, and Calliptamus barbarus are three species that are widely found in northern China. These species are polyphagous, feeding on a variety of wild plants as well as crops, particularly legumes. The genome sizes, phylogenetic position, and transcriptome analysis of the genus Calliptamus were already known previous to this research. The repeatome analysis of these species was missing, which is directly linked to the larger genome sizes of the grasshoppers. Here, we classified repetitive DNA sequences at the level of superfamilies and sub-families, and found that LINE, TcMar-Tc1 and Ty3-gypsy LTR retrotransposons dominated the repeatomes of all genomes, accounting for 16–34% of the total genomes of these species. Satellite DNA dynamic evolutionary changes in all three genomes played a role in genome size evolution. This study would be a valuable source for future genome assemblies.

Abstract

Transposable elements (TEs) play a significant role in both eukaryotes and prokaryotes genome size evolution, structural changes, duplication, and functional variabilities. However, the large number of different repetitive DNA has hindered the process of assembling reference genomes, and the genus level TEs diversification of the grasshopper massive genomes is still under investigation. The genus Calliptamus diverged from Peripolus around 17 mya and its species divergence dated back about 8.5 mya, but their genome size shows rather large differences. Here, we used low-coverage Illumina unassembled short reads to investigate the effects of evolutionary dynamics of satDNAs and TEs on genome size variations. The Repeatexplorer2 analysis with 0.5X data resulted in 52%, 56%, and 55% as repetitive elements in the genomes of Calliptamus barbarus, Calliptamus italicus, and Calliptamus abbreviatus, respectively. The LINE and Ty3-gypsy LTR retrotransposons and TcMar-Tc1 dominated the repeatomes of all genomes, accounting for 16–35% of the total genomes of these species. Comparative analysis unveiled that most of the transposable elements (TEs) except satDNAs were highly conserved across three genomes in the genus Calliptamus grasshoppers. Out of a total of 20 satDNA families, 17 satDNA families were commonly shared with minor variations in abundance and divergence between three genomes, and 3 were Calliptamus barbarus specific. Our findings suggest that there is a significant amplification or contraction of satDNAs at genus phylogeny which is the main cause that made genome size different.

Cytogenomic analysis unveils mixed molecular evolution and recurrent chromosomal rearrangements shaping the multigene families on Schistocerca grasshopper genomes

Multigene families are essential components of eukaryotic genomes and play key roles either structurally and functionally. Their modes of evolution remain elusive even in the era of genomics, because multiple multigene family sequences coexist in genomes, particularly in large repetitive genomes. Here, we investigate how the multigene families 18S rDNA, U2 snDNA, and H3 histone evolved in 10 species of Schistocerca grasshoppers with very large and repeat-enriched genomes. Using sequenced genomes and fluorescence in situ hybridization mapping, we find substantial differences between species, including the number of chromosomal clusters, changes in sequence abundance and nucleotide composition, pseudogenization, and association with transposable elements (TEs). The intragenomic analysis of Schistocerca gregaria using long-read sequencing and genome assembly unveils conservation for H3 histone and recurrent pseudogenization for 18S rDNA and U2 snDNA, likely promoted by association with TEs and sequence truncation. Remarkably, TEs were frequently associated with truncated copies, were also among the most abundant in the genome, and revealed signatures of recent activity. Our findings suggest a combined effect of concerted and birth-and-death models driving the evolution of multigene families in Schistocerca over the last 8 million years, and the occurrence of intra- and interchromosomal rearrangements shaping their chromosomal distribution. Despite the conserved karyotype in Schistocerca, our analysis highlights the extensive reorganization of repetitive DNAs in Schistocerca, contributing to the advance of comparative genomics for this important grasshopper genus.

The B Chromosomes of Prochilodus lineatus (Teleostei, Characiformes) Are Highly Enriched in Satellite DNAs

B or supernumerary chromosomes are dispensable elements that are widely present in numerous eukaryotes. Due to their non-recombining nature, there is an evident tendency for repetitive DNA accumulation in these elements. Thus, satellite DNA plays an important role in the evolution and diversification of B chromosomes and can provide clues regarding their origin. The characiform Prochilodus lineatus was one of the first discovered fish species bearing B chromosomes, with all populations analyzed so far showing one to nine micro-B chromosomes and exhibiting at least three morphological variants (Ba, Bsm, and Bm). To date, a single satellite DNA is known to be located on the B chromosomes of this species, but no information regarding the differentiation of the proposed B-types is available. Here, we characterized the satellitome of P. lineatus and mapped 35 satellite DNAs against the chromosomes of P. lineatus, of which six were equally located on all B-types and this indicates a similar genomic content. In addition, we describe, for the first time, an entire population without B chromosomes.

A Long-Term Conserved Satellite DNA That Remains Unexpanded in Several Genomes of Characiformes Fish Is Actively Transcribed

Eukaryotic genomes contain large amounts of repetitive DNA sequences, such as tandemly repeated satellite DNAs (satDNAs). These sequences are highly dynamic and tend to be genus- or species-specific due to their particular evolutionary pathways, although there are few unusual cases of conserved satDNAs over long periods of time. Here, we used multiple approaches to reveal that an satDNA named CharSat01-52 originated in the last common ancestor of Characoidei fish, a superfamily within the Characiformes order, ∼140-78 Ma, whereas its nucleotide composition has remained considerably conserved in several taxa. We show that 14 distantly related species within Characoidei share the presence of this satDNA, which is highly amplified and clustered in subtelomeric regions in a single species (Characidium gomesi), while remained organized as small clusters in all the other species. Defying predictions of the molecular drive of satellite evolution, CharSat01-52 shows similar values of intra- and interspecific divergence. Although we did not provide evidence for a specific functional role of CharSat01-52, its transcriptional activity was demonstrated in different species. In addition, we identified short tandem arrays of CharSat01-52 embedded within single-molecule real-time long reads of Astyanax paranae (536 bp-3.1 kb) and A. mexicanus (501 bp-3.9 kb). Such arrays consisted of head-to-tail repeats and could be found interspersed with other sequences, inverted sequences, or neighbored by other satellites. Our results provide a detailed characterization of an old and conserved satDNA, challenging general predictions of satDNA evolution.

Draft Genome Assemblies and Annotations of Agrypnia vestita Walker, and Hesperophylax magnus Banks Reveal Substantial Repetitive Element Expansion in Tube Case-Making Caddisflies (Insecta: Trichoptera)

Trichoptera (caddisflies) play an essential role in freshwater ecosystems; for instance, larvae process organic material from the water and are food for a variety of predators. Knowledge on the genomic diversity of caddisflies can facilitate comparative and phylogenetic studies thereby allowing scientists to better understand the evolutionary history of caddisflies. Although Trichoptera are the most diverse aquatic insect order, they remain poorly represented in terms of genomic resources. To date, all long-read based genomes have been sequenced from individuals in the retreat-making suborder, Annulipalpia, leaving ∼275 Ma of evolution without high-quality genomic resources. Here, we report the first long-read based de novo genome assemblies of two tube case-making Trichoptera from the suborder Integripalpia, Agrypnia vestita Walker and Hesperophylax magnus Banks. We find that these tube case-making caddisflies have genome sizes that are at least 3-fold larger than those of currently sequenced annulipalpian genomes and that this pattern is at least partly driven by major expansion of repetitive elements. In H. magnus, long interspersed nuclear elements alone exceed the entire genome size of some annulipalpian counterparts suggesting that caddisflies have high potential as a model for understanding genome size evolution in diverse insect lineages.