De novo identification of satellite DNAs in the sequenced genomes of Drosophila virilis and D. americana using the RepeatExplorer and TAREAN pipelines

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.

High-Quality Genome Assemblies Reveal Evolutionary Dynamics of Repetitive DNA and Structural Rearrangements in the Drosophila virilis Subgroup

High-quality genome assemblies across a range of nontraditional model organisms can accelerate the discovery of novel aspects of genome evolution. The Drosophila virilis group has several attributes that distinguish it from more highly studied species in the Drosophila genus, such as an unusual abundance of repetitive elements and extensive karyotype evolution, in addition to being an attractive model for speciation genetics. Here, we used long-read sequencing to assemble five genomes of three virilis group species and characterized sequence and structural divergence and repetitive DNA evolution. We find that our contiguous genome assemblies allow characterization of chromosomal arrangements with ease and can facilitate analysis of inversion breakpoints. We also leverage a small panel of resequenced strains to explore the genomic pattern of divergence and polymorphism in this species and show that known demographic histories largely predicts the extent of genome-wide segregating polymorphism. We further find that a neo-X chromosome in Drosophila americana displays X-like levels of nucleotide diversity. We also found that unusual repetitive elements were responsible for much of the divergence in genome composition among species. Helitron-derived tandem repeats tripled in abundance on the Y chromosome in D. americana compared to Drosophila novamexicana, accounting for most of the difference in repeat content between these sister species. Repeats with characteristics of both transposable elements and satellite DNAs expanded by 3-fold, mostly in euchromatin, in both D. americana and D. novamexicana compared to D. virilis. Our results represent a major advance in our understanding of genome biology in this emerging model clade.

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.

High quality genome assemblies reveal evolutionary dynamics of repetitive DNA and structural rearrangements in the Drosophila virilis sub-group

High-quality genome assemblies across a range of non-traditional model organisms can accelerate the discovery of novel aspects of genome evolution. The Drosophila virilis group has several attributes that distinguish it from more highly studied species in the Drosophila genus, such as an unusual abundance of repetitive elements and extensive karyotype evolution, in addition to being an attractive model for speciation genetics. Here we used long-read sequencing to assemble five genomes of three virilis group species and characterized sequence and structural divergence and repetitive DNA evolution. We find that our contiguous genome assemblies allow characterization of chromosomal arrangements with ease and can facilitate analysis of inversion breakpoints. We also leverage a small panel of resequenced strains to explore the genomic pattern of divergence and polymorphism in this species and show that known demographic histories largely predicts the extent of genome-wide segregating polymorphism. We further find that a neo-X chromosome in D. americana displays X-like levels of nucleotide diversity. We also found that unusual repetitive elements were responsible for much of the divergence in genome composition among species. Helitron-derived tandem repeats tripled in abundance on the Y chromosome in D. americana compared to D. novamexicana, accounting for most of the difference in repeat content between these sister species. Repeats with characteristics of both transposable elements and satellite DNAs expanded by three-fold, mostly in euchromatin, in both D. americana and D. novamexicana compared to D. virilis. Our results represent a major advance in our understanding of genome biology in this emerging model clade.

Comparison of the evolutionary patterns of DNA repeats in ancient and young invertebrate species flocks of Lake Baikal

DNA repeat composition of low coverage (0.1-0.5) genomic libraries of four amphipods species endemic to Lake Baikal (East Siberia) and four endemic gastropod species of the fam. Baicaliidae have been compared to each other. In order to do so, a neighbor joining tree was inferred for each quartet of species (amphipods and mollusks) based on the ratio of repeat classes shared in each pair of species. The topology of this tree was compared to the phylogenies inferred for the same species from the concatenated protein-coding mitochondrial nucleotide sequences. In all species analyzed, the fraction of DNA repeats involved circa half of the genome. In relatively more ancient amphipods (most recent common ancestor, MRCA, existed approximately sixty millions years ago), the most abundant were species-specific repeats, while in much younger Baicaliidae (MRCA equal to ca. three millions years) most of the DNA repeats were shared among all four species. If the presence/absence of a repeat is regarded as a separate independent trait, and the ratio of shared to total numbers of repeats in a species pair is used as the measure of distance, the topology of the NJ tree is the same as the quartet phylogeny inferred for the mitogenomes protein coding nucleotide sequences. Meanwhile, in each group of species, a substantial number of repeats were detected pointing to the possibility of non-neutral evolution or a horizontal transfer between species occupying the same biotope. These repeats were shared by non-sister groups while being absent in the sister genomes. On the other hand, in such cases some traits of ecological significance were also shared.

Structure, Organization, and Evolution of Satellite DNAs: Insights from the Drosophila repleta and D. virilis Species Groups

The fact that satellite DNAs (satDNAs) in eukaryotes are abundant genomic components, can perform functional roles, but can also change rapidly across species while being homogenous within a species, makes them an intriguing and fascinating genomic component to study. It is also becoming clear that satDNAs represent an important piece in genome architecture and that changes in their structure, organization, and abundance can affect the evolution of genomes and species in many ways. Since the discovery of satDNAs more than 50 years ago, species from the Drosophila genus have continuously been used as models to study several aspects of satDNA biology. These studies have been largely concentrated in D. melanogaster and closely related species from the Sophophora subgenus, even though the vast majority of all Drosophila species belong to the Drosophila subgenus. This chapter highlights some studies on the satDNA structure, organization, and evolution in two species groups from the Drosophila subgenus: the repleta and virilis groups. We also discuss and review the classification of other abundant tandem repeats found in these species in the light of the current information available.

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.

Identification and characterization of repetitive DNA in the genus Didelphis Linnaeus, 1758 (Didelphimorphia, Didelphidae) and the use of satellite DNAs as phylogenetic markers

Didelphis species have been shown to exhibit very conservative karyotypes, which mainly differ in their constitutive heterochromatin, known to be mostly composed by repetitive DNAs. In this study, we used genome skimming data combined with computational pipelines to identify the most abundant repetitive DNA families of Lutreolina crassicaudata and all six Didelphis species. We found that transposable elements (TEs), particularly LINE-1, endogenous retroviruses, and SINEs, are the most abundant mobile elements in the studied species. Despite overall similar TE proportions, we report that species of the D. albiventris group consistently present a less diverse TE composition and smaller proportions of LINEs and LTRs in their genomes than other studied species. We also identified four new putative satDNAs (sat206, sat907, sat1430 and sat2324) in the genomes of Didelphis species, which show differences in abundance and nucleotide composition. Phylogenies based on satDNA sequences showed well supported relationships at the species (sat1430) and groups of species (sat206) level, recovering topologies congruent with previous studies. Our study is one of the first attempts to present a characterization of the most abundant families of repetitive DNAs of Lutreolina and Didelphis species providing insights into the repetitive DNA composition in the genome landscape of American marsupials.

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

Tandem repeats are important parts of eukaryotic genomes being crucial e.g., for centromere and telomere function and chromatin modulation. In Lepidoptera, knowledge of tandem repeats is very limited despite the growing number of sequenced genomes. Here we introduce seven new satellite DNAs (satDNAs), which more than doubles the number of currently known lepidopteran satDNAs. The satDNAs were identified in genomes of three species of Crambidae moths, namely Ostrinia nubilalis, Cydalima perspectalis, and Diatraea postlineella, using graph-based computational pipeline RepeatExplorer. These repeats varied in their abundance and showed high variability within and between species, although some degree of conservation was noted. The satDNAs showed a scattered distribution, often on both autosomes and sex chromosomes, with the exception of both satellites in D. postlineella, in which the satDNAs were located at a single autosomal locus. Three satDNAs were abundant on the W chromosomes of O. nubilalis and C. perspectalis, thus contributing to their differentiation from the Z chromosomes. To provide background for the in situ localization of the satDNAs, we performed a detailed cytogenetic analysis of the karyotypes of all three species. This comparative analysis revealed differences in chromosome number, number and location of rDNA clusters, and molecular differentiation of sex chromosomes.

Sequence Composition Underlying Centromeric and Heterochromatic Genome Compartments of the Pacific Oyster Crassostrea gigas

Segments of the genome enriched in repetitive sequences still present a challenge and are omitted in genome assemblies. For that reason, the exact composition of DNA sequences underlying the heterochromatic regions and the active centromeres are still unexplored for many organisms. The centromere is a crucial region of eukaryotic chromosomes responsible for the accurate segregation of genetic material. The typical landmark of centromere chromatin is the rapidly-evolving variant of the histone H3, CenH3, while DNA sequences packed in constitutive heterochromatin are associated with H3K9me3-modified histones. In the Pacific oyster Crassostrea gigas we identified its centromere histone variant, Cg-CenH3, that shows stage-specific distribution in gonadal cells. In order to investigate the DNA composition of genomic regions associated with the two specific chromatin types, we employed chromatin immunoprecipitation followed by high-throughput next-generation sequencing of the Cg-CenH3- and H3K9me3-associated sequences. CenH3-associated sequences were assigned to six groups of repetitive elements, while H3K9me3-associated-ones were assigned only to three. Those associated with CenH3 indicate the lack of uniformity in the chromosomal distribution of sequences building the centromeres, being also in the same time dispersed throughout the genome. The heterochromatin of C. gigas exhibited general paucity and limited chromosomal localization as predicted, with H3K9me3-associated sequences being predominantly constituted of DNA transposons.